资源描述：
《线性相位fir数字均衡器设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、汕头大学工学院三级项目报告课程名称：数字信号处理项目题目：线性相位FIR数字均衡器设计指导教师：李旭涛系别：电子工程系专业：电子信息工程学号：姓名：合作者完成时间：2011年11月25日至11月26日成绩：评阅人：李旭涛线性相位FIR数字均衡器的设计千千静听、酷狗等常用音乐播放软件都有均衡器，界面如下面所示：电脑音量输出的均衡器界面：设计思路：人们能听到的声音的频谱范围大约为20Hz—20kHz，由于20KHz远大于20Hz，为了方便，假设人们能听到的声音的频谱范围大约为0—20kHz。根据酷狗均衡器的界面，我们把这20kHz分成8个频段，调节八个频
2、段的信号强度就能达到均衡器的效果。这八个频段的分频点的值为fa=2.5kHz,fb=5kHz,fc=7.5kHz,fd=10kHz,fe=12.5kHz,ff=15.0kHz,fg=17.5kHz;采样频率为fs=40kHz。则8个频段的频率分别为[0,2500Hz]频带1[2500,5000Hz]频带2[5000,7500Hz]频带3[7500,10000Hz]频带4[10000,12500Hz]频带5[12500,15000Hz]频带6[15000,17500Hz]频带7[17500,20000kHz]频带8各个频段加凯泽窗后所对应的单位冲激响应
3、为其中M=(N-1)/2;选取为凯泽窗；8个不同频段的单位冲激响应对应8个不同的滤波器，经过滤波后得到各频段的幅频响应
4、H(f)
5、,把这8个频段组合起来作为输出，即得到8频段的均衡器。Matlab代码如下：%8频段均衡器clearNwin=205;%凯泽窗的阶数dealt=10;%凯泽窗的形状参数Nf=512;fs=40000;M=(Nwin-1)/2;n=0:Nwin-1;fa=2500;fb=5000;fc=7500;fd=10000;fe=12500;ff=15000;fg=17500;%各频率点频率值wa=2*pi*fa/fs;wb=2*pi
6、*fb/fs;wc=2*pi*fc/fs;wd=2*pi*fd/fs;we=2*pi*fe/fs;%频率归一化wf=2*pi*ff/fs;wg=2*pi*fg/fs;dap=sin(wa*(n-M))./(pi*(n-M)+eps);%频带1滤波器的单位冲激响应dbp=(sin(wb*(n-M))-sin(wa*(n-M)))./(pi*(n-M)+eps);%频带2滤波器的单位冲激响应dcp=(sin(wc*(n-M))-sin(wb*(n-M)))./(pi*(n-M)+eps);%频带3滤波器的单位冲激响应ddp=(sin(wd*(n-M))-
7、sin(wc*(n-M)))./(pi*(n-M)+eps);%频带4滤波器的单位冲激响应dep=(sin(we*(n-M))-sin(wd*(n-M)))./(pi*(n-M)+eps);%频带5滤波器的单位冲激响应dfp=(sin(wf*(n-M))-sin(we*(n-M)))./(pi*(n-M)+eps);%频带6滤波器的单位冲激响应dgp=(sin(wg*(n-M))-sin(wf*(n-M)))./(pi*(n-M)+eps);%频带7滤波器的单位冲激响应dhp=[(n-M)==0]-sin(wg*(n-M))./(pi*(n-M)+e
8、ps);%频带8滤波器的单位冲激响应w_kaiser=kaiser(Nwin,dealt);%凯泽窗wn=w_kaiser';hap=wn.*dap;%对单位冲激响应加凯泽窗hbp=wn.*dbp;hcp=wn.*dcp;hdp=wn.*ddp;hep=wn.*dep;hfp=wn.*dfp;hgp=wn.*dgp;hhp=wn.*dhp;%求冲激响应的傅里叶变换[y1,f]=freqz(hap,1,Nf,fs);[y2,f]=freqz(hbp,1,Nf,fs);[y3,f]=freqz(hcp,1,Nf,fs);[y4,f]=freqz(hdp,
9、1,Nf,fs);[y5,f]=freqz(hep,1,Nf,fs);[y6,f]=freqz(hfp,1,Nf,fs);[y7,f]=freqz(hgp,1,Nf,fs);[y8,f]=freqz(hhp,1,Nf,fs);mag1=abs(y1);mag2=abs(y2);mag3=abs(y3);mag4=abs(y4);mag5=abs(y5);mag6=abs(y6);mag7=abs(y7);mag8=abs(y8);%各个波段的幅频图subplot(4,2,1),plot(f,20*log(mag1));xlabel('f/Hz');y
10、label('
11、H(f)
12、/dB');gridon;title('频带1滤波器的幅频图');subplot(
