V. Homomorphic Signal Processing.ppt

* X = fft(x); plot(abs(X)) fft 橫軸 轉換的方法 (1) Using normalized frequency F: F = m / N. (2) Using frequency f, f = F ? fs = m ? (fs / N). B. 繪出頻譜 (請參考附錄四) * (1) wavplay(x): 將 x 以 11025Hz 的頻率播放 (時間間隔 = 1/11025 = 9.07 ? 10?5 秒) (2) sound(x): 將 x 以 8192Hz 的頻率播放 (3) wavplay(x, fs) 或 sound(x, fs): 將 x 以 fs Hz 的頻率播放 Note: (1)~(3) 中 x 必需是1 個column (或2個 columns)，且 x 的值應該 介於 ?1 和 +1 之間 (4) soundsc(x, fs): 自動把 x 的值調到 ?1 和 +1 之間 再播放 * C. 聲音的播放 wavwrite(x, fs, waveFile) 將數據 x 變成一個 *.wav 檔，取樣速率為 fs Hz ? x 必需是1 個column (或2個 columns) ? x 值應該 介於 ?1 和 +1 之間 ? 若沒有設定fs，則預設的fs 為 8000Hz * D. 用 Matlab 製作 *.wav 檔： wavwrite E. 用 Matlab 錄音的方法 錄音之前，要先將電腦接上麥克風，且確定電腦有音效卡 (部分的 notebooks 不需裝麥克風即可錄音) 範例程式： Sec = 3; Fs = 8000; recorder = audiorecorder(Fs, 16, 1); recordblocking(recorder, Sec); audioarray = getaudiodata(recorder); 執行以上的程式，即可錄音。 錄音的時間為三秒，sampling frequency 為 8000 Hz 錄音結果為 audioarray，是一個 column vector (如果是雙聲道，則是兩個 column vectors) * wavplay(audioarray, Fs); % 播放錄音的結果 t = [0:length(audioarray)-1]./Fs; plot (t, audioarray‘); % 將錄音的結果用圖畫出來 xlabel(sec,FontSize,16); wavwrite(audioarray, Fs, ‘test.wav’) % 將錄音的結果存成 *.wav 檔 範例程式 (續)： * recordblocking(recorder, Sec); (錄音的指令) recorder: the parameters obtained by the command “audiorecorder” Sec: the time length for recording audioarray = getaudiodata(recorder); (將錄音的結果，變成 audioarray 這個 column vector，如果是雙聲道，則 audioarray 是兩個 column vectors) 以上這三個指令，要並用，才可以錄音 指令說明： recorder = audiorecorder(Fs, nb, nch); Fs: sampling frequency, nb: using nb bits to record each data nch: number of channels (1 or 2) (提供錄音相關的參數) * Video 檔讀取: aviread Image 檔讀取: imread Image 檔製作: imwrite Image 檔顯示: imshow, image, imagesc F：影像檔的處理 基本概念：灰階影像在 Matlab 當中是一個矩陣 彩色影像在 Matlab 當中是三個矩陣，分別代表 Red, Green, Blue * *.bmp: 沒有經過任何壓縮處理的圖