fftfilter.m

来自「在matlab中实现数字信号处理的一些源程序」· M 代码 · 共 10 行

t=linspace(0,2*pi,2^8);
y=exp(-(cos(t).^2)).*(sin(2*t)+2*cos(4*t)+0.4*sin(t).*sin(50*t));
figure(1);
plot(t,y);
fy=fft(y);
filterfy=[fy(1:6) zeros(1,2^8-11) fy(2^8-4:2^8)];
filtery=ifft(filterfy);
figure(2);
subplot 211; plot(t,real(filtery));
subplot 212; plot(t,imag(filtery));