Discussion :

[Audpidou911]

Bonjour, je dois tracer une courbe avec les valeurs d'une variable placée dans une boucle, mais quand je plot, il ne me retourne que la dernière valeur calculée.
Comment faire pour garder en mémoire toutes valeurs précédentes ?

Voici le code:

Code :

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

1
2
3
4
 
        xm_nvo=(sum (x_v*w))/sum(w);
        xp=xm_nvo;
        plot(xm, 'r')

[duf42]

Bonjour,

Le code que tu nous montre ne nous aide pas beaucoup puisque tu traces la variable 'xm' et on ne voit pas comment elle est calculée...

Duf

[Audpidou911]

Oups pardon !

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
VT=0.5
N=100
x_v est un vecteur 1*3
w est un vecteur 3*1
 
for i=1:N;
    if d(i)<VT;
        xm=(sum (x_v*w))/sum(w);
        xp=xm;
    else
        fprintf ('output (%d) failed\n',i);
        w(i,:)=[];
        x_v(:,i)=[];
        N=N-1;
        xm_nvo=(sum (x_v*w))/sum(w);
        xp=xm_nvo;
    end
end
 
plot (xp, 'r')

Voilà je crois que j'ai tout dit !

[oesterlj]

Place à l'aide d'une boucle, toutes tes valeurs dans un vecteur, tu pourras ensuite tracer correctement ta courbe

[Audpidou911]

C'est pour faire ça que j'ai tout mis dans xp, sauf que les valeurs sont écrasées à chaque itération.

[Ptinéwik]

Salut,

C'est pour faire ça que j'ai tout mis dans xp, sauf que les valeurs sont écrasées à chaque itération.

Oui effectivement tu ecrases xp à chaque fois, puisque tu redéfinis xp à chaque itération.
Pour incrémenter xp au fur et à mesure de ta boucle, il faut que tu précise que le nouveau xp se mets à la ligne correspondant à ton i:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
VT=0.5
N=100
x_v est un vecteur 1*3
w est un vecteur 3*1
 
cp=0;
for i=1:N;
cp=cp+1;
    if d(i)<VT;
        xm=(sum (x_v*w))/sum(w);
        xp(cp,:)=xm;
    else
        fprintf ('output (%d) failed\n',i);
        w(i,:)=[];
        x_v(:,i)=[];
        N=N-1;
        xm_nvo=(sum (x_v*w))/sum(w);
        xp(cp,:)=xm_nvo;
    end
end
 
plot (xp, 'r')

[Audpidou911]

J'ai testé ce code, et je ne sors rien du tout !
Et le graphe est vide...

[duf42]

Que vaut la variable xp à la fin de ton code?

[Merel]

Bonjour,

Une autre question en passant concernant ton code. Comment peux-tu vider un vecteur qui n'existe pas ?

w et x_v sont des vecteurs 3*1 et 1*3. Or dans ta boucle, tu vas de i=1 à i=100.

Donc

Code :

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

1
2
w(i,:)=[];
x_v(:,i)=[];

n'ont aucun sens.


Autre point : ton premier sum ne sert à rien. La multiplication d'un vecteur ligne avec un vecteur colonne donne un scalaire.

Code :

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

xm=x_v*w/sum(w);

[Audpidou911]

Mon xp devrait être un vecteur formé par tous les résultats issus du calcul xm_nvo= (...).

