impinvar0.m

单个脉冲序列生成子程序 & 用脉冲响应不变法把模拟滤波器系数变换为数字滤波器系数的函数

function [bd,ad]=impinvar0(ba,aa,Fs)
%《数字信号处理教程——MATLAB释义与实现》第八章子程序impinvar
% 用脉冲响应不变法把模拟滤波器系数变换为数字滤波器系数的函数
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% [bd,ad]=impinvar0(ba,aa,Fs)
% ----------------------------------------------
% ba 为模拟系统函数的分子系数向量
% aa 为模拟系统函数的分母系数向量
% bd 为对应数字系统函数的分子系数向量
% ad 为对应数字系统函数的分母系数向量
% Fs 为采样频率(Hz)
%
[Ra,pa, Ca]=residue(ba, aa);	% 将模拟滤波器系数向量变为模拟极点和留数
T=1/Fs; pd=exp(pa*T);			% 将模拟极点变为数字（z平面）极点pd
[bd,ad]=residuez(T*Ra, pd, Ca);	% 用原留数Ra和数字极点pd求得数字滤波器系数向量
bd=real(bd); ad=real(ad);		% 去除运算误差所产生的微小虚数