Discussion :

[leontief]

Bonjour,

J'ai une question concernant un problème que je n'arrive à résoudre moi même, je me permets donc de solliciter votre aide, si le temps vous le permet.
Je dispose de deux paramètres contribuant au calcul d'un troisième, le premier variant de 1 à 10 par pas de 0.01, le deuxième de 0 à 1 par pas de 0.01.
et j'aimerais effectuer la surface de la matrice définie comme cela :


Cependant je ne vois pas comment obtenir cette matrice.

En vous remerciant d'avance.

[Jerome Briot]

Utilise la fonction meshgrid :

Code :

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

1
2
3
[x,y] = meshgrid(1:1:10,0:0.1:1);
 
[x(:) y(:) x(:).*y(:)]

A toi d'adapter ce code avec les bonnes valeurs correspondant à ton besoin

[leontief]

Je vous remercie,
Néanmoins, en me relisant je pense que je me suis mal exprimé.
voici mon code

Code :

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

1
2
3
4
5
6
7
8
 Pi = delta * Omega * weq';
    Sigma = zeros (2,2); 
    Sigma(1,1)=0.0009;
    Sigma(2,2)=0.0031;
    omegabis = tau*Omega;
    Er = Pi + omegabis * P' * inv(P * omegabis * P' + Sigma) * (Q - P * Pi);
    I = ones(size(Omega,1),1);
    w = ((inv(Omega)*I)/(I'*inv(Omega)*I)) + inv(delta * Omega)*(Er-(I'*inv(Omega)*Er)/(I'*inv(Omega)*I)*I);

Il me faudrait donc faire une boucle faisant le calcul pour les couples (1,0),(1,0.01) etc.. et qu'ils ressortent en tant que vecteurs distincts pour chaque paramètre, et que cela me renvoie la matrice w calculée à partir de ces couples. ( quant à elle composé de 6 lignes ) je n'aurais juste qu'à retirer la première ligne et en faire un vecteur.
ainsi en utilisant reshape, je pourrais utiliser mesh(x,y,z), au final je pourrais disposer de 6 surfaces.
La section "reshape" de ce lien explique bien mon idée. http://fr.mathworks.com/help/matlab/...a-surface.html
Je ne sais pas si je suis très clair, c'est assez compliqué à expliquer et ça m'ennuie énormément, j'ai réussi à faire une boucle mais qui ne me renvoie pas chaque couple possible.

P.S : Les paramètres a bouger étant delta et tau.

[Beltharion]

A priori Dut a raison. La fonction meshgrid te permet de faire une grille de calcul. Si tu boucles sur les éléments de cette grille tu pourras obtenir w pour tout les couples possibles.