Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: phase-space versus controlling parameter surface

  • To: mathgroup at smc.vnet.net
  • Subject: [mg91733] Re: phase-space versus controlling parameter surface
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Sun, 7 Sep 2008 05:32:28 -0400 (EDT)
  • References: <g9t6r5$jd9$1@smc.vnet.net>

Hi,

your input is a bit corrupted but, and a bundel of lines
is not a surface

ode = {x'[t] == v[t],
    v'[t] == -x[t] 3 - k v[t] + B Cos[t]};
init = {x[0] == 1, v[0] == 0};

DuffingLine[kk_?NumericQ, Bb_?NumericQ, {t_, t0_, t1_}] :=
  Module[{sol},
   sol = NDSolve[Join[ode, init] /. {k -> kk, B -> Bb},
             {x[t], v[t]},  {t, t0, t1}];
   ParametricPlot3D[{t, x[t], v[t]} /. sol, {t, t0, t1}][[1]]
   ]

alllines =
   Table[DuffingLine[0.1, 1, {t, 0, 4 Pi}], {k, 0.0001, 1, 0.01}, {b,
     0.0, 2, 0.05}];

alllines =
   Table[DuffingLine[0.1, 1, {t, 0, 4 Pi}], {k, 0.0001, 1, 0.01}, {b,
     0.0, 2, 0.05}];

Graphics3D[alllines]

BTW MathGL3d line illumination will give you the illusion of a surface
with

Get["MathGL3d`"]

MVShow3D[Graphics3D[alllines /. _Hue :> Sequence[]],
  MVNewScene -> True]


Regards
   Jens

Luca Petrone wrote:
> Dear All,
> 
> I am interested in plotting a surface of the phase-space versus a controlling parameter, e.g. the B in a Duffing's equation
> 
> x'[t] == v[t]v'[t] == - x[t]^3 - k v[t] + B Cos[t]
> 
> that is=2C in the space {x[t]=2C v[t]=2C B} for a known k and B varying from Bmin to BmaxI tryed something like :
> 
> ParametricPlot3D[ Evaluate[{x[t]=2C v[t]=2C B} /.    NDSolve[{v'[t] == - x[t]^3 - k v[t] + B Cos[t]=2C            x'[t] == v[t]=2C x[0] ===
>  1=2C v[0] == 0}=2C {x=2C v}=2C {t=2C 0=2C 2000}=2C     MaxSteps -> Infinity]  ]=2C {t=2C 1950=2C 1950 + 4 Pi}=2C {B=2C 0.2=2C 0.6}]
> 
> but without success.Is there any way to get it ?
> Thank you very much for your help.
> 
> Yours=2C
> Luca P.Milano - Italy
> _________________________________________________________________
> 


  • Prev by Date: How to do efficient NDSolve Error control
  • Next by Date: Re: An Argument Problem
  • Previous by thread: Re: Re: phase-space versus controlling parameter surface
  • Next by thread: Re: phase-space versus controlling parameter surface