%《數字信號處理教程——MATLAB釋義與實現》第三章例3.3.1程序hc331
% 一階離散系統的頻率響應和輸入輸出
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
b=[0.5 0.8]; 		% 分子多項式系數向量
a=[1 -0.9];		% 分母多項式系數向量
[H,w]=freqz(b,a);	% 求出頻率響應（0到pi分成500點）
figure(1)
subplot(2,1,1),plot(w,abs(H)),grid on
xlabel('頻率（1/樣本）'),ylabel('幅度dB')
subplot(2,1,2),plot(w,angle(H)),grid on
xlabel('頻率（1/樣本）'),ylabel('相角(度)')
set(gcf,'color','w')
w0=0.1*pi; 
nw0=floor(w0/pi*512)	% 求出頻率w0對應的下標nw0
Hw0=H(nw0);         % 求出頻率w0處的頻率響應
Aw0=abs(Hw0),       % 求出頻率w0處的幅特性
thetaw0=angle(Hw0)  % 求出頻率w0處的相移
pause
n=0:30;x=cos(0.1*pi*n);         % 給定輸入序列
y=Aw0*cos(0.1*pi*n-thetaw0);    % 算出系統輸出
figure(2),
subplot(2,1,1),stem(n,x,'.')    % 畫輸入序列
ylabel('輸入序列')
line([0,30],[0,0])
set(gcf,'color','w')              % 置圖形背景色為白
subplot(2,1,2),stem(n,y,'.')   % 畫輸出序列
ylabel('輸出序列')
line([0,30],[0,0])