mybilinear.m

按照给定指标采用模拟低通滤波器（巴特沃思低通）去逼近理想模拟低通滤波器

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