ex31a.m

来自「数字信号分析处理方面的好教材,有DSP的MATLAB源程序.」· M 代码 · 共 14 行

%
x = [1,1,1,1];
subplot(1,1,1);
%
% a) 付里叶变换
w = [0:1:500]*2*pi/500;  %取Δt=1/500
[H] = freqz(x,1,w);
magH = abs(H); phaH = angle(H); phaH(126)=-47.5841*pi/180;
subplot(2,1,1); plot(w/pi,magH); grid
xlabel('');
ylabel('|X|'); title(' 频谱的幅度')
subplot(2,1,2); plot(w/pi,phaH/pi*180); grid
xlabel('以pi为单位的频率');
ylabel('度'); title('频谱的相角')