Discussion :

[arlyk]

Bonjour à tous,

J'ai besoin de comparer 2 vecteurs de taille différentes (ex: size(A,2)=10 et size(B,2)=15), c'est à dire qu'il me faut trouver les valeurs de A qui sont inférieur aux valeurs de B+500.
Lorsque j'utilise la fonction "find"

Code :

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

find (A<=B)

j'obtiens le message d'erreur classique "Matrix dimensions must agree."

Merci d'avance

[Jerome Briot]

Tes explications manquent de précisons.

Que représentent les valeurs contenues dans ces vecteurs ?

[arlyk]

Merci d'avoir regardé aussi rapidement mon message.

Il s'agit de différents temps, A correspondant au moment où le sujet a répondu à l'expérience

Code :

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

1
2
3
4
5
6
7
8
9
10
11
A = Real_Time_KeyIds_left =
 
  1.0e+004 *
 
  Columns 1 through 11
 
    0.2090    0.2670    0.9550    1.0379    1.2235    1.4855    1.6065    1.8734    2.0740    2.1353    2.2563
 
  Columns 12 through 14
 
    2.4138    2.5746    2.9509

et B correspond aux moments des changements de stimuli auxquels devait répondre le sujet.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
B=Real_Time_Cross_change =
 
  1.0e+004 *
 
  Columns 1 through 11
 
    0.1742    0.2322    0.9168    0.9749    0.9981    1.1837    1.4274    1.5551    1.8336    2.0425    2.0773
 
  Columns 12 through 15
 
    2.2281    2.5299    2.9244    2.9940

[Newenda]

c'est à dire qu'il me faut trouver les valeurs de A qui sont inférieur aux valeurs de B+500.

Ceci n'est pas clair.

Trouver plusieurs valeurs d'un vecteur par rapport à plusieurs valeurs d'un autre vecteur n'a pas de sens sans plus de précision.

Tu veux comparer quelles valeurs de A par rapport à quelles valeurs de B ? Car si tu cherches simplement à trouver toutes les valeurs de A inférieurs à B+500, il te suffit de trouver toutes les valeurs de A qui sont plus petites que la plus petite valeur de B+500.

PS: de plus les dimensions que tu donnes au début de ton message de A et B ne correspondent pas à l'exemple que tu donnes de A et B, qui sont simplement des vecteur de 14 et 15 (et non 10 et 15). N'aurais tu pas oublié une valeur dans A ?

[arlyk]

Le sujet doit répondre le plus rapidement possible à des changements de stimuli.
A est un vecteur temps: les différents moments où le sujet a répondu.
B est aussi un vecteur temps: les différents moments où les stimuli ont changé.

A est donc la réponse à B et je veux comparer toutes les valeurs de A à B.

Par exemple la différence A(1)-B(1) = 348.1621ms
Le sujet à donc mis 348 ms pour répondre au changement B(1).

Par contre il y a des distracteurs qui font que le sujet peut signaler un changement alors qu'il n'y en a pas eu et inversement lui faire manquer un changement, c'est pourquoi mes vecteurs A et B ne font pas forcément la même taille.

Ce que je cherche à faire c'est connaitre les différents temps de réaction.
Mon problème si je fais des boucles for est que A(x) ne correspond pas forcément à B(x). Par exemple si le sujet a manqué le 6eme changement de simulus correspndant à B(6), la 6eme réponse du sujet A(6) correspondra donc à B(7) et ainsi de suite.
D'où l'idée de faire B+500 pour vérifier que la réponse A(x) correspondrait bien au changement B(x).

Les dimension que j'ai données au début du message sont aléatoire comme elles peuvent changer d'un sujet à l'autre.

J'espère être un peu plus clair qu'avant et merci de m'avoir répondu .

[Jerome Briot]

Si j'ai bien compris, tu cherches à faire ceci :

Code :

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

1
2
3
4
5
for n = 1:numel(A)
 
   idx(n) = find((B+500)>A(n),1,'first');
 
end

C'est ça ?

[arlyk]

Salut, désolé pour l'absence.
Dut, je ne sais pas trop si c'est ça que je cherche à faire, je ne comprend pas ce que ta ligne de code fait

J'ai écrit ça:

