用MATLAB求周期矩形脉冲幅频谱图

用Matlab编程，绘制出周期矩形脉冲信号的幅频谱图，其中周期矩形脉冲的周期为T=1ms，脉冲宽度为τ=0.1ms，如下图所示TT220t)(tx1T2T2解：用MATLAB编写的程序如下：fs=150;N=200;%采样频率为150Hz，点数为200点Dt=1/fs;%数字角频率t=(0:N-1)*Dt;A=1;P=2*pi;y=A*square(P*(t+0.05),10);%周期脉冲函数subplot(2,1,1)plot(t,y);%周期脉冲波形axis([01.501.5])xlabel('时间')ylabel('幅值')title('周期矩形脉冲信号')Y=fft(y);%快速傅里叶变换fy=(abs(Y)*2)/N;%频谱subplot(2,1,2)stem(fy);axis([03002])xlabel('频率')ylabel('幅值')title('周期矩形脉冲信号的幅频谱图')00.511.500.511.5时间幅值周期矩形脉冲信号05101520253000.511.52频率幅值周期矩形脉冲信号的幅频谱图