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