非均匀滤波器的设计.doc

N=4000;Fs=8000;T=1/Fs;Tp=N*T;

t=0:T:(N-1)*T;

fc=Fs/10;f0=fc/10;

xt=sin(2*pi*15*t)+sin(2*pi*40*t)+sin(2*pi*100*t)+sin(2*pi*200*t)+sin(2*pi*400*t)+sin(2*pi*800*t)+sin(2*pi*1500*t)+sin(2*pi*3000*t);

fst=fft(xt,N);k=0:N-1;f=k/Tp;

figure(9)

subplot(2,1,1);plot(t,xt);grid;xlabel(t/s);ylabel(x(t));

axis([0