鍩轰簬MATLAB鐨勬暟瀛楁护娉㈠櫒鐨勮�捐�＄▼搴?- 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

IIR 

浣庨�氭护娉㈠櫒鐨勮�捐�＄▼搴忎负锛?/p>

 

 

Ft=8000;

 

 

Fp=1000;

 

 

Fs=1200;

 

 

As=100 ;

 

 

Ap=1;

 

 

wp=2*pi*Fp/Ft;

 

 

ws=2*pi*Fs/Ft;

 

 

fp=2*Fp*tan(wp/2);

 

 

fs=2*Fs*tan(ws/2);

 

 

[n11,wn11]=buttord(wp,ws,1,50,'s');

 

 

[b11,a11]=butter(n11,wn11,'s');

 

 

[num11,den11]=bilinear(b11,a11,0.5);

 

 

[h,w]=freqz(num11,den11);

 

 

axes(handles.axes1);

 

 

plot(w*8000*0.5/pi,abs(h));

 

 

xlabel('Frequency/Hz');

 

 

ylabel('Magnitude');

 

 

title('

宸寸壒娌冩柉鏁板瓧浣庨�氭护娉㈠櫒

'); 

 

宸寸壒娌冩柉甯﹂�氭护娉㈠櫒璁捐�＄▼搴忎�?/p>

: 

Ft=8000;

 

 

Fp1=1200;

 

 

Fp2=3000;

 

 

Fs1=1000;

 

 

Fs2=3200;

 

 

As=100;

 

 

Ap=1;

 

 

wp1=tan(pi*Fp1/Ft);

 

 

wp2=tan(pi*Fp2/Ft);

 

 

ws1=tan(pi*Fs1/Ft);

 

 

ws2=tan(pi*Fs2/Ft);

 

 

w=wp1*wp2/ws2;