方波濾波matlab的實(shí)現(xiàn)

1、a)對(duì)周期方波信號(hào)進(jìn)行濾波（ 1）生成一個(gè)基頻為10Hz 的周期方波信號(hào)；（ 2）設(shè)計(jì)一個(gè)濾波器，濾去該周期信號(hào)中40Hz 以后的頻率成分，觀察濾波前后的信號(hào)波形和頻譜。%采樣頻率取200Hz ，可以改變看看。輸入、輸出信號(hào)頻譜同時(shí)繪制在一起對(duì)比采用帶通濾波器，只保留40Hz 頻率成分clearfs=200;t=0:1/fs:1;x=square(2*pi*10*t);wp=35 45*2/fs;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=

2、filter(b,a,x);figure(1)plot(t,x,r-,t,y)grid onaxis(0 1.2 -1.2 1.2)title( 紅色代表原信號(hào)，藍(lán)色代表只保留40Hz 頻率成分 )figure(2)H,W=freqz(b,a);k=0:511;plot(fs/2)/512*k,abs(H);grid ontitle( 濾波器頻率響應(yīng))T=1/fs;N=4*(fs/10);n=0:N-1;xn=square(2*pi*10*n*T);X=fftshift(fft(xn,512);xk=1/N*X;Y=fftshift(fft(y,512);yk=1/N*Y;figure(3)p

3、lot(-fs/2+fs/512*k,abs(xk),b,-fs/2+fs/512*k,abs(yk),r)grid onlegend(原信號(hào)的頻譜 ,濾波后信號(hào)的頻譜)紅 色 代 表 原 信 號(hào) ， 藍(lán) 色 代 表 只 保 留 40Hz 頻 率 成 分10.80.60.40.20-0.2-0.4-0.6-0.8-100.20.40.60.811.2濾波器頻率響應(yīng)1.41.210.80.60.40.2001020304050607080901000.7原信號(hào)的頻譜濾波后信號(hào)的頻譜0.60.50.40.30.20.10-80-60-40-20020406080100-100%采樣頻率取200Hz

4、 ，可以改變看看。輸入、輸出信號(hào)頻譜同時(shí)繪制在一起對(duì)比%采用低通濾波器，保留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,a=butter(N,wn,low)y=filter(b,a,x);figure(1)plot(t,x,r-,t,y)%grid onaxis(0 1.2 -2 2)title( 紅色代表原信號(hào)，藍(lán)色代表只保留40Hz 以下頻率成分 )figure(2)H,W=freqz(b,a)

5、;k=0:511;plot(fs/2)/512*k,abs(H);grid ontitle( 濾波器頻率響應(yīng))T=1/fs;N=4*(fs/10);n=0:N-1;xn=square(2*pi*10*n*T);X=fftshift(fft(xn,512);xk=1/N*X;Y=fftshift(fft(y,512);yk=1/N*Y;figure(3)plot(-fs/2+fs/512*k,abs(xk)grid onlegend(方波信號(hào)的頻譜)figure(4)plot(-fs/2+fs/512*k,abs(xk),b,-fs/2+fs/512*k,abs(yk),r)grid onlegend(原信號(hào)的頻譜 ,濾波后信號(hào)的頻譜)紅 色 代 表 原 信 號(hào) ， 藍(lán) 色 代 表 只 保 留 40Hz 以 下 頻 率 成 分21.510.50-0.5-1-1.5-200.20.40.60.811.2濾波器頻率響應(yīng)1.41.210.80.60.40.2001020304050607080901000.7方波信號(hào)的頻譜0.60.50.4