Numerical solution to coupled diff. eq.

mrandersdk

I have three equations

[tex] (\frac{\partial}{\partial t} + c \frac{\partial}{\partial z}) E_p(z,t) = i N_r(z) \sigma_{ba}(z,t)[/tex]

[tex] \frac{1}{i E_c^*(z,t)}(\frac{\partial}{\partial t} + \Gamma_{bc}) \sigma_{bc}(z,t) = \sigma_{ba}(z,t)[/tex]

[tex] -\frac{E_p(z,t)}{E_c(z,t)} + \frac{1}{i E_c(z,t)}(\frac{\partial}{\partial t} + \Gamma_{ba}) \sigma_{ba}(z,t) = \sigma_{bc}(z,t)[/tex]

where the functions [tex]E_p(z,t)[/tex] and [tex]N_r(z)[/tex] is know, and all other things that do not depend on time or position are known constants.

How can I solve this numericaly with some give initial conditions. I have matlab 7.0 or mathematica to my disposal but can't make anything work.

Chris Hillman

Did you try finding an exact solution? And are these complex functions or what? Is * Hermitian transpose?

mrandersdk

yes the function can be complex, but are of real variable, so not need for complex function theory.

the 'i' in the first equation is the complex i.

Do you think it is possible to find and exact solution, that would of cause be great.

Rainbow Child

It is not really too hard to find an exact solution if [tex]E_p(z,t)[/tex] and [tex]N_r(z)[/tex] are known.

From the first equation you have [tex] \sigma_{ba}(z,t) [/tex] in terms of the known funtions [tex]E_p(z,t)[/tex] and [tex]N_r(z)[/tex].

From the 3nd equation you have [tex] \sigma_{bc}(z,t) [/tex] in terms of the known funtions [tex]E_p(z,t),\,N_r(z)[/tex] and the unknown function [tex] E_c(z,t) [/tex].

Plugging these informations into the 2nd equation you have a PDE which involves [tex] E_c(z,t) \, and \, E^*_c(z,t) [/tex], say it (A).

• If you are looking for real solutions, i.e. [tex] E_c(z,t) =E^*_c(z,t) [/tex] then equation (A) is just an ODE with respect to [tex]E_c(z,t) [/tex] since it involves only the derivative [tex]\partial_t E_c(z,t) [/tex]. It looks like
  
  
  [tex] \partial_t E_c(z,t) \sim \alpha_1\,E_c(z,t) +\alpha_2\,E^3_c(z,t) [/tex]
  
  where [tex] \alpha_i [/tex] are known funtions of [tex] (z,t) [/tex]. It can be full integrated either by hand or with the help of Mathematica.
  
• If [tex]E_c(z,t) [/tex] is a complex funtion you have to split every term in (A) at it's Real and Imaginary part, in order to end up with two [I]messy[/] DE. Fortunately, in this case too you can integrate the resulting system, by imposing the integrabity conditions on [tex] Re[E_c(z,t)]\, and\, Im[E_c(z,t)] [/tex]

I hope that was useful for you!

momentum_waves

Wave phenomena are predicted. Thankfully, at this point, they look to be linear (if gamma & c are constant).

If so, then you may be reasonably safe in using simple numeric strategies.