q729.asv

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

%《MATLAB及在电子信息课程中的应用》第七章例7.29程序q729
% 用脉冲响应不变法和双线性变换法将模拟滤波器离散化
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
%脉冲响应不变法实现模拟到数字滤波器的转换
clear;close all
b=1000;a=[1,1000];
w=[0:1000*2*pi]; 			% 设定模拟频率
[hf,w]=freqs(b,a,w);		% 计算模拟滤波器频响函数
subplot(2,3,1)				% 画出模拟滤波器幅频特性
plot(w/2/pi,abs(hf));
grid on;title('（a)|Ha(jf)|')
set(gcf,'color','w')        % 设置图形背景色为白色
xlabel('f (Hz)');ylabel('幅度');
title('（a）模拟滤波器频响特性');
Fs0=[1000,500];
for m=1:2 
   Fs=Fs0(m)                % T=0.001s及T=0.002s
   [d,c]=impinvar(b,a,Fs)	% 用impinvar 函数实现离散化
   [f,e]=bilinear(b,a,Fs)	% 用bilinear 函数实现离散化
   wd=[0:512]*pi/512;		% 设定数字归一化频率
   hw1=freqz(d,c,wd);		% 计算数字滤波器频响函数
   hw2=freqz(f,e,wd); 
   % 画出数字滤波器幅频特性
   subplot(2,3,2);plot(wd/pi,abs(hw1)/abs(hw1(1)));
   title('（b）脉冲响应不变法');
   hold on,grid on
   subplot(2,3,3);plot(wd/pi,abs(hw2)/abs(hw2(1)));
   title('（C）双线性变换法');
   hold on;grid on
end