%プログラム MATLAB2.3 の補足 %（離散時間フーリエ変換の振幅スペクトルと位相スペクトル） w = -pi:0.01:pi; %周波数の範囲 x = 1/8*[1 1 1 1 1 1 1 1]; %離散時間信号 X = freqz(x,1,w); %離散時間フーリエ変換 figure(1) plot(w,abs(X)); %振幅スペクトルの図示 axis([-pi pi 0 1]); grid; xlabel('Frequency \omega [rad]'); ylabel('|X(e^{j\omega})|'); figure(2) plot(w,angle(X)); %位相スペクトルの図示 axis([-pi pi -pi pi]); grid; xlabel('Frequency \omega [rad]'); ylabel('\angle X(e^{j\omega})');
