Que fait cette ligne de code ?

Version imprimable

15/07/2008, 13h05
ANOVA

Que fait cette ligne de code ?

Bonjour,

je dois traduire un algorithme MATLAB en R. Seulement, je ne connais absolument pas Matlab et je bloque sur cette instruction :

Citation:

V(:,iter+1)=max(aux,2)

Que fait elle ? Pouvez vous m'expliquer plus particulièrement max(aux,2), car j'ai trouvé ce que signifie V(:,iter+1) (on prend le iter+1 ième vecteur colonne de la matrice V).

V est une matrice
aux est une matrice
iter est une variable qui compte les itérations d'une boucle for

Le commentaire à côté de cette ligne de code est :

Citation:

% take the max value over vij ignoring the missing values

Merci de votre aide.

PS : pour info, cet algorithme vient de la page 11 de ce document : matlab
15/07/2008, 13h22
Caro-Line

La version de MATLAB utilisée est-elle indiquée ?

Car c'est étrange dans ma version, max(aux,2) renverra une matrice de la taille de aux, avec les valeurs inférieures à 2 remplacées par 2 (en gros pour chque élément il prend le max entre aux et 2).
Donc là tu essayerais de faire rentrer une matrice dans un vecteur colonne (puisque ton interprétation est exacte).
Donc pour moi ce code ne devrait pas fonctionner...

Et ce ne serait pas plus simple d'avoir l'algo en algo plutôt que de passer par du MATLAB ?
15/07/2008, 13h26
phryte
Salut.

Citation:

max(aux,2) renverra une matrice de la taille de aux, avec les valeurs inférieures à 2 remplacées par 2 (en gros pour chque élément il prend le max entre aux et 2).

Caro a raison. Exemple :
Code:

1 2 3 4 V=zeros(5,5); iter=1; aux=[0 1 2 3 4]'; V(:,iter+1)=max(aux,2)
Résultats :
Code:

1 2 3 4 5 6 7 V = 0 2 0 0 0 0 2 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 0 0 0
15/07/2008, 13h38
ANOVA

Citation:

Envoyé par caro95470

Et ce ne serait pas plus simple d'avoir l'algo en algo plutôt que de passer par du MATLAB ?

Si je pouvais avoir l'algo plutôt que du code MATLAB, ça m'arrangerait, seulement, le seul document dont je dispose est celui ci :

matlab

et l'algo d'où vient la ligne de code qui pose problème est page 11.

Citation:

Envoyé par caro95470

La version de MATLAB utilisée est-elle indiquée ?

Non, la version de Matlab n'est pas indiquée, mais le document date de 2003...en plus, j'ai même pas matlab, donc je ne peux pas faire de tests moi même...

Je me permet de faire remonter ce sujet pour demander un petit service à quelqu'un qui aurait matlab. Je ne dispose pas de Matlab (malheureusement). Je ne peux donc pas faire les tests seul. Est ce que quelqu'un pourrait exécuter ce code et me dire s'il fonctionne et ce qu'il fait ?

Code:

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
29
30
31
32
33
34
35
36
37
38
39
clear
MaxIter=30;
beta=0.9;
K=0:0.01:1;
[rk,dimK]=size(K);
V=zero(dimK,MaxIter);
 
for iter=1:MaxIter
    iter
    aux=zeros(dimK,dimK)+NaN;
 
    for ik=1: dimK
        for ik2=1:(ik-1)
            aux(ik,ik2)=log(K(ik)-K(ik2))+beta*V(ik2,iter);
        end
    end
 
V(:,iter+1)=max(aux,2);
 
end
[Val,Ind]=max(aux,2);
 
optK=K(Ind);
optK=optK+Val*0;
optC=K'-optK';
 
