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)');