【文章內(nèi)容簡介】 置頻譜的間隔，分辨率 Subplot(212)。 plot([fs/2:derta_fs:fs/2derta_fs],abs(z))。%畫出延時后的信號頻譜圖 title(39。延時后的信號頻譜圖39。)。grid on。 運行結(jié)果： 0 2 4 6 8 10 12 14x 104101 個個個個個個個個個3 2 1 0 1 2 3x 1040500100015002022 個個個個個個個個個 圖 33 延時后語音信號的時域波形與頻譜圖結(jié)果分析：由運行結(jié)果可以看出，時域波形向右平移，頻譜無變化。 混響 調(diào)整原始采樣信號的長度，將其與延時信號相加，就得到混響。畫出混響后語音信號的時域波形與頻譜圖，并回放延時后的語音信號。武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20程序： %混響后信號 [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 z=[zeros(2022,1)。x]。 %對語音信號進行延時 x=[x。zeros(2022,1)]。 %使語音信號與延時后信號同等長度 y1=x+z。 %信號的混響 sound(y1,fs)。 %回放混響后的信號 figure(3)。subplot(211)。plot(y1)。 title(39?；祉懙臅r域圖39。)。 %混響時域圖 Y1=fft(y1)。 %做length(y1)點的FFT Y1=fftshift(Y1)。 %平移，是頻率中心為0 derta_fs = fs/length(y1)。 %設(shè)置頻譜的間隔，分辨率 Subplot(212)。 plot([fs/2:derta_fs:fs/2derta_fs],abs(Y1))。%畫出混響頻譜圖 title(39?；祉懙念l譜圖39。)。grid on。運行結(jié)果： 0 2 4 6 8 10 12 14x 104101 個個個個個個3 2 1 0 1 2 3x 1040100020223000 個個個個個個 圖 34 混響后語音信號的時域波形與頻譜圖武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20結(jié)果分析：由運行結(jié)果可以看出，原始信號的延時聲陣列非常密集、復(fù)雜，混響后時域波形音量增大了，各個脈沖變寬，語音信號的頻譜展寬了，幅值增大，回放的語音聲音變得融合，能聽到回音。 單回聲濾波器 用單回聲濾波器對原始語音信號濾波,畫出濾波后語音信號的時域波形與頻譜圖，并回濾波時后的語音信號。程序： %單回聲濾波器 [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 a=。 %a取小于等于1 z=[zeros(2022,1)。x]。 %對語音信號進行延時 x=[x。zeros(2022,1)]。 %使語音信號與延時后信號同等長度 y2=x+z*a。 %信號經(jīng)單回聲濾波 sound(y2,fs)。 %回放經(jīng)單回聲濾波器后的語音信號 figure。subplot(211)。plot(y2)。title(39。單回聲濾波的信號時域圖39。)。 Y2=fft(y2)。 %做length(y2)點的FFT Y2=fftshift(Y2)。 %平移，是頻率中心為0 derta_fs = fs/length(y2)。 %設(shè)置頻譜的間隔，分辨率 Subplot(212)。 plot([fs/2:derta_fs:fs/2derta_fs],abs(Y2))。%畫出單回聲信號頻譜圖 title(39。單回聲濾波的信號頻譜圖39。)。grid on。武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20運行結(jié)果： 0 2 4 6 8 10 12 14x 104101 個個個個個個個個個個個3 2 1 0 1 2 3x 1040500100015002022 個個個個個個個個個個個 圖 35 單回聲濾波器濾波后信號的時域波形與頻譜圖結(jié)果分析：由運行結(jié)果可以看出，時域波形的幅值變大，毛刺變多，頻譜的中心頻率兩邊的頻率分量展寬、幅值變大，回放的聲音立體感變強。 多回聲濾波器 用多回聲濾波器對原始語音信號濾波,畫出濾波后語音信號的時域波形與頻譜圖，并回濾波時后的語音信號。程序：%多回聲濾波器 [x,fs,bits]=wavread(39。D:\39。)。 x1=x(:,1)。 %變換為單聲道語音信號 delay_t=。 %回聲延遲時間 a=。 %衰減系數(shù) N=3。 %N階延時 Bz=[1。zeros(fs*delay_t,1)*N。a^N]。 Az=[1。zeros(fs*delay_t,1)。a]。武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20 y=filter(Bz, Az ,x1 )。 sound(y,fs,bits)。 %回放經(jīng)多回聲濾波器后的語音信號 figure。subplot(211)。plot(y)。title(39。多回聲濾波的信號時域圖39。)。 Y2=fft(y)。 %做length(y)點的FFT Y2=fftshift(Y2)。 %平移，是頻率中心為0 derta_fs = fs/length(y)。 %設(shè)置頻譜的間隔，分辨率 Subplot(212)。 plot([fs/2:derta_fs:fs/2derta_fs],abs(Y2))。 %畫出多回聲信號頻譜圖 title(39。多回聲濾波的信號頻譜圖39。)。grid on。運行結(jié)果： 0 2 4 6 8 10 12 14x 104101 個個個個個個個個個個個3 2 1 0 1 2 3x 1040500100015002022 個個個個個個個個個個個圖 36 單回聲濾波器濾波后信號的時域波形與頻譜圖結(jié)果分析：由運行結(jié)果可以看出，時域波形的幅值變大，毛刺變得更多，頻譜的中心頻率幅值變大，回放的聲音立體感變的更強。 無限回聲濾波器 用無限回聲濾波器對原始語音信號濾波,畫出濾波后語音信號的時域波形與武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20頻譜圖，并回濾波時后的語音信號。程序： %無限回聲濾波器 [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 a=。 %a取小于等于1 Bz=[0,0,0,0,0,0,0,0,0,0,1]。 %分子的系數(shù) Az=[1,0,0,0,0,0,0,0,0,0,a]。 %分母的系數(shù) yy1=filter(Bz,Az,x)。 %FIR濾波器進行濾波 sound(yy1,fs,bits)。 %回放經(jīng)無限回聲濾波器后的語音信號 figure(5)。subplot(211)。plot(yy1)。 %無限回聲信號時域波形 title(39。無限個回聲濾波的信號時域圖39。)。 YY1=fft(yy1)。 %經(jīng)無限回聲濾波器后的信號做length(yy1) 點的FFT YY1=fftshift(YY1)。 %平移，是頻率中心為0 derta_fs = fs/length(yy1)。 %設(shè)置頻譜的間隔，分辨率 Subplot(212)。plot([fs/2:derta_fs:fs/2derta_fs],abs(YY1))。%無限回聲頻譜圖 title(39。無限個回聲濾波的信號頻譜圖39。)。grid on。 運行結(jié)果：武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20 0 2 4 6 8 10 12 14x 104101 個個個個個個個個個個個個個3 2 1 0 1 2 3x 10401000202230004000 個個個個個個個個個個個個個圖 37 無限回聲濾波器濾波后信號的時域波形與頻譜圖結(jié)果分析：由運行結(jié)果可以看出，時域波形的幅值變大，頻譜的中心頻率不變，中心頻率兩邊的頻率分量變得密集、幅值變小，回放的聲音變得清晰，音量變大。 全通濾波器 用全通濾波器對原始語音信號濾波,畫出濾波后語音信號的時域波形與頻譜圖，并回濾波時后的語音信號。程序： %全通結(jié)構(gòu)的混響器 [x,fs,bits]=wavread(39。D:\39。)。 %讀取語音信號 x=x(:,1)。 %只取單聲道 a=。 %a取小于等于1 Bz1=[a,0,0,0,0,0,0,0,0,0,1]。 %分子的系數(shù) Az1=[1,0,0,0,0,0,0,0,0,0,a]。 %分母的系數(shù) yy2=filter(Bz1,Az1,x)。 %FIR濾波器進行濾波 sound(yy2,fs,bits)。 %回放經(jīng)全通結(jié)構(gòu)的混響器后的語音信號 figure(6)。 subplot(211)。plot(yy2)。 title(39。全通濾波信號的時域圖39。)武漢理工大學(xué)《信息群處理課程設(shè)計》課程設(shè)計說明書20 %全通結(jié)構(gòu)的混響器時域波形 YY2=fft(yy2)。 %經(jīng)全通濾波器后的信號做length(yy2)點的FFT YY2=fftshift(YY2)。 %平移，是頻率中心為0 derta_fs = fs/length(yy2)。 %設(shè)置頻譜的間隔，分辨率 Subplot(212)。 plot([fs/2:derta_fs:fs/2derta_fs],abs(YY2))。 %全通濾波信號的頻譜圖 title(39。全通濾波信號的頻譜圖39。)。grid on。運行結(jié)果： 0 2 4 6 8 10 12 14x 104101 個個個個個個個個個個3 2 1 0 1 2 3x 1040500100015002022 個個個個個個個個個個圖 38 無限回聲濾波器濾波后信號的時域波形與頻譜圖結(jié)果分析：由運行結(jié)果可以看出，時域波形和頻譜基本不變，回放的聲音跟