MathGroup Archive 2000

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

Search the Archive

Re: Mathematica won't solve simple diff. eqn. system

  • To: mathgroup at smc.vnet.net
  • Subject: [mg24727] Re: Mathematica won't solve simple diff. eqn. system
  • From: changpoohbear at my-deja.com
  • Date: Wed, 9 Aug 2000 02:31:26 -0400 (EDT)
  • References: <8mdlt6$5jt@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com


Christopher,

It appears that the general linear system that you are trying to solve
is in the form of

a i1'[t] + b i2'[t] + c i1[t] + d i2[t] ==0
e i1'[t] + f i2'[t] + g i1[t] + h i2[t] ==0
i1[0] == alpha
i2[0] == beta

Since this is an autonomous system, it also can be represented in
matrix form as

A i'[t] = B i[t]

where

A={{a,b},{e,f}}; B=-{{c,d},{g,h}}; i[t]={{i1[t]},{i2[t]}};

*IF* A is invertible, the solution can be given as

i[t] = MatrixExp[Inverse[A].B].i[0]

where

i[0] = {{i1[0]},{i2[0]}}

In your example given below, I don't believe the solutions given are
correct ({{-4,8},{-4,8}} also is not invertible!).  For your specific
example, observe that

Eliminate[{-4 i1'[t] + 8i2'[t] - 25i1[t] + 20 i2[t] ==
      0,-4i1'[t] + 8 i2'[t] - 10 i1[t] + 40 i2[t] == 0},
      {i1'[t],i2'[t]}]

yields

3 i1[t]==-4 i2[t]

Using this fact, together with the fact that i1[0]=i2[0]=0, one should
obtain

i1[t]=i2[t]=0

I don't know why my version of Mathematica 3.0.x cannot realise this,
and just returns its input as its output when using DSolve.  <Shrug>

I hope this helps!

Mike

In article <8mdlt6$5jt at smc.vnet.net>,
  "Christopher R. Carlen" <crcarle at sandia.gov> wrote:
> Mathematica 4.0 and linear constant coefficient differential
equations:
>
> I have the following system:
>
> -4 i1'[t] + 8 i2'[t] - 25 i1[t] + 20 i2[t] == 0
> -4 i1'[t] + 8 i2'[t] - 10 i1[t] + 40 i2[t] == 0
> i1[0]==0
> i2[0]==0
>
> Which when I try to solve with DSolve, it fails.
>
> It seems any system of the form:
>
> y1'[t] + y2'[t] + C1 y1[t] + C2 y2[t] == 0
> y1'[t] + y2'[t] + C3 y1[t] + C4 y2[t] == 0
>
> can't be solved.  If the coefficients on y1' and y2' are not the same
> between the two equations, then it can be solved.
>
> The problem is that there is a solution to the above system, which I
> have verified.  That solution is:
>
> i1[t_] = 4 + 64 E^(-5 t) - 68 E^(-4 t)
> i2[t_] = 1 - 52 E^(-5 t) + 51 E^(-4 t)
>
> So the question is:  If there is a solution (and not a very difficult
> one) why can't Mathematica find it???  Is there some way to coerce
Mathematica to
> produce the equation, in both the symbolic and numerical situations?
>
> These types of systems arise frequently in the study of electronic
> circuits.  Numerical solvers like SPICE solve them without any
> difficulty.  I have struggled with getting Mathematica to solve them
for a long
> time.  Sometimes I force a numerical solution by perturbing the
> coefficients a bit, as long as the error is acceptible.
>
> But I would like to understand better what the hangup is.  I have had
a
> diff. eqns. course, but haven't gone into systems yet.
>
> Thanks.
> --
> _______________________
> Christopher R. Carlen
> Sr. Laser/Optical Tech.
> Sandia National Labs
>
>


Sent via Deja.com http://www.deja.com/
Before you buy.


  • Prev by Date: Re: Package naming
  • Next by Date: Re: Mathematica won't solve simple diff. eqn. system
  • Previous by thread: Re: Mathematica won't solve simple diff. eqn. system
  • Next by thread: Re: Mathematica won't solve simple diff. eqn. system