[工学]数字信号处理第10章-45.pptVIP

　　2． 实验原理　　设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法）， 应用最广泛的是双线性变换法。 其基本设计过程是： ① 先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； ② 设计过渡模拟滤波器； ③ 将过渡模拟滤波器的系统函数转换成数字滤波器的系统函数。 MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。 　　3. 实验内容及步骤　　（1） 调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号s(t)，该函数还会自动绘图显示s(t)的时域波形和幅频特性曲线， 如图10.4.1所示。 由图可见， 三路信号时域混叠无法在时域分离。 但频域是分离的， 所以可以通过滤波的方法在频域分离，这就是本实验的目的。 　　stem(f， abs(fxt)/max(abs(fxt))， ′.′)； 　　grid； title(′(b) s(t)的频谱′)　　axis(［0， Fs/5， 0， 1.2］)； 　　xlabel(′f/Hz′)； ylabel(′幅度′)　　5． 实验程序框图　　实验程序框图如图10.4.2所示， 供读者参考。 　　6． 思考题　　（1） 请阅读信号产生函数mstg， 确定三路调幅信号的载波频率和调制信号频率。 　　（2） 信号产生函数mstg中采样点数N=800， 对st进行N点FFT可以得到6根理想谱线。 如果取N=1000， 可否得到6根理想谱线？为什么？N=2000呢？请改变函数mstg中采样点数N的值， 观察频谱图， 验证您的判断是否正确。 　　（3） 修改信号产生函数mstg， 给每路调幅信号加入载波成分， 产生调幅（AM）信号，重复本实验， 观察AM信号与抑制载波调幅信号的时域波形及其频谱的差别。 　　提示： AM信号表示式： 　　　　　s(t)=［1+cos(2πf0t)］cos(2πfct) 10.4.2 滤波器参数及实验程序清单　　1. 滤波器参数选取　　观察图10.4.1可知， 三路调幅信号的载波频率分别为250 Hz、 500 Hz、 1000 Hz。 带宽（也可以由信号产生函数mstg清单看出）分别为50 Hz、 100 Hz、 200 Hz。 所以， 分离混合信号st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、 带通滤波器、 高通滤波器）的指标参数选取如下：  对载波频率为250 Hz的调幅信号， 可以用低通滤波器分离，其指标为　　通带截止频率fp=280 Hz， 通带最大衰减ap=0.1 dB；　　阻带截止频率fs=450 Hz， 阻带最小衰减as=60 dB。 椭圆数字滤波器设计 MATLAB语法： [b,a]=ellip(n,Rp,Rs,wn) (5) [b,a]=ellip(n, Rp,Rs, wn,’ftype’) (6) [z,p,k]=ellip(…) (7) [a,b,c,d]=ellip(…) (8) 语句(5)~ (8)都可以设计一个阶数为n，截止频率为wn,通带波纹最大衰减为Rp、阻带波纹最小衰减为Rs的 数字低通滤波器。对于语句(6)，若参数ftype取‘high’或‘stop’，则可以来设计相应的高通或带阻滤波器。 返回值a和b为系统函数的分子和分母的系数。返回 10.4.4 简答思考题　　思考题（1）： 已经在10.4.2节解答。 思考题（3）很简单， 请读者按照该题的提示修改程序， 运行观察。 　　思考题（2）： 因为信号s(t)是周期序列， 谱分析时要求观察时间为整数倍周期。 所以，本题的一般解答方法是， 先确定信号s(t)的周期，再判断所给采样点数N对应的观察时间Tp=NT是否为s(t)的整数个周期。 但信号产生函数mstg产生的信号s(t)共有6个频率成分， 求其周期比较麻烦， 故采用下面的方法解答。 10.5 实验五： FIR数字滤波器设计 与软件实现10.5.1 实验指导　　1． 实验目的　　（1） 掌握用窗函数法设计FIR数字滤波器的原理