function [b,a]=imp_invr(c,d,T)
%用脈沖響應不變法實現模擬濾波器到數字濾波器的變換
%b數字濾波器的分子多項式系數；a數字濾波器的分母多項式系數；
%c模擬濾波器的分子多項式系數；d模擬濾波器的分母多項式系數；T采樣參數
%
[R,p,k]=residue(c,d);
p=exp(p*T);
[b,a]=residuez(R,p,k);
b=real(b);a=real(a);