Discussion :

[azassma]

Bonjour,

S'il vous plait, j'ai besoin de calculer la vitesse de phase en utilisant cette formule

Code :

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

 V = w2L/(Phase1 - Phase2)

Phase 1 et Phase 2 sont des vecteurs de taille m. Le problème est que je dois travailler avec des spectres de phase continues, pour cela j'aurai besoin de dérouler les phases sous forme continue mais je ne sais pas comment faire.
Est ce que vous pouvez m'aider s'il vous plait?

Merci

[jyber]

unwrap ?

[azassma]

Les deux phases indiquées dans la formule de la vitesse de phase correspondent aux phases des deux échos de notre signal temporel.
Or, la transformée de Fourrier du premier écho isolé au moyen du programme de filtrage temporel, donne les parties réelle R₁(w) et imaginaire I₁(w) ce qui permet de calculer la phase selon la relation suivante :

Code :

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

Phase1=arctg (I<sub>1</sub>(w)/R<sub>1</sub>(w))

De la même manière on calcule la phase 2 du deuxième écho.

Code :

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

Phase2=arctg (I<sub>2</sub>(w)/R<sub>2</sub>(w))

Dans le calcul de ces deux phases, on obtient deux spectres dont la phase varie de -Pi/2 et Pi/2;
Alors qu'il est nécessaire d'aboutir à des spectres de phase continues pour le calcul de la vitesse de phase et je ne sais pas comment les rendre continues.

Merci beaucoup

[jyber]

help unwrap
-- Function File: B = unwrap (X)
-- Function File: B = unwrap (X, TOL)
-- Function File: B = unwrap (X, TOL, DIM)
Unwrap radian phases by adding multiples of 2*pi as appropriate to
remove jumps greater than TOL. TOL defaults to pi.

[azassma]

Merci beaucoup pour votre aide
Effectivement le problème est résolu avec la fonction matlab "unwrap"