use.m

信号处理－滤波器设计（基于matlab和Mathmatica的设计方法）中关于无乘法器椭圆IIR滤波器设计的原创程序。

Fp=0.135;Fs=0.2;Ap=0.2;As=30;
[num,den]=Low_Mufilter(Fp,Fs,Ap,As)
%总的函数执行过程如下：
%>> use
%最小阶数为n=5
%滤波器的阶数（必须为奇数）5
%a值的范围：-0.661312<α<-0.309017
%α＝-0.505971(参考值)
%选择： α=-0.5
%βmax值的范围：0.787095<βmax<0.800091
%选择：βmax=0.7891
%输入量化的顺序，其中Seq=1表示α1,Seq=2,3,...表示β2,β3,...
%  Seq=[3 1]
%β3值的范围：0.342948<β3<0.352782
%选择：β3=0.3438
%α1值的范围：-0.268877<α1<-0.256615
%选择：α1=-0.2656
% 
%Transfer function from input 1 to output:
%-0.2679 z + 1
%-------------
% z - 0.2679
%
%Transfer function from input 2 to output:
%0.7891 z^2 - 0.8945 z + 1
%-------------------------
% z^2 - 0.8945 z + 0.7891
% 
%Transfer function from input 3 to output:
%0.3468 z^2 - 0.6734 z + 1
%-------------------------
% z^2 - 0.6734 z + 0.3468
% 
%Sampling time: unspecified
%总的表达式为：
 
%Transfer function:
%0.06766 z^5 + 0.04734 z^4 + 0.1055 z^3 + 0.1055 z^2 + 0.04734 z + 0.06766
%-------------------------------------------------------------------------
%      z^5 - 1.836 z^4 + 2.158 z^3 - 1.307 z^2 + 0.4991 z - 0.07332
% 
%Sampling time: unspecified
%
%num =
%
%    0.0677    0.0473    0.1055    0.1055    0.0473    0.0677
%
%
%den =
%
%    1.0000   -1.8359    2.1584   -1.3073    0.4991   -0.0733