基于matlab在聲音信號(hào)采集與處理中的應(yīng)用畢業(yè)論文-在線瀏覽
2024-08-04 15:46本頁(yè)面
【正文】 速傅里葉變換,得到信號(hào)的頻譜圖。我的電腦自帶麥克風(fēng)聲卡濾波采樣進(jìn)行A/D轉(zhuǎn)換得到wave格式文件（其中從聲卡到A/D轉(zhuǎn)換是Windows自帶的錄音機(jī)）。噪聲信號(hào)通常為隨機(jī)序列,在本設(shè)計(jì)中用隨機(jī)函數(shù)干擾信號(hào)構(gòu)建命令函數(shù)為Noise=*randn(n,2)。再對(duì)噪音信號(hào)進(jìn)行頻譜變換得到其頻譜圖,從圖中可以看出干擾信號(hào)，在4000Hz和 6000Hz頻點(diǎn)處有一高峰 ,其中 5500Hz 正是本設(shè)計(jì)所要利用的。在濾波器的選取時(shí), 由于設(shè)計(jì)方法的側(cè)重點(diǎn)不同, 作出比較是困難的。用雙線性變換法設(shè)計(jì)低通濾波器。在 MATLAB中 ,可以利用函數(shù) butterworth,設(shè)計(jì)FIR濾波器,利用 MATLAB 中的函數(shù)freqz 畫出各濾波器的頻率響應(yīng)。在一個(gè)窗口同時(shí)畫出濾波前后的波形及頻譜。利用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行分析和處理，采集語(yǔ)音信號(hào)后，利用MATLAB軟件平臺(tái)進(jìn)行頻譜分析；并對(duì)所采集的語(yǔ)音信號(hào)加入干擾噪聲，對(duì)加入噪聲的信號(hào)進(jìn)行頻譜分析，設(shè)計(jì)合適的濾波器濾除噪聲，恢復(fù)原信號(hào)。C:\Program Files\MATLAB704\work\39。l=length(y)。 %采樣時(shí)間%現(xiàn)在y：l*2雙聲道，f=,b=16%需要20e3Hz的采樣頻率,得到新數(shù)據(jù)y2f2=20e3。 B=l/l2。 %重新采樣endsubplot(211)plot(y)。wavwrite(y2,20e3,39。)。 %聲音信號(hào)采樣頻率為22050x1=wavread(39。)。 %播放聲音信號(hào)y1=fft(x1,1024)。 %將0到511，步長(zhǎng)為1的序列的值與fs相乘并除以1024的值，賦值給ffigure(1) %創(chuàng)建圖形窗1plot(x1) %做原始聲音信號(hào)的時(shí)域圖形title(39。)。time n39。 %x軸的名字是“time n”ylabel(39。)。頻率響應(yīng)圖39。 %創(chuàng)建兩行一列繪圖區(qū)間的第1個(gè)繪圖區(qū)間 plot(abs(y1(1:512))) %做原始聲音信號(hào)的FFT頻譜圖title(39。)subplot(2,1,2)。 %abs是絕對(duì)值，plot是直角坐標(biāo)下線性刻度曲線title(39。)xlabel(39。)。幅值39。 圖1 原始聲音信號(hào)圖2 聲音信號(hào)頻率響應(yīng)圖 圖3 原始聲音信號(hào)fft與信號(hào)頻譜（2）. 給原始的聲音信號(hào)加上一個(gè)隨機(jī)函數(shù)噪聲。程序：fs=22050。C:\Program Files\MATLAB704\work\39。 %讀取聲音信號(hào)的數(shù)據(jù)，賦給變量x1f=fs*(0:511)/1024。 %將0到x1的長(zhǎng)度減1后的值除以fs的值，且步長(zhǎng)為1/fs的值，的序列的值，賦予tn=length(x1)。 %所加的噪聲是隨機(jī)函數(shù)噪聲x2=x1+Noise。 %播放聲音信號(hào)y1=fft(x1,1024)。 %對(duì)信號(hào)y2做1024點(diǎn)FFT變換figure(1)。 %做加噪后的信號(hào)時(shí)域圖形（繪出t對(duì)x2的線性圖）title(39。)。time n39。 %x軸的名字是“time n”ylabel(39。)。 % 創(chuàng)建兩行一列繪圖區(qū)間的第1個(gè)繪圖區(qū)間plot(f,abs(y1(1:512)))。原始語(yǔ)音信號(hào)頻譜39。xlabel(39。)。fuzhi39。 % y軸的名字是“fuzhi”subplot(2,1
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1