Tek Boyutlu Bir İşaretin Zaman Frekans Düzleminde Görüntülenmesi


Tek boyutlu zamana bağlı bir fonksiyonun zaman-frekans düzleminde gözlemlenmesi bir çok kullanıcıyı ilgilendiren bir konudur. Matlabdaki sürekli zaman dalgacık dönüşümünü gerçekleştiren CWT komutu bu işlemi gerçekleştirmektedir. Ancak CWT komutu zaman-ölçek (scale) bilgisini vermekte olup frekans sonucu için ilave eklentiye ihtiyaç duyulmaktadır. cwt_signal fonksiyonu girilen bir vektör için Fs örnekleme frekansında (Hz), F_ilk başlangıç frekansı ve F_son bitiş frekansı aralığında kalan zaman-frekans düzleminde delta_F çözünürlükte sonucu vermektedir. Baz fonksiyonu olarak çoğu uygulama en tercih edilen ve en sağlıklı sonucu veren baz fonksiyonu olduğu için komleks Morlet seçilmiş ve aşağıdaki kodda 1.5 Hzlik merkez frekansı ve 1.5 Hzlik genişlik anlamına gelecek şekilde "cmor1.5-1.5" olarak girilmiştir. İstenirse bu kullanıcı tarafından fonksiyonun içine girilip değiştirilebilir.



function coef=cwt_signal(isaret,Fs,F_ilk,F_son,delta_F);
delta=1/Fs;
Fa=F_ilk:delta_F:F_son;
Fc=centfrq('cmor1.5-1.5'); %mother wavelet complex morlet center frequency 1.5 Hz %Width 1.5 Hz
scale=Fc./(Fa*delta);
coef=cwt(isaret,scale,'cmor1.5-1.5');
mesh(abs(coef))


Comments

Popular posts from this blog

Latex'te Denklem İçerisine Ufak Boşluklar Koymak

LaTeX'te Sunum Hazırlamak

Octave'da Grafik Çizdirme