Discussion :

[MrPingouin]

Bonjour,

Existe t-il un équivalent à cette opération sous une forme vectorielle et donc plus rapide pour matlab..

Code :

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

1
2
3
4
5
 
o=zeros(1,length(M))
for i=1:length(mex)
    o(mex(:,1))=mex(:,2)
end

mex est une matrice avec 2 colonnes dont la premiere represente la position des pics d'un signal et la seconde la valeur de ces pics.

L'opération ci dessus a pour but de creer un vecteur avec les valeurs des pics à leur position afin de garder une taille de vecteur equivalente à M le signal original.

Merci d 'avance

[christophe_halgand]

Bonjour,

Ta boucle ne sert à rien puisque tu n'utilises pas l'index de ta boucle...

Est ce que

Code :

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

length(M) = max(mex(:,1))

?

Christophe

[Newenda]

L'opération ci dessus a pour but de creer un vecteur avec les valeurs des pics à leur position afin de garder une taille de vecteur equivalente à M le signal original.

Comment peux tu avoir un vecteur avec uniquement certains max (tes pics) qui serait de la même taille que ton vecteur original?

[christophe_halgand]

Je pense qu'il veut créer un vecteur de même longueur que celui de base mais seule la valeur des pics serait différent de zéro.

Tu mets tout à zéros et là où tu as des pics, tu insères la valeur de ce pic.

Christophe

[Newenda]

Ha oui ok mais dans tous les cas il est nécessaire de définir ce qu'est un max (bcp de discussion sur ce sujet) :

-Un max dans une fenêtre de position donnée ?
-là ou l'on observe un changement de signe de la dérivée ou bien là ou la dérivée est nulle (=tous les "pics") ?
-Un max qui dépasse un seuil donné ?

Ce n'est qu'après qu'on pourra commencer à discuter matlab.

[Caro-Line]

Tout ça m'a l'air bien embrouillé.

@Newenda : je pense que tu pars un peu trop loin par rapport à la question posée

@MrPingouin :
Christophe a raison sur le fait que la boucle ne sert à rien.
Par contre il faut que tu sois sur que length(M) soit supérieur au maximum des mex(:,1).
Sinon ton initialisation ne sert à rien.

Je te conseille fortement de changer le nom de ta variable mex, MEX étant une fonction de MATLAB.

[Newenda]

@Newenda : je pense que tu pars un peu trop loin par rapport à la question posée

Bah je sais pas trop mais il parait tout de même nécessaire de faire une condition pour trouver ces max non? comment peut on définir un max avec une simple boucle? sauf si on travail justement sur une valeur max dans un offset donné. Une condition simple serait un truc du genre

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
for i=2:length(X)-1
 
if X(i-1) > X(i) & X(i+1) > X(i) 
X(i)=pic bas
elseif  X(i-1) < X(i) & X(i+1) < X(i) 
X(i)=pic haut
else
X(i)=pas un pic
end
 
end

[MrPingouin]

Bon la réponse était tout simplement dans le premier post.
La boucle ne sert a rien
Désolé d avoir pris de votre temps pour ça.

Par contre ,vos réponses m'ont permis de simplifier pas mal de choses dans cette partie du code et je vous en remercie.

Pour Newenda,je n'ai pas parlé de max mais de pics... Leur détection se fait quelques lignes en amont de ce que je vous ait présenté.
Elle utilise une methode assez classique qui utilise la différence entre un point et ces proches voisins en "lisant" la courbe dans un sens puis dans l'autre puis en faisant la somme des deux lectures.Il n'y a donc pas besoin de seuil dans ce cas.

Merci encore de vos réponses.

[Newenda]

Rha effectivement, j'ai encore lu trop vite