1.仿真實驗中將提供一個待分析處理的音頻信號(.wav文件)或者也可以現(xiàn)用麥克錄段你的聲音。
2.實驗步驟
⑴ 選擇一個wav文件作為分析的對象,當然也可以現(xiàn)用麥克錄段你的聲音。
⑵ 讀取音頻信號的數(shù)據(jù),并畫出其時域波形和頻譜圖。例如:
⑶ 反演它的時域值。當然可以使用富氏反變換,但這里我們利用剛才分析的幾個主要的沖擊峰值來重構。并聽一聽其聲音的變化。
⑷ 給出仿真結論及不足之處。
1.幾個要用到的函數(shù):
⑴ [x,fs,bits]=wavread('filename')這是一個matlab中讀取wav文件的數(shù)據(jù)的函數(shù)。其中的x表示一長串的數(shù)據(jù),一般是兩列(立體聲);fs是該wav文件在采集時用的采樣頻率;bits是指在進行A/D轉化時用的量化位長(一般是8bits或16bits)。
⑵ [d]=fft(w,l)這是matlab中fft函數(shù)的一種輸入輸出形式。w是一列波形數(shù)據(jù);l是指示用多少點的fft,我們應該選擇2的乘方的數(shù)(如16,128,1024等),因為這樣就可以使用優(yōu)化的蝶形算法;d是頻域的輸出。由于fft(dft)的對稱性,又輸入的是實數(shù),fft的結果的復數(shù)序列是共軛反對稱的,所以它們的模的大小對稱,一般來說只用取一半的數(shù)據(jù)就可以了。
⑶ sound(w,fs,bits)和前面的wavread一樣的參數(shù)表示,它將數(shù)列的數(shù)據(jù)通過聲卡轉化為聲音。
你自己做
本站僅提供存儲服務,所有內容均由用戶發(fā)布,如發(fā)現(xiàn)有害或侵權內容,請
點擊舉報。