matlab第10章MATLAB在信号.ppt

【例10.43】设计一个带阻数字滤波器，要求 利用椭圆滤波器原型。 解： %例10.43的程序：eg10_43.m %设计IIR数字带阻滤波器 %数字带阻滤波器的技术指标 wl=0.19*pi;wsl=0.198*pi;ws2=0.202*pi;wu=0.21*pi; Rp=3;Rs=13; %设置周期为1 T=1; * %用双线性变换法转换为模拟带阻滤波器的指标 Omegau=2/T*tan(wu/2);Omegal=2/T*tan(wl/2); Omegas1=2/T*tan(ws1/2);Omegas2=2/T*tan(ws2/2); Omega0=sqrt(Omegau.*Omegal); B=Omegau-Omegal; %转换为归一化模拟低通的指标 yitas2=Omegas2/B;yita0=Omega0/B; lanmdas=yitas2/(yitas2.^2-yita0.^2);lanmdap=1; %设计归一化模拟低通滤波器 [N,Wn]=ellipord(lanmdap,lanmdas,Rp,Rs,s) [z0,p0,k0]=ellipap(N,Rp,Rs); p=p0*Wn;z=z0*Wn; k=k0*Wn^N; ba=k*real(poly(z)); aa=real(poly(p)); * %转换为模拟带阻 w0=(Omegau+Omegal)/2; %模拟带阻的中心频率 [bt,at] = lp2bs (ba,aa,w0,B) ; %双线性不变法转换为数字带阻 [bd,ad]=bilinear(bt,at,1/T); %求数字带阻的频率响应 [h,w]=freqz (bd,ad); %画幅度响应图 plot(w/pi,20*log10(abs(h)/max(abs(h)))); axis([0.15,0.25,-20,0]); grid; title(数字带阻滤波器的幅度响应); xlabel(w（pi rad));ylabel(分贝数); * 所设计滤波器幅度特性如上图所示，由图可见，设计结果满足指标要求。 * 六、应用举例 【例10.44】已知由三路抑制载波调幅信号相加构成的复合信号st，设计合适的滤波器，分离出这三路调幅信号。要求滤波器的通带最大衰减为0.1dB，阻带最小衰减为60dB。 * （1）由s(t)的时域波形可见在时域根本无法分离这三个正弦分量的信号，我们需要画出信号的频谱，从频域观察3个分量的特征。 这三路信号分布在不同的频带内，可以完全分离开 只要我们设计出合适的滤波器，就可以将这三路信 号在频域分离。 * （2）设计一个低通滤波器，取出第1路调幅信号。由频谱图可见，该低通滤波器的通带截止频率可取为 ，阻带截止频率可取为 ，按照题目要求 取为0.1dB， 取为60dB。为了使所设计的滤波器具有较低的阶数，我们选择椭圆滤波器。 * （3）调用滤波器实现函数filter，用第（2）步所设计的低通滤波器对复合信号st进行滤波，分离出第1路调幅信号。 * （4）要分离出第2、3路信号，只要设计出合适的带通和高通滤波器即可。其中带通滤波器的通带截止频率为fpl=400，fpu=600，阻带截止频率为fsl=350，fsu=650。 * 分离出第2路信号的时域波形 分离出第2路信号的频谱 * 为分离出第3路调幅信号所设计高通滤波器的指 标为fp=800，fs=700。 * 分离出第3路信号的时域波形 分离出第3路信号的频谱 * 七、IIR和FIR数字滤波器的比较 （1）在相同的技术指标下，IIR滤波器由于存在输入对输出的反馈，所以可用比FIR滤波器较少的阶数来满足指标的要求，所用的存储单元少，运算次数少，较为经济。一般在相同的指标下，FIR滤波器的阶数比IIR滤波器的阶数高5~10倍左右。 （2）FIR滤波器可得到严格的线性相位，而IIR滤波器的选择性越好，其相位的非线性越严重。要得到具有线性相位的IIR滤波器需要通过全通网络进行相位校正，这会增加滤波器的阶数。从这一点上，FIR滤波器优于IIR滤波器。 * （3）FIR滤波器主要采用非递归结构，系统总是稳定的，有限精度运算的误差也较小；IIR滤波器必须采用递归结构，极点在z平面的单位圆内时系统才稳定，而运算中的四舍五入处理有时会引起寄生振荡。 （4）FIR滤波器的脉冲响应为有限长，因而可以用快速傅立叶变换算法提高运算速度；IIR滤波器则不能。 * （5）从设计上看，IIR滤波器可以利用模拟滤波器设计的