信号分析与处理第一次上机作业
例:矩形脉冲周期信号频谱分析的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)');