hc543.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第四章例5.4.3程序hc543
% 非单位圆上螺旋线Chirp的频谱分析
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
p1=0.9*exp(j*pi/3);p2=0.9*exp(-j*pi/3);
%a=conv([1,-p1],[1,-p2]);b=1;
%[r,p,k]=residuez(b,a)
M=64,n=0:M-1;dw=pi/M;                         % 规定点数，求频率分辨率
x=0.9.^n.*(cos(pi*n/3)+2*0.2887*sin(pi*n/3)); % 给定x序列
subplot(3,2,1),stem(n,x,'.')                  % 画序列曲线
H=freqz(x,1,M);                               % 求序列频率响应
subplot(3,2,2),plot([0:M-1]*dw,abs(H))        % 画H的幅频特性
x1=x.*(0.9).^-n; H1=freqz(x1,1,64);           % 求r=0.9的中心圆上的H(z)，先修改x
subplot(3,2,3),plot([0:M-1]*dw,abs(H1))       % 画H1的幅频特性
A=1;W0=0.9^(-3/M),W=W0*exp(-j*dw);            % 确定螺旋线轨迹参数A和W
z = A * W.^(-(0:M-1));                        % 写出螺旋线轨迹方程
subplot(3,2,4),plot(z),grid on                % 画出螺旋线轨迹
axis([-1.2,1.2,-1.2,1.2]),axis equal
hold on,plot(p1,'x')                        % 把序列极点画在螺旋线轨迹图上
G=czt(x,M,W,A);                               % 调用czt函数
subplot(3,2,5),plot([0:M-1]*dw,abs(G(1:M)))   % 画出czt函数的幅特性
subplot(3,2,6),plot([0:M/2-1]*dw,abs(G(1:M/2)))  % 画前半段czt的幅特性益3显示峰值
set(gcf,'color','w');                           % 置图形背景色为白