1、并联系统的系统函数为级联系统的系统函数为宁可瑞滤波器(Linkwitz-Riley),由两个巴特沃斯滤波器级联而成。N阶巴特沃夫滤波器等效宁可瑞滤波器的设计为了使设计的IIR滤波器方便在DSP上实现,常将滤波器转换为二阶节级联的形式。设计好分频器后,为验证分频后的信号是否正确,可用白噪声信号作为输入信号,然后对分频后的信号进行频谱分析。三、仪器设备1.硬件:计算机一台,耳机。2.软件:MATLAB R2010b 四、实验步骤任意选择两种类型的IIR数字滤波器,设计一个二分频的数字分频器,已知系统的采样率为48000Hz。(1)分频点为2000Hz;(2)要求给出类似图8.3的幅频特性图,分频器
2、的幅频响应平坦,在分频点处最多不能超过3dB的偏差;(3)滤波器必须是二阶节形式;(4)给出相位特性图;(5)用频谱分析的方法验证设计好的分频器;(6)对选用的两种类型的滤波器效果进行对比。滤波器设计的基本步骤:五、数据记录我选择要设计的合成滤波器为ButterWorth IIR滤波器和Linkwitz-Riley IIR滤波器。1.设计程序设计程序如下:(以4阶巴特沃斯滤波器、宁可瑞滤波器设计的分频器程序为例(分频器阶数为8阶)%设计分频器clear;clc;fs = 48000;%采样频率为48000Hzfc = 2000;%分频点为2000Hzwc = 2 * fc / fs;N = 4
3、; %滤波器阶数,分频器阶数为2*NBL,AL = butter(N,wc); %计算巴特沃思低通滤波器系统函数B,A系数BH,AH = butter(N,wc,high); %计算巴特沃思高通滤波器系统函数B,A系数magHH,w=freqz(BH,AH);%计算巴特沃思高通滤波器幅频特性magHH=20*log10(abs(magHH);f=w*fs/(2*pi);%把数字频率w转换为模拟频率fBL1,AL1 = butter(N/2,wc);BH1,AH1 = butter(N/2,wc, BL1=conv(BL1,BL1); %计算宁可瑞低通滤波器系统函数B,A系数AL1 = conv
4、(AL1,AL1);BH1=conv(BH1,BH1);%计算宁可瑞高通滤波器系统函数B,A系数AH1 = conv(AH1,AH1);magHH1,w1=freqz(BH1,AH1);%计算宁可瑞高通滤波器幅频特性magHH1=20*log10(abs(magHH1);f1=w1*fs/(2*pi);semilogx(f,magHH,-.r,f1,magHH1,bhold on;magHL,w=freqz(BL,AL);%计算巴特沃思低通滤波器幅频特性magHL=20*log10(abs(magHL);magHL1,w1=freqz(BL1,AL1);%计算宁可瑞低通滤波器幅频特性magHL
5、1=20*log10(abs(magHL1);semilogx(f,magHL,f1,magHL1,B=conv(BL,AH)+conv(BH,AL); %计算巴特沃思滤波器并联系统的系统函数A=conv(AL,AH);magH,w=freqz(B,A); %计算巴特沃思滤波器并联系统幅频特性magH=20*log10(abs(magH);B1=conv(BL1,AH1)+conv(BH1,AL1); %计算宁可瑞滤波器并联系统的系统函数A1=conv(AL1,AH1);magH1,w1=freqz(B1,A1); %计算宁可瑞滤波器并联系统幅频特性magH1=20*log10(abs(mag
6、H1);semilogx(f,magH,f1,magH1,legend(巴特沃斯滤波器,宁可瑞滤波器title(IIR分频器的幅度特性axis(100 20000 -40 10);hold ongrid on%分析巴特沃斯滤波器及其设计的分频器的幅频特性、零极点分布%巴特沃斯低通subplot(2,2,1);zplane(BL,AL);巴特沃斯低通滤波器的零极点分布)HL,wL=freqz(BL,AL);subplot(2,2,2);plot(wL/pi,abs(HL);巴特沃斯低通滤波器的幅度特性xlabel(omega/piylabel(|H(ejomega)|subplot(2,2,4)
7、;plot(wL/pi,angle(HL);phi(omega)巴特沃斯低通滤波器的相频特性%巴特沃斯高通zplane(BH,AH);巴特沃斯高通滤波器的零极点分布HH,wH=freqz(BH,AH);plot(wH/pi,abs(HH);巴特沃斯高通滤波器的幅度特性plot(wH/pi,angle(HH);巴特沃斯高通滤波器的相频特性%设计的分频器zplane(B,A);分频器的零极点分布H,w=freqz(B,A);plot(w/pi,abs(H);分频器的幅度特性plot(w/pi,angle(H);分频器的相频特性%分析宁可瑞滤波器及其设计的分频器的幅频特性、零极点分布%宁可瑞低通zp
8、lane(BL1,AL1);宁可瑞低通滤波器的零极点分布HL1,wL1=freqz(BL1,AL1);plot(wL1/pi,abs(HL1);宁可瑞低通滤波器的幅度特性plot(wL1/pi,angle(HL);宁可瑞低通滤波器的相频特性%宁可瑞高通zplane(BH1,AH1);宁可瑞高通滤波器的零极点分布HH1,wH1=freqz(BH1,AH1);plot(wH1/pi,abs(HH1);宁可瑞高通滤波器的幅度特性plot(wH1/pi,angle(HH1);宁可瑞高通滤波器的相频特性zplane(B1,A1);H1,w1=freqz(B1,A1);plot(w1/pi,abs(H1)
9、;plot(w1/pi,angle(H);%滤波效果验证%巴特沃斯设计的分频器滤波效果hB,g=tf2sos(B,A) %调用函数tf2sos,将巴特沃斯滤波器设计的分频器转换成二阶节形式xB,fs,bits=wavread(E:white.wavX=fft(xB,1024);for i=1:size(hB)xB=filter(hB(i,1:3),hB(i,4:6),xB);%二阶节级联形式对白噪声进行滤波处理endwavwrite(xB,fs,bits,e:巴特沃斯设计的分频器滤波后信号.wav)%将滤波后的噪声保存 YB=fft(xB,1024);k=0:1023;N=1024;wk=2*
10、k/N;subplot(211);plot(wk,abs(X); title(原始白噪声信号频谱subplot(212);plot(wk,abs(YB);巴特沃斯设计的分频器滤波后信号频谱hL,gL=tf2sos(B1,A1) %调用函数tf2sos,将宁可瑞滤波器设计的分频器转换成二阶节形式xL,fs,bits=wavread(X=fft(xL,1024);size(hL)xL=filter(hL(i,1:3),hL(i,4:6),xL); %二阶节级联形式对白噪声进行滤波处理wavwrite(xL,fs,bits,宁可瑞设计的分频器滤波后信号.wavYL=fft(xL,1024);plot
11、(wk,abs(YL);宁可瑞设计的分频器滤波后信号频谱2.图形记录figure1两种滤波器设计的分频器的幅度特性曲线:figure2巴特沃思低通滤波器的零极点分布和幅频特性:figure3巴特沃思高通滤波器的零极点分布和幅频特性:figure4巴特沃思滤波器设计的分频器的零极点分布和幅频特性:figure5宁可瑞低通滤波器的零极点分布和幅频特性:figure6宁可瑞高通滤波器的零极点分布和幅频特性:figure7宁可瑞滤波器设计的分频器的零极点分布和幅频特性:figure8巴特沃思滤波器设计的分频器的滤波效果的频谱分析:figure9宁可瑞滤波器设计的分频器的滤波效果的频谱分析:figure
12、10将白噪声音频文件通过分频器处理后保存为wav文件:3.数据记录hB,g=tf2sos(B,A) %将巴特沃斯滤波器设计的分频器转换成二阶节形式得到数据:hL,gL=tf2sos(B1,A1) %将宁可瑞滤波器设计的分频器转换成二阶节形式六、实验结果讨论1.对巴特沃思滤波器设计的分频器及滤波效果进行讨论根据调用tf2sos函数得到的数据可以得出二阶节级联形式的分频器的系统函数: 从figure1分频器的幅度特性曲线可以看出巴特沃思滤波器设计的分频器整体较为平整,高低通并联而成的分频器系统在分频点2000hz的地方有3dB左右的偏差。粗略可以看出设计的数字分频器满足设计要求;而从figure2
13、到figure4这三幅零极点分布和幅频特性图可以看出巴特沃思高低通滤波器和分频器的极点都在单位圆内,所以设计的分频器系统是稳定的,分频器的幅度比较平缓,在约0.1的地方有提升;将figure10中巴特沃思滤波器设计的分频器处理后生成的wav文件打开,和原白噪声wav文件进行对比,听觉上未发现有差别,基本一致;将该分频器处理后的效果进行频谱分析,如figure8所示,可以看到处理得到的频谱在0.1左右的地方与原来的有所差别,可以与figure4显示的分频器在该处的幅值达到1.4(大约3dB的偏差)相符合,设计结果在可接受范围内。2.对宁可瑞滤波器设计的分频器及滤波效果进行讨论 从figure1分
14、频器的幅度特性曲线可以看出宁可瑞滤波器设计的分频器幅度特性曲线很平坦,高低通并联而成的分频器系统在分频点2000hz的地方幅度也基本没有偏差,满足设计要求;而从figure5到figure7这三幅零极点分布和幅频特性图可以看出宁可瑞高低通滤波器和分频器的极点都在单位圆内,所以设计的分频器系统是稳定的,分频器的幅度比较平缓,在约0.1的地方(分频点)有小波动,不过稳定保持在1;将figure10中宁可瑞滤波器设计的分频器处理后生成的wav文件打开,和原白噪声wav文件进行对比,听觉上未发现有差别,基本一致;将该分频器处理后的效果进行频谱分析,如figure9所示,可以看到处理得到的频谱基本跟原音
15、频的频谱图一致,符合设计要求。七、结论: 本实验设计内容是二分频的数字分频器,输入信号的低通和高通分量通过设计的并联系统分两路输出,实际中常采用分频器将音频信号进行处理,将高频部分输出在小音箱,低频部分输出在大音箱。而实验中设计的较为理想的分频器相当于全通滤波器,通过分频器处理得到的wav音频文件已经将低频高频部分叠加在一起,所以处理后的效果应该与原信号一致。根据实验结果可以得出结论:设计的二分频数字分频器符合设计要求,实现的分频效果较好。而宁可瑞滤波器设计的分频器的效果比巴特沃思滤波器设计的分频器的效果要好,基本实现信号的全通,与原信号基本一致。八、参考文献1高西全、丁玉美、阔永红 数字信号处理原理、实现与应用(第2版) 电子工业出版社 2012年5月2董长虹、高成、金涛 Matlab信号处理与应用 国防工业出版社 20053张志涌、杨祖樱 MATLAB教程 北京航空航天大学出版社 2011年7月
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1