Hilbert Huang transform

**puzzle1989** · 18/08/2016, 16h01

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

**johnette** · 18/08/2016, 16h11

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

**puzzle1989** · 18/08/2016, 16h16

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

**puzzle1989** · 18/08/2016, 16h24

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

y'a t-il une solution pour que j'applique hilbert-huang transform sur tt mon signal?
merci d'avance

**johnette** · 18/08/2016, 16h58

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.

**puzzle1989** · 18/08/2016, 18h34

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

.
merci

**johnette** · 19/08/2016, 08h51

Ca donnerait qqchose comme ça dans un premier temps

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**puzzle1989** · 19/08/2016, 15h04

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

.
avez-vous une idée pour concaténer tous ses échantillons afin d'avoir les imfs de tout le signal?.

merci