Kalau Anda belum tahu cerita saya tentang Transformasi Fourier silahkan klik disini. Jika sudah saya punya pertanyaan…

Baiklah pakai contoh saja dengan Matlab…

Kita buat sebuah sinyal sebagai berikut, masih sama seperti sebelumnya, namun kali ini kedua frekuensi, yaitu 100 Hz dan 200Hz tidak muncul bersamaan tetapi bergantian, apakah Transformasi Fourier mampu melihat kedua frekuensi ini?

fs = 1000;
t = 0:1/fs:0.5;
tx = [t t+t(length(t))];
y1 = sin(2*pi*100*t);
y2 = sin(2*pi*200*t);
y = [y1 y2];
plot(tx,y);
title(‘Sinyal dengan kandungan 2 frekuensi’)
xlabel(‘waktu (detik)’);

Hasilnya, gambaran dari sinyal yang saya tanyakan adalah sebagai berikut…

Kemudian kita lakukan FFT…

Y = fft(y);

Kita gambarkan hasilnya dengan perintah-perintah berikut…

f = fs*(0:(length(Y)-1)/2)/length(Y);
figure;
plot(f,abs(Y(1:(length(Y)+1)/2)));
title(‘Kandungan frekuensi sinyal y (gambar 1 sisi)’);
xlabel(‘frekuensi (Hz)’);

hasilnya sebagai berikut…

Mmmm bisa khan? Sama saja hasilnya… bisa diperoleh 2 frekuensi sesuai dengan dugaan kita, bagaimana jika ditambahkan derau kemudian di-FFT…

ya= y + 2*randn(size(tx));
figure;
plot(tx,ya);
title(‘Sinyal apakah ini….??’);
xlabel(‘waktu (detik)’)
;

YA = fft(ya);
f = fs*(0:(length(YA)-1)/2)/length(YA);
figure;
plot(f,abs(YA(1:(length(YA)+1)/2)));
title(‘Kandungan frekuensi sinyal apa ini….??’);
xlabel(‘frekuensi (Hz)’);

Nah hasilnya…

Sama seperti artikel saya yang lalu khan? Luar biasa Transformasi Fourier ini…

Sekarang pertanyaan saya, lantas bisakah kita tahu bahwasanya ke-2 frekuensi tidak bersamaan munculnya? Ya tidak bisa-lah… hanya kandungan frekuensi saja, sedangkan kapan dan lama waktu muncul masing-masing freuensi itu kita tidak tahu… lantas kalo mo tahu? Ya pake lainnya donk… apaan tuch? Pake STFT (Short Time Fourier Trasnform)…

Wah apa lagi nich…

Ceritanya begini, jika TF bekerja untuk seluruh sinyal, tapi STFT hanya bekerja pada sebuah jendela yang kecil yang kemudian digeser-geser mulai dari awal hingga akhir untuk mendapatkan interpretasi data keseluruhan secara waktu dan frekuensi atau istilahnya time-frequency domain… di Matlab pake perintah specgram()

figure;
specgram(y,256,1000);

Keterangan:
256 sebagai panjang jendela, sedangkan 1000 merupakan fs-nya

Hasilnya…

Nah tuch… kelihatan bahwa kedua frekuensi muncul secara tidak bersamaan, lebih tepat berturutan, hanya saja tidak terlalu jelas dimana tepatnya frekuensi mulai bergantian… Baik sekarang Anda perhatikan masing-masing perintah dan hasil gambarnya sebagai berikut:

figure;
specgram(y,64,1000);
figure;
specgram(y,128,1000);
figure;
specgram(y,256,1000);

Hasilnya secara berturutan…

Mm menarik hasilnya, dengan semakin besar ukuran jendela, semakin akurat resolusi frekuensinya, tapi semakin gak jelas resolusi waktunya. Demikian juga sebaliknya, semakin kecil ukuran jendelanya, semakin bagus resolusi waktunya, tapi resolusi frekuensi-nya makin jelek…

