【正文】 39。, untitled_OutputF, ... 39。, [] , ... 39。, [])。amp。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。 end function untitled_OpeningF(hObject, eventdata, handles, varargin) = hObject。 function varargout = untitled_OutputF(hObject, eventdata, handles) varargout{1} = 。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。%讀取原始音頻信號 y=y(:,1)。%播放原始信號 Y=fft(y)。plot(y)。原始信號時域波形 39。 subplot(2,2,3)。title(39。)。plot(angle(Y))。原始信號相頻 39。 （ 3） 采樣后語音信號的時域波形和頻譜 function pushbutton2_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。,[50000 500000])。 Y=fft(y,6001)。plot(y)。采樣信號時域波形 39。 subplot(2,2,3)。title(39。)。plot(angle(Y))。采樣信號相頻 39。 sound(y,fs,bits)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。 y=y(:,1)。y]。 subplot(2,2,1:2)。title(39。)。plot(abs(Z))。延時后幅頻 39。 subplot(2,2,4)。title(39。)。%回放延遲信號 （ 5） 對采樣后的信號混響 function pushbutton4_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。,[50000 150000])。 z=[zeros(3000,1)。 y=[y。%補 3000 y1=y+z。 subplot(2,2,1:2)。title(39。)。plot(abs(Y1))。混響的幅頻 39。 subplot(2,2,4)。title(39。)。%回放混響信號 （ 6） 單回聲濾波器程序設(shè)計 function pushbutton5_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。,[50000 500000])。 z=[zeros(3000,1)。 y=[y。 a=。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[1,zeros(1,R1),a]。%單回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。 [h,w]=freqz(Bz1,Az1)。plot(abs(h))。單回聲濾波器幅頻響應(yīng) 39。 subplot(3,2,2)。 title(39。)。plot(y2)。單回聲濾波器時域圖 39。 subplot(3,2,5)。title(39。)。plot(angle(Y2))。單回聲濾波器相頻 39。 sound(y2,fs,bits)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。 y=y(:,1)。y]。zeros(3000,1)]。 N=5。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[1,zeros(1,N*R1),a^N]。%多重回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。 [h,w]=freqz(Bz1,Az1)。plot(abs(h))。多重回聲濾波器幅頻響應(yīng) 39。 subplot(3,2,2)。 title(39。)。plot(y2)。多重回聲濾波器時域圖 39。 subplot(3,2,5)。title(39。)。plot(angle(Y2))。多重回聲濾波器相頻 39。 sound(y2,fs,bits)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。 y=y(:,1)。y]。zeros(3000,1)]。 R=5000。%無限個回聲濾波器系統(tǒng)函數(shù)分子 Az=[1,zeros(1,R1),a]。 Y2=fft(y2,6001)。 subplot(3,2,1)。 title(39。)。 plot(angle(h))。無限個回聲濾波 器相頻響應(yīng) 39。 subplot(3,2,3:4)。title(39。)。plot(abs(Y2))。無限個回聲濾波器信號幅頻 39。 subplot(3,2,6)。title(39。)。 (9) 全通結(jié)構(gòu)濾波器程序設(shè)計 function pushbutton8_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。,[50000 500000])。 z=[zeros(3000,1)。 y=[y。 a=。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[a,zeros(1,R1),1]。%無限個回聲濾波器系統(tǒng)函數(shù)分母 [h,w]=freqz(Bz1,Az1)。 YY2=fft(yy2,6001)。plot(abs(h))。全通濾波器幅頻響應(yīng) 39。 subplot(3,2,2)。 title(39。)。plot(yy2)。全通結(jié)構(gòu)時域圖 39。 subplot(3,2,5)。title(39。)。plot(angle(YY2))。全通結(jié)構(gòu)相頻 39。 sound(yy2,fs,bits)。 混響：在時域圖看前面幾乎沒變化，后面的幅度有一定變化。 時域上形 狀大致相同，幅度有微小變化。 時域上濾波前的幅度比濾波后的數(shù)值大，形狀上濾波后是連續(xù)的，濾波前是間斷的 頻譜上變化兩邊上下跳的動變化，中間變化抵消 時域上波形變得稀疏了，幅度相對變化較大。 7 總結(jié)及心得體會 選擇這段音頻的原因：這段音頻雖然只有 15秒，但是它無伴奏音樂，因此處理起來比 較得心應(yīng)手，且聲音優(yōu)美清晰 處理目的：基于 matlab實現(xiàn)語音信號的延遲與混響 處理方式： 處理結(jié)果：利用 matlab成功實現(xiàn)了語音信號的延遲與混響 通過此次課程設(shè)計，使我更加扎實的掌握了有關(guān)語音信號處理方面的知識，在設(shè)計過程中雖然遇到了一些問題，但經(jīng)過一次又一次的思考，一遍又一遍的檢查終于找出了原因所在，也暴露出了前期我在這方面的知識欠缺和經(jīng)驗不足。 過而能改，善莫大焉。這次課程設(shè)計終于順利完成了，在設(shè)計中遇到了很多問題，最后我們大家的不懈努力下，終于游逆而解。同時，設(shè)計讓我感觸很深。通過這次課程設(shè)計，我掌握了一些 matlab編程的常用指令，同時也對 matlab運用得更加得心應(yīng)手，同時也和課上老師講解的內(nèi)容互相印證，能把書本上的知識轉(zhuǎn)變?yōu)榉抡娼Y(jié)果，著實讓我們很高興。通過這次課程設(shè)計使我懂得了理論與實際相結(jié)合是很重要的，只有理論知識是遠遠不夠的，只有把所學(xué)的理論知識與實踐相結(jié)合起來，從理論中得出結(jié)論，才能真正為社會服務(wù)，從而提高自己的實際動手能力和獨立思考的能力。 實驗過程中，也對團隊精神的進行了考察，讓我們在合作起來更加默契，在成功后一起體會喜悅的心情。 此次設(shè)計也讓我明白了思路即出路，有什么不懂不明白的地方要及時請教或 上網(wǎng)查詢，只要認(rèn)真鉆研，動腦思考，動手實踐，就沒有弄不懂的知識，收獲頗豐。CALLBACK39。Property39。Value39。s Tools menu. Choose GUI allows only one % instance to run (singleton). % % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help untitled % Last Modified by GUIDE 21Dec2020 21:01:58 % Begin initialization code DO NOT EDIT gui_Singleton = 1。gui_Name39。gui_Singleton39。gui_OpeningF39。gui_OutputF39。gui_LayoutF39。gui_Callback39。 if nargin amp。 ischar(varargin{1}) = str2func(varargin{1})。 else gui_mainf(gui_State, varargin{:})。 % Update handles structure guidata(hObject, handles)。 % Outputs from this function are returned to the mand line. function varargout = untitled_OutputF(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT)。 % Executes on button press in pushbutton1. function pushbutton1_Callback(hObject, eventdata, handles) % hObject handle to pushbutton1 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [y,fs,bits]=wavread(39。)。%只取單信道 sound(y,fs,bits)。%對原始信號取傅里葉變換 subplot(2,2,1:2)。title(39。)。plot(abs(Y))。原始信號幅頻 39。 subplot(2,2,4)。title(39。)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。%只取原信號的一部分，即時間區(qū)間 [50000 150000] y=y(:,1)。%抽取 6001點 subplot(2,2,1:2)。title(39。)。plot(abs(Y))。采樣信號幅頻 39。 subplot(2,2,4)。title(39。)。%回放采樣信號 % Executes on button press in pushbutton3. function pushbutton3_Callback(hObject, eventdata, handles) % hObject handle to pushbutton3 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [y,fs,bits]=wavread(39。,[50000 500000])。 z=[zeros(20000,1)。%延遲 3000 Z=fft(z,6001)。plot(z)。延時后時域波形 39。 subplot(2,2,3)。title(39。)。plot(angle(Z))。延時后相頻 39。 sound(y,fs,bits)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。 y=y(:,1)。y]。zeros(3000,1)]。%將原始信號與延遲信號相加，產(chǎn)生混響信號 Y1=fft(y1,6001)。plot(y1)。混響的時域波形 39。 subplot(2,2,3)。title(39。)。plot(angle(Y1))?；祉懙南囝l 39。 sound(y,fs,bits)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。 y=y(:,1)。y]。zeros(3000,1)]。 R=5000。%單回聲濾波器系統(tǒng)函數(shù)分子 Az1=(1)。%單回聲濾波器系統(tǒng)函數(shù) Y2=fft(y2,6001)。%求設(shè)計的濾波器頻譜 subplot(3,2,1)。 title(39。)。 plot(angle(h))。單回聲濾波器相頻響應(yīng) 39。 subplot(3,2,3:4)。title(39。)。plot(abs(Y2))。單回聲濾波器幅頻 39。 subplot(3,2,6)。title(39。)。 % Executes on button press in pushbutton6. function pushbutton6_Callback(hObject, eventdat