数字信号处理 第10章.pptVIP

　　2． 实验原理 　　设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是： ① 将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； ② 设计过渡模拟滤波器；③ 将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第6章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2以及椭圆模拟与数字滤波器。本实验要求读者调用如上函数直接设计IIR数字滤波器。 　　本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数 filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。 　　3. 实验内容及步骤 　　（1） 调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，该函数还会自动绘图显示st的时域波形和幅频特性曲线，如图10.4.1所示。由图可见，三路信号时域混叠无法在时域分离。但频域是分离的，所以可以通过滤波器的方法在频域分离，这就是本实验的目的。 图10.4.1 三路调幅信号st(即s(t))的时域波形和幅频特性曲线 　　（2） 要求将st中三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1 dB, 阻带最小衰减为60 dB。 　　提示： 抑制载波单频调幅信号的数学表示式为 其中，cos(2πfct)称为载波，fc为载波频率，cos(2πf0t)称为单频调制信号，f0为调制正弦波信号频率，且满足fcf0。由上式可见，所谓抑制载波单频调幅信号，就是两个正弦信号相乘，它有2个频率成分： 和频fc+f0、差频fc－f0，这两个频率成分关于载波频率fc对称。所以，1路抑制载波单频调幅信号的频谱图是关于载波频率fc对称的两根谱线。容易看出，图10.4.1中三路调幅信号的载波频率分别为250 Hz、500 Hz、1000 Hz。有关调幅(AM)和抑制载波调幅(SCAM)的一般原理与概念，请参考通信原理教材。 　　（3） 编程序调用MATLAB滤波器设计函数ellipord和ellip分别设计这三个椭圆滤波器，并绘图显示其损耗函数曲线。  　　（4） 调用滤波器实现函数filter，用三个滤波器分别对信号产生函数mstg产生的信号st进行滤波，分离出st中的三路不同载波频率的调幅信号y1(n)、y2(n)和y3(n)， 并绘图显示y1(n)、y2(n)和y3(n)的时域波形，观察分离效果。 　　4． 信号产生函数mstg清单 　　function st=mstg 　　%产生信号序列向量st, 并显示st的时域波形和频谱 　　%st=mstg 返回三路调幅信号相加形成的混合信号，长度N=800 　　N=800 %N为信号st的长度 　　Fs=10000;T=1/Fs;Tp=N*T; 　%采样频率Fs=10 kHz，Tp为采样时间 　　t=0:T:(N－1)*T;k=0:N－1;f=k/Tp; 　　fc1=Fs/10;　%第1路调幅信号的载波频率fc1=1000 Hz 　　fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100 Hz 　　fc2=Fs/20;　 %第2路调幅信号的载波频率fc2=500 Hz 　　fm2=fc2/10;　 %第2路调幅信号的调制信号频率fm2=50 Hz 　　fc3=Fs/40;　 %第3路调幅信号的载波频率fc3=250 Hz  　　fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25 Hz 　　xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t);%产生第1路调幅信号 　　xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号 　　xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号 　　st=xt1+xt2+xt3;　%三路调幅信号相加 　　fxt=fft(st, N); 　　%计算信号st的频谱 　　%= = = =以下为绘图部分，绘制st的时域波形和幅频特性曲线= = = = = = 　　subplot(3, 1, 1) 　　plot(t, st);grid;xlabel(t/s);ylabel(s(t)); 　　axis(［