电子科大《数字信号处理(DSP)》第6章数字滤波器.pptVIP

全通滤波器：相位均衡器 零点设置在单位圆外，与单位圆内的极点互为倒数根 全通滤波器：相位均衡器 a序列和b序列正好相反排列； 幅频特性恒定为1；相频特性可以变动。 b=[-0.4 0.6 -0.9 1 ];a=[1 -0.9 0.6 -0.4]; [h,w]=freqz(b,a);subplot(2,1,1),plot(w,abs(h),w,angle(h),:); 全通滤波器：相位均衡器 最小相位系统 所有零点都设置在单位圆内。 稳定因果的可逆系统； 任何稳定因果系统都可由一个最小相位系统和一个稳定因果的全通系统级联构成； 在所有具有相同幅频特性的系统中，最小相位系统具有最小群延迟。 最小相位系统 最小相位系统的设计 对于要求的幅频特性： 将幅频特性平方表达为： 进行代换： 得到： 取 在单位圆内的全部零极点，就可以构成最小相位系统。 二阶带通滤波器的仿真结果 极点幅度与-3dB带宽的关系 对窄带滤波器，可近似表达为： 二阶带通滤波器设计实例 设计带通滤波器分离音频信号： 考虑音频范围： 进行频率范围的归一化 二阶带通滤波器设计实例 设置带宽范围： 决定极点位置： 代入数值得出系统函数： 二阶带通滤波器设计仿真 n=(0:1/40000:3); x=sin(1200*n)-cos(1600*n)+sin(2000*n); subplot(3,1,1),plot(n,x);axis([0.1,0.4,-3.2,3.2]); pause;sound(x,40000);pause; b1=[0.001 0 -0.001];a1=[1 -1.9971 0.998];[h1,w1]=freqz(b1,a1,500000); y1=filter(b1,a1,x);subplot(3,3,4),plot(w1,abs(h1));axis([0,0.1,-0.2,1.2]); subplot(3,3,7),plot(n,y1);axis([0.2,0.25,-1.2,1.2]);pause;sound(y1,40000);pause; b2=[0.001 0 -0.001];a2=[1 -1.9964 0.998];[h2,w2]=freqz(b2,a2,500000); y2=filter(b2,a2,x);subplot(3,3,5),plot(w2,abs(h2));axis([0,0.1,-0.2,1.2]); subplot(3,3,8),plot(n,y2);axis([0.2,0.25,-1.2,1.2]);pause;sound(y2,40000);pause; b3=[0.001 0 -0.001];a3=[1 -1.9955 0.998];[h3,w3]=freqz(b3,a3,500000); y3=filter(b3,a3,x);subplot(3,3,6),plot(w3,abs(h3));axis([0,0.1,-0.2,1.2]); subplot(3,3,9),plot(n,y3);axis([0.2,0.25,-1.2,1.2]);pause;sound(y3,40000); 二阶带通滤波器设计仿真结果 二阶带阻滤波器：阻塞滤波器 将一对共轭零点设置在单位圆特定频率点上，对该频率形成阻塞；极点靠近零点。 二阶带阻滤波器：阻塞滤波器 极点与零点相位相同，决定阻塞频率； 极点的幅度决定带宽； 二阶带阻滤波器的仿真 for i=1:8, a0=0.9;w0=i*pi/10; b=[1 -2*cos(w0) 1]; a=[1 -2*a0*cos(w0) a0*a0]; [h,w]=freqz(b,a);hm=max(abs(h)); subplot(1,2,1),zplane(b,a); subplot(1,2,2),plot(w,abs(h)/hm); pause;end; 二阶带阻滤波器的仿真结果 根据计算机音频信号的特点，设置采样周期： 二阶带阻滤波器设计实例 60Hz交流信号经常形成对微弱音频信号的干扰，利用阻塞滤波器可以有效消除这种干扰，提取时域音频微弱信号。 二阶带阻滤波器设计实例 归一化阻塞频率为 设置阻带宽度 得出极点幅度 二阶带阻滤波器设计仿真 b=[1 -1.999707696 1];a=[1 -1.997707989 0.998001]; [h,w]=freqz(b,a,50000); subplot(2,1,1),plot(w,abs(h));axis([0,0.1,-0.2,1.2]); pause;n=0:1/22050:5; x0=[zeros(1,200) 0.1 0.1 0.1 zeros(1,200) -