Code :

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

1
2
3
4
5
6
7
8
9
10
Correct = 0;
Incorrect = 0;
for c = 1:size(A,2);
    for cc = 1:size(B,2);
               if A(c)<= B(cc)+500 && A(c)> B(cc)
            Correct = Correct+1
        else Incorrect = Incorrect+1;
        end;
    end;
end;

C'est pas encore bon mais je me rapproche...

[tachmou]

Salut !

Ton problème est intéressant je trouve. Par contre il me semble un peu complexe, et je ne pense pas que passer directement à de la programmation soit le plus avisé...

Je te conseille avant tout de faire une liste de ce que tu cherches. Dis moi si je me trompe, mais je crois que tu cherches :

1) Les éléments de A qui viennent juste après B ("juste après" étant à définir)
2) Les éléments de B qui n'ont pas de suite immédiate dans A ("immédiate" étant à définir)
3) Les éléments de A qui n'ont pas d’antécédent immédiat dans B ("immédiat" étant à définir)

Si c'est à peu près ça... excuse moi, mais le "trouver les valeurs de A qui sont inférieur aux valeurs de B+500" me semble un peu à côté de la plaque, ou alors c'est loin d'être complet ! De plus, même cette décomposition que je viens de te donner me paraît simpliste, mais tout dépend de tes enregistrements. Par ailleurs, si ton enregistrement est "sale", il faudra peut-être penser à faire un "nettoyage"... (j'ai déjà effectué des traitements de ce genre pour analyser le comportement de nuées de robots imitateurs donc je parle en connaissance de cause)

Cordialement,

ps: je te mets +1 pour les deux messages que tu as posté qui donnent vraiment des informations intéressantes pour t'aider. Le reste va trop vite mais c'est bien d'essayer
ps2: la liste que je viens de donner est peut-être - voire c'est presque certain - à modifier/rectifier/améliorer etc.

[arlyk]

Salut,
Merci tachmou, la liste des points que tu m'as donnée est tout à fait juste et peut-être me clarifie aussi ce que je veux faire avec Matlab (je devrais travailler plus sur papier avant de m'attaquer à Matlab).


PS1: je suis flatté que tu trouves mon problème intéressant

PS2: c'est pour quand le robo néo qui fait le ménage ?

[tachmou]

Ok, dans ce cas, tu peux commencer à réfléchir au premier point en le décortiquant. Tiens nous au courant de ton avancée (tu peux passer à de l'algorithmique voire à de la programmation maintenant, mais pose tranquillement les choses, en donnant des noms logiques à tes indices, en créant des petites fonctions qui s'occuperont de faire des tâches basiques etc.)

[shaiHulud]

Bonjour,
Pour comparer (via <) tous les éléments de A à tous les éléments de B

Code :

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

1
2
 
cmp = bsxfun(@le,A(:),B(:)')

renvoie une matrice de comparaison cmp t.q

cmp(i,j) == true <=> A(i) < B(j)

Cela évite (et est équivalent à) la boucle. Ensuite il faut traiter cmp. En général on trouve tout ce qu'on veut en faisant des somme ou des find dans la matrice (mais il faut savoir ce que l'on cherche). Elle permettra de répondre au 3 questions.

[tachmou]

Je ne connaissais pas bsxfun, elle peut être intéressante dans ton cas effectivement, mais elle ne peut pas répondre à ton problème tel quel... il faut tout de même faire quelques traitements supplémentaires... Tu peux essayer de l'utiliser, mais si tu maîtrises mal le calcul vectoriel sous matlab, je te conseille de rester à des fonctionnalités basiques. Il sera toujours temps d'optimiser ton code plus tard.

[arlyk]

Salut,
Je dois passer à autre chose voici pour ceux que ça pourrait aider ma version finale (qui n'est pas très stricte mais suffisante pour ce que je veux en faire)

Code :

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

1
2
3
4
5
6
7
8
9
10
Correct = 0;
Incorrect = 0;
for c = 1:size(A,2);
    for cc = 1:size(B,2);
               if A(c)<= B(cc)+500 && A(c)> B(cc)+100
            Correct = Correct+1
        else Incorrect = Incorrect+1;
        end;
    end;
end;

Merci