实验五:双线性变换法的设计IIR 数字滤波器

xlabel('t');

ylabel('sin(2*pi*100*t)'); title('时域信号');

Y=fft(x,N);%对信号进行FFT变换 magY=abs(Y);%求得FFT变换后的幅度 f=n*fs/N;%频率序列 subplot(3,1,2);

plot(f(1:N/2),magY(1:N/2));%画出幅频响应 xlabel('f'); ylabel('幅度'); title('N=2000'); grid; %滤波器

fp=90; %通带截止频率 fs=110; %阻带起始频率 Rp=3; %通带最大衰减 Rs=25;%阻带最小衰减

Wp=fp/(fn/2);%计算归一化角频率 Ws=fs/(fn/2);

[n,Wn]=buttord(Wp,Ws,Rp,Rs);%计算阶数和截止频率 [b,a]=butter(n,Wn);%计算H(z)分子、分母多项式系数

[H,F]=freqz(b,a,1000,2000);%计算H(z)的幅频响应,freqz(b,a,计算点数,采样速率) subplot(3,1,3)

plot(F,20*log10(abs(H))) %画出幅频特性图 axis([0,1000,-200,3]); y=filter(b,a,x); figure; plot(t,y);

11

12

>>闂傚倸鍊搁崐宄懊归崶顒夋晪鐟滃繘骞戦姀銈呯婵°倐鍋撶痪鎯ь煼閺岋綁骞囬锝嗏挅濠电偛妯婃禍婊堝礃閳ь剙顪冮妶鍡楀Ё缂傚秴妫楀玻鍧楁偄閸忓皷鎷虹紒缁㈠幖閹冲繗銇愯缁辨帡鎮╅崘鑼患缂備緡鍠栭悧鎾崇暦閹烘鍊烽悗鐢登圭敮鎯р攽閻樺灚鏆╁┑顔碱嚟閹广垹螣娓氼垳鈧埖銇勯弴妤€浜鹃梺鍝勭焿缁查箖骞嗛弮鍫晬婵炴垶锕╂导锟�<<
12@gma联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4