程式设计入门篇音讯读写.pptVIP

MATLAB 程式設計入門篇：音訊讀寫、錄製與播放 音訊讀寫、錄製與播放 張智星 (Roger Jang) jang@ /jang 清大資工系 多媒體檢索實驗室 20-1 音訊的基本介紹 聲音訊號（Audio Signal）簡稱音訊，泛指由人耳聽到的各種聲音的訊號。音訊有些基本的特質，可說明如下： 音量（Volume）：聲音的大小稱為音量，又稱為力度、強度（Intensity）或是能量（Energy）。音量越大，代表音訊波形的震幅越大。 音高（Pitch）：聲音的基本頻率（Fundamental Frequency）越高，代表音高越高（例如女高音的歌聲）；反之，聲音的基本頻率越低，代表音高越低（例如男低音的歌聲）。（有關基本頻率的說明，將在本章其後各小節說明。） 音色（Timber）：音訊波形在每個週期內的變化，就形成了此音訊的音色。不同的音色即代表不同的音訊內容，例如不同的字有不同的發音，或是不同的歌手有不同的特色，這些都是由於音色不同而產生。 20-2 WAV檔案的讀取範例之一 例如，若要讀取光碟中的檔案 welcome.wav，畫出音訊的波形並播放出此音訊，可使用下列程式： 範例20-1： readWave01.m [y, fs]=wavread(welcome.wav); sound(y, fs); % 播放此音訊 time=(1:length(y))/fs; % 時間軸的向量 plot(time, y); % 畫出時間軸上的波形 20-2 WAV檔案的讀取範例之一 執行結果:歡迎光臨波形圖 20-2 WAV檔案的讀取範例之二 若要知道 welcome.wav 的取樣點是由多少個位元來表示，可使用 [y, fs, nbits] = wavread(‘welcome.wav’)。若要知道音訊長度，則可使用 length(y)/fs。以下範例可以印出音訊檔 welcome.wav 的各種相關資訊： 範例20-2： readWave02.m fileName=welcome.wav; [y, fs, nbits]=wavread(fileName); fprintf(音訊檔案 %s 的資訊：\n, fileName); fprintf(音訊長度 = %g 秒\n, length(y)/fs); fprintf(取樣頻率 = %g 取樣點/秒\n, fs); fprintf(解析度 = %g 位元/取樣點\n, nbits); 20-2 WAV檔案的讀取範例之二 範例20-2執行結果: 音訊檔案 welcome.wav 的資訊： 音訊長度 = 1.45134 秒 取樣頻率 = 11025 取樣點/秒 解析度 = 8 位元/取樣點 20-2 WAV檔案的讀取範例之三 wav 檔案的 8 位元是以 unsigned integer的方式來儲存，因此所能表示的數值是介於 0 和 255 之間，MATLAB 再將此值設定至變數 y 時，會自動將其數值調整至介於 –1 和 1 之間，因此若要將 MATLAB 讀出之數值轉回原先 8 位元所表示之數值，只要將變數 y 乘以 128，再加上 128，就可以得到原先的整數值，例如： 範例20-3： readWave03.m difference = 0 fileName=welcome.wav; [y, fs, nbits]=wavread(fileName); y0=y*(2^nbits/2)+(2^nbits/2); % y0 是原先儲存在音訊檔案中的值 difference=sum(abs(y0-round(y0))) 20-2 WAV檔案的讀取範例之四 wavread 也可以讀取雙聲道或立體聲（Stereo）的音訊檔案，此時傳回的變數為具有兩直行的陣列，每一直行代表一個聲道的音訊，例如： 範例20-4： readWave04.m fileName=flanger.wav; [y, fs]=wavread(fileName); % 讀取音訊檔 sound(y, fs); % 播放音訊 left=y(:,1); % 左聲道音訊 right=y(:,2); % 右聲道音訊 subplot(2,1,1), plot((1:length(left))/fs, left); subplot(2,1,2), plot((1:length(right))/fs, right); 20-2 WAV檔案的讀取範例之四 此範例會讀取雙聲道的音訊檔 flanger.wav，播放此雙聲道的音訊，並畫出兩個聲道的音訊波形如下： 20-2 WAV檔案的讀取範例之五 如果音訊檔案很大，無法一次讀入記憶體，我們也可以使用 wavread