Dalgacık dönüşümü büyüklüğüne karşı frekansı nasıl çizebilirim?

Morlet sürekli dalgacık dönüşümü yapıyorum. wscalogramSinyalim var ve şimdi aşağıdaki resimdeki gibi frekans büyüklüğünü çizmek istiyorum, ancak nasıl yapılacağını bilmiyorum:

scal2freqÖlçekleri sahte frekanslara dönüştürmek için MATLAB işlevini kullandım . Ayrıca sinyalimde büyük bir sönümleme oranına (% 4) sahip bazı frekanslar var, bu yüzden arsada iyi görünmüyorlar. Bu yüksek sönümlü modları nasıl abartırım?

MATLAB kullanıyorum, işte kodlarım:

% Import the text4.txt to matlab workspace. and save it under name "data"
t=linspace(0,30,301);
Fs=ceil(inv(t(2)-t(1)));
x=data(:,4); % use x=data(:,3),x=data(:,5) too. first column is time,second is refrence
wname = 'morl';
scales = 1:1:256;
coefs = cwt(x,scales,wname,'lvlabs');
freq = scal2frq(scales,wname,1/Fs);
surf(t,freq,abs(coefs));shading('interp');
axis tight; xlabel('Seconds'); ylabel('Pseudo-Frequency (Hz)');
axis([0 30 0 1 0 60])
xlabel('Time'); ylabel('Frequency');
figure;
sc=wscalogram('image',coefs,'scales',freq,'ydata',x);
hold on 

Şu anki planım şöyle:

Kodunuzu biraz değiştiriyoruz, ancak büyük değişiklikler yok ve her iki şekilde de doğru sonuçlar alıyorum. Bu şablon kodunu burada kullanın ve herhangi bir sorun görmemeniz gerekir. Doğru sonuçları alıyorum.

clear all; 
t=linspace(0,30,301); 
Fs = (inv(t(2)-t(1))); 
x=randn(100,1);  
wname = 'morl'; 
scales = 1:1:256; 
chefs = cwt(x,scales,wname,'lvlabs');

freq = scal2frq(scales,wname,1/Fs);

figure; 
coefsSquared = abs(coefs).^2; 
imagesc(coefsSquared); 
grid off;

%Pick one of the columsn to plot: 
figure; 
plot(coefsSquared(:,47))

İlk görüntünün sağ tarafındaki rakamı elde etmek için, bu sadece zaman dilimi spektrumudur. Başka bir deyişle, frekans-zaman matrisinizin bir sütunu.

Örneğin, frekansa göre zamana göre aşağıdaki kare katsayılı matrisi alabiliriz:

Diyelim ki sütun 47'deki zaman dilimi spektrumuna bakmak istiyoruz. Sonra, bu burada verilmiştir:

Burada, 47. sütuna karşılık gelen zaman dilimi için katsayıların frekansa karşı gücünü görürsünüz.

Sürekli dalgacık dönüşümü (CWT), zaman ölçeği analizi için bir yöntemdir . Evet, doğru okudunuz , sıklığı değil ölçeği . Ancak, ölçekleri frekanslara ve hatta oldukça kolay bir şekilde eşleştirmek mümkündür. Bir MATLAB kullanıcısı olduğunuzdan, muhtemelen aşağıdakileri yapan bu işlevi kullanmak isteyeceksiniz :

F = scal2frq (A, 'wname', DELTA), A tarafından verilen ölçeklere karşılık gelen sözde frekansları, dalgacık fonksiyonu 'wname' (daha fazla bilgi için dalga kaynağına bakın) ve DELTA örnekleme periyodunu döndürür.

Okuduğunuz gibi, frekanslar sinyaldeki gerçek frekanslar değil, yaklaşık değerlerdir. Gerçek frekansları istiyorsanız, kısa süreli Fourier dönüşümünü (STFT) kullanmalısınız.

Peki, CWT ne zaman kullanılmalıdır? Genel bir ilke olarak, sinyaller geçici olduğunda (yani hızlı değişimler olduğunda) dalgacıklar iyi çalışır .

EDIT: haha, burada ne yaptığımı bilmiyorum. özür dilerim: - Sorunuzun noktasını tamamen kaçırdım. Sadece başlığı oku ....

Hala yardıma ihtiyacınız olup olmadığından emin değilim ...

Frekans ve büyüklük istiyorsanız, katsayıları çizmek için "mesh / surf" komutunu kullanın.

Sonra arsa yönünü değiştirin, böylece aradığınız sonucu elde edersiniz. Bunun için "görünüm" komutuna bakın.