【文章內(nèi)容簡介】 。混響的時域波形39。)。subplot(2,2,3)。plot(abs(Y1))。title(39?；祉懙姆l39。)。subplot(2,2,4)。plot(angle(Y1))。title(39?；祉懙南囝l39。)。sound(y,fs,bits)。%回放混響信號 單回聲濾波器程序設(shè)計[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=[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è)計的濾波器頻譜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)。 多重回聲濾波器程序設(shè)計[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)。 無限個回聲濾波器程序設(shè)計[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)。 全通結(jié)構(gòu)濾波器程序設(shè)計[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)。 4 結(jié)果分析與體會運行結(jié)果：對信號采樣后，信號時域波形出現(xiàn)了離散化的分布，不像原始圖像那樣密集分布，頻譜圖也是同樣的變換，幅度和頻率都是根據(jù)采樣區(qū)間變化的從時域圖出延時是讓原來的波形向右移動3000，頻域的幅度發(fā)生了變化。混響：在時域圖看前面幾乎沒變化，后面的幅度有一定變化。在頻譜圖幅度變化很大，中間的上下相互抵消，兩邊的跳動很大。時域上形狀大致相同，幅度有微小變化。頻譜上也是濾波前后兩邊跳動變換較大，幅度變化較大。時域上濾波前的幅度比濾波后的數(shù)值大，形狀上濾波后是連續(xù)的，濾波前是間斷的頻譜上變化兩邊上下跳的動變化，中間變化抵消時域上波形變得稀疏了，幅度相對變化較大。頻譜上變化兩邊上下跳的動變化，中間變化抵消全通結(jié)構(gòu)時域大致一樣，濾波前比濾波后幅度有變化6000以后上下跳動抵消頻域上變化的不為明顯，幅度發(fā)生了變化。5 設(shè)計心得體會 課程設(shè)計是培養(yǎng)我們綜合運用所學(xué)知識,發(fā)現(xiàn),提出,分析和解決實際問題,鍛煉實踐能力的重要環(huán)節(jié),是對我們實際工作能力的具體訓(xùn)練和考察過程，通過課程設(shè)計我們能夠比較系統(tǒng)的了解理論知識，把理論和實踐相結(jié)合，并且用到生活當(dāng)中。在做設(shè)計的過程中總會出現(xiàn)各種問題，在這種情況下我們都會努力尋求最佳路徑解決問題，無形間提高了我們的動手，動腦能力，并且同學(xué)之間還能相互探討問題，研究解決方案，增進大家的團隊意識。 此次課程設(shè)計綜合了數(shù)字信號處理的理論知識和實踐經(jīng)驗，經(jīng)過這段時間查閱的相關(guān)資料并結(jié)合以前學(xué)過的數(shù)字信號處理課程和實驗，我如期完成了設(shè)計任務(wù)要求，并加深了對知識的理解，提高了自己的實踐動手能力。拿到課程設(shè)計任務(wù)的時候，不知道從哪兒入手，選定題目后，上網(wǎng)查閱了相關(guān)資源，去圖書館搜集資料，確定了設(shè)計的基本方案，步驟，開始用軟件進行實現(xiàn)，實現(xiàn)過程中出現(xiàn)了很多問題，上網(wǎng)查閱資料，在專業(yè)的論壇里面去看，下載MATLAB的視頻，慢慢的開始有了進一步的理解。在以前的實驗中沒有接觸過GUI界面的問題，根本不知道從哪入