Ini-lah yang dimaksudkan dengan Ketidak-pastian Heisenberg… ada semacam trade-off antara resolusi waktu dan frekuensi, tapi minimal sudah kita peroleh ranah waktu-frekuensi, alhamdulillah…

Ada komentar saudara-saudari sekalian?

Tags: , ,

19 Responses to “Transformasi Fourier – Kapan muncul frekuensi-nya?”

  1. tau cra buat simulasi TF pake simulink g???

  2. Waduh, saya lagi pusing bikin spectrum analyser sederhana untuk freq 1-100KHz. kira2 yang cocok pake mikrokontroler jenis apa ya?

    Please, FFT-nya diporting ke mikro aja mas, dibuat spectrum analyzer.
    Frequency samplernya kecil gapapa yang penting bisa buat belajar. Kalo ga ya flowchartnya aja.

  3. @Fitria:
    Sebaiknya untuk pemilihan uC menyesuaikan frekuensi sampling yang digunakan. Untuk 100 kHz maks berarti frekuensi samplingnya minimal harus dua kali (Teorema Nyquist) atau 200 kHz.

    Untuk uC sekarang sudah dalam frekuensi kerja MHz, jadi pilih saja yang murah dan enak dipake (he he he kayak baju aja)…

    Terima kasih, kapan2, Insya Alloh saya coba mengulas implementasi DSP menggunakan uC…

  4. mas mo mnt tolong nih…
    sya da tgs dari dosen,bkn paper aplikasi dari transformasi fourier

  5. loh ko sama mas QQ,,jngn2 Qta kul di kmps yang sama,heee
    tapi aq jg blm dpt bahannya,,kl dah dpt krim emeil y,, Yuddy21tea@yahoo.com
    y nih Ag bth bagt paper tentang aplikasi transformasi fourier

  6. permisi mo tanya.. untuk magnitude dalam frekuensi audio itu satuannya apa Pak..? trims

  7. sip sekali pak tulisannya.. kebetulan sangat berguna sekali buat saya yang tidak punya dasar apa2 dibidang sinyal.. kalo bisa, diperbanyak tentang persinyalan pak..hehehe..

    terima kasih.

  8. Pak Agfi, saya lagi ngerjain tugas akhir nih
    tentang speech recognition untuk aplikasi command and control
    yang saya mau tanya pak, apakah cukup parameter yang diambil dari sinyal suara hanya dari magnitude FFTnya saja?
    ide ini saya dapat dari tulisan bapak yang berjudul
    “Pemisahan derau dan suara – Speech Noise Reduction”
    Terimakasih sebelumnya

  9. mohon petunjuk yang lebih jelas pak…
    terima kasih

  10. Sekarang pertanyaan saya, lantas bisakah kita tahu bahwasanya ke-2 frekuensi tidak bersamaan munculnya?

    Bisa pak..pake filter bank dan juga ada teori Natural Observation Method…
    Kalo pake FFT ya pake window….
    Tks

  11. untuk mencari kapan munculnya – bisa menggunakan STFT (dengan resolusi sama rata) atau Wavelet (resolusi tidak sama rata, tradeoff)…

  12. terimakasih telah berbagi pak

  13. thank you for sharing it’s so helpful

  14. terimakasih sangat membantu

  15. thank you for sharing it’s so helpful

  16. HASIL DATA & LIVE DRAW TOGEL VEGAS
    Paito Warna Terpopuler
    live draw sydney
    Paito Warna singapore
    Paito Warna Singapore

  17. The most favorite website I’ve ever read with a collection of high-knowledge words. Everything is very opeen with a really clear description of the issues.

  18. A very interesting place and full of information that is very supportive among young people.
    Syair Sgp Mbah Semar

Trackbacks/Pingbacks

  1. Rangkuman artikel Pemrosesan Sinyal Digital sepanjang 2008-2012 | DSP & Embedded Electronics

Leave a Reply

You can use these tags: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>