用MATLAB求周期矩形脉冲幅频谱图
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
A=1;P=2*pi;
y=A*square(P*(t+0.05),10);%周期脉冲函数
subplot(2,1,1)
plot(t,y);
axis([0 1.5 0 1.5])
xlabel('时间')
ylabel('幅值')
title('周期矩形脉冲信号')
Y=fft(y);%快速傅里叶变换
fy=(abs(Y)*2)/N;%频谱
subplot(2,1,2)
stem(fy);
axis([0 30 0 2])
xlabel('频率')
ylabel('幅值')
title('周期矩形脉冲信号的幅频谱图')%周期脉冲波形
周期矩形脉冲信号
1.5
1
幅
值
0.5
0
00.5ቤተ መጻሕፍቲ ባይዱ
时间
周期矩形脉冲信号的幅频谱图11.52
1.5
幅
值
1
0.5
0
051015
用Matlab编程，绘制出周期矩形脉冲信号的幅频谱图，其中周期矩形脉冲的周期为T=1ms，脉冲宽度为τ=0.1ms，如下图所示
1x(t)
2T
解：T
0
22T2Tt用MATLAB编写的程序如下：
fs=150;N=200;%采样频率为150Hz，点数为200点
Dt=1/fs;%数字角频率
t=(0:N-1)*Dt;
频率202530
y=A*square(P*(t+0.05),10);%周期脉冲函数
subplot(2,1,1)
plot(t,y);
axis([0 1.5 0 1.5])
xlabel('时间')
ylabel('幅值')
title('周期矩形脉冲信号')
Y=fft(y);%快速傅里叶变换
fy=(abs(Y)*2)/N;%频谱
subplot(2,1,2)
stem(fy);
axis([0 30 0 2])
xlabel('频率')
ylabel('幅值')
title('周期矩形脉冲信号的幅频谱图')%周期脉冲波形
周期矩形脉冲信号
1.5
1
幅
值
0.5
0
00.5ቤተ መጻሕፍቲ ባይዱ
时间
周期矩形脉冲信号的幅频谱图11.52
1.5
幅
值
1
0.5
0
051015
用Matlab编程，绘制出周期矩形脉冲信号的幅频谱图，其中周期矩形脉冲的周期为T=1ms，脉冲宽度为τ=0.1ms，如下图所示
1x(t)
2T
解：T
0
22T2Tt用MATLAB编写的程序如下：
fs=150;N=200;%采样频率为150Hz，点数为200点
Dt=1/fs;%数字角频率
t=(0:N-1)*Dt;
频率202530