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);		% 去除運算誤差所產生的微小虛數