【導(dǎo)讀】①較全面了解常用的數(shù)據(jù)分析與處理原理及方法；②能夠運(yùn)用相關(guān)軟件進(jìn)行模擬分析；③掌握基本的文獻(xiàn)檢索和文獻(xiàn)閱讀的方法；④提高正確地撰寫(xiě)論文的基本能力。，數(shù)字信號(hào)處理原理與實(shí)現(xiàn)（第二版），電子工業(yè)出版社，2020；2．張雄偉，DSP集成開(kāi)發(fā)與應(yīng)用實(shí)例，第1版，電子工業(yè)出版社，2020；

　　

【正文】 2,2,4)。plot(angle(Y1))。title(39?；祉懙南囝l 39。)。 sound(y,fs,bits)。%回放混響信號(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。C:\Users\Administrator\Desktop\ 音頻信號(hào) .wav39。,[50000 150000])。 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è)計(jì)的濾波器頻譜 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。單回聲濾波器時(shí)域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。單回聲濾波器幅頻 39。)。 30 subplot(3,2,6)。plot(angle(Y2))。title(39。單回聲濾波器相頻 39。)。 sound(y2,fs,bits)。 % Executes on button press in pushbutton6. function pushbutton6_Callback(hObject, eventdata, handles) % hObject handle to pushbutton6 (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。C:\Users\Administrator\Desktop\ 音頻信號(hào) .wav39。,[50000 150000])。 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。多重回聲濾波器時(shí)域圖 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)。 31 % 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。C:\Users\Administrator\Desktop\ 音頻信號(hào) .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz=[0,zeros(1,R1),1]。%無(wú)限個(gè)回聲濾波器系統(tǒng)函數(shù)分子 Az=[1,zeros(1,R1),a]。%無(wú)限個(gè)回聲濾波器系統(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。無(wú)限個(gè)回聲濾波器幅頻響應(yīng) 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。無(wú)限個(gè)回聲濾波器相頻響應(yīng) 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。無(wú)限個(gè)回聲濾波器時(shí)域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。無(wú)限個(gè)回聲濾波器信號(hào)幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。無(wú)限個(gè)回聲濾波器信號(hào)相頻 39。)。 sound(y2,fs,bits)。 % Executes on button press in pushbutton8. function pushbutton8_Callback(hObject, eventdata, handles) 32 % hObject handle to pushbutton8 (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。C:\Users\Administrator\Desktop\ 音頻信號(hào) .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯 Bz1=[a,zeros(1,R1),1]。%無(wú)限個(gè)回聲濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。%無(wú)限個(gè)回聲濾波器系統(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)時(shí)域圖 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)。 33 本科生課程設(shè)計(jì)成績(jī)?cè)u(píng)定表 姓 名 陳哲 性 別 男 專業(yè)班級(jí) 通信 GJ1001 題 目 ： 基于 MATLAB 的聲音延遲與混響 答辯或質(zhì)疑記錄： 1. 什么是延時(shí)？ 答： 延時(shí)就是將音源延遲一段時(shí)間后，再欲播放的效果處理。依其延遲時(shí)間的不同，可分別產(chǎn)生合唱、鑲邊、回音等效果 。 2. 如何產(chǎn)生延時(shí)？ 答： 當(dāng)延遲時(shí)間在 3~35ms之間時(shí)人耳感覺(jué)不到滯后音的存在，并且他與原音源疊加后，會(huì)因其相位干涉而產(chǎn)生 梳狀濾波 效應(yīng)，這就是鑲邊效果。如果延遲時(shí)間在50ms以上時(shí)，其延遲音就清晰可辨，此時(shí)的處理效果才是回音?；匾籼幚硪话愣际怯糜诋a(chǎn)生簡(jiǎn)單的混響效果。 成績(jī)?cè)u(píng)定依據(jù) ： 最終評(píng)定成績(jī)（以優(yōu)、良、中、 及格、不及格評(píng)定） 指導(dǎo)教師簽字： __________________ 年 月 日