bandpiir.m

双音多频编解码系统

%bandpIIR.m函数文件，通过调用butter()函数，重新设计数字带通巴特沃斯滤波器
function [Bn,An,N]=bandpIIR(alfa,omega,fs);

%预畸变求模拟带通滤波器参数
ome=2*tan(omega/2);
omer=ome(2)-ome(1);
lam=ome/omer;
lam0=sqrt(lam(1)*lam(2));
ome2=lam-lam0^2*lam.^(-1);

%确定模拟低通原型阶数和截止频率
[N0,omgc]=buttord(ome2(2),ome2(4),alfa(1),alfa(2),'s');
%确定满足指标要求的Hdn(z)
[Bn,An]=butter(N0,[omega(1)/pi omega(2)/pi]);
N=2*N0;