Discussion :

[bendesarts]

Bonjour,

J'ai deux fonctions qui dépendent du temps.
Pour la fonction 1, je dispose :
- du vecteur t de 83 elements
- du vecteur x1(t) de 83 élements

Pour la fonction 2, je dispose:
- du vecteur t de 130 élements
- du vecteur x2(t) de 130 éléments.

J'aimerais tracer la fonction qui a :
- en abscisse : t de dimension 83 élements
- en ordonnée : abs(x1(t)-x2(t)) de dimension 83 élements.

Comme j'ai des vecteurs de tailles différentes, pouvez m'indiquer la manip à faire ?

Je vous remercie pour votre aide.

[Gooby]

Bonjour,

Souhaites-tu effectuer cette opération sur les 83 éléments de l'un et les 83 autres premiers éléments de l'autre?

Si oui, il te suffit de faire ceci:

Code :

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

abs(x1(:)-x2(1:83))

sinon, précise ta pensée.

[bendesarts]

Alors j'essayer de préciser ma pensée...

J'aimerais tracer la différence x1(t) - x2(t).

Par contre, je n'ai pas le même nombre de points pour la variable temps pour les deux fonctions.

Pour que ma différence est un sens, il faudrait faire x1(t)-x2(t) mais avec la même variable du temps. Donc si je prends par exemple la variable temps associée à la fonction x1(t), il faut que je calcule x2(t) pour tous les instants définies par le vecteur temps associées à la fonction 1 avant de pouvoir faire la différence.

Est-ce que quelqu'un aurait une idée de la manip à faire ?

Je pense qu'il doit s'agir d'une opération assez classique en matlab ?

Merci pour votre aide

[Gooby]

ça me paraît pas si classique que ça. Est-ce que tu as les échelles de temps associer à tes signaux histoire qu'on ait un ordre de grandeur?

[thecrazydonut]

Est-ce que quelqu'un aurait une idée de la manip à faire ?

Bonjour,

Si j'ai bien compris, le début et la fin de tes vecteurs x1 et x2 sont les mêmes. La solution la plus simple, selon moi, est d'interpoler tes données contenues dans ton vecteur le plus court. J'ai posté une question sur l'interpolation il y a déjà quelque temps là et la solution de Winjerome est très facile à mettre en oeuvre. En revanche, je ne sais pas si tu peux te permettre d'interpoler tes données dans ton cas, c'est à toi de voir .

[bendesarts]

Exactement, cette piste m'a l'air intéressante.

Rappel:
Ma fonction 1 se compose deux vecteurs (de taille 83 éléments):
t1
x1(t1)

Ma fonction 2 se compose deux vecteurs (de taille 130 éléments):
t2
x2(t2)

Voici ce que j'ai fait :

Pour ma fonction 2, j'ai crée le vecteur x2interp qui calcule x2 aux points du vecteur t1 par interpolation de la manière suivante:

Code :

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

x2interp=interp1(t2,x2,t1)

Code :

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

plot(t1,abs(x1-x2interp))

Par contre, effectivement, dans mon cas, l'interpolation doit être bien maîtrisée car je souhaite mesurer des écarts très petits en deux courbes.

Il faut donc que j'arrive à faire que les erreurs dues à l'interpolation soit inférieur au moins 10% des erreurs que je souhaite mesurer entre mes deux courbes.

Aussi, il faut que je règle bien:
- soit la grille d'échantillonnage (pas discrétisation) de mon interpolation
- soit la méthode utilisée pour l'interpolation.

Avez-vous des conseils à ce sujet afin que je puisse limiter les erreurs liées à l'interpolation ?

Merci d'avance pour vos conseils