MathGroup Archive 2010

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

Search the Archive

Re: NDSolve and hybrid dynamics (Differential Algebraic Equation DAE)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg113556] Re: NDSolve and hybrid dynamics (Differential Algebraic Equation DAE)
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Tue, 2 Nov 2010 05:02:50 -0500 (EST)

The equations are not the same if you change y[t] to y[x[t]]. Further, if you want the break in y[t] when x[t] is 0.1 then that must be the boundary condition on x rather than x[-Pi] == 0.5

eqns = {x'[t] == y[t] + Sin[t],
   y[t] == Piecewise[{{1, t < t1}}, 2],
   x[t1] == 1/10};

sol = {x[t], y[t]} /.
   DSolve[eqns, {x[t], y[t]}, t][[1]] //
  Simplify

{1/10 - 2*t1 + Cos[t1] + Piecewise[{{t - Cos[t], t < t1}}, 2*t - Cos[t]], 
   Piecewise[{{1, t < t1}}, 2]}

Simplify[sol, t == t1]

{1/10, 2}

Solve[(Simplify[sol[[1]], t < t1] /. t -> t1) == 1/10, t1][[1]]

{t1 -> 0}

sol = Simplify[sol /. t1 -> 0]

{11/10 + Piecewise[{{t - Cos[t], t < 0}}, 2*t - Cos[t]], 
 Piecewise[{{1, t < 0}}, 2]}

Plot[sol, {t, -5, 5},
 Frame -> True,
 Axes -> False,
 Exclusions -> 0,
 ExclusionsStyle -> Red]


Bob Hanlon

---- Ravi Balasubramanian <ravi.balasubramanian at yale.edu> wrote: 

=============

Thank you, Bob.  I guess I oversimplified the original problem.  The 
variable y is state-dependent such as
y[x] := Piecewise[{{1, x < 1/10}}, 2]

I tried the following system of equations, but it failed:

eqns = {x'[t] == y[t] + Sin[t],
    y[t] == Piecewise[{{1, x[t] < 1/10}}, 2], x[-Pi] == 1/2};
sol = {x[t], y[t]} /. NDSolve[eqns, {x[t], y[t]}, {t, -5, 5}]

NDSolve::ndsz: At t == -3.33251, step size is effectively zero; 
singularity or stiff system suspected.

I presume this is what you meant by breaking the problem into two 
regions, but that transition is dependent on x.

With your suggestion of using Piecewise, I was able to run the same 
system by creating a function y

y//Clear
y[xVar_?NumericQ] := Piecewise[{{1, xVar < 1/10}}, 2]
eqns = {x'[t] == y[x[t]] + Sin[t], x[-Pi] == 1/2};

xSoln = NDSolve[eqns, {x[t]}, {t, -5, 5}][[1, 1, 2]]
Plot[xSoln, {t, -5, 5}, Frame -> True, Axes -> False,
  PlotStyle -> Thick]

Plot[y[xSoln], {t, -5, 5}, Frame -> True, Axes -> False,
  PlotStyle -> Thick]

Is there a way that I can get NDSolve to manage the y variable for me as 
part of the state because my system will have multiple such non-smooth 
variables?

Ravi

Bob Hanlon wrote:
> If is a programming construct not a mathematical function. Use Piecewise
> 
> eqns = {x'[t] == y[t] + Sin[t],
>    y[t] == Piecewise[{{1, t < 1/10}}, 2],
>    x[-Pi] == 1/2};
> 
> sol = {x[t], y[t]} /. DSolve[eqns, {x[t], y[t]}, t][[1]]
> 
> {(1/2)*(-1 + 2*Pi + 2*Piecewise[{{t - Cos[t], t <= 1/10}}, 
>             -(1/10) + 2*t - Cos[t]]), Piecewise[{{1, t < 1/10}}, 2]}
> 
> Plot[sol, {t, -5, 1},
>  Frame -> True,
>  Axes -> False,
>  Exclusions -> 1/10,
>  ExclusionsStyle -> Red]
> 
> For NDSolve break the problem into two regions
> 
> 
> Bob Hanlon
> 
> ---- Ravi Balasubramanian <ravi.balasubramanian at yale.edu> wrote: 
> 
> =============
> Friends,
> 
> I am trying to solve a hybrid dynamics problem using NDSolve.
> 
> I am aware that I could use the event locator methods in NDSolve and 
> manually take care of the transitions (set initial conditions etc).  But 
> my system is complicated, and I would have to do a lot of book-keeping.
> 
> The following represents a simplified version of the problem.
> 
> eqns = {x'[t] == y[t] + Sin[t], y[t] == If[t < 0.1, 1, 2],
>     x[-Pi] == 1/2};
> sol = NDSolve[eqns, {x[t], y[t]}, {t, -5, 5}]
> 
> It is a first order differential equation in x[t], but it includes a 
> y[t] variable that is discontinuous.  There should be solution for this 
> problem, but I get the following error message from NDSolve.
> 
> LinearSolve::"nosol" : "Linear equation encountered that has no solution.
> 
> The error is probably arising from the "If" condition.
> 
> FYI, in my real system the y[t] variable is a Lagrange multiplier that 
> is enforcing a constraint.  Under certain conditions, that constraint is 
> broken and the system transitions into another state.  It will be great 
> if this problem can be solved with NDSolve keeping track of the y[t] 
> variable for me.  Any suggestions greatly appreciated.
> 
> Thanks!
> 
> Ravi Balasubramanian
> Yale University
> 



  • Prev by Date: Re: Balance point of a solid
  • Next by Date: Re: how to plot nminimized result
  • Previous by thread: Re: NDSolve and hybrid dynamics (Differential Algebraic Equation DAE)
  • Next by thread: Re: NDSolve and hybrid dynamics (Differential Algebraic Equation DAE)