hc542.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第四章图5.4.2程序hcf542
% 小于单位圆圆周上的频率特性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
n=0:63;p=0.91;
x=p.^n;                         % 输入序列
subplot(3,2,1),stem(n,x,'.')    % 画输入序列图
title('(a) x=(0.91)^n')
xlabel('n'),ylabel('x(n)')
H=freqz(x,1,64);                % 求序列的（单位圆）频率特性
subplot(3,2,3),plot(n*pi/64,abs(H))   % 画幅频特性图
title('(c) 单位圆上的幅频特性')
xlabel('\omega'),ylabel('|X(\omega)|')
subplot(3,2,5),plot(n*pi/64,angle(H))   % 画相频特性图
title('(e) 单位圆上的相频特性')
xlabel('\omega'),ylabel('angle(X)')
x1=x.*(0.9).^-n;                % 求出对非单位圆折算的输入序列
H1=freqz(x1,1,64);      % 求折算序列的（单位圆）频率特性
subplot(3,2,2),stem(n,x1,'.')
title('(b) x1=x.*(0.9)^-n')
xlabel('n'),ylabel('x1(n)')
subplot(3,2,4),plot(n*pi/64,abs(H1))   % 画非单位圆幅频特性图
title('(d) r=0.9的圆上的幅频特性')
xlabel('\omega'),ylabel('|X1(\omega)|')
subplot(3,2,6),plot(n*pi/64,angle(H1))   % 画非单位圆相频特性图
title('(f) r=0.9的圆上的相频特性')
xlabel('\omega'),ylabel('angle(X1)')
set(gcf,'color','w');                % 置图形背景色为白