Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Numerical solution to coupled diff. eq.

  1. Nov 28, 2007 #1
    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.
  2. jcsd
  3. Nov 28, 2007 #2

    Chris Hillman

    User Avatar
    Science Advisor

    Why Numerical?

    Did you try finding an exact solution? And are these complex functions or what? Is * Hermitian transpose?
  4. Nov 29, 2007 #3
    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.
  5. Dec 24, 2007 #4
    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 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! :smile:
  6. Dec 24, 2007 #5
    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.
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook