【正文】 t, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 500000])。%只取原信號的 一部分，即時間區(qū)間 [50000 150000] y=y(:,1)。 Y=fft(y,6001)。%抽取 6001點(diǎn) subplot(2,2,1:2)。plot(y)。title(39。采樣信號時域波形 39。)。 subplot(2,2,3)。plot(abs(Y))。title(39。采樣信號幅頻 39。)。 subplot(2,2,4)。plot(angle(Y))。title(39。采樣信號相頻 39。)。 sound(y,fs,bits)。%回放采樣信號 （ 4） 對采樣后的信號延時 function pushbutton3_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 500000])。 y=y(:,1)。 z=[zeros(20000,1)。y]。%延遲 3000 Z=fft(z,6001)。 subplot(2,2,1:2)。plot(z)。title(39。延時后時域波形 39。)。 subplot(2,2,3)。plot(abs(Z))。title(39。延時后幅頻 39。)。 subplot(2,2,4)。plot(angle(Z))。title(39。延時后相頻 39。)。 sound(y,fs,bits)。%回放延遲信號 （ 5） 對采樣后的信號混響 function pushbutton4_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。%補(bǔ) 3000 y1=y+z。%將原始信號與延遲信號相加，產(chǎn)生混響信號 Y1=fft(y1,6001)。 subplot(2,2,1:2)。plot(y1)。title(39?；祉懙臅r域波形 39。)。 subplot(2,2,3)。plot(abs(Y1))。title(39?；祉懙姆l 39。)。 subplot(2,2,4)。plot(angle(Y1))。title(39?；祉懙南囝l 39。)。 sound(y,fs,bits)。%回放混響信號 （ 6） 單回聲濾波器程序設(shè)計 function pushbutton5_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 500000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[1,zeros(1,R1),a]。%單回聲濾波器系統(tǒng)函數(shù)分子 Az1=(1)。%單回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。%單回聲濾波器系統(tǒng)函數(shù) Y2=fft(y2,6001)。 [h,w]=freqz(Bz1,Az1)。%求設(shè)計的濾波器頻譜 subplot(3,2,1)。plot(abs(h))。 title(39。單回聲濾波器幅頻響應(yīng) 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。單回聲濾波器相頻響應(yīng) 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。單回聲濾波器時域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。單回聲濾波器幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。單回聲濾波器相頻 39。)。 sound(y2,fs,bits)。 （ 7） 多重回聲濾波器程序設(shè)計 function pushbutton6_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 500000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 N=5。 R=9000。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[1,zeros(1,N*R1),a^N]。%多重回聲濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。%多重回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。 Y2=fft(y2,6001)。 [h,w]=freqz(Bz1,Az1)。 subplot(3,2,1)。plot(abs(h))。 title(39。多重回聲濾波器幅頻響應(yīng) 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。多重回聲濾波器相頻響應(yīng) 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。多重回聲濾波器時域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。多重回聲濾波器幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。多重回聲濾波器相頻 39。)。 sound(y2,fs,bits)。 （ 8） 無限個回聲濾波器程序設(shè)計 function pushbutton7_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 500000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz=[0,zeros(1,R1),1]。%無限個回聲濾波器系統(tǒng)函數(shù)分子 Az=[1,zeros(1,R1),a]。%無限個回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz,Az,y)。 Y2=fft(y2,6001)。 [h,w]=freqz(Bz,Az)。 subplot(3,2,1)。plot(abs(h))。 title(39。無限個回聲濾波器幅頻響應(yīng) 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。無限個回聲濾波 器相頻響應(yīng) 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。無限個回聲濾波器時域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。無限個回聲濾波器信號幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。無限個回聲濾波器信號相頻 39。)。 sound(y2,fs,bits)。 (9) 全通結(jié)構(gòu)濾波器程序設(shè)計 function pushbutton8_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。,[50000 500000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[a,zeros(1,R1),1]。%無限個回聲濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。%無限個回聲濾波器系統(tǒng)函數(shù)分母 [h,w]=freqz(Bz1,Az1)。 yy2=filter(Bz1,Az1,y)。 YY2=fft(yy2,6001)。 subplot(3,2,1)。plot(abs(h))。 title(39。全通濾波器幅頻響應(yīng) 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。全通濾波器相頻響應(yīng) 39。)。 subplot(3,2,3:4)。plot(yy2)。title(39。全通結(jié)構(gòu)時域圖 39。)。 subplot(3,2,5)。plot(abs(YY2))。title(39。全通結(jié)構(gòu)幅頻 39。)。 subplot(3,2,6)。plot(angle(YY2))。title(39。全通結(jié)構(gòu)相頻 39。)。 sound(yy2,fs,bits)。 6 結(jié)果分析 GUI 設(shè)計界面 對信號采樣后，信號時域波形出現(xiàn)了離散化的分布，不像原始圖像那樣密集分布，頻譜圖也是同樣的變換，幅度和頻率都是根據(jù)采樣區(qū)間變化的 從時域圖出延時是讓原來的波形向右移動 3000，頻域的幅度發(fā)生了變化。 混響：在時域圖看前面幾乎沒變化，后面的幅度有一定變化。 在頻譜圖幅度變化很大，中間的上下相互抵消，兩邊的跳動很大。 時域上形 狀大致相同，幅度有微小變化。 頻譜上也是濾波前后兩邊跳動變換較大，幅度變化較大。 時域上濾波前的幅度比濾波后的數(shù)值大，形狀上濾波后是連續(xù)的，濾波前是間斷的 頻譜上變化兩邊上下跳的動變化，中間變化抵消 時域上波形變得稀疏了，幅度相對變化較大。 頻譜上變化兩邊上下跳的動變化，中間變化抵消 全通結(jié)構(gòu)時域大致一樣，濾波前比濾波后幅度有變化 6000 以后上下跳動抵消 頻域上變化的不為明顯，幅度發(fā)生了變化。 7 總結(jié)及心得體會 選擇這段音頻的原因：這段音頻雖然只有 15秒，但是它無伴奏音樂，因此處理起來比 較得心應(yīng)手，且聲音優(yōu)美清晰 處理目的：基于 matlab實(shí)現(xiàn)語音信號的延遲與混響 處理方式： 處理結(jié)果：利用 matlab成功實(shí)現(xiàn)了語音信號的延遲與混響 通過此次課程設(shè)計，使我更加扎實(shí)的掌握了有關(guān)語音信號處理方面的知識，在設(shè)計過程中雖然遇到了一些問題，但經(jīng)過一次又一次的思考，一遍又一遍的檢查終于找出了原因所在，也暴露出了前期我在這方面的知識欠缺和經(jīng)驗(yàn)不足。實(shí)踐出真知，通過親自動手制作，使我們掌握的知識不再是紙上談兵。 過而能改，善莫大焉。在課程設(shè)計過程中，我們不斷發(fā)現(xiàn)錯誤，不斷改正，不斷領(lǐng)悟 ，不斷獲取，最終的檢測調(diào)試環(huán)節(jié)，本身就是在踐行“過而能改，善莫大焉”的知行觀。這次課程設(shè)計終于順利完成了，在設(shè)計中遇到了很多問題，最后我們大家的不懈努力下，終于游逆而解。在今后社會的發(fā)展和學(xué)習(xí)實(shí)踐過程中，一定要不懈努力，不能遇到問題就想到要退縮，一定要不厭其煩的發(fā)現(xiàn)問題所在，然后一一進(jìn)行解決，只有這樣，才能成功的做成想做的事，才能在今后的道路上劈荊斬棘，而不是知難而退，那樣永遠(yuǎn)不可能收獲成功，收獲喜悅，也永遠(yuǎn)不可能得到社會及他人對你的認(rèn)可 ! 課程設(shè)計誠然是一門專業(yè)課，給我很多專業(yè)知識以及專業(yè)技能上的提升 ，同時又是一門講道課，一門辯思課，給了我許多道，給了我很多思，給了我莫大的空間。同時，設(shè)計讓我感觸很深。使我對抽象的理論有了具體的認(rèn)識。通過這次課程設(shè)計，我掌握了一些 matlab編程的常用指令，同時也對 matlab運(yùn)用得更加得心應(yīng)手，同時也和課上老師講解的內(nèi)容互相印證，能把書本上的知識轉(zhuǎn)變?yōu)榉抡娼Y(jié)果，著實(shí)讓我們很高興。 回顧起此課程設(shè)計，至今我仍感慨頗多，從理論到實(shí)踐，在這段日子里，可以說得是苦多于甜，但是可以學(xué)到很多很多的東西，同時不僅可以鞏固了以前所學(xué)過的知識，而且學(xué)到了很多在書本上所沒有學(xué)到過的知 識。通過這次課程設(shè)計使我懂得了理論與實(shí)際相結(jié)合是很重要的，只有理論知識是遠(yuǎn)遠(yuǎn)不夠的，只有把所學(xué)的理論知識與實(shí)踐相結(jié)合起來，從理論中得出結(jié)論，才能真正為社會服務(wù)，從而提高自己的實(shí)際動手能力和獨(dú)立思考的能力。在設(shè)計的過程中遇到問題，可以說得是困難重重，但可喜的是最終都得到了解決。 實(shí)驗(yàn)過程中，也對團(tuán)隊(duì)精神的進(jìn)行了考察，讓我們在合作起來更加默契，在成功后一起體會喜悅的心情。果然是團(tuán)結(jié)就是力量，只有互相之間默契融洽的配合才能換來最終完美的結(jié)果。 此次設(shè)計也讓我明白了思路即出路，有什么不懂不明白的地方要及時請教或 上網(wǎng)查詢，只要認(rèn)真鉆研，動腦思考，動手實(shí)踐，就沒有弄不懂的知識，收獲頗豐。 8 參考資料 1 彭啟琮等， DSP 技術(shù)的發(fā)展與應(yīng)用，第二版，高等教育出版社， 2020 2 張雄偉， DSP 集成開發(fā)與應(yīng)用實(shí)例，第 1 版，電子工業(yè)出版社， 2002 3 楊述斌等，數(shù)字信號處理實(shí)踐教程，第一版，華中科技大學(xué)出版社， 2020 9 源代碼 function varargout = untitled(varargin) % UNTITLED Mfile for % UNTITLED, by itself, creates a new UNTITLED or raises th