举个例子,希望有所帮助.\x0d代码clc; clear all; close all;\x0dwp=5000*2*pi;\x0dws=15000*2*pi;ap=1;as=70;Fs=15000;Fp=9000;%选择滤波器的最小阶数\x0d[N,Wc]=buttord (wp,ws,ap,as,'s');\x0d%创建Butterworth低通滤波器原型\x0d[Z,P,K]=buttap(N);\x0d%零极点增益模型转换为状态空间模型\x0d[A,B,C,D]=zp2ss(Z,P,K);\x0d%实现低通向低通的转变\x0d[AT,BT,CT,DT]=lp2lp(A,B,C,D,Wc);\x0d%状态空间模型转换为传递函数模型\x0d[num1,den1]=ss2tf(AT,BT,CT,DT);\x0d%运用双线性变换法把模拟滤波器转换成数字滤波器\x0d[num2,den2]=bilinear(num1,den1,35000);\x0d%求频率响应\x0d[H,W]=freqz(num2,den2);\x0d%绘出频率响应曲线figure;plot(W*Fs/(2*pi),abs(H));grid;\x0dxlabel('频率/Hz');ylabel('幅值');结果
