发布时间 : 星期三 文章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)的波形');
.