mybilinear.m

来自「按照给定指标采用模拟低通滤波器（巴特沃思低通）去逼近理想模拟低通滤波器」· M 代码 · 共 9 行

function [b,a]=mybilinear(z,p,k,Fs)   
%Fs=2*采样频率
zk=(1+p/Fs)./(1-p/Fs);          %通过双线性Z变换求得对应极点zk
zd=(1+z/Fs)./(1-z/Fs);
KZ=(k*prod(Fs-z)./prod(Fs-p));
zd=[zd;-ones(length(p)-length(z),1)];
[b,a]=zp2tf(zd,zk,KZ);