IIR数字滤波器设计及软件实现[1].docVIP

. .. 实验四 function st=mstg N=800; Fs=10000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; fm1=fc1/10; fc2=Fs/20; fm2=fc2/10; fc3=Fs/40; fm3=fc3/10; xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); st=xt1+xt2+xt3; fxt=fft(st,N); subplot(2,1,1); plot(t,st);grid;xlabel(t(s));ylabel(s(t)); axis([0 Tp/4 min(st) max(st)]);title((a) s(t)的波形); subplot(2,1,2); stem(f,abs(fxt)/max(abs(fxt)),.);grid;title((b) s(t)的频谱); axis([0 Fs/5 0 1.2]); xlabel(f(Hz));ylabel(幅度); function st=mstg1 N=800; Fs=10000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T;k=0