Probleme pour générer un graph sur Matlab

[invitee0cc8045]

Bonjour,

Je tente de générer un graphique sur matlab, le but est de montrer l'instabilité des prix via une fonction:

j'en arrive à écrire ça:

X=[0, 2.5];
P(:,1)=0.1;
N=length(X);

for n=1:N
    for t=1:10
    P(n,t+1)=X(n)*P(n,t)*[1-P(n,t)];
    end
end
plot(t,P)

La machine me répond quant à elle ceci :

??? Attempted to access P(2,1); index out of bounds because
size(P)=[1,101].

Error in ==> test at 7
    P(n,t+1)=X(n)*P(n,t)*[1-P(n,t)];

Ce que biensur je ne saisis pas, à la fin j'aimerais tendre vers ça lien image: http://hpics.li/57b346a

Si vous avez des conseils je vous en remercierai !

Thomas
-----

[invite15928b85]

Bonjour.

Déclarez les dimensions de vos matrices complètement avant emploi. Cela se passera mieux.

Exemple : P = zeros[2,101]; // avant les boucles for

Matlab peut redimensionner dynamiquement mais il y a des limites.

[invitee0cc8045]

Bonsoir,

Merci pour cette réponse,

Finalement, j'ai en effet compris mon erreur, j'ai procédé comme ceci :

X=[1.5, 3.1];
P = 0.5 * ones(N,101);
N=length(X);

for n=1:N
    for t=1:100
    P(n,t+1)=X(n)*P(n,t)*[1-P(n,t)];
    end
end
plot(1:101,P);

Merci en tout cas, j'espère que d'autres pourront apprendre de mes erreurs.

A tout hasard, auriez vous une idée de comment manipuler le Subplot ?

Cordialement.

Thomas

[invitea29b3af3]

Salut

Dans ton code la déclaration de N doit bien évidemment venir avant celle de P
A noter que l'usage des [ ] au lieu des ( ) dans ta boucle for ne génère pas d'erreur (en réalité du crées une matrice de taille 1x1... donc un scalaire), mais des parenthèses suffisent. Et si tu veux pouvoir un jour réutiliser ce code mais avec un signal P plus long, ce serait bien que tu n'aies qu'à changer le 101 dans sa déclaration, et pas devoir ajuster partout. Donc tu peux utiliser size(P,2) à la place. C'est juste des petites remarques, fais-en ce que tu veux

X=[1.5, 3.1];
N=length(X);
P = 0.5 * ones(N,101);

for n=1:N
    for t=1:size(P,2)-1
        P(n,t+1)=X(n)*P(n,t)*(1-P(n,t));
    end
end
plot(1:size(P,2),P);

Pour subplot, j'imagine que tu veux séparer les 2 courbes ?
Par exemple si tu voulais afficher 6 graphes, 3 par ligne (donc 2 "colonnes de graphes") dans une fenêtre, genre ça:

1 graphe là            1 graphe là            1 graphe là
1 graphe là            1 graphe là            1 graphe là

tu utilises: subplot(NOMBRES_DE_LIGNES, NOMBRE_DE_COLONNES, NUMERO_DU_GRAPHE). Dans ce cas tu veux NOMBRES_DE_LIGNES = 2, NOMBRES_DE_COLONNES = 3, et NUMERO_DU_GRAPHE va de 1 à 6, puisque tu as 6 graphes. Donc :

figure
subplot(2,3,1), plot(signal_1)
subplot(2,3,2), plot(signal_2)
subplot(2,3,3), plot(signal_3)
subplot(2,3,4), plot(signal_4)
subplot(2,3,5), plot(signal_5)
subplot(2,3,6), plot(signal_6)

L'ordre des graphes se fait par ligne, donc 1 à 3 sur la 1ère ligne et 4 à 6 sur la 2ème ligne dans cet exemple.

Dans ton cas, disons qu'on veut 2 lignes et 1 colonne.

figure
subplot(2,1,1), plot(P(1,:)), axis tight
subplot(2,1,2), plot(P(2,:)), axis tight

Le "figure" sert juste à ce qu'il t'ouvre une nouvelle figure pour faire ce plot, au cas où t'avais déjà un autre graphe d'affiché.
axis tight sert à resserré les axes sur 1:101. Enfin... essaie de l'enlever et tu verras ce que je veux dire
A noter que tu peux omettre la base de temps 1:size(P,2) dans "plot" quand l'incrément est 1 (puisque c'est la valeur par défaut)

[A voir en vidéo sur Futura]