资源描述：
《matlabIIR滤波器的设计7》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验七IIR滤波器的设计一.实验内容1、调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st,画出信号的时域波形和幅度频谱特性(以Hz为横坐标单位)。确定三路调幅信号的载波频率和调制信号频率。(这个怎么可能是抑制载波信号？？？)function[st]=msgt(N,Fs);■T=l/Fs;t=O:T:(N-l)*T;k=O:N-l;■fcl=Fs/10;fml=fcl/10;fc2=Fs/20;fm2=fc2/10■■fc3=Fs/40;fm3=fc3/l0;■xtl二cos(2*pi*fml*t).*cos(2*
2、pi*fcl*t)+cos(2*pi*fcl*t);■xt2二cos(2*pi*fm2*t)・*cos(2*pi*fc2*t)+cos(2*pi*fc2*t);■xt3二cos(2*pi*fm3*t)・*cos(2*pi*fc3*t)+cos(2*pi*fc3*t);■st=xtl+xt2+xt3;2、通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0・1dB，阻带最小衰减为60dBo3、编程序调用MATL
3、AB滤波器设计函数ellipord和ellip,分别设计这三个椭圆滤波器，并绘图显示其衰减函数曲线。4、调用滤波器输出函数filter,用三个滤波器分别对信号产生函数mstg产生的信号st进行滤波，分离出st中的三路不同载波频率的调幅信号yl(n)、y2(n)和y3(n),并绘图显示yl(n)>y2(n)和y3(n)的时域波形，观察分离效果。5、信号产生函数mstg中采样点数N=1600,对st进行N点FFT可以得到6根理想谱线。如果取N=1800,可否得到6根理想谱线？为什么？N=2000呢？画图验证结论。K6、修改信号产生函数mst
4、g,给每路调幅信号加入载波成分，产生调幅（AM）信号，观察AM信号与抑制载波调幅信号的时域波形及其频谱的差别（ppt有错误）编程原理、思路和公式三.程序脚本并注释Msgt子程序function[st]=msgt(N,Fs);T=l/Fs;t=0:T:(N-l)*T;k=0:N-l;fcl=Fs/10;fml=fcl/10;fc2=Fs/20;fm2=fc2/10;fc3=Fs/40;fm3=fc3/10;xtl=cos(2*pi*fml*t)•Pos(2*pi*fcl*t);xt2=cos(2*pi*fm2*t)•*cos(2*pi*f
5、c2*t);xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t);st=xtl+xt2+xt3;function[st]=msgtl(N,Fs);T=l/Fs;t=0:T:(N-l)*T;k=0:N-l;fcl=Fs/10;fml=fcl/10;fc2=Fs/20;fm2=fc2/10;fc3=Fs/40;fm3=fc3/10;xtl=cos(2*pi*fm•*cos(2*pi*fcl*t)+cos(2*pi*fcl壮);xt2=cos(2*pi*fm2*t)•Pos(2*pi*fc2*t)+cos(2*pi*fc2
6、*t);xt3=cos(2*pi*fm3*t)•*cos(2*pi*fc3*t)+cos(2*pi*fc3*t);st=xtl+xt2+xt3;closeall;clearall;N=1600;Fs=10000;T=1/Fs;Tp=N*T;%采样频率Fs=10kHz,Tp为采样时间t=O:T:(N-1)*T;k=O:N-1;f=k/Tp;st=msgt(N,Fs);fxt=fft(st,N);%计算信号st的频谱subplot(2J,l)plot(t,st);grid;xlabel('t/s');ylabel('s(ty);axis([
7、O,Tp/&min(st),max(st)]);title('(a)s(t)的波形')subplot(2,l,2)stem(f,abs(fxt)/max(abs(fxt)),,.');grid;title('(b)s(t)的频谱')axis([0,Fs/5,0,1.2]);xlabel('f/Hz');ylabel('幅度');％三路调幅信号的载波频率分别为250Hz、500Hz、1000Hz。%低通滤波fp=280;fs=450;%fp必须大于250,fs必须小于500wp=2*fp/Fs;ws=2*fs/Fs;rp=O.1;rs=6
8、0;%DF指标(低通滤波器的通、阻带边界频)[N,wpO]=ellipord(wp,ws,rp,rs);%调用ellipord计算椭圆DF阶数N和通带截止频率wp[B,A]=ellip(N,rp,rs,wp
