參與評(píng)論

評(píng)論

登錄后參與評(píng)論

全部評(píng)論(1條)

• 

  我昱云南的人民 2018-04-12 18:21:30

  fp = 30; fs = 35; Fs = 100; wp = 2*pi*fp/Fs; ws = 2*pi*fs/Fs; wp = tan(wp/2); ws = tan(ws/2); % 通帶Z大衰減為0.5dB，阻帶Z小衰減為40dB [N， wn] = buttord(wp， ws， 0.5， 40， 's'); % 模擬低通濾波器極零點(diǎn) [z， p， k] = buttap(N); % 由極零點(diǎn)獲得轉(zhuǎn)移函數(shù)參數(shù) [b， a] = zp2tf(z， p， k); % 由原型濾波器獲得實(shí)際低通濾波器 [B， A] = lp2lp(b， a， wp); [bz， az] = bilinear(B， A， .5); [h， w] = freqz(bz， az， 256， Fs); figure plot(w， abs(h)) grid on 圖1 巴特沃斯數(shù)字低通濾波器 1-2基于Butterworth模擬濾波器原型，使用雙線性狀換設(shè)計(jì)數(shù)字濾波器：各參數(shù)值為：通帶截止頻率Omega=0.2*pi，阻帶截止頻率Omega=0.3*pi，通帶波動(dòng)值Rp=1dB，阻帶波動(dòng)值Rs=15dB，設(shè)Fs=4000Hz。 代碼： wp=0.2*pi;ws=0.3*pi; Fs=4000;T=1/Fs; OmegaP=(2/T)*tan(wp/2); OmegaS=(2/T)*tan(ws/2); rp=1;rs=15;as=15; ripple=10^(-rp/20);attn=10^(-rs/20); [n，wn]=buttord(OmegaP，OmegaS，rp，rs，'s'); [z，p，k]=Buttap(n); [b，a]=zp2tf(z，p，k); [bt，at]=lp2lp(b，a，wn); [b，a]=bilinear(bt，at，F(xiàn)s); [db，mag，pha，grd，w]=freqz_m(b，a); % %下面繪出各條曲線 subplot(2，2，1);plot(w/pi，mag);title('Magnitude Frequency幅頻特性'); xlabel('w(/pi)');ylabel('|H(jw)|'); axis([0，1，0，1.1]); set(gca，'XTickMode'，'manual'，'XTick'，[0 0.2 0.3 1]); set(gca，'YTickMode'，'manual'，'YTick'，[0 attn ripple 1]);grid subplot(2，2，2);plot(w/pi，db);title('Magnitude Frequency幅頻特性(db)'); xlabel('w(/pi)');ylabel('dB'); axis([0，1，-30，5]); set(gca，'XTickMode'，'manual'，'XTick'，[0 0.2 0.3 1]); set(gca，'YTickMode'，'manual'，'YTick'，[-60 -as -rp 0]);grid subplot(2，2，3);plot(w/pi，pha/pi);title('Phase Frequency相頻特性'); xlabel('w(/pi)');ylabel('pha(/pi)'); axis([0，1，-1，1]); subplot(2，2，4);plot(w/pi，grd);title('Group Delay群延時(shí)'); xlabel('w(/pi)');ylabel('Sample'); axis([0，1，0，15]); set(gca，'XTickMode'，'manual'，'XTick'，[0 0.2 0.3 1]);grid 運(yùn)行結(jié)果： 圖2巴特沃思數(shù)字低通濾波器幅頻-相頻特性 1-3設(shè)計(jì)一巴特沃斯高通數(shù)字濾波器，要求通帶截止頻率0.6*pi，通帶衰減不大于1dB，阻帶衰減15DB，抽樣T=1。 Wp=0.6*pi; Ws=0.4*pi; Ap=1; As=15; [N，wn]=buttord(Wp/pi，Ws/pi，Ap，As); %計(jì)算巴特沃斯濾波器階次和截止頻率 %頻率變換法設(shè)計(jì)巴特沃斯高通濾波器 [db，mag，pha，grd，w]=freqz_m(b，a); %數(shù)字濾波器響應(yīng) plot(w，mag); title('數(shù)字濾波器幅頻響應(yīng)|H(ej\Omega)|') 圖3巴特沃斯數(shù)字高通濾波器 2-1用窗函數(shù)法設(shè)計(jì)一個(gè)線性相位FIR低通濾波器，并滿足性能指標(biāo):通帶邊界頻率 Wp=0.5*pi，阻帶邊界頻率Ws=0.66*pi，阻帶衰減不小于40dB，通帶波紋不大于3dB。選擇漢寧窗。 代碼： wp =0.5*pi; ws=0.66*pi; wdelta =ws-wp; N= ceil(8*pi/wdelta) if rem(N，2)==0 N=N+1; end );

  贊(7)

  回復(fù)(0)

  評(píng)論

  評(píng)論

  登錄后參與評(píng)論