【正文】 data, handles, varargin) % This function has no output args, see OutputF. % hObject handle to figure % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % varargin mand line arguments to untitled (see VARARGIN) % Choose default mand line output for untitled = hObject。, [] , ... 39。, mfilename, ... 39。,hObject,eventData,handles,...) calls the local % function named CALLBACK in with the given input arguments. % % UNTITLED(39。 回顧起此課程設(shè)計，至今我仍感慨頗多，從理論到實踐，在這段日子里，可以說得是苦多于甜，但是可以學(xué)到很多很多的東西，同時不僅可以鞏固了以前所學(xué)過的知識，而且學(xué)到了很多在書本上所沒有學(xué)到過的知 識。實踐出真知，通過親自動手制作，使我們掌握的知識不再是紙上談兵。 6 結(jié)果分析 GUI 設(shè)計界面 對信號采樣后，信號時域波形出現(xiàn)了離散化的分布，不像原始圖像那樣密集分布，頻譜圖也是同樣的變換，幅度和頻率都是根據(jù)采樣區(qū)間變化的 從時域圖出延時是讓原來的波形向右移動 3000，頻域的幅度發(fā)生了變化。全通結(jié)構(gòu)幅頻 39。 subplot(3,2,3:4)。 title(39。 R=5000。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。)。plot(y2)。無限個回聲濾波器幅頻響應(yīng) 39。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz=[0,zeros(1,R1),1]。,[50000 500000])。 subplot(3,2,6)。title(39。)。%多重回聲濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。 z=[zeros(3000,1)。title(39。)。 plot(angle(h))。%單回聲濾波器系統(tǒng)函數(shù) Y2=fft(y2,6001)。y]。混響的相頻 39。 subplot(2,2,3)。zeros(3000,1)]。 sound(y,fs,bits)。title(39。%延遲 3000 Z=fft(z,6001)。)。plot(abs(Y))。%只取原信號的 一部分，即時間區(qū)間 [50000 150000] y=y(:,1)。 subplot(2,2,4)。title(39。 （ 2） 原始語音信號時域波形和頻譜 function pushbutton1_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。 if nargin amp。gui_OpeningF39。 GUI 設(shè)計面板是 GUI 設(shè)計工具應(yīng)用的平臺 ,面板上部提供了菜單和常用工具按鈕 ,左邊提供了多種 GUI 控件 ,如按鈕、單選按鈕、復(fù)選框、文本框等。 （ 3）無限個回聲濾波器的系統(tǒng)函數(shù) : RRazzzH???? 1)( , a 1 （ 4）全通結(jié)構(gòu)的混響器的系統(tǒng)函數(shù) : 全通濾波器的傳遞函數(shù)公式為 H(Z)=y(Z)/X(Z)=(K+Z^(m))/(1K*Z^(m)) 其中 m 為回聲延時取樣， k 為反饋系數(shù)。 *均衡頻率（ EQ F），這里的頻率均衡用于音色調(diào)整，此為均衡的中點頻率選擇?；匾籼幚硪话愣际怯糜诋a(chǎn)生簡單的混響效果。 聲陣密度：此項參數(shù)可調(diào)整聲陣的密度，其值較高時，混響效果較為溫暖，但有明顯的聲染色；其值較低時，混響效果較深邃，切聲染色也較弱。 混響與延時 (1)混響效果主 要是用于增加音源的融合感。 wav 文件是一種數(shù)字聲音文件格式 ,本課程設(shè)計基于 Matlab 分析了 wav 聲音文件頻譜與聲音的關(guān)系 。用 MATLAB 提供的 GUI界面開發(fā)環(huán)境設(shè)計用戶圖形界面使布局編程簡化。 數(shù)字信號處理主要是研究用數(shù)字或符號序列表示和處理信號。 (2)將音頻文件加載到 matlab 中，用合適的方式對信號經(jīng)行處理 (3)回答以下四個問題：為什么要選擇這段語音信號？你怎么處理這段語音信號？為什么你選擇用這種方式進行處理？處理結(jié)果是什么？ 3 設(shè)計的總體方案與原理 基本思路 根據(jù)設(shè)計要求，我們商定的處理方式是基于 matlab 來進行聲音的延遲與混響 設(shè)計流程 錄制一段語音信號，時長為 15s，并對錄制的信號進行采樣 語音信號的頻譜分析，畫出采樣后語音信號的時域波形和頻譜圖 將信號加入延時和混響，再分析其頻譜，并與原始信號頻譜進行比較 設(shè)計幾種特殊類型的濾波器：單回聲濾波器，多重回聲濾波器，無限個回聲濾波器，全通結(jié)構(gòu)的混響器，并畫出濾波器的頻域響應(yīng) 用自己設(shè)計的濾波器對采集的語音信號進行濾波 分析得到信號的頻譜，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進行對比，分析信號的變化 回放語音信號 設(shè)計的詳細原理 信號采樣 （ 1）采樣頻率 采樣頻率是指計算機每秒鐘采集多少個聲音樣本，采樣頻率越高，即 采樣的間隔時間越短，則在單位時間內(nèi)計算機得到的聲音樣本數(shù)據(jù)就越多，對聲音波形的表示也越精確。 高頻滾降：此項參數(shù)用于模擬自然混響當中，空氣對高頻的吸收效應(yīng)，以產(chǎn)生較為自然的混響效果。 調(diào)治深度：指上述調(diào)頻電路的調(diào)治深度。 *反饋高頻比（ Hi Ratio），即反饋回路上的高頻衰減控制。其中向量 b和 a分別表示系統(tǒng)函數(shù)的分子、分母多項式的系數(shù)，若 a＝ 1，此時表示 FIR濾波器，否則就是 IIR濾波器。 GUI 界面 圖形用戶界面 GUI (Graphics User Interface) 是由各種圖形對象 ,如圖形窗口、圖軸、菜單、按鈕、文本框等構(gòu)建的用戶界面 ,是人機交流信息的工具和方法。 gui_State = struct(39。, untitled_OutputF, ... 39。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。%讀取原始音頻信號 y=y(:,1)。 subplot(2,2,3)。原始信號相頻 39。plot(y)。)。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。title(39。 subplot(2,2,4)。,[50000 150000])。 subplot(2,2,1:2)?；祉懙姆l 39。%回放混響信號 （ 6） 單回聲濾波器程序設(shè)計 function pushbutton5_Callback(hObject, eventdata, handles) [y,fs,bits]=wavread(39。 a=。plot(abs(h))。)。title(39。 sound(y2,fs,bits)。zeros(3000,1)]。 [h,w]=freqz(Bz1,Az1)。 title(39。 subplot(3,2,5)。多重回聲濾波器相頻 39。y]。 Y2=fft(y2,6001)。 plot(angle(h))。)。title(39。 z=[zeros(3000,1)。%無限個回聲濾波器系統(tǒng)函數(shù)分母 [h,w]=freqz(Bz1,Az1)。 subplot(3,2,2)。全通結(jié)構(gòu)時域圖 39。plot(angle(YY2))。 時域上形 狀大致相同，幅度有微小變化。這次課程設(shè)計終于順利完成了，在設(shè)計中遇到了很多問題，最后我們大家的不懈努力下，終于游逆而解。 實驗過程中，也對團隊精神的進行了考察，讓我們在合作起來更加默契，在成功后一起體會喜悅的心情。Value39。gui_OpeningF39。 if nargin amp。 % 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)。%對原始信號取傅里葉變換 subplot(2,2,1:2)。原始信號幅頻 39。C:\Users\Administrator\Desktop\ 音 頻 信號 .wav39。)。title(39。 z=[zeros(20000,1)。 subplot(2,2,3)。延時后相頻 39。y]?；祉懙臅r域波形 39。plot(angle(Y1))。 y=y(:,1)。%單回聲濾波器系統(tǒng)函數(shù)分子 Az1=(1)。)。title(39。 subplot(3,2,6)。,[50000 500000])。 R=9000。 title(39。 subplot(3,2,3:4)。多重回聲濾波器幅頻 39。 % Executes on button press in pushbutton7. function pushbutton7_Callback(hObject, eventdata, handles) % hObject handle to pushbutton7 (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。 a=。plot(abs(h))。)。title(39。 sound(y2,fs,bits)。zeros(3000,1)]。 subplot(3,2,1)。全通濾波器相頻響應(yīng) 39。plot(abs(YY2))。)。 gTXRm 6X4NGpP$vSTTamp。gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm6X4NGpP$vSTTamp。 UE9aQGn8xp$Ramp。 gTXRm6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 UE9aQGn8xp$Ramp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am UE9aQGn8xp$Ramp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpaza dNuKNamp。 ksv*3t nGK8! z89Am YWv*3t nGK8!z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8! z89Am v^$UE9wEwZQcUE% amp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz849Gx^Gj qv^$UE9wEwZQcUE% amp。 ksv*3tnGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWv*3t nGK8!z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkumamp。 gTXRm6X4NGpP$vSTTa