【正文】 1,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)。 % 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\39。,[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。多重回聲濾波器時域圖 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)。 % 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\39。,[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]。%無限個回聲濾波器系統(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)。 % Executes on button press in pushbutton8.function pushbutton8_Callback(hObject, eventdata, handles)% 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\39。,[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]。%無限個回聲濾波器系統(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)。 本科生課程設(shè)計成績評定表姓 名陳哲性 別男專業(yè)班級通信GJ1001題 目：基于MATLAB的聲音延遲與混響 答辯或質(zhì)疑記錄：1. 什么是延時？答：延時就是將音源延遲一段時間后，再欲播放的效果處理。依其延遲時間的不同，可分別產(chǎn)生合唱、鑲邊、回音等效果。2. 如何產(chǎn)生延時？答：當(dāng)延遲時間在3~35ms之間時人耳感覺不到滯后音的存在，并且他與原音源疊加后，會因其相位干涉而產(chǎn)生梳狀濾波效應(yīng)，這就是鑲邊效果。如果延遲時間在50ms以上時，其延遲音就清晰可辨，此時的處理效果才是回音?；匾籼幚硪话愣际怯糜诋a(chǎn)生簡單的混響效果。成績評定依據(jù)：最終評定成績（以優(yōu)、良、中、及格、不及格評定）指導(dǎo)教師簽字：__________________ 年 月