矩形周期脉冲信号MATLAB实现

信号分析与处理第一次上机作业例：矩形脉冲周期信号频谱分析的MATLAB实现。①矩形脉冲周期信号的时域波形；②矩形脉冲周期信号的频谱：实部和虚部；③矩形脉冲周期信号的频谱：幅值和相位。MATLAB程序：t=-10:0.01:10;y=0.5*(square(0.4*pi*(t+0.5),20)+1);plot(t,y);grid;axis([-10,10,-0.1,1.2]);title('矩形脉冲周期信号'),xlabel('t'),ylabel('f(t)');n=-30:30;e=1;tao=2;zq=5;w=(2*pi)/zq;xr=(e*tao/zq).*sinc(n.*tao./zq);xi=zeros(61,1);figure(2)subplot(2,1,1),stem(n,xr,'.');grid;xlabel('k'),ylabel('RealPartofX(k)');subplot(2,1,2),stem(n,xi,'.');grid;xlabel('k'),ylabel('ImaginaryPartofX(k)');n=-30:30;e=1;tao=2;zq=5;w=(2*pi)/zq;x=abs((e*tao/zq).*sinc(n.*tao./zq));y=atan2(0,(e*tao/zq).*sinc(n.*tao./zq));figure(3)subplot(2,1,1),stem(n,x,'.');grid;xlabel('k'),ylabel('MagnitudePartofX(k)');subplot(2,1,2),stem(n,y,'.');grid;xlabel('k'),ylabel('PhasePartofX(k)');①矩形脉冲周期信号的时域波形；-10-8-6-4-2024681000.20.40.60.81矩形脉冲周期信号tf(t)②矩形脉冲周期信号的频谱：实部和虚部；-30-20-100102030-0.200.20.40.6kRealPartofX(k)-30-20-100102030-1-0.500.51kRealPartofX(k)③矩形脉冲周期信号的频谱：幅值和相位。-30-20-10010203000.10.20.30.4kMagnitudePartofX(k)-30-20-10010203001234kPhasePartofX(k)