【正文】 )。,39。 set(,39。on39。Enable39。)。,39。 set(,39。,toc))。,sprintf(39。 end set(,39。Value39。)。,39。 27 % 濾波按鈕 set(,39。on39。Enable39。 guidata(hObject,handles)。 =b。 end [H f]=freqz(b,a,512,2)。)。 end else errordlg(39。)。 otherwise errordlg(39。bandstop39。)。 case 3 [b a]=cheby2(order,20,[W1 W2],39。bandstop39。)。) case 1 [b a]=butter(order,[W1 W2],39。)==1 switch get(,39。 end elseif get(,39。)。 otherwise errordlg(39。 case 3 [b a]=cheby2(order,20,[W1 W2])。) case 1 [b a]=butter(order,[W1 W2])。)==1 switch get(,39。 end elseif get(,39。)。 otherwise errordlg(39。high39。)。 case 3 [b a]=cheby2(order,20,W1,39。high39。)。) case 1 [b a]=butter(order,W1,39。)==1 switch get(,39。 end elseif get(,39。)。 otherwise errordlg(39。low39。)。 25 case 3 [b a]=cheby2(order,20,W1,39。low39。)。) case 1 [b a]=butter(order,W1,39。)==1 switch get(,39。 end if get(,39。 if W11 W1=2*W1/。String39。))。 W1=str2double(get(,39。String39。)。,39。 set(,39。on39。Enable39。)。,39。 set(,39。on39。Enable39。 end guidata(hObject,handles)。 =+y。,i))。)。Value39。k=get(%d,39。)。 else errordlg(39。Value39。)==1 [b a]=cheby2(order,20,2*fband(i:i+1)/)。 elseif get(,39。Value39。)==1 [b a]=butter(order,2*fband(i:i+1)/)。 for i=1:num if get(,39。 =zeros(size())。().*ones(1,num+1)])。%8 filters totally fband=。))。 order=str2double(get(,39。)。,39。 set(,39。off39。Enable39。))。,get(0,39。 else set(hObject,39。white39。BackgroundColor39。))。,get(0,39。 else set(hObject,39。BackgroundColor39。、 % 設(shè)置滑動條的屬性 usewhitebg = 1。on39。Enable39。 to=to+n。 ylim(,[ymin ymax])。 ylim(,[0 60])。 for i=1:num yb(i)=60+20*log10(mean(Pyy(floor(fband(i)*n/):floor(fband(i+1)*n/))))。 loglog(,f,Pyy(1:n/2))。 tic P=fft(y,n)。async39。 pause(steptoc*)。 tic while to=N。 22 fband=。 ymin=min()。 to=n。 % fs 是采樣速率（ Hz）。 n=4096。 % step=。%播放 WAV文件 % 播放按鈕 . N=size()。async39。% 是采樣數(shù)據(jù) guidata(hObject,handles)。, fullfile(pathname, filename)])%文件名和路徑 [ bits]=wavread(filename)。) else disp([39。%生成打開文件對話框 if isequal(filename,0) disp(39。Pick an WAVEfile39。*.wav39。 最后感謝基信學(xué)院和我的母校 — 云南農(nóng)業(yè)大學(xué)四年來對我的教導(dǎo)與栽培。同時還要感謝大學(xué)四年來所有的老師，為我們打下電信專業(yè)知識的基礎(chǔ)；也還要感謝所有的同學(xué)們，正是因為有了你們的支持和鼓勵。 在這里首先要感謝我的指導(dǎo)老師 趙紅偉老師。 2020. [8]． MATLAB語音處理及合成 [M]．北京 ： 清華大學(xué) 出版社。王玲 精通 MATLAB GUI 設(shè)計 .[M] 電子工業(yè)出版社 .. [6]胡昌華．基于 MATLAB的系統(tǒng)分析與設(shè)計 [M]．西安：電子科技大學(xué)出版社。帶通濾波器頻率交界處的 誤差可以通過優(yōu)化設(shè)計的方法來將其減小，即各帶通濾波器不完全以給定的頻帶為自己的生成參數(shù)，而是以附近的某個頻率值作為其生成參數(shù)，這樣就需要用到優(yōu)化的理論。最后給出了軟件的運行結(jié)果，并對軟件進行了誤差分析。界面上采用 MATLAB提供的 GUIDE工具進行設(shè)計，使得軟件開發(fā)效率較高。 18 5．總結(jié)與展望 總結(jié) 本文首先對數(shù)字濾波器的原理進行了分析，并 選擇了一種實現(xiàn)數(shù)字均衡器的方案，在軟件的設(shè)計上采用軟件工程的理論，使用面向數(shù)據(jù)流的設(shè)計方法進行軟件設(shè)計，并且采用了模塊劃分的方法。故提高精度關(guān)鍵在于提高濾波器的設(shè)計精度。對于 IIR 數(shù)字濾波器，階數(shù)越高，則 濾波器的精度越高，即越接近于理想的情況，但計算機運行速度就越慢了，所以實際中應(yīng)該權(quán)衡運行速度與濾波器的精度。圖 9 為選擇高通濾波器時的運行界面，它和低通剛好相反，只允許高頻的通過。如圖7所示，濾波器類型為低通濾波器，其條形圖和頻譜圖反映只允許低頻通過，其高頻基本都被濾去，故實現(xiàn)了低通濾波器的功能。但是并不能完全濾掉某個頻段的信號。為了使得播放流暢，需要與 pause 的使用結(jié)合，具體方法為用 tic 和 toc 計算出其它代碼運行時間，然后用播放時間減去其它代碼運行時間，從而得到需要暫停運行的時間。 to=to+n。 ylim(,[ymin ymax])。 axis(,[fband(1) fband(num+1) 1e3 ymax])。 Pyy =2*abs(P)/n。)。 wavplay(y,39。 y=(from:to,:)。 num=length(fband)1。 f=lins