Discussion :

[souf308]

Bonjour a tous, j'aurai besoin d'un peu d'aide pour l'utilisation de la fonction mscohere, je n'y comprend rien.

j'ai 2 signaux de frequence d'echantillonage 1250 et les vecteurs sont de longueur 12500*2 (les deux).
Quand j'utilise la fonction mscohere de la sorte

Code :

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

Coh = mscohere(x,y,1000,0,1000,1000)

ca me renvoie l'erreur suivante

Code :

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

1
2
3
4
5
6
7
8
9
??? Error using ==> welchparse at 32
The length of the two input vectors must be equal to calculate the cross
spectral density.
 
Error in ==> welch at 40
[x,M,isreal_x,y,Ly,win,winName,winParam,noverlap,k,L,options] = ...
 
Error in ==> mscohere at 89
[varargout{1:nargout}] = welch({x,y},esttype,varargin{:});

En fait cette erreur reviens a chanque fois sachant que la longueur de LFPA et LFPB est bien la meme. Je comprend pas ou est le probleme ayant essayer de modifier les autres options, je me retrouve toujours a cette meme erreur "The length of the two input vectors must be equal to calculate the cross
spectral density."

Une idée?
Merci

[Jerome Briot]

les vecteurs sont de longueur 12500*2 (les deux).

Je ne connais pas cette fonction mais les vecteur ne devraient-ils pas avoir comme dimension 12500x1 ?

[souf308]

Merci pour ta reponse, j'avais essayer ca avant mais ca n'avais pas marcher, mais en reessayant 12500*1 encore avec differents parametre ca marche mieux j'ai plus le message d'erreur.

J'ai un probleme de parametrage cependant car le resultat me done une coherence de 1 partout ce qui est impossible.

Comment doit on choisir les parametres s'il te plait (fenetre etc..)? J'ai eu des cours de traitement du signal mais passer du théorique a la pratique j'ai un peu de mal la..

mscohere(x,y,window,noverlap,nfft,fs)

-pour windows je prend une hanning
-noverlap je sais pas trop quoi mettre
-nfft non plus..

Y'a t'il des calculs a faire pour definirs les bon paramettres?
Merci

[souf308]

Bon j'ai bossé un peu et je comprend mieu comment choisir le paraammetre lenght of FFT.

(Dites moi svp si j'ai faux..)

lenght of FFT permet d'avoir une meilleur resolution plus il tend vers Fs car
resolution=Fs/(FFT size)

Il me rest encore le nooverlaping et quel fennetre choisir. Sur quel criteres on se base?
merci

[tachmou]

Salut,

J'ai déjà rencontré ce genre de problème, et je n'ai trouvé aucun document expliquant clairement comment paramétrer cette fonction... j'ai l'impression que le flair et les tests est le mieux que tu puisses faire malheureusement. De plus, ça doit dépendre de ton cas d'étude et des bandes que tu étudies. De mes propres tests et de ce que je me souviens, tu risque d'être limité à deux décades d'étude mais il vaut mieux étendre l'exécution de la fonction à 4 décades puis tronquer à deux... ce qui fait que la fenêtre et la précision que tu obtiens à la fin est pour ainsi dire... merdique. Mais malgré tout, c'est comme ça que j'ai pu obtenir les meilleurs résultats.

Bon courage, et si tu as/trouves de bons documents / bonnes explications : je suis preneur (et je ne serai pas le seul je pense)

Cordialement,