Butterworth滤波器
wp=100*2*pi;ws=200*2*pi;
ap=2;
as=15;
Fs=500;
=buttord (wp,ws,ap,as,'s'); %选择滤波器的最小阶数
=buttap(N); %创建Butterworth低通滤波器原型
=zp2ss(Z,P,K); %零极点增益模型转换为状态空间模型
=lp2lp(A,B,C,D,Wc); %实现低通向低通的转变
=ss2tf(AT,BT,CT,DT); %状态空间模型转换为传递函数模型
%运用双线性变换法把模拟滤波器转换成数字滤波器
=bilinear(num1,den1,500);
=freqz(num2,den2); %求频率响应
plot(W*Fs/(2*pi),abs(H));grid; %绘出频率响应曲线
xlabel('频率/Hz');ylabel('幅值')
其中:
=zp2ss(Z,P,K); %零极点增益模型转换为状态空间模型
=lp2lp(A,B,C,D,Wc); %实现低通向低通的转变
=ss2tf(AT,BT,CT,DT); %状态空间模型转换为传递函数模型
这三行虽然有注释,但是我还是看不懂,原因应该是对Butterworth设计的原理思路不懂,请高手帮忙啊
[ 本帖最后由 eight 于 2007-7-16 11:15 编辑 ] 滤波器设计的例子论坛很多,请搜索论坛。原理不懂,就看一下书了。缺啥补啥了,反正是为了学知识。如果理论不扎实,看英文的help文档会很费劲了。
页:
[1]