mat_lvbo.m

来自「使用DSP和MATLAB设计的一个低通FIR滤波器」· M 代码 · 共 27 行

clear;
t=0:0.0001:0.01;
sig=cos(2*pi*1000*t)+cos(2*pi*2500*t)/2;
subplot(2,2,1);
plot(t,sig);
title('原始信号(时域)');

y=fft(sig,512);
py=y.*conj(y)/512;
f=10000*(0:155)/512;
subplot(2,2,2);
plot(f,py(1:156));
title('原始信号(频域)');
　　
c=kaiserord([1000 1500],[1 0],[0.01 0.1],8000,'cell')
[b,a]=fir1(c{:})
newsig=filter(b,a,sig);
subplot(2,2,3);
plot(t,newsig);
title('输出信号(时域)');

y=fft(newsig,512);
py=y.*conj(y)/512;
f=10000*(0:155)/512;
subplot(2,2,4);
plot(f,py(1:156))
title('输出信号(频域)');