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

对周期方波信号进行滤波

) 1、生成一个基频为10Hz的周期方波信号；

) 2、设计一个滤波器，滤去该周期信号中40Hz以后的频率成分，观察滤波前后的信号波形和频谱。

)3、如果该信号x(t)淹没在噪声中(随机噪声用randn(1,N)生成，N为其样点数)，试滤去噪声。

%采样频率取200Hz，可以改变看看。输入、输出信号频谱同时绘制在一起对比

%采用低通滤波器，保留40Hz以下频率成分

clear

fs=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 on

axis([0 1.2 -2 2])

title('红色代表原信号，蓝色代表只保留40Hz以下频率成分')

figure(2)

[H,W]=freqz(b,a);

k=0:511;

plot((fs/2)/512*k,abs(H));

grid on

title('滤波器频率响应')

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;