最新數(shù)字信號處理課程設(shè)計報告-資料下載頁
2025-08-01 19:29本頁面
【正文】 繪制截取后的語音信號title(39。原始語音信號39。)。xlabel(39。時間(單位：s)39。)。ylabel(39。幅度39。)。axis tightsubplot(2,2,3)。plot(t,y) % 繪制截取后的語音信號title(39。加入單頻干擾后的語音信號39。)。xlabel(39。時間(單位：s)39。)。ylabel(39。幅度39。)。axis tightX=abs(fft(x))。 Y=abs(fft(y))。 % 對原始信號和加噪信號進(jìn)行fft變換,取幅度譜X=X(1:N/2)。 Y=Y(1:N/2)。 % 截取前半部分deltaf=fs/N。 % 計算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。 % 計算頻譜頻率范圍subplot(2,2,2)。plot(f,abs(X)) % 繪制頻譜圖title(39。語音信號幅度譜圖39。)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。) 。axis tightsubplot(2,2,4)。plot(f,abs(Y)) % 繪制頻譜圖title(39。加入干擾后的語音信號幅度譜圖39。)。xlabel(39。頻率（單位：Hz）39。)。ylabel(39。幅度譜39。) 。axis tightfp=2000。fc=2150。 fs=8000 %定義通帶和阻帶截止頻率Rp=1。As=50。 % 定義通帶波紋和阻帶衰減wp=fp/fs*2*pi。ws=fc/fs*2*pi。 %計算對應(yīng)的數(shù)字頻率T=2。 %定義采樣間隔Omegap=2/T*tan(wp/2)。Omegas=2/T*tan(ws/2)。 %截止頻率預(yù)畸變[cs,ds]=afd_chb1(Omegap,Omegas,Rp,As) % 計算濾波器系統(tǒng)函數(shù)分子分母系數(shù)[b,a]=bilinear(cs,ds,1/T) % 雙線性變換得到數(shù)字濾波器系統(tǒng)函數(shù)分子分母系數(shù)[db,mag,pha,grd,w]=freqz_m(b,a)。 %驗證濾波器是否達(dá)到指定性能delta=[1,zeros(1,99)]。ha=filter(b,a,delta)。Subplot(2,2,1)。plot(w/pi,db,39。LineWidth39。,2)。axis([0,1,80,10])。grid on。% 繪制相對幅度響應(yīng)曲線圖，限制橫縱坐標(biāo)范圍，添加網(wǎng)格X_l=[0,0,wp/pi,ws/pi。1,1,wp/pi,ws/pi]。Y_l=[As,Rp,80,80。As,Rp,10,10]。 % 在wp，ws，Rp，As處畫線以更直觀判斷設(shè)計是否達(dá)標(biāo)，每列參數(shù)是每個線條的端點(diǎn)坐標(biāo)line(X_l,Y_l,39。Color39。,39。r39。,39。LineWidth39。,2,39。LineStyle39。,39。39。) % 添加線條，紅色，線寬為2title(39。濾波器幅度響應(yīng)圖39。)。xlabel(39。w/pi39。)。ylabel(39。dB39。)。subplot(2,2,2)。plot(w/pi,mag,39。LineWidth39。,2)。axis([0,1,0,])。grid on。title(39。濾波器幅度響應(yīng)圖39。)。xlabel(39。w/pi39。)。ylabel(39。幅度mag39。)。% 繪制絕對幅度響應(yīng)曲線圖subplot(2,2,3)。plot(w/pi,pha)。axis tight。grid on。title(39。濾波器相位響應(yīng)圖39。)。xlabel(39。w/pi39。)。ylabel(39。相位pha39。)。% 繪制相位特性曲線圖subplot(2,2,4)。plot(ha,39。LineWidth39。,2)。 axis([0,100,])。grid on。title(39。濾波器脈沖響應(yīng)圖39。)。xlabel(39。n39。)。ylabel(39。ha39。)。% 繪制脈沖響應(yīng)圖y_fil=filter(b,a,y)。 % 對信號y進(jìn)行濾波處理，其中b,a為上面設(shè)計好的濾波器參數(shù)Y_fil=abs(fft(y_fil))。Y_fil=Y_fil(1:N/2)。 % 計算頻譜取前一半sound (y_fil,fs,bits)figure。subplot(3,2,1)。plot(t,x)。axis tight。grid on。title(39。原始語音信號時間x39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制原始錄制的語音信號時間波形圖subplot(3,2,3)。plot(t,y)。axis tight。grid on。title(39。加干擾語音信號時間y39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制加噪語音信號時間波形圖subplot(3,2,5)。plot(t,y_fil)。axis tight。grid on。title(39。濾波后語音信號時間y_fil39。)。xlabel(39。時間t39。)。ylabel(39。幅度39。)。% 繪制去噪語音信號時間波形圖subplot(3,2,2)。plot(f,abs(X))。axis([0 3000 0 60])。grid on。title(39。原始語音信號幅度譜X39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制原始錄制的語音信號頻譜圖subplot(3,2,4)。plot(f,abs(Y))。axis([0 3000 0 60])。grid on。title(39。加干擾語音信號幅度譜Y39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制加噪語音信號頻譜圖subplot(3,2,6)。plot(f,abs(Y_fil))。axis([0 3000 0 60])。grid on。title(39。濾波后語音信號幅度譜Y_fil39。)。xlabel(39。頻率f39。)。ylabel(39。幅度39。)。% 繪制去噪語音信號頻譜圖[C,B,A]=dir2par(b,a)
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1