IIR数字滤波器的设计及软件实现

.

subplot(3,1,2); tplot(y1t,T,yt); subplot(3,1,3);

stem(f,abs(fyt)/max(abs(fyt)),'.'); grid;

title('(c) s(t)的频谱'); axis([0,Fs/5,0,1.2]);

xlabel('f/Hz');ylabel('幅度'); clc;clear

%子程序

%产生信号程序 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;%第一路调幅信号的载波频率fc1=1000HZ

fm1=fc1/10;%第一路调幅信号的调制信号频率为fm1=100hz. fc2=FS/20;%第二路调幅信号的载波频率fc2=500HZ

fm2=fc2/10;%第二路调幅信号的调制信号频率为fm2=50hz. fc3=FS/40;%第三路调幅信号的载波频率fc3=250HZ

fm3=fc3/10;%第三路调幅信号的调制信号频率为fm3=25hz. 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);%计算信号st的频谱.%绘图 subplot(2,1,1) plot(t,st); grid on;

xlabel('t/s');ylabel('s(t)');

axis([0,TP/8,min(st),max(st)]); title('(a)s(t)的波形') subplot(2,1,2)

stem(f,abs(fxt)/max(abs(fxt)),'.'); grid on;title('(b) s(t)的频谱') axis([0,FS/5,0,1.2]);

xlabel('f/HZ');ylabel('幅度');

%损耗输出波形

function myplot(B,A) [H,W]=freqz(B,A,1000); m=abs(H);

.

.

plot(W/pi,20*log10(m/max(m))); grid on;

xlabel('\\omega/\\pi'); ylabel('幅度（dB）'); axis([0,1,-80,5]);

title('(a) s(t)损耗函数曲线');

%滤波器输出波形 function tplot(xn,T,yn) n=0:length(xn)-1; t=n*T; plot(t,xn); grid on; xlabel('t/s'); ylabel('y(n)');

axis([0,t(end),min(xn),1.2*max(xn)]) title('(b) s(t)的波形');

.