q616.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第六章例6.16程序q616
% 二阶巴特沃斯低通滤波器的频率响应
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
% 二阶巴特沃斯低通滤波器的离散系统函数为：
% H(z)=(z^2+2z+1)/((2+sqrt(2))z^2+(2+sqrt(2));
% 求其频率响应可将z用exp(i*w)代入

b=[1,2,1]; a=[2+sqrt(2),0,2-sqrt(2)];	% 给出滤波器分子分母系数
N=input('取频率数组的点数N= ');		
w=[0:N-1]*pi/N;							% 给出0到pi之间的频率数组
H=polyval(b,exp(i*w))./polyval(a,exp(i*w));	% 求频率响应
figure(1)							    % 在线性坐标内画频率特性
set(gcf,'color','w')                    % 设置图形背景色为白色
subplot(2,1,1),plot(w,abs(H)),grid on
subplot(2,1,2),plot(w,unwrap(angle(H)))
grid on
figure(2)							    % 在对数坐标内画频率特性
set(gcf,'color','w')                    % 设置图形背景色为白色
subplot(2,1,1),
semilogx(w,20*log10(abs(H))),grid on    % 画幅频特性
subplot(2,1,2)
semilogx(w,unwrap(angle(H))),grid on    % 画相频特性