fgp821.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第八章图8.2.1生成程序fgp821
% 巴特沃思滤波器的极点和平方幅度响应图生成
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
wc=1;w=-2:0.01:2;                   % 画频谱的自变量范围
M=1./(1+(ones(3,1)*w).^(2*[1;2;10]*ones(size(w))));% 巴特沃思滤波器的幅特性
figure(1),plot(w,M)                 % 绘幅特性图
axis([0,2,0,1]),grid on
ylabel('J=|Ha(j\Omega)|^2'),xlabel('\Omega/\Omegac')    % 标注
set(gcf,'color','w')
figure(2),
N=3;k=0:2*N-1;wc=1;                 % 设定巴特沃思滤波器的参数
subplot(1,2,1),
polar((k+0.5*mod(N+1,2))*pi/N,wc*ones(1,2*N),'*')   % 画N=3的极点图
title('N=3')
subplot(1,2,2),
N=N+1;k=0:2*N-1;
polar((k+0.5*mod(N+1,2))*pi/N,wc*ones(1,2*N),'*')   % 画N=4的极点图
title('N=4')
set(gcf,'color','w')                        % 置图形背景色为白