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

MathGroup Archive 2008

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

Search the Archive

Solving algebraic equations with some conditions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg92007] Solving algebraic equations with some conditions
  • From: Kay-Michael Voit <kay at voits.net>
  • Date: Tue, 16 Sep 2008 19:26:36 -0400 (EDT)

Hello,
I am trying to solve a system of 3 PDEs using a Fourier ansatz:
nL[x_, t_] =
  Sum[cn[p, l] Exp[I (p k x + l \[Omega] t)], {p, -1, 1}, {l, -1, 1}];
EscL[x_, t_] =
  Sum[cE[p, l] Exp[I (p k x + l \[Omega] t)], {p, -1, 1}, {l, -1, 1}];
iL[t_] = Sum[ci[l] Exp[I l \[Omega] t], {l, -1, 1}];

So I have a system of up to 55 algebraic equations (Not every Exp[I (p k
x + l \[Omega] t)] occures after inserting the ansatz, thus some will be
trivial) for the 55 Fourier coefficient. The equations are not linear,
but linear in every coefficient. I was told that this would take some
time but should be possible.

Additionally, I know that the functions are all real, thus I get more
equations:
K[p,l]=Conjugate[K[p,l]]

Unfortunately, I have the impression that Mathematica cannot handle
these in Solve[], as simple examples produce things like Conjugate^-1[I].
A possibility I considered is constructing the coefficients with real
and imaginary part:

iL[t_]
  = Sum[(ci[l, re] + ci[l, im] I) Exp[I l \[Omega] t], {l, -1, 1}];

Then I had to give Mathematica assumtions that these are all real so
that an expression like

-H Wges+(I cn[0,0,im]+cn[0,0,re])/\[Tau]==0

will be correktly solved to

{cn[0,0,im] -> 0, cn[0,0,re]->\[Tau] H Wges}
(also assuming that \[Tau], H, Wges are real)
I am not sure how to achieve that.

Indeed I have another condition for the coefficients:
K[p,l][\[omega]]=K[p,-l]][-\[omega]]
In this case, I have completely no idea how to model this in Mathematica.

Can anybody give me a hint (on the problem or on literature) how to
solve this problem?

Thanks in advance
Kay-Michael Voit


  • Prev by Date: Re: Advanced plotting
  • Next by Date: Re: Apparent error integrating product of DiracDelta's
  • Previous by thread: Re: Problem in PL spectrum analysis
  • Next by thread: Re: Solving algebraic equations with some conditions