对周期方波信号进行滤波matlab的实现

1、对周期方波信号进行滤波1、生成一个基频为10Hz的周期方波信号；2、设计一个滤波器，滤去该周期信号中40Hz以后的频率成分，观察滤波前后的信号波形和 频谱。3、如果该信号x淹没在噪声中(随机噪声用randn(1,N)生成，N为其样点数)，试滤去噪 声。%采样频率取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滤波器频率响应)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滤波器频率响应)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