帶通濾波器濾波

2023-01-21 21:39:05 字數 1471 閱讀 9021

物件1:乙個含有10hz、30hz和60hz的混合正弦波訊號x=sin(2*pi*t*10)+ sin(2*pi*t*30)+ sin(2*pi*t*60),通過乙個設計的iir帶通橢圓濾波器(elliptic)後,使輸入的混合訊號中頻率為30hz的正弦波訊號通過,而將頻率為10hz和60hz的正弦波訊號大大衰減。

產生乙個含有10hz、30hz和60hz的混合正弦波訊號x

matlab程式實現如下:

fs=200;

t=(1:200)/fs;

x1=sin(2*pi*t*10);

x2=sin(2*pi*t*30);

x3=sin(2*pi*t*60);

x= x1+ x2+ x3;

plot(t,x);

title('混合正弦波訊號x(t)-濾波前');

grid;

產生乙個4階iir帶通橢圓濾波器

matlab程式實現如下:

wp1=20;

wp2=40;

fs=200;

rp=0.2;

rs=40;

wp1=2*pi*wp1;

wp2=2*pi*wp2;

bw=wp2-wp1;

wo=sqrt(wp2*wp1);

[z,p,k]=ellipap(4,rp,rs建立模擬低通濾波器原型

[a,b,c,d]=zp2ss(z,p,k); %把模擬低通濾波器原型轉換成模擬低通濾波器

[at,bt,ct,dt]=lp2bp(a,b,c,d,wo,bw); %把模擬低通濾波器轉換成帶通濾波器

[at1,bt1,ct1,dt1]=bilinear(at,bt,ct,dt,fs); %雙線性變換法把模擬濾波

[num,den]=ss2tf(at1,bt1,ct1,dt1器轉換成數字濾波器

[h,w]=freqz(num,den);

plot(w*fs/(2*pi),abs(h));

title('頻率響應特性');

grid;

xlabel('頻率/hz');

ylabel('幅值');

對混合正弦波訊號x進行濾波

matlab程式實現如下

y=filter(num,den,x);

plot(t,y);

title('濾波後訊號-y');

gri;

axis([0 1 -1 1]);

繪出訊號濾波前、後的幅頻圖

matlab程式實現如下

s=fft(x,512);

sf=fft(y,512);

w=(0:255)/256*(fs/2);

plot(w,abs(s(1:256)),'-',w,abs(sf(1:256)),'k:*');

title('訊號濾波前、後的幅頻圖');

xlabel('頻率/hz');

ylabel('幅值');

legend('濾波前的幅頻', '濾波後的幅頻');

濾波器例項

根據語音頻號選擇合適的濾波器,程式1 用matlab對原始語音頻號進行分析,畫出它的時域波形和頻譜 fs 22050 語音頻號取樣頻率為22050 x1 w read h 課程設計2 讀取語音頻號的資料,賦給變數x1 sound x1,22050 語音頻號 y1 fft x1,1024 對訊號做10...

無源濾波器和有源濾波器的區別

工作原理無源濾波器由lc等被動元件組成,將其設計為某頻率下極低阻抗,對相應頻率諧波電流進行分流,其行為模式為提供被動式諧波電流旁路通道 而有源濾波器由電力電子元件和dsp等構成的電能變換裝置,檢測負載諧波電流並主動提供對應的補償電流,補償後的源電流幾乎為純正弦波,其行為模式為主動式電流源輸出。諧波處...

FTTX光濾波器

光濾波器簡介 光濾波器 optical filter 是用來進行波長選擇的儀器,它可以從眾多的波長中挑選出所需的波長,而除此波長以外的光將會被拒絕通過。廣泛應用於dwdm系統內動態波長選擇,dwdm訊號分離,光效能監測,可調諧光雜訊過濾及光放大器雜訊抑制等領域。光復用器 耦合器 是把不同波長的光復用...