最新語音信號處理及matlab實現(xiàn)課程設(shè)計-資料下載頁
2025-06-29 03:30本頁面
【正文】 E:\xinhao139。)。t=(0:length(x0)1)/16000。f=fs*(0:10239)/20480。 y1=fft(x0,20480)。%畫語音信號時域圖figure(1)subplot(1,1,1)。 plot(t,x0)grid on。 axis tight。 title(39。原始語音信號時域39。)。 xlabel(39。time(s)39。)。 ylabel(39。幅度39。)。 %畫語音信號頻域圖 figure(2) subplot(1,1,1)。 plot(f,abs(y1(1:10240)))。 grid on。 axis([0,8000,0,20])。 title(39。原始語音信號頻譜39。)。 xlabel(39。Hz39。)。 ylabel(39。幅度39。)。 %加噪后語音信號變化 d=[*cos(2*pi*4000*t)]39。 %加噪 x1=x0+d。 y2=fft(x1,20480)。%加噪后語音信號的時域圖 figure(3) subplot(1,1,1)。 plot(t,x0) grid on。ylim([,])。 axis tight。 title(39。加噪后的信號時域39。)。 xlabel(39。time(s)39。)。 ylabel(39。幅度39。)。 %加噪后語音信號的頻域圖 figure(4) subplot(1,1,1)。 plot(f,abs(y2(1:10240)))。 grid on。 axis([0,8000,0,20])。 title(39。加噪語音信號頻譜39。)。 xlabel(39。Hz39。)。 ylabel(39。幅度39。)。 %設(shè)計帶阻濾波器wlp=2850*pi*2/fs。wls=3400*pi*2/fs。wup=4900*pi*2/fs。wus=4600*pi*2/fs。B=(wuswls)。M=ceil(12*pi/B)1。wp=[(wls+wlp)/2/pi,(wus+wup)/2/pi]。[bz,az]=fir1(M,wp,39。stop39。,blackman(M+2))。hn=fir1(M,wp,39。stop39。,blackman(M+2))。%帶阻濾波器窗函數(shù)圖像figure(5)subplot(1,1,1)。freqz(hn)。grid on。axis tight。title(39。窗函數(shù)圖像39。)。xlabel(39。Hz39。)。ylabel(39。幅度39。)。%畫濾波后的圖figure(6) freqz(bz,az)。 grid on。 axis tight。y3=filter(bz,az,x1)。X=fft(y3,20480)。%濾波后時域圖subplot(2,1,1)。plot(y3)。 grid on。ylim([ ,])。title(39。濾波后信號時域波形39。)。 xlabel(39。time(s)39。)。 ylabel(39。幅度39。)。 %濾波后頻譜subplot(2,1,2)。plot(abs(X(1:10240)))。 title(39。濾波后信號頻譜39。)。grid on。ylim([0,20])。山東建筑大學(xué)信電學(xué)院課程設(shè)計說明書第 20 頁xlabel(39。Hz39。)。ylabel(39。幅度39。)。 d1=[*cos(2*pi*200*t)]39。 %加噪 x4=x0+d1。 y4=fft(x4,20480)。%加噪后語音信號的時域圖 figure(7) subplot(1,1,1)。 plot(t,x4) grid on。 axis tight。 title(39。加噪后的信號時域39。)。 xlabel(39。time(s)39。)。 ylabel(39。幅度39。)。 %加噪后語音信號的頻域圖 figure(8) subplot(1,1,1)。 plot(f,abs(y4(1:10240)))。 grid on。 axis([0,8000,0,20])。 title(39。加噪語音信號頻譜39。)。 xlabel(39。Hz39。)。 ylabel(39。幅度39。)。 %加高通濾波 wp1=100*pi*2/fs。 ws1=60*pi*2/fs。 Rp1=1。 As1=100。 [N1,wso]=ellipord(wp1,ws1,Rp1,As1)。 [B1,A1]=ellip(N1,Rp1,As1,wso,39。high39。)。 figure(9) freqz(B1,A1)。 grid on。 ylim([0,1])。y4=filter(B1,A1,x4)。X1=fft(y4,20480)。%濾波后時域圖subplot(2,1,1)。plot(y4)。 grid on。%axis tight。ylim([,])。title(39。濾波后信號時域波形39。)。 xlabel(39。time(s)39。)。 ylabel(39。幅度39。)。 %濾波后頻譜subplot(2,1,2)。plot(abs(X1(1:10240)))。 title(39。濾波后信號頻譜39。)。grid on。%ylim([0,20])。 axis tight。xlabel(39。Hz39。)。ylabel(39。幅度39。)。
公司管理相關(guān)推薦
鄂ICP備17016276號-1