Come posso tracciare la frequenza rispetto all'entità della trasformazione wavelet?

Sto eseguendo la trasformazione wavelet continua di Morlet. Ho il wscalogramsegnale e ora voglio tracciare la grandezza di freq come nella figura seguente, ma non so come farlo:

Ho usato la scal2freqfunzione MATLAB per convertire le scale in pseudo frequenze. Inoltre ho alcune frequenze nel mio segnale che hanno un ampio rapporto di smorzamento (4%), quindi non sono ben visibili nella trama. Come posso esagerare con queste modalità altamente smorzate?

Sto usando MATLAB, ecco i miei codici:

% 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 

La mia trama attuale si presenta così:

Modificando un po 'il tuo codice, ma senza grandi cambiamenti, e ottengo risultati corretti in entrambi i modi. Usa questo codice modello qui e non dovresti vedere alcun problema. Ottengo i risultati corretti.

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))

Per ottenere la figura sul lato destro della prima immagine che hai messo, questo è semplicemente lo spettro della fascia oraria. In altre parole, una colonna della tua matrice frequenza-tempo.

Ad esempio, potremmo ottenere la seguente matrice a coefficienti quadrati, in frequenza rispetto al tempo:

Supponiamo quindi che desideriamo esaminare lo spettro per la fascia oraria nella colonna 47. Quindi, questo è dato qui:

Qui, quindi, viene visualizzata la potenza dei coefficienti, rispetto alla frequenza, per l'intervallo di tempo corrispondente alla colonna 47.

La trasformazione wavelet continua (CWT) è un metodo per l' analisi della scala temporale . Sì, l'hai letto correttamente, scala , non frequenza . Tuttavia, è possibile mappare le scale alle frequenze e persino abbastanza facilmente. Dato che sei un utente MATLAB, probabilmente vorrai utilizzare questa funzione , che procede come segue:

F = scal2frq (A, 'wname', DELTA) restituisce le pseudo-frequenze corrispondenti alle scale date da A, la funzione wavelet 'wname' (vedi wavefun per maggiori informazioni) e il periodo di campionamento DELTA.

Come puoi leggere, le frequenze non sono le frequenze reali nel segnale, ma le approssimazioni. Se vuoi le frequenze reali, dovresti usare la trasformata di Fourier di breve durata (STFT).

Quindi, quando si dovrebbe usare il CWT? Come principio generale, le wavelet funzionano bene quando i segnali sono transitori (cioè hanno rapidi cambiamenti).

EDIT: haha, non so cosa sto facendo qui. scusa MrGreen completamente perso il punto della tua domanda. Leggi il titolo ....

Non sono sicuro se hai ancora bisogno di aiuto ...

Se si desidera la frequenza rispetto all'ampiezza, utilizzare un comando "mesh / surf" per tracciare i coefficienti.

Quindi cambia l'orientamento della trama in modo da avere il risultato che cerchi. Per questo vedere il comando "visualizza".