使用Matlab软件画出频谱图
Matlab参考程序如下所示:
M=4; N=2*M+1; T=0.5; n=-4*M:4*M;
x = [zeros(1,3*M),ones(1,N),zeros(1,3*M)]; % 给出输入序列
w = [-15:0.1:15] + 1e-10;
X = sin(0.5*N*w*T)./sin(0.5*w*T); % 给出频谱序列
subplot(1,3,1); stem(n,x,'.'); % 画出输入序列
axis([-20,20,-0.1,1.1]), grid on
xlabel('n'); title('(a) 序列幅度');
subplot(1,3,2); plot(w,X),grid on % 画出频谱序列
xlabel('\Omega'); title('(b) 幅频特性');
subplot(1,3,3); plot(w,X),grid on % 改变横轴比例,画出频谱序列
v = axis;axis([-pi/T,pi/T,v(3),v