Discussion :

[kizaa]

Bonjour,
Besoin d'aide svp..

fonction contour en 2D, j'ai trouvé

Code matlab :

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

1
2
3
4
[x1grid,x2grid] = meshgrid(0:0.05:5);
K=c*x1grid+b*x2grid-d*log(x1grid)-a*log(x2grid); % solution du système..équations diff
v = -5:0.2:5;
contour(x1grid,x2grid,K,v,'-'); %diagramme de phase

j'ai un système 3D, dont K la solution
donc normalement x1grid,x2grid et x3grid, je veux faire la même chose!

Code matlab :

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

contour3 (x1grid,x2grid,x3grid....

merci d'avance

[Kangourou]

essaie isosurface.
A+

[kizaa]

Bonjour,

Code matlab :

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

1
2
3
4
5
 
[x1grid,x2grid] = meshgrid(0:0.05:5)
x3grid=x1grid;
K= (c-x2grid)*x1grid - 0.3*d*x1grid.^3 - 0.5*x2grid.^2 - (val-x3grid +b*x1grid.^2-aa*x1grid.^3)*x2grid + rr*s*(0.5*x1grid.^2-x0*x1grid+(x1grid-x0)*x2grid) - (c+(b-d)*x1grid.^2-aa*x1grid.^3+val+rr*(x1grid+x2grid))*x3grid + 0.5*x3grid.^2;
isosurface(x1grid,x2grid,x3grid,K)

??? Error using ==> isosurface at 74
V must be a 3D array.

qd je mets

Code matlab :

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

1
2
3
 
[x1grid,x2grid,x3grid] = meshgrid(0:0.05:5);
K= (c-x2grid)*x1grid - 0.3*d*x1grid.^3 - 0.5*x2grid.^2 - (val-x3grid +b*x1grid.^2-aa*x1grid.^3)*x2grid + rr*s*(0.5*x1grid.^2-x0*x1grid+(x1grid-x0)*x2grid) - (c+(b-d)*x1grid.^2-aa*x1grid.^3+val+rr*(x1grid+x2grid))*x3grid + 0.5*x3grid.^2;

j'obtiens l'erreur

??? Error using ==> mtimes
Input arguments must be 2-D.

vu que
size(x1grid)= 101 101 101; idem pour x2grid et x3grid

[Jean Dumoncel]

L'opérateur '*' effectue une multiplication matricielle alors qe tu veux certainement faire une multiplication terme a terme. Pour cela utilise l'opérateur '.*'

[kizaa]

résolu
mercii