fir_demo.m

来自「很多matlab的源代码」· M 代码 · 共 30 行

f1=100;
f2=6000;
fs=16000;
fstop=1000;
npass=fstop/fs;
omg1=2*pi*f1/fs;
omg2=2*pi*f2/fs;
N=256;
n=0:(N-1);
x=sin(n*omg1)+sin(n*omg2);
M=51;
m=0:(M-1);
h=2*npass*sinc(2*npass*(m-(M-1)/2));
y=conv(x,h);
subplot(3,2,1);
plot(n/fs,x);
X=fft(x,N);
subplot(3,2,2);
plot(n/N*fs,abs(X));
subplot(3,2,3);
stem(m,h);
subplot(3,2,4);
H=fft(h,N);
plot(n/N*fs,20*log10(abs(H)));
subplot(3,2,5);
plot((0:(length(y)-1))/fs,y);
Y=fft(y,N);
subplot(3,2,6);
plot(n/N*fs,abs(Y));