基于matlab的語音信號(hào)的分析與處理(編輯修改稿)
2025-07-15 17:07 本頁面
【文章內(nèi)容簡介】 figure(2)freqz(x1) %繪制原始語音信號(hào)的頻率響應(yīng)圖title(39。頻率響應(yīng)圖39。)plot(f,abs(y1(1:512)))。figure(3) title(39。原始語音信號(hào)FFT頻譜39。)xlabel(39。Hz39。)。ylabel(39。幅值39。)；圖31 語音信號(hào)的時(shí)域波形圖32 語音信號(hào)的頻譜圖%對(duì)噪音信號(hào)進(jìn)行時(shí)域和頻域分析。fs=11025。 x1=wavread(39。39。)。 %讀取語音信號(hào)的數(shù)據(jù)，賦給變量x1f=fs*(0:511)/1024。 %將0到511，步長為1的序列的值與fs相乘并除以1024的值，賦值給ft=0:1/fs:(length(x1)1)/fs。 %將0到x1的長度減1后的值除以fs的值，且步長為1/fs的值，的序列的值，賦予tAu=。 %噪聲幅值d=[Au*sin(2*pi*5500*t)]39。 %所加的噪聲是正弦信號(hào)x2=x1+d。 %將正弦信號(hào)噪聲加在語音信號(hào)上sound(x2,22050)。 %播放語音信號(hào)y1=fft(x1,1024)。 %對(duì)信號(hào)y1做1024點(diǎn)FFT變換y2=fft(x2,1024)。 %對(duì)信號(hào)y2做1024點(diǎn)FFT變換figure(1)。 %創(chuàng)建圖形窗1plot(t,x2)。 %做加噪后的信號(hào)時(shí)域圖形（繪出t對(duì)x2的線性圖）title(39。加噪后的信號(hào)39。)。xlabel(39。時(shí)間39。)。 %x軸的名字是“時(shí)間”ylabel(39。幅值39。)。 %y軸的名字是“幅值”figure(2) %創(chuàng)建圖形窗2subplot(2,1,1)。 % 創(chuàng)建兩行一列繪圖區(qū)間的第1個(gè)繪圖區(qū)間plot(f,abs(y1(1:512)))。 %做原始語音信號(hào)的頻譜圖title(39。原始語音信號(hào)頻譜39。)。xlabel(39。Hz39。)。 %x軸的名字是“HZ”ylabel(39。幅值39。)。 % y軸的名字是“幅值”subplot(2,1,2)。 %創(chuàng)建兩行一列繪圖區(qū)間的第2個(gè)繪圖區(qū)間plot(f,abs(y2(1:512)))。 %做加噪后的語音信號(hào)的頻譜圖（abs是絕對(duì)值，plot是直角坐標(biāo)下線性刻度曲線）title(39。加噪后的信號(hào)頻譜39。)。xlabel(39。Hz39。)。 %x軸的名字是“HZ”ylabel(39。幅值39。)。 % y軸的名字是“幅值”圖33 加噪后的語音信號(hào)時(shí)域波形圖34原始語音信號(hào)與加噪后的語音信號(hào)的頻譜比較圖%butter設(shè)計(jì)模擬低通濾波器再用雙線性變換法設(shè)計(jì)IIR數(shù)字低通濾波器fp=1000。fc=1200。FS=8000。rp=1。rs=100。 wpz=2*pi*1000/FS。wsz=2*pi*1200/FS。wp=2*tan(wpz/2)*FS。ws=2*tan(wsz/2)*FS。 %預(yù)畸校正轉(zhuǎn)換指標(biāo)[N,wc]=buttord(wp,ws,rp,rs,39。s39。)。[B,A]=butter(N,wc,39。s39。)。[Bz,Az]=bilinear(B,A,FS)。w=0:*pi:pi。[h,w]=freqz(Bz,Az,w)。figure(1)。plot(w/pi,20*log(abs(h)),39。k39。)。axis([0,1,800,100])。xlabel(39。w/pi39。)。ylabel(39。幅度/dB39。)。grid。title(39。butter數(shù)字低通濾波器39。)。[y,fn,nbits]=wavread(39。39。)。 Y=fft(y)。y1=filter(Bz,Az,y)。Y1=fft(y1)。n=0:length(y)1。figure(2)。subplot(221)。plot(y)。title(39。未濾波語音波形39。)。subplot(222)。plot(y1)。title(39。濾波后語音波形39。)。subplot(223)。plot(n,Y)。title(39。未濾波語音頻譜39。)。subplot(224)。plot(n,Y1)。title(39。濾波后語音頻譜39。)。sound(y1,fn,nbits)。 %濾波后語音回放圖35 設(shè)計(jì)butter數(shù)字低通濾波器圖36濾波前與濾波后語音信號(hào)波形、頻譜的比較%用凱瑟窗設(shè)計(jì)FIR低通濾波器fp=1000。fs=1200。rs=100。Fs=8000。 %kaise
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1