【正文】 Dsp 課程設(shè)計(jì) 基于 Matlab 的聲音的延遲與混響 摘要 數(shù)字信號(hào)處理技術(shù)自誕生以來(lái)，有了快速的發(fā)展， 主要是研究用數(shù)字或符號(hào)序列表示和處理信號(hào)，被廣泛應(yīng)用于各個(gè)領(lǐng)域。 聲音信號(hào)是一維連續(xù)信號(hào) ,而計(jì)算機(jī)只能處理離散信號(hào)。聲卡的位客觀地反映了數(shù)字聲音信號(hào)對(duì)輸入聲音信號(hào)描述的準(zhǔn)確程度。 擴(kuò)散度：此項(xiàng)參數(shù)可調(diào)整混響聲陣密度 的增長(zhǎng)速度，其可調(diào)范圍為 0~10，其值較高時(shí)，混響效果比較豐厚、溫暖；其值較低時(shí)，混響效果則較空曠、冷僻。 當(dāng)延遲時(shí)間在 3~35ms之間時(shí)人耳感覺不到滯后音的存在，并且他與原音源疊加后，會(huì)因其相位干涉而產(chǎn)生 梳狀濾波 效應(yīng)，這就是鑲邊效果。 *高頻增益（ HF），指高頻均衡控制。直達(dá)聲和在 R 抽樣周 期后出現(xiàn)的一種單個(gè)回音，可以用 FIR 濾波器產(chǎn)生，微分方程為： y[n]=x[n]+α x[nR] |α |1 傳輸函數(shù)為： H(z)=1+α zR 傳遞函數(shù)的幅頻響應(yīng)形狀象梳子，這種濾波器又叫梳狀濾波器。調(diào)用 GUIDE 的方法有 2 種 ,在 MATLAB 命令窗口中輸入 guide 命令 ,或在 MATLAB 主菜單中點(diǎn)擊 File → New → GUI 即可打開一個(gè)可編輯的新窗口。gui_Singleton39。gui_Callback39。 guidata(hObject, handles)。%對(duì)原始信號(hào)取傅里葉變換 subplot(2,2,1:2)。原始信號(hào)幅頻 39。C:\Users\Administrator\Desktop\ 音 頻 信號(hào) .wav39。)。title(39。 z=[zeros(20000,1)。 subplot(2,2,3)。延時(shí)后相頻 39。y]?；祉懙臅r(shí)域波形 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。 （ 8） 無(wú)限個(gè)回聲濾波器程序設(shè)計(jì) function pushbutton7_Callback(hObject, eventdata, handles) [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))。)。 頻譜上變化兩邊上下跳的動(dòng)變化，中間變化抵消 全通結(jié)構(gòu)時(shí)域大致一樣，濾波前比濾波后幅度有變化 6000 以后上下跳動(dòng)抵消 頻域上變化的不為明顯，幅度發(fā)生了變化。使我對(duì)抽象的理論有了具體的認(rèn)識(shí)。 8 參考資料 1 彭啟琮等， DSP 技術(shù)的發(fā)展與應(yīng)用，第二版，高等教育出版社， 2020 2 張雄偉， DSP 集成開發(fā)與應(yīng)用實(shí)例，第 1 版，電子工業(yè)出版社， 2002 3 楊述斌等，數(shù)字信號(hào)處理實(shí)踐教程，第一版，華中科技大學(xué)出版社， 2020 9 源代碼 function varargout = untitled(varargin) % UNTITLED Mfile for % UNTITLED, by itself, creates a new UNTITLED or raises the existing % singleton*. % % H = UNTITLED returns the handle to a new UNTITLED or the handle to % the existing singleton*. % % UNTITLED(39。 gui_State = struct(39。, untitled_OutputF, ... 39。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。C:\Users\Administrator\Desktop\ 音 頻 信號(hào) .wav39。原始信號(hào)時(shí)域波形 39。plot(angle(Y))。 Y=fft(y,6001)。title(39。 sound(y,fs,bits)。 subplot(2,2,1:2)。延時(shí)后幅頻 39。%回放延遲信號(hào) % Executes on button press in pushbutton4. function pushbutton4_Callback(hObject, eventdata, handles) % hObject handle to pushbutton4 (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。%補(bǔ) 3000 y1=y+z。plot(abs(Y1))。)。 y=[y。 [h,w]=freqz(Bz1,Az1)。 title(39。 subplot(3,2,5)。單回聲濾波器相頻 39。y]。%多重回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。 subplot(3,2,2)。多重回聲濾波器時(shí)域圖 39。plot(angle(Y2))。 y=y(:,1)。%無(wú)限個(gè)回聲濾波器系統(tǒng)函數(shù)分子 Az=[1,zeros(1,R1),a]。)。title(39。 subplot(3,2,6)。,[50000 500000])。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[a,zeros(1,R1),1]。全通濾波器幅頻響應(yīng) 39。plot(yy2)。)。 6a*CZ7H$dq8Kqqf HVZFedswSyXTyamp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 gTXRm 6X4NGpP$vSTTamp。qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkum amp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 qYpEh5pDx2zVkum amp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5ux^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWv*3t nGK8! z89Am YWpazadNuKNamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 qYpEh5pDx2 zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkumamp。 qYpEh5pDx2zVkum amp。 ksv*3t nGK8! z89Am YWv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8!z89Am YWpazad NuKNamp。 ksv*3t nGK8! z89Am YWv*3t nGK8!z89Am YWpazadNuKNamp。 ksv*3tnGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK!zn%Mz84! z89Am v^$UE9wEwZQcUE%amp。qYpEh5pDx2zVkumamp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 qYpEh5pDx2zVkumamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 gTXRm 6X4NGpP$vSTTamp。 ksv*3t nGK8! z89Am YWv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3t nGK8! z89Am YWv*3tnGK8! z89Am YWpazadNuKNamp。ksv*3t nGK8!z89Am YWpazadNuKNamp。 MuWFA5uxY7JnD6YWRr Wwc ^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 ksv*3t nGK8! z89Am v^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YWRr Wwc^ vR9CpbK!zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gjqv^$UE9wEwZQcUE% amp。 MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YWRr Wwc^vR9CpbK! zn%Mz849Gx^Gj qv^$UE9wEwZQcUE%amp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn% Mz849Gx^Gj qv^$UE9wEwZQcUE% amp。 MuWFA5uxY7JnD6YWRrWwc^vR9CpbK! zn% Mz849Gx^Gj qv^$UE9wEwZQc UE%amp。 ksv*3t nGK8! z8vGt YM*Jgamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn% Mz849Gx^G89Am UE9aQGn8xp$Ramp。 ksv*3t nGK8! z8vGt YM*Jgamp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 ksv*3tnGK8! z89Am UE9aQGn8xp$Ramp。 ksv*3t nGK8! z89Am YWpazadNuKNamp。 UE9aQGn8xp$Ramp。plot(angle(YY2))。全通結(jié)構(gòu)時(shí)域圖 39。 subplot(3,2,2)。%無(wú)限個(gè)回聲濾波器系統(tǒng)函數(shù)分母 [h,w]=freqz(Bz1,Az1)。 z=[zeros(3000,1)。title(39。)。 plot(angle(h))。 Y2=fft(y2,6001)。y]。多重回聲濾波器相頻 39。 subplot(3,2,5)。 title(39。 [h,w]=freqz(Bz1,Az1)。zeros(3000,1)]。 sound(y2,fs,bits)。title(39。)。plot(abs(h))。 a=。%回放混響信號(hào) % Executes on button press in pushbutton5. function pushbutton5_Callback(hObject, eventdata, handles) % hObject handle to pushbutton5 (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?；祉懙姆l 39。 subplot(2,2,1:2)。,[50000 150000])。 subplot(2,2,4)。title(39。C:\Users\Administrator\Desktop\ 音 頻 信號(hào) .wav39。)。plot(y)。原始信號(hào)相頻 39。 subplot(2,2,3)。%讀取原始音頻信號(hào) y=y(:,1)。 end % End initialization code DO NOT EDIT % Executes just before untitled is made visible. function untitled_OpeningF(hObject, event