imp_invr.m

来自「数字信号处理的matlab实现代码 对电子信息专业同学的学习很有帮助」· M 代码 · 共 15 行

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');