Discussion :

[Gazelle Arctique]

Re-bonjour à tous,

J'ai un graphique issu d'un test expérimental. Je voudrais calculer l'intégrale de ce signal entre deux points. Pour cela, j'utilise la fonction integral qui nécessite une fonction handle. Comme lu dans la FAQ, j'utilise l'arobase.

Voici un extrait de mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
clear all; close all;
 
[EDA_L1,date_L1]=xlsread('EDA\EDAL1.xlsx');
 
%[...]
 
T0=0; T1=15; T2=16.07;T3=25;T4=27.38;T5=30;T6=31.4;T7=40;T8=41.6;
 
%[...]
handle=@EDA_L1;
AUC_Zone1 = integral(handle, T0, T1)

J'obtiens l'erreur matlab suivante:

"EDA_L1" was previously used as a variable, conflicting with its use here as the name of a function or command.

Après quelques recherches google, j'ai vu qu'une solution serait d'écrire directement l'équation de la courbe mais dans mon cas, comme il s'agit plus d'un signal que d'une simple droite ou courbe, ce n'est pas si évident.

Comment pourrais-je faire pour contrer cette erreur Matlab ?

Je vous remercie pour votre aide.

[Invité]

Bonjour,

La fonction integral utilise une fonction dont tu connais l'expression (ex. y = f(x), tu prendrais alors le handle de cette fonction @f), or ici tu disposes d'un vecteur de données numériques.
C'est la fonction trapz à utiliser.

[Gazelle Arctique]

Merci pour votre réactivité !

En fait, j'avais déjà essayé la fonction trapz :

Code :

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

1
2
3
4
5
[EDA_L1,date_L1]=xlsread('EDA\EDAL1.xlsx');
 
T = [0 15 16.07 25 27.38 30 31.4 40 41.6];
 
AUC_Zone1 = trapz([T(1) T(2)], EDA_L1)

J'obtiens l'erreur suivante :
"LENGTH(X) must equal the length of Y in dim 1."
Si je cherche à calculer l'intégrale entre les temps T3 et T4 par exemple, là Matlab capte que j'ai un souci de longueur entre X et Y...

Vous pourriez m'apporter vos lumières ?

Je vous remercie vraiment de prendre le temps de m'aider

[Invité]

EDA_L1 contient un certain nombre de valeurs, qui correspondent plus ou moins j'imagine aux données temporelles de ton vecteur T.
Mais comme MATLAB, je me demande comment est cette répartition.

Si par exemple je te donne deux vecteurs :

Code :

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

1
2
values = rand(1,100);
T = [1 30 56];

Sauras-tu me dire (sans plus de précision) quelles sont les valeurs de values contenues dans chaque intervalle de temps ? (Rien ne dit qu'elles sont équitablement réparties, ce serait trop facile )

La fonction trapz nécessite que tu lui donnes deux vecteurs :
Celui de tes données bien sûr, mais aussi celui du vecteur temps de même longueur auquel correspond chaque valeur, afin d'avoir les couples :

Code :

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

1
2
3
4
T(1)     EDA_L1(1)
T(2)     EDA_L1(2)
...
T(end)   EDA_L1(end)

[Gazelle Arctique]

Oh non ! Il pourrait faire un effort le p'tit Matlab...

Donc si j'ai bien compris, je dois créer un nouveau vecteur avec les 9 points en ordonnées correspondant à mon vecteur temporel T; c'est-à-dire :

Code :

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

1
2
T = [ 0 1 2 3 4 5 6 7 8 ]
EDA_L1 = [ 3 2 4 9 6 8 4 5 7 ]

et la ligne de code pour calculer l'aide sous la courbe entre deux points s'écrirait :

Code :

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

AUC_Zone1 = trapz([T(1) T(2)], [EDA_L1(1) EDA_L1(2)])

?

Je te remercie beaucoup pour ton explication, très claire !

[Invité]

Ici oui, tu calculerais l'aire en gris :

Mais si tu étends ta méthode à

Code :

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

AUC_Zone1 = trapz([T(i) T(j)], [EDA_L1(i) EDA_L1(j)])

Cela ne marchera pas, tu calculeras seulement (ici i = 3; j = 7;) cette aire cantonnée aux seuls deux points que tu fournis à trapz :

Si tu veux prendre en compte le reste, il faudra aussi le fournir

Code :

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

AUC_Zone1 = trapz(T(i : j), EDA_L1(i : j));

Pour au final bien avoir cette aire :

[Gazelle Arctique]

Aaah ok !

Merci pour la nuance ! Effectivement, ça change tout.

Je vais m'atteler à récupérer mes valeurs et essayer ça.

Merci beaucoup !!!

[tachmou]

Pour la peine, je mets un groooos +1 à Winjerome qui a passé du temps à t'expliquer ça avec des graphiques !