hc831.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第八章例8.3.1程序hc831
% 脉冲响应不变法模拟系统与数字系统的比较
% 电子工业出版社出版  陈怀琛编著 2004年8月
%
clear,
ba=[1,1];aa=[1,5,6];        % 系统分子分母系数向量
T=0.1;Fs=1/T;
[Ra,pa,ha]=residue(ba, aa);	% 将模拟滤波器系数向量变为模拟极点和留数
pd=exp(pa*T);			% 将模拟极点变为数字（z平面）极点pd
[bd,ad]=residuez(T*Ra, pd, ha)	% 用原留数Ra和数字极点pd求得数字滤波器系数
t=0:0.1:3;
subplot(2,1,1)
set(gcf,'color','w')              % 置图形背景色为白
ha=impulse(ba,aa,t);     % 计算模拟系统的脉冲响应，画曲线
plot(t,ha*T),hold on,grid on
pause
hd=impz(bd,ad,31);      % 数字系统的脉冲响应，以同样时间轴画点列
stem(0.1*[0:30],hd,'.')
legend('原模拟滤波器的脉冲响应ha(t)乘以T','等价数字滤波器的脉冲响应he(n)')
subplot(2,1,2)
w=[0:0.1:10]*2*pi;
H=freqs(ba,aa,w);       % 计算模拟系统的频率特性，画曲线
plot(w,abs(H)),hold on
pause
Hd=freqz(bd,ad,w/Fs); 
plot(w/Fs,abs(Hd),'.')	% 计算数字系统的频率特性，画曲线
plot(w,abs(Hd),'-.')				% 将数字频率线性放大与模拟频率的幅特性比较
legend('原模拟滤波器的频率响应|Ha(\Omega)|','等价数字滤波器的频率响应|He(\omega)|','频率轴放大Fs倍后数字响应')
axis([0.,50,0,0.3]),grid on
[bd1,ad1] = impinvar(ba,aa,1/T) % 与MATLAB中的脉冲响应不变法函数比较