2012-01-18 57 views

Respuesta

17

El comando que necesita (desde la versión 7) es VectorPlot. Hay buenos ejemplos en la documentación.

Creo que el caso que le interesa es una ecuación diferencial

y'[x] == f[x, y[x]] 

En el caso que dio en su pregunta,

f[x_, y_] := y 

que integra a la exponencial

In[]:= sol = DSolve[y'[x] == f[x, y[x]], y, x] 
Out[]= {{y -> Function[{x}, E^x c]}} 

Podemos trazar el campo de pendiente (ver wikibooks:ODE:Graphing) usando

VectorPlot[{1, f[x, y]}, {x, -2, 2}, {y, -2, 2}] 

y

Esto pueden representarse gráficamente con las soluciones a la DE de usar algo como

Show[VectorPlot[{1, f[x, y]}, {x, -2, 2}, {y, -2, 8}, 
    VectorStyle -> Arrowheads[0.03]], 
Plot[Evaluate[Table[y[x] /. sol, {c, -10, 10, 1}]], {x, -2, 2}, 
    PlotRange -> All]] 

y again

Tal vez un ejemplo más interesante es la de Gauss

In[]:= f[x_, y_] := -x y 

In[]:= sol = DSolve[y'[x] == f[x, y[x]], y, x] /. C[1] -> c 
Out[]= {{y -> Function[{x}, E^(-(x^2/2)) c]}} 

Show[VectorPlot[{1, f[x, y]}, {x, -2, 2}, {y, -2, 8}, 
    VectorStyle -> Arrowheads[0.026]], 
Plot[Evaluate[Table[y[x] /. sol, {c, -10, 10, 1}]], {x, -2, 2}, 
    PlotRange -> All]] 

-xy


Por último, no es un concepto relacionado del gradiente de campo, donde se mira el (derivado del vector) gradiente de una función:

In[]:= f[x_, y_] := Sin[x y] 
     D[f[x, y], {{x, y}}] 
     VectorPlot[%, {x, -2, 2}, {y, -2, 2}] 

Out[]= {y Cos[x y], x Cos[x y]} 

Sin[x y]

0

Al parecer, de la demostración se conectó que toma una función f (x, y), pero hay una serie de diferencias. Sin embargo, sabiendo que f(x,y)=y(x)', puede simplemente usar f(x,y)=C*E^x donde x=t. Mis diferenciales pueden estar un poco oxidados, pero estoy bastante seguro de que es correcto.

+0

alguien tiene un trazador de líneas 1 para los campos de pendiente? – user968102

+0

1 liner? Si necesita saber cómo colocarlo, puede consultar http://www.physicsforums.com/showthread.php?t=152157. Es un poco más de una línea, ya que Mathematica aparentemente necesita una biblioteca importada, entonces la llamada es la trama. la función en sí misma es un poco más. –

Cuestiones relacionadas