基于matlab的語音信號采集與處理課程設(shè)計(jì)報(bào)告(編輯修改稿)
2024-12-16 03:34 本頁面
【文章內(nèi)容簡介】 %繪制原始語音信號的頻率響應(yīng)圖 title(39。頻率響應(yīng)圖 39。) figure(3) subplot(2,1,1)。 %創(chuàng)建 兩行一列繪圖區(qū)間的第 1個繪圖區(qū)間 plot(abs(y1(1:512))) %做原始語音信號的 FFT頻譜圖 4 title(39。原始語音信號 FFT頻譜 39。) subplot(2,1,2)。 plot(f,abs(y1(1:512)))。 %abs是 絕對值 ， plot是 直角坐標(biāo)下線性刻度曲線 title(39。原始語音信號頻譜 39。) xlabel(39。Hz39。)。 ylabel(39。幅值 39。)。 時域波形和頻譜 ： 圖 1 原始語音信號 5 圖 2 語音信號頻率響應(yīng)圖 圖 3 原始語音信號 FFT 與信號頻譜 2. 給原始的 語音 信號加上一個高頻余弦 噪聲 ，頻率為 5500hz。 對 加噪后的 語音進(jìn)行分析，并畫出其 信號時域和頻譜圖 。 6 程序： fs=22050。 x1=wavread(39。C:\Documents and Settings\Administrator\桌面 \聲音文件 .wav39。)。 %讀取語音信號的數(shù)據(jù)，賦給變量 x1 f=fs*(0:511)/1024。 %將 0到 511，步長為 1的序列的值與 fs相乘并除以 1024的值，賦值給 f t=0:1/fs:(length(x1)1)/fs。 %將 0 到 x1 的長度減 1 后的值除以 fs 的值，且步長為1/fs的值，的序列的值，賦予 t Au=。 %噪聲幅值 d=[Au*sin(2*pi*5500*t)]39。 %所加的噪聲是正弦信號 x2=x1+d。 %將正弦信號噪聲加在語音信號上 sound(x2,22050)。 %播放語音信號 y1=fft(x1,1024)。 %對信號 y1做 1024點(diǎn) FFT變換 y2=fft(x2,1024)。 %對信號 y2做 1024點(diǎn) FFT變換 figure(1)。 %創(chuàng)建圖形窗 1 plot(t,x2)。 %做加噪后的信號時域圖形（繪出 t對 x2的線性圖） title(39。加噪后的信號 39。)。 xlabel(39。time n39。)。 %x軸的名字是“ time n” ylabel(39。fuzhi n39。)。 %y軸的名字是“ fuzhi n” figure(2) %創(chuàng)建圖形窗 2 subplot(2,1,1)。 % 創(chuàng)建兩行一列繪圖區(qū)間的第 1個繪圖區(qū)間 plot(f,abs(y1(1:512)))。 %做原始語音信號的頻譜圖 title(39。原始語音信號頻譜 39。)。 xlabel(39。Hz39。)。 %x軸的名字是“ time n” ylabel(39。fuzhi39。)。 % y軸的名字是“ fuzhi” subplot(2,1,2)。 %創(chuàng)建兩行一列繪圖區(qū)間的第 2個繪圖區(qū)間 plot(f,abs(y2(1:512)))。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1