figure(1)
plot(K,V);
xlabel('Size of Cake');
ylabel('Value Function');
 
 
% plot optimal consumption rule as a function of cake size
figure(2)
plot(K,[optC K'],'LineWidth',2)
xlabel('Size of Cake');
ylabel('Optimal Consumption');
text(0.4,0.65,'45 deg.line','FontSize',18)
text(0.4,0.13,'Optimal Consumption','FontSize',18)

Merci.

16/07/2008, 09h57
Caro-Line

Ben...il plante :aie:

Ligne 6 : il faut remplacer zero par zeros :

Code:

V=zeros(dimK,MaxIter);

Et (comme je l'avais prédit, je vais me reconvertir dans la voyance :roll:), là :

Code:

V(:,iter+1)=max(aux,2);

donne :

Citation:

??? Subscripted assignment dimension mismatch.

Pour les explications voir mon post #2.

Je suppose (avec un degré d'incertitude maximum cependant) que cela vient d'une mauvaise utilisation de la fonction MAX.
Il faudrait voir la documentation de cette fonction dans la version de MATLAB utilisée par l'auteur (2003 si on regarde ici :Quelles sont les différentes versions de MATLAB ?, on peut supposer que c'est la r13..., ayant seulement accès à la documentation de la R13SP2, ben l'aide de MAX est la même)

Ou alors comprendre ce qu'il a voulu faire avec son max.
16/07/2008, 15h56
ANOVA

Effectivement, cela vient peut être d'un problème de version. En allant sur ce site (http://www.mathworks.com/access/help...x+matlab&meta=), on trouve :

Citation:

C = max(A,[],dim) returns the largest elements along the dimension of A specified by scalar dim. For example, max(A,[],1) produces the maximum values along the first dimension (the rows) of A.

Donc je pense que ce qu'il veut faire avec cette instruction, c'est récupérer la valeur maximale par rapport aux colonnes. Du coup, avec la dernière version de matlab, la syntaxe serait la suivante :

Code:

V(:,iter+1)=max(aux,[],2);

Est ce que vous pouvez tester le code avec cette modification et me dire si ça marche svp ?

Merci beaucoup à caro95470 pour sa rapidité.:)
16/07/2008, 16h28
Caro-Line

J'avais déjà essayé mais ça plante à cette ligne là :

Code:

optK=optK+Val*0;

Citation:

??? Error using ==> plus
Matrix dimensions must agree.

Error in ==> Base2 at 24
optK=optK+Val*0;

car optK est de dimension 1*101 et Val 101*1...
Donc je modifie ainsi :

Code:

optK=optK+Val'*0;

Les graphiques obtenus sont en pièces jointes.
17/07/2008, 09h05
ANOVA

Merci pour les graphiques, ça me donne une idée de ce à quoi je dois aboutir. Dernières questions (probablement) : que signifie cette ligne de code ?

Code:

optK=optK+Val'*0;

Le prime, c'est pour transposer le vecteur, mais la multiplication par 0, j'imagine que ça ne multiplie pas par 0 sinon, je ne vois pas l'intérêt de l'instruction.

Et que fait celle ci ?

Code:

[Val, Ind]=max(aux,[],2)

(là, j'ai aucune idée...)
17/07/2008, 09h20
Caro-Line

Oui c'est un peu tordu comme code, je n'avais même pas essayé de comprendre en fait.

Alors je commence par ta dernière question :

Code:

[Val,Ind]=max(aux,[],2);

Comme pour la ligne avec V(:,iter+1) on récupère les max (dans Val) et en plus on récupère les indices de ces maximums.
C'est-à-dire qu'on a pour un i donné Val(i)=aux(i,Ind(i)).

Ensuite :

Code:

optK=optK+Val'*0;

Sisi c'est bien ce que tu penses sauf que la fourberie ici vient du fait que dans Val on a des NaN.
Or : 0*qqch=0 sauf :
- 0*NaN=NaN (NaN est plus fort que 0 !!!)
- 0*(+/-)inf = NaN (aussi même si ici ça ne t'intéresse pas)

Du coup ce que ça fait c'est que on conserve optK sauf qu'aux indices où Val=NaN, on met NaN aussi.