對周期方波信號進行濾波matlab的實現(xiàn)

1、對周期方波信號進行濾波1、生成一個基頻為10Hz的周期方波信號；2、設(shè)計一個濾波器，濾去該周期信號中40Hz以后的頻率成分，觀察濾波前后的信號波形和 頻譜。3、如果該信號x淹沒在噪聲中(隨機噪聲用randn(1,N)生成，N為其樣點數(shù))，試濾去噪 聲。%采樣頻率取200Hz,可以改變看看。輸入、輸出信號頻譜同時繪制在一起比照%采用低通濾波器，保存40Hz以下頻率成分clearfs=200;t=0:1/fs:1;x二square(2*pi*10*t);wp=40*2/fs;ws=45*2/fs;Rp=3;Rs=45;Nn=128;N, wn=buttord (wp, ws, Rp, Rs)b,

2、a=butter (N, wn, low，)y=filter(b, a, x);plot (t, x, r- , t, y)%grid onaxis(0 1.2 -2 2)title(，紅色代表原信號，藍色代表只保存40Hz以下頻率成分)H, W=freqz (b, a);k=0:511;plot (fs/2)/512*k, abs(H);grid ontitleC濾波器頻率響應(yīng))T=l/fs;N=4* (fs/10);n=0:N-l;xn=square(2*pi*10*n*T);X=fftshift(fft(xn, 512);xk=l/N*X;Y=fftshift(fft (y,512);y

3、k=l/N*Y;plot (-fs/2+fs/512*k, abs(xk)grid onlegend (方波信號的頻譜)plot (-fs/2+fs/512*k, abs (xk),J b，, -fs/2+fs/512*k, abs (yk), r)grid onlegend (原信號的頻譜，濾波后信號的頻譜)喊哪厥煙%燃96炮溉96蹴%保存40Hz頻率成分溉厥厥%9概%采樣頻率取200Hz,可以改變看看。輸入、輸出信號頻譜同時繪制在一起比照%采用帶通濾波器，只保存40Hz頻率成分clearfs=200;t=0:1/fs:1;x二square(2*pi*10*t);wp=35 45*2/fs;

4、ws=30 50*2/fs;Rp=3;Rs=45;Nn=128;N, wn=buttord (wp, ws, Rp, Rs)b, a=butter (N, wn, bandpass，) %只保存 40Hz 頻率成分y=filter (b, a, x);plot (t, x,，r-，t, y)grid onaxis(0 1.2 -1.2 1.2)titleC紅色代表原信號，藍色代表只保存40Hz頻率成分)H, W=freqz (b, a);k=0:511;plot (fs/2)/512*k, abs(H);grid ontitleC濾波器頻率響應(yīng))T=1/fs;N=4*(fs/10);n=0:N-l;xn二square (2*pi*10*n*T);X=fftshift(fft (xn,512);xk=l/N*X;Y=fftshift(fft(y,512);yk=l/N*Y;plot (-fs/2+fs/512*k, ab