matlab IIR滤波器的设计7.docVIP

实验七IIR滤波器的设计 实验内容 1、调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，画出信号的时域波形和幅度频谱特性（以Hz为横坐标单位）。确定三路调幅信号的载波频率和调制信号频率。 (这个怎么可能是抑制载波信号？？？) function [st]=msgt(N,Fs); T=1/Fs;t=0:T:(N-1)*T;k=0:N-1; 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)+cos(2*pi*fc1*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=xt1+xt2+xt3; 2、通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1 dB, 阻带最小衰减为60 dB。 3、编程序调用MATLAB滤波器设计函数ellipord和ellip，分别设计这三个椭圆滤波器，并绘图显示其衰减函数曲线。 4、调用滤波器输出函数filter，用三个滤波器分别对信号产生函数mstg产生的信号st进行滤波，分离出st中的三路不同载波频率的调幅信号y1(n)、y2(n)和y3(n)，并绘图显示y1(n)、y2(n)和y3(n)的时域波形，观察分离效果。 5、 信号产生函数mstg中采样点数N=1600，对st进行N点FFT可以得到6根理想谱线。如果取N=1800，可否得到6根理想谱线？为什么？N=2000呢？画图验证结论。 6、修改信号产生函数mstg，给每路调幅信号加入载波成分，产生调幅（AM）信号，观察AM信号与抑制载波调幅信号的时域波形及其频谱的差别 (ppt有错误） 编程原理、思路和公式 程序脚本并注释 Msgt 子程序 function [st]=msgt(N,Fs); T=1/Fs;t=0:T:(N-1)*T;k=0:N-1; 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; function [st]=msgt1(N,Fs); T=1/Fs;t=0:T:(N-1)*T;k=0:N-1; 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)+cos(2*pi*fc1*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=xt1+xt2+xt3; close all; clear all; N=1600; Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; st=msgt(N,Fs); fxt=fft(st,N); %计算信号st的频谱 subplot(2,1,1) plot(t,st);grid;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;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必须小于500 wp=2*fp/Fs;ws=2*fs/Fs; rp=0.1;rs=60;