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