【正文】 1 [b a]=butter(order,W1,39。 end elseif get(,39。 otherwise errordlg(39。)。low39。) case 1 [b a]=butter(order,W1,39。 end if get(,39。String39。 W1=str2double(get(,39。)。 set(,39。Enable39。,39。on39。 end guidata(hObject,handles)。,i))。Value39。)。Value39。 elseif get(,39。)==1 [b a]=butter(order,2*fband(i:i+1)/)。 =zeros(size())。%8 filters totally fband=。 order=str2double(get(,39。,39。off39。))。 else set(hObject,39。BackgroundColor39。,get(0,39。BackgroundColor39。on39。 to=to+n。 ylim(,[0 60])。 loglog(,f,Pyy(1:n/2))。async39。 tic while to=N。 ymin=min()。 % fs 是采樣速率（ Hz）。 % step=。async39。, fullfile(pathname, filename)])%文件名和路徑 [ bits]=wavread(filename)。%生成打開(kāi)文件對(duì)話框 if isequal(filename,0) disp(39。*.wav39。同時(shí)還要感謝大學(xué)四年來(lái)所有的老師，為我們打下電信專業(yè)知識(shí)的基礎(chǔ)；也還要感謝所有的同學(xué)們，正是因?yàn)橛辛四銈兊闹С趾凸膭?lì)。 2020. [8]． MATLAB語(yǔ)音處理及合成 [M]．北京 ： 清華大學(xué) 出版社。帶通濾波器頻率交界處的 誤差可以通過(guò)優(yōu)化設(shè)計(jì)的方法來(lái)將其減小，即各帶通濾波器不完全以給定的頻帶為自己的生成參數(shù)，而是以附近的某個(gè)頻率值作為其生成參數(shù)，這樣就需要用到優(yōu)化的理論。界面上采用 MATLAB提供的 GUIDE工具進(jìn)行設(shè)計(jì)，使得軟件開(kāi)發(fā)效率較高。故提高精度關(guān)鍵在于提高濾波器的設(shè)計(jì)精度。圖 9 為選擇高通濾波器時(shí)的運(yùn)行界面，它和低通剛好相反，只允許高頻的通過(guò)。但是并不能完全濾掉某個(gè)頻段的信號(hào)。 to=to+n。 axis(,[fband(1) fband(num+1) 1e3 ymax])。)。 y=(from:to,:)。 f=linspace(0,n/2)。 from=1。即先取一批信號(hào)，將其播放出來(lái)，同時(shí)進(jìn)行 FFT，并繪出時(shí)域波形與頻譜分析結(jié)果，然后取下一批信號(hào)??由于各批信號(hào)間的時(shí)間間隔很小，從而感覺(jué)像連續(xù)播放，而時(shí)域波形與頻譜則是實(shí)時(shí)地變化。輸出序列的長(zhǎng)度與輸入序列的長(zhǎng)度相同。 %保存數(shù)據(jù) end 在程序中有一句 y=(9^k1)/8*qfilter(b,a,)，作用是加權(quán)以便下一句的求和，其中 k是從滾動(dòng)條中得到的參數(shù)，從 0到 2，但需要的權(quán)值是從 0 到 10，令 w 為需要的權(quán)，即 ，由于實(shí)際 需要 w 為指數(shù)變化，故設(shè) (7) 代入 w=0, k=0。,i))。Value39。)。Value39。 elseif get(,39。)==1 %巴特沃斯濾波器 [b a]=butter(order,2*fband(i:i+1)/)。 nn=length()。))。 均衡器模塊的實(shí)現(xiàn) 均衡器模塊功能是生成帶通濾波器并對(duì)原信號(hào)濾波?！盀V波”用于對(duì)輸入的音頻信號(hào)用生成的濾波器濾波。左下角是條 形圖，它是許多音頻軟件使用的頻譜輸11 出方式。“原聲播放”用于播放輸入的原信頻信號(hào)，區(qū)別于濾波后的信號(hào)，用于對(duì)比濾波前后的信號(hào)。故本文采用MATLAB 作為主要編程語(yǔ)言實(shí)現(xiàn)數(shù)字均衡器。濾波模塊實(shí)現(xiàn) 的功能是濾波器的生成和實(shí)時(shí) 圖形顯示模塊則實(shí)現(xiàn)波形的實(shí)時(shí)顯示和頻譜的實(shí)時(shí)顯示。 在模塊劃分時(shí)應(yīng)遵循如下規(guī)則：改進(jìn)軟件結(jié)構(gòu)提高模塊獨(dú)立性；模塊規(guī)模應(yīng)該適中；深度、寬度、扇出和扇入都應(yīng)適當(dāng) ；模塊的作用域應(yīng)該在控制域之內(nèi)；力爭(zhēng)降低模塊接口的復(fù)雜程度；設(shè)計(jì)單入口單出口的模塊；模塊功能應(yīng)該可以預(yù)測(cè)。結(jié)合上一節(jié)的內(nèi)容，可以得出軟件的數(shù)據(jù)流圖如圖 2所示。 軟件設(shè)計(jì) 數(shù)據(jù)流圖 數(shù)據(jù)流圖 (DFD)是一種圖形化技術(shù)，它描繪信息流和數(shù)據(jù)從輸入移動(dòng)到輸出的過(guò)程中所經(jīng)受的變換。上述方法得到的是低通濾波器，利用變換公式作頻率變換得到高通、帶通濾波器。下面介紹 IIR 濾波器的設(shè)計(jì)方法。 FIR 濾波器沒(méi)有現(xiàn)在 的設(shè)計(jì)公式，對(duì)計(jì)算工具要求較高，需要借助計(jì)算機(jī)來(lái)設(shè)計(jì)。 FIR 濾波器可得到嚴(yán)格的相位，而 IIR 濾波器不能得到。 (2) FIR 濾波器 FIR 濾波器的特點(diǎn)是：系統(tǒng)的單位沖激響應(yīng) h(n)是個(gè)有限長(zhǎng)序列；系統(tǒng)函數(shù) |H(z)|在|z|0 處收斂，極點(diǎn)全部在 z=0 處 (即 FIR 一定為穩(wěn)定系統(tǒng) )；結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒(méi)有輸出到輸入反饋?？紤]到軟件的靈活性及易于實(shí)現(xiàn)，本文采用軟件方式實(shí)現(xiàn)數(shù)字濾波器。 頻段分得越細(xì)，調(diào)節(jié)的峰越尖銳，即 Q 值（品質(zhì)因數(shù)）越高， 調(diào)節(jié)時(shí)補(bǔ)償?shù)迷郊?xì)致，頻段分的越粗則調(diào)節(jié)的峰就比較寬，當(dāng)聲場(chǎng)傳輸頻率特性曲線比較復(fù)雜時(shí)較難補(bǔ)償。這樣人們根據(jù)不同的要求分別選擇不同段數(shù)的頻率均衡器。頻響多峰谷時(shí)：聲音粗糙混濁，高音刺耳發(fā)毛，無(wú)層次感擴(kuò)聲易發(fā)生反饋嘯叫。 4K~8Kz 清脆 6800Hz 形成尖嘯，銳利的感覺(jué)。 100~200Hz 豐滿 200~500Hz 力度 易引起嗡嗡聲的煩悶心理。通過(guò)調(diào)節(jié) 均衡器可以得 到最佳音樂(lè)狀態(tài)。Digital filter； 3 1 前言 隨著數(shù)字化技術(shù)的快速、深入發(fā)展 ,人們對(duì)數(shù)字化電子產(chǎn)品所產(chǎn)生的圖像、圖形以及聲音等質(zhì)量的要求越來(lái)越高。附加的工具箱（單獨(dú)提供的專用 MATLAB 函數(shù)集）擴(kuò)展了 MATLAB 環(huán)境， 以解決這些應(yīng)用領(lǐng)域內(nèi) 特定類型的問(wèn)題 。 18 展望 14 運(yùn)行實(shí)例 12 4運(yùn)行實(shí)例與實(shí)驗(yàn) 分析 12 實(shí)時(shí)顯示的實(shí)現(xiàn) 9 3軟件實(shí)現(xiàn) 7 軟件設(shè)計(jì) 7 IIR 數(shù)字濾波器的設(shè)計(jì)方法 1 ABSTRACT 3 2設(shè)計(jì)原理 5 均衡器分類 5 數(shù)字濾波器 6 數(shù)字濾波器的原理簡(jiǎn)介 6 FIR 與 IIR 濾波器的比較與選擇 7 均衡器的原理 8 模塊劃分 10 均衡器模塊的實(shí)現(xiàn) 21 1 摘 要 MATLAB 的應(yīng)用范圍非常廣，包括信號(hào)和圖像處理、通訊、控制系統(tǒng)設(shè)計(jì)、測(cè)試和測(cè)量、財(cái)務(wù)建模和分析以及計(jì)算生物學(xué)等眾多應(yīng)用領(lǐng)域。MATLAB?，F(xiàn)有的許多播放器都具有均衡器的功能，如千千靜聽(tīng) (如圖 1 所示 )等。 60~100Hz 沉重 80Hz 附近能產(chǎn)生極強(qiáng)的“重感”效果，響度很高也不會(huì)給人舒服的感覺(jué)，可給人以強(qiáng)烈的刺激作用。 2K~4Kz 尖銳 3400Hz 易引起聽(tīng)覺(jué)疲勞。 整個(gè)頻響特性平直時(shí)：聲音自然豐滿而有彈性，層次清晰 圓 順悅耳。常用的專業(yè)圖示均衡器則是將 20Hz~20kHz 的信號(hào)分成 8 段、 10 段、 15 段、 27 段、 31 段來(lái)進(jìn)行調(diào)節(jié)。 3．房間均衡器，用于調(diào)整房間內(nèi)的頻率響應(yīng)特性曲線的均衡器，由于裝飾材料對(duì)不同頻率的吸收（或反射