音讯处理二.PPTVIP

音訊處理簡介傅立葉轉換與濾波器之應用 清大資工系 多媒體資訊檢索實驗室 傅立葉轉換 將訊號轉換成不同頻率的正(餘)弦波疊加的方法 假設x(n)是原訊號，X(k)是轉換後的結果： 可使用Matlab的fft或SAP toolbox內的fftOneSide進行轉換 轉換的結果，稱之為頻譜 X為長度512的向量，分別平均對應到0~16000的頻率 t = linspace(0, 1, 16000); % fs is 16000 x = sin(2*pi*440.*t); % 440 Hz 的正弦波 X = fft(x(1:512)); plot(abs(X), .-); 傅立葉轉換 藉由頻譜的不同，可以辨別出不同的聲音 [y1 fs nbits] = wavread(clarinet.wav); [y2 fs nbits] = wavread(violin.wav); frameSize = 512; Y1 = fft(y1(8000:8000+frameSize-1)); Y2 = fft(y2(8000:8000+frameSize-1)); freqStep = fs/frameSize; freq = freqStep*(0:frameSize/2); figure; subplot(2,1,1); plot(freq, abs(Y1(1:257)), .-); subplot(2,1,2); plot(freq, abs(Y2(1:257)), .-); 傅立葉轉換 將頻譜並排起來，就是頻譜圖： [y fs] = wavread(ㄚ、ㄧ、ㄨ、ㄝ、ㄛ.wav); frameMat = buffer2(y, 512, 0); magSpec = fftOneSide(frameMat, fs); powerSpecInDb = 20*log(magSpec); subplot(2,1,1); x = (0:length(y)-1)./fs; plot(x, y); axis tight subplot(2,1,2); imagesc(powerSpecInDb); axis xy; 傅立葉轉換 執行結果： 濾波器 一個濾波器，可以用兩個向量a和b來代表 a = [1, a2, a3, ..., ap] b = [b1, b2, b3, ..., bq]; 如果應用到訊號裡面，則可表示如下 y[n] = b1x[n] + b2x[n-1] + ... + bqx[n-q+1] - a2y[n-1] - a3y[n-2] - ... - apy[n-p+1] 可使用filter函式，來計算濾波後的結果 濾波器 如果： a = [1]; b = [0.2 0.2 0.2 0.2 0.2]; 則： y[n] = 0.2x[n] + 0.2x[n-1] + ... + 0.2x[n-4] 相當於五點求平均 較尖銳、變化劇烈的部分會變平滑 高頻部分會被壓低→低通濾波器 濾波器應用 範例：低通濾波器 [x fs nbits] = wavread(movie.wav); a = [1]; b = [0.2 0.2 0.2 0.2 0.2]; y = filter(b, a, x); wavwrite(y, fs, nbits, result.wav); time = (1:length(x))/fs; subplot(2,1,1); plot(time, x); title(Original signal); subplot(2,1,2); plot(time, y); title(Output signal); 濾波器應用 範例：高通濾波器 [x fs nbits] = wavread(movie.wav); a = [1]; b = [1 -1]; % y[n] = x[n] - x[n-1] % 變化劇烈的訊號（即高頻訊號），其差異會被突顯出來 y = filter(b, a, x); wavwrite(y, fs, nbits, result.wav); time = (1:length(x))/fs; subplot(2,1,1); plot(time, x); title(Original signal); subplot(2,1,2); plot(time, y); title(Output signal); 濾波器應用 範例：單次回音器 [x fs nbits] = wavread(movie.wav); delay = 0.2; gain = 0.7; a = [1]; b = [1 zeros(1, delay*fs-1) gain]; % y[n] = x[n] + gain*