Buka – buka file lama di folder kuliah S1 Teknik Fisika, saya menemukan file script matlab yang saat itu dimintai tolong seorang teman untuk mengerjakan Tugas Akhirnya. Script ini untuk memfilter sinyal (suara) pada frekuensi tertentu dengan menggunakan bandpass filter, yakni gabungan antara low pass filter dan high-pass filter.
Misalnya kita mempunyai sinyal suara dengan frekuensi 0 – 20000 Hz, namun karena suatu hal (misalnya noise) kita cuma butuh sinyal di frekuensi 700 – 1400 Hz saja, maka kita bisa menggunakan script matlab / octave ini. Script ini berjalan lancar baik pada Matlab maupun GNU Octave.
Gambar spectrogram di bawah ini membandingkan spectrogram sinyal sebelum difilter dan setelah difilter. Script untuk GNU / Octave ada di bawahnya.
|
Sinyal Sebelum di filter |
|
|
Sinyal Setelah di filter |
|
Berikut script kode Octave / Matlab (link Github disini).
% demo program untuk memfilter suara pada range fekuensi oktaf tertentu
% by bagustris bagustris@linuxmail.org
% adopted from somewhere (forget aka lupa)
% set general variables
clear all; close all; clc;
sf = 22050; % sample frequency
nf = sf / 2; % nyquist frequency
d = 1.0; % duration, ganti untuk merubah durasi suara
n = sf * d; % number of samples
nh = n / 2; % half number of samples
% set variables for filter
lf = 700; % lowest frequency, ganti dengan batas bawah bandpass
hf = 1400; % highest frequency, ganti dengan batas atas bandpass
lp = lf * d; % ls point in frequency domain
hp = hf * d; % hf point in frequency domain
% band pass filter!
filter = zeros(1, n); % initializaiton by 0
filter(1, lp : hp) = 1; % filter design in real number
filter(1, n - hp : n - lp) = 1; % filter design in imaginary number
% make noise
rand('state',sum(100 * clock)); % initialize random seed
noise = randn(1, n); % Gausian noise
noise = noise / max(abs(noise)); % -1 to 1 normalization
% do filter
s = fft(noise); % FFT
s = s .* filter; % filtering
s = ifft(s); % inverse FFT
s = real(s);
%untuk mendengarkan suara asli
sound(noise,sf);
pause(d+0.5);
%untuk mendengarkan suara terfilter (bandpass);
sound(s,sf);
pause(d+0.5);
%untuk menyimpan file suara dalam format .wav
%ketik perintah berikut dalam command window (tanpa tanda persen%)
% wavwrite(s,sf 'namafile.wav');
Untuk mengecek hasilnya (spectrum yang telah terfilter) bisa menggunakan spectrogram sbb:
figure(1); specgram(noise, 512, sf);
figure(2); specgram(s, 512, sf);
Silakan dicoba. Please note the correctness and accurateness is not guaranteed.