Hilbert Huang transform

Version imprimable

18/08/2016, 16h01
puzzle1989

Hilbert Huang transform

bonjour à tous,

je veux utilisé les codes plot_hht.m et emd.m fournis par ce lien:
https://www.mathworks.com/matlabcent...ent/plot_hht.m

tout d'abord mon signal il est de 2500000 échantillons et de période d'échantillonage(Ts=4e-9).
quand j'applique ce code sous matlab 7.10.0(R2010a), mon ordinateur de 8 Go de ram se plante :( sans donner aucun résultat.
aidez moi S.V.P.

merci d'avance :)
18/08/2016, 16h11
johnette

Essaye de découper ton signal en plusieurs parties.
Tu verras si c'est Matlab qui plante.
18/08/2016, 16h16
puzzle1989

Citation:

Envoyé par johnette

Essaye de découper ton signal en plusieurs parties.
Tu verras si c'est Matlab qui plante.

merci je vais essayer tout de suite
18/08/2016, 16h24
puzzle1989

Citation:

Envoyé par puzzle1989

merci je vais essayer tout de suite

j'ai essayé et pour une partie de 10000 échantillon ça marche et matlab ne se plante pas :mouarf::mouarf:
y'a t-il une solution pour que j'applique hilbert-huang transform sur tt mon signal?
merci d'avance :D
18/08/2016, 16h58
johnette

Tu peux modifier la fonction plot_hht, afin quelle renvoie les variables à ploter.

Ensuite tu boucles la fonction modifiée en plusieurs fois, puis tu concatènes les multiples variables dans une seule.

Pour finir tu peux essayer de l'afficher, si c'est possible !!!

Sinon tu plot 1 point sur 2, 3, 4 ou plus si ça marche.
18/08/2016, 18h34
puzzle1989

Citation:

Envoyé par johnette

Tu peux modifier la fonction plot_hht, afin quelle renvoie les variables à ploter.

Ensuite tu boucles la fonction modifiée en plusieurs fois, puis tu concatènes les multiples variables dans une seule.

Pour finir tu peux essayer de l'afficher, si c'est possible !!!

Sinon tu plot 1 point sur 2, 3, 4 ou plus si ça marche.

tout d'abord,je vous remercie infiniment pour ta réponse;
si c'est possible, pouvez-vous m'expliquer plus votre solution :oops:.
merci

Ca donnerait qqchose comme ça dans un premier temps

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
 
Ts=4e-9
s=signal
for i =1:100
x=s(1+(length(s)/100)*(i-1):i*length(s)/100)
[M(:,i),k1(:,i),k2(:,i),imf(:,i)] = function plot_hht_modif(x,Ts)
 
end
 
 
 
 [M,k1,k2,imf] = function plot_hht_modif(x,Ts)
% Plot the HHT.
% plot_hht(x,Ts)
% 
% :: Syntax
%    The array x is the input signal and Ts is the sampling period.
%    Example on use: [x,Fs] = wavread('Hum.wav');
%                    plot_hht(x(1:6000),1/Fs);
% Func : emd
 
% Get HHT.
imf = emd(x);
for k = 1:length(imf)
   b(k) = sum(imf{k}.*imf{k});
   th   = angle(hilbert(imf{k}));
   d{k} = diff(th)/Ts/(2*pi);
end
[u,v] = sort(-b);
b     = 1-b/max(b);
 
% Set time-frequency plots.
N = length(x);
c = linspace(0,(N-2)*Ts,N-1);
for k = v(1:2)
   figure, plot(c,d{k},'k.','Color',b([k k k]),'MarkerSize',3);
   set(gca,'FontSize',8,'XLim',[0 c(end)],'YLim',[0 1/2/Ts]); xlabel('Time'), ylabel('Frequency');
end
 
% Set IMF plots.
M = length(imf);
N = length(x);
c = linspace(0,(N-1)*Ts,N);
 
end

19/08/2016, 15h04
puzzle1989

bonjour,

j'arrive à extraire les IMFS de chaque échantillon mais concernant les imfs de tout mon signal , je me reste bloqué.
je fais une boucle for sur plot_hht en ajoutant des "hold on" mais ça donne les courbes de chaque echantillon qui sont plus que mille figure:oops:.
avez-vous une idée pour concaténer tous ses échantillons afin d'avoir les imfs de tout le signal?.

merci ;)
22/08/2016, 11h04
johnette

Tu as la fonction reshape pour concaténer ta matrice
22/08/2016, 15h23
puzzle1989

Citation:

Envoyé par johnette

Tu as la fonction reshape pour concaténer ta matrice

bonjour jobnette,

avez-vous une idée s'il est possible d'utiliser le code trouvé sur Mathworks pour des signaux longs?

merci
22/08/2016, 15h38
johnette

Aucune idée ...
Il faudrait faire des tests pour ça.
22/08/2016, 15h41
puzzle1989

Citation:

Envoyé par johnette

Aucune idée ...
Il faudrait faire des tests pour ça.

ok merci jobnette ;););)