Je vous envoie le programme en entier, ce sera peut-être plus compréhensible que des explications en toutes lettres (enfin j'espère )

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
[x]= xlsread ('Classeur1.xls');
 
w= [1; 1; 0.85];
VT=2;
M=size(x,1);
t=0:0.04:10;
 
for k=1:M; t=1:0.1:100;
x_v=x(k,:);
d=sort(abs(diff([x_v x_v(:,1)],1,2)));
N=size(d,2);
end
 
for i=1:N;
    if d(i)<VT;
        xm=(sum (x_v*w))/sum(w);
        xp=xm;
    else
        fprintf ('output (%d) failed\n',i);
        w(i,:)=[];
        x_v(:,i)=[];
        N=N-1;
        xm_nvo=(sum (x_v*w))/sum(w);
        xp=xm_nvo;
    end
end
 
plot (xp, 'r')

Classeur1 est un tableau (M*N)
M=100
N=3

En fait, ce que je veux faire, c'est une moyenne pondérée des coefficients pour chaque ligne de ma matrice.

w représente le poids de chaque sortie.

[Merel]

Le programme en entier ne nous aide pas beaucoup si tu nous ne donne pas plus d'information.

Qui a-t-il dans ta variable x ? Quelle est sa taille ? ordre de grandeur des valeurs ?

si tu peux nous mettre quelques lignes de ce fichier ca serait pas mal.

Et surtout que cherches-tui à faire. Un code c'est bien beau mais si on ne sait pas ce que tu veux faire, ca nous aide pas. Enfin si on voit ce que tu veux faire mais ca serait bien que tu le précises avant qu'on ne passe quelques minutes à chercher ;-)

PS: utilise la balise code sur ton/tes messages et indente ton code ca t'aidera...

edit : t'es devin t'as répondu à mes questions avant mon post !

[Merel]

Heu c'est quoi ton histoire d'output failed et de sort ?

Si tu veux faire la moyenne pondéré de tes 3 coef, t'as pas besoin de tout ca ou alors j'ai pas compris ce que tu cherches à faire.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
%x=fix(rand(100,3)*10);
%ou
x= xlsread ('Classeur1.xls');
 
w=[1 1 0.85]';
 
for i=1:100
    xp(i)=sum(x(i,:)*w)/sum(w);
end
 
plot(xp)

[Audpidou911]

En fait les sorites doivent toutes tourner autour d'une certaine valeur, ici 12.
Donc dès qu'une valeur s'en éloigne trop, je dois déclarer que la sortie est erronée et qu'on l'exclut pour le reste des calculs.
Par exemple, si à l'instant t=3 (donc 3e ligne), la sortie 2 a pour valeur 15, l'algorithme doit la déclarer défaillante et dire que désormais, le poids de la sortie est 0 et ne plus considérer que les sorties 1 et 3.

[Merel]

Ha ben voilà, c'est un peu plus clair maintenant.

tu te compliques la vie avec ton algorithme.

Il suffit de regarder si un des éléments de ta ligne est plus grand que 12. Si c'est le cas, tu passes w à 0 et non pas à [] qui est l'élément vide.

Si il n'y a pas de valeur éronnée, tu gardes le tout.

Voici un exemple avec des valeurs bidons.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
x=[5 7 9;10 9 7;15 9 8;9 15 7; 6 8 15];
M=size(x,1);
 
for i=1:M
 w=[1 1 0.85]';
 
        if sum(x(i,:)>12)<0
 
        xp(i)=sum(x(i,:)*w)/sum(w);
 
        else
        [num ind]=find(x(i,:)>12);
 
        w(ind)=0;
        xp(i)=sum(x(i,:)*w)/sum(w);
 
        end
end

[Audpidou911]

Tu es magique !!
Merci pour tout !!

[Merel]

Pas de quoi mais essaye de bien comprendre le code ;-)
C'est plus important que tout...

[Jerome Briot]

Voici un exemple avec des valeurs bidons.

La même chose en une ligne :

Code :

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

1
2
3
4
x = [5 7 9;10 9 7;15 9 8;9 15 7; 6 8 15];
w = [1;1;0.85];
 
xp = x.*(x<=12)*w./sum((x<=12)*w,2);

[Merel]

Le dernier sum est même inutile puisqu'il s'agit d'un calcul matriciel.

Code :

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

1
2
3
4
x = [5 7 9;10 9 7;15 9 8;9 15 7; 6 8 15];
w = [1;1;0.85];
 
xp = x.*(x<=12)*w./((x<=12)*w);

[Jerome Briot]

Le dernier sum est même inutile puisqu'il s'agit d'un calcul matriciel.

C'était juste pour voir si tu suivais