利用matlab实现数字滤波的一个小实验
关键词: matlab 数字滤波
%程序利用低通滤波器实现将信号y=2*sin(10*pi*t)+5*cos(40*pi*t)变成y=2*sin(10*pi*t)信号
N=128; %
fs=80; %信号的采用频率
n=0:N-1;
x=2*sin(10*pi*n/fs)+5*cos(40*pi*n/fs);
subplot(2,2,1)
plot(n,x);
y=fft(x,N);
q=(0:N/2-1)*fs/N;
mx=abs(y(1:N/2));
subplot(2,2,2)
plot(q,mx);
[b,a]=cheby1(12,0.5,100/500); %设计一个切必雪夫滤波器
y=filter(b,a,x); %滤波处理
subplot(2,2,3)
plot(n,y);
yy=fft(y,N);
mx2=abs(yy(1:N/2));
subplot(2,2,4)
plot(q,mx2);
实验结果
说明:图1:未滤波前信号波形
图2:未滤波前信号频谱
图3 :滤波后信号波形
图4 :滤波后信号频谱