hc331.m

某本matlab书籍中的代码

%《数字信号处理教程——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])