hc834.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第八章例8.3.4程序hc834
% 双线性变换法模拟系统与数字系统的比较
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
ba=[1,1]; aa=[1,5,6];                 % 系统分子分母系数向量
T=0.1; Fs=1/T; 
[bd,ad] = bilinear0(ba,aa,Fs);           % 双线性变换求等价数字滤波器系数
%[bd,ad] = bilinear(ba,aa,Fs);          % 调用信号处理工具箱中的双线性变换函数
t=0:0.1:3;
subplot(2,1,1)
ha=impulse(ba,aa,t);plot(t,ha*T),hold on		% 计算模拟系统的脉冲响应，画曲线
h=impz(bd,ad,31);stem(0.1*[0:30],h,'.') 	% 数字系统的脉冲响应，以同样时间轴画点列
legend('原模拟滤波器的脉冲响应ha(t)乘以T','等价数字滤波器的脉冲响应he(n)')
w=[0:0.1:10]*2*pi;
subplot(2,1,2)
Ha=freqs(ba,aa,w); plot(w,abs(Ha)),hold on	% 计算模拟系统的频率特性，画曲线
H=freqz(bd,ad,w/Fs); plot(w/Fs,abs(H),'.')	 
plot(w,abs(H),'-.')	   % 计算数字系统的频率特性，画曲线
legend('原模拟滤波器的频率响应|Ha(\Omega)|','等价数字滤波器的频率响应|He(\omega)|','频率轴放大Fs倍后数字响应')
set(gcf,'color','w')              % 置图形背景色为白