矩形周期脉冲信号MATLAB实现

下面是小编为大家整理的矩形周期脉冲信号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("Real Part of X(k)");

　subplot(2,1,2),stem(n,xi,".");grid;

　xlabel("k"),ylabel("Imaginary Part of X(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("Magnitude Part of X(k)");

　subplot(2,1,2),stem(n,y,".");grid;

　xlabel("k"),ylabel("Phase Part of X(k)");

　 ①矩形脉冲周期信号的时域波形； -10 -8 -6 -4 -2 0 2 4 6 8 1000.20.40.60.81矩 形 脉 冲 周 期 信 号tf(t) ②矩形脉冲周期信号的频谱：实部和虚部； -30 -20 -10 0 10 20 30-0.200.20.40.6kReal Part of X(k)-30 -20 -10 0 10 20 30-1-0.500.51kReal Part of X(k) ③矩形脉冲周期信号的频谱：幅值和相位。

　-30 -20 -10 0 10 20 3000.10.20.30.4kMagnitude Part of X(k)-30 -20 -10 0 10 20 3001234kPhase Part of X(k)