【正文】 UIDATA) [y,fs,bits]=wavread(39。 Y=fft(y,6001)。采樣信號時域波形 39。title(39。plot(angle(Y))。 28 sound(y,fs,bits)。 y=y(:,1)。 subplot(2,2,1:2)。)。延時后幅頻 39。title(39。%回放延遲信號 % 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。 z=[zeros(3000,1)。%補 3000 y1=y+z。title(39。plot(abs(Y1))。 subplot(2,2,4)。)。,[50000 150000])。 y=[y。%濾波器階數(shù)設置，其值越高，回聲越明顯 Bz1=[1,zeros(1,R1),a]。 [h,w]=freqz(Bz1,Az1)。單回聲濾波器幅頻響應 39。 title(39。plot(y2)。 subplot(3,2,5)。)。單回聲濾波器相頻 39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。y]。 N=5。%多重回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。plot(abs(h))。 subplot(3,2,2)。)。多重回聲濾波器時域圖 39。title(39。plot(angle(Y2))。 sound(y2,fs,bits)。 y=y(:,1)。zeros(3000,1)]。%無限個回聲濾波器系統(tǒng)函數(shù)分子 Az=[1,zeros(1,R1),a]。 subplot(3,2,1)。)。無限個回聲濾波器相頻響應 39。title(39。plot(abs(Y2))。 subplot(3,2,6)。)。,[50000 150000])。 y=[y。%濾波器階數(shù)設置，其值越高，回聲越明顯 Bz1=[a,zeros(1,R1),1]。 YY2=fft(yy2,6001)。全通濾波器幅頻響應 39。 title(39。plot(yy2)。 subplot(3,2,5)。)。全通結(jié)構(gòu)相頻 39。依其延遲時間的不同，可分別產(chǎn)生合唱、鑲邊、回音等效果 。 成績評定依據(jù) ： 最終評定成績（以優(yōu)、良、中、 及格、不及格評定） 指導教師簽字： __________________ 年 月 日 。如果延遲時間在50ms以上時，其延遲音就清晰可辨，此時的處理效果才是回音。 sound(yy2,fs,bits)。plot(angle(YY2))。title(39。全通結(jié)構(gòu)時域圖 39。)。 subplot(3,2,2)。plot(abs(h))。%無限個回聲濾波器系統(tǒng)函數(shù)分母 [h,w]=freqz(Bz1,Az1)。 a=。 z=[zeros(3000,1)。 % 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。title(39。無限個回聲濾波器信號幅頻 39。)。 subplot(3,2,3:4)。 plot(angle(h))。 title(39。 Y2=fft(y2,6001)。 R=5000。y]。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。多重回聲濾波器相頻 39。)。 subplot(3,2,5)。plot(y2)。 title(39。多重回聲濾波器幅頻響應 39。 [h,w]=freqz(Bz1,Az1)。%濾波器階數(shù)設置，其值越高，回聲越明顯 Bz1=[1,zeros(1,N*R1),a^N]。zeros(3000,1)]。 y=y(:,1)。 sound(y2,fs,bits)。plot(angle(Y2))。title(39。單回聲濾波器時域圖 39。)。 subplot(3,2,2)。plot(abs(h))。%單回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。 a=。 z=[zeros(3000,1)。%回放混響信號 % 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。title(39?；祉懙姆l 39。)。 29 subplot(2,2,1:2)。 y=[y。,[50000 150000])。)。 subplot(2,2,4)。plot(abs(Z))。title(39。y]。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。采樣信號相頻 39。)。 subplot(2,2,3)。plot(y)。,[50000 150000])。原始信號相頻 39。)。 subplot(2,2,3)。plot(y)。%讀取原始音頻信號 y=y(:,1)。我們還能從設計中檢驗我們 所學的理論知識到底有多少，鞏固我們已經(jīng)學會的，不斷學習我們所遺漏的新知識，把這門課學的扎實。 經(jīng)過此次課程設計，使我對數(shù)字信號處理的理論知識有了更深了的了解，更近一步的理解了用 MATLAB語言對信號進行采樣，濾波等處理。在做設計的過程中總會出現(xiàn)各種問題，在這種情況下我們都會努力尋求最佳路徑解決問題，無形間提高了我們的動手，動腦能力，并且同學之間還能相互探討問題，研究解決方案，增進大家的團隊意識。 頻譜上也是濾波前后兩邊跳動變換較大，幅度變化較大。 4 結(jié)果分析與體會 運行結(jié)果： 19 對信號采樣后，信號時域波形出現(xiàn)了離散化的分布，不像原始圖像那樣密集分布，頻譜圖也是同樣的變換，幅度和頻率都是根據(jù)采樣區(qū)間變化的 20 從時域圖出延時是讓原來的波形向右移動 3000，頻域的幅度發(fā)生了變化。title(39。全通結(jié)構(gòu)幅頻 39。)。 subplot(3,2,3:4)。 plot(angle(h))。 18 title(39。 yy2=filter(Bz1,Az1,y)。 R=5000。y]。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。無限個回聲濾波器信 號相頻 39。)。 subplot(3,2,5)。plot(y2)。 title(39。無限個回聲濾波器幅頻響應 39。 [h,w]=freqz(Bz,Az)。%濾波器階數(shù)設置，其值越高，回聲越明顯 17 Bz=[0,zeros(1,R1),1]。 y=[y。,[50000 150000])。)。 subplot(3,2,6)。plot(abs(Y2))。title(39。多重回聲濾波器相頻響應 39。)。 subplot(3,2,1)。%多重回聲濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。 a=。 16 z=[zeros(3000,1)。 多重回聲濾波器程序設計 [y,fs,bits]=wavread(39。title(39。單回聲濾波器幅頻 39。)。 subplot(3,2,3:4)。 plot(angle(h))。 title(39。%單回聲濾波器系統(tǒng)函數(shù) Y2=fft(y2,6001)。 R=5000。y]。C:\Users\Administrator\Desktop\ 音頻信號 .wav39?；祉懙南囝l 39。)。 subplot(2,2,3)。plot(y1)。zeros(3000,1)]。 y=y(:,1)。 sound(y,fs,bits)。plot(angle(Z))。title(39。延時后時域波形 39。%延遲 3000 Z=fft(z,6001)。,[50000 150000])。)。 subplot(2,2,4)。plot(abs(Y))。title(39。%只取原信號的一部分，即時間區(qū)間 [50000 150000]