MathGroup Archive 2011

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

Search the Archive

Help with non-linear system of equations (Trigonometric and Polynomial)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg115578] Help with non-linear system of equations (Trigonometric and Polynomial)
  • From: Manjinder Benning <manjbenning at gmail.com>
  • Date: Sat, 15 Jan 2011 04:45:29 -0500 (EST)

Hi i am trying to do some 3D modeling by finding a 6 unknown 6 equation
simultaneous system. I am using the equation:

Cos(a) = A.B/|A||B|
where A and B are lines (vectors) in 3 space.
a is the angle between the lines.

When i run the calculation, Mathematica just runs forever and doesn't come
to a solution. It would be great if a closed form solution could be found.
I am understanding that it is difficult to solve systems with polynomial and
trig parts?

Description of Problem:

So i have 3 points in 3-space.
Pl (Xl, Yl, Zl) - A constant point in the air
Pc(Xc, Yc, Zc) - A constant point on the ground (Z=0)
Pe's (Xe_n, Ye_n, Ze_n) - A series of points that lay on the edges of an
ellipse co-planar to Z=0 (which are known experimentally)
A - The angle between the line (Pl -> Pc) and each (Pl -> Pe_n)

Due to the geometry of the problem (a ellipse) angle A is constant.


Are there any methods i can use to make it easier to solve (ie bounding the
problem) or am i dreaming that i will find a closed form solution?
...any help would be appreciated:

Heres the code:

eq1 =Cos[A] ==  ((Xc - Xl)*(Xe1-Xl) + (Yc-Yl)*(Ye1-Yl) + (-Zl)*(-Zl))/Sqrt[(
(Xc - Xl)^2 + (Yc-Yl)^2+(-Zl)^2 ) * ((Xe1-Xl)^2 +(Ye1-Yl)^2+(-Zl)^2)]

eq2 =Cos[A] ==  ((Xc - Xl)*(Xe2-Xl) + (Yc-Yl)*(Ye2-Yl) + (-Zl)*(-Zl))/Sqrt[(
(Xc - Xl)^2 + (Yc-Yl)^2+(-Zl)^2 ) * ((Xe2-Xl)^2 +(Ye2-Yl)^2+(-Zl)^2)]

eq3 =Cos[A] ==  ((Xc - Xl)*(Xe3-Xl) + (Yc-Yl)*(Ye3-Yl) + (-Zl)*(-Zl))/Sqrt[(
(Xc - Xl)^2 + (Yc-Yl)^2+(-Zl)^2 ) * ((Xe3-Xl)^2 +(Ye3-Yl)^2+(-Zl)^2)]

eq4 =Cos[A] ==  ((Xc - Xl)*(Xe4-Xl) + (Yc-Yl)*(Ye4-Yl) + (-Zl)*(-Zl))/Sqrt[(
(Xc - Xl)^2 + (Yc-Yl)^2+(-Zl)^2 ) * ((Xe4-Xl)^2 +(Ye4-Yl)^2+(-Zl)^2)]

eq5 =Cos[A] ==  ((Xc - Xl)*(Xe5-Xl) + (Yc-Yl)*(Ye5-Yl) + (-Zl)*(-Zl))/Sqrt[(
(Xc - Xl)^2 + (Yc-Yl)^2+(-Zl)^2 ) * ((Xe5-Xl)^2 +(Ye5-Yl)^2+(-Zl)^2)]

eq6 =Cos[A] ==  ((Xc - Xl)*(Xe6-Xl) + (Yc-Yl)*(Ye6-Yl) + (-Zl)*(-Zl))/Sqrt[(
(Xc - Xl)^2 + (Yc-Yl)^2+(-Zl)^2 ) * ((Xe6-Xl)^2 +(Ye6-Yl)^2+(-Zl)^2)]

Solve[{eq1,eq2,eq3,eq4,eq5,
eq6},{A,Xl,Yl,Zl,Xc,Yc}]

much gratitude to all things Mathematica!
Manjinder.


-- 
www.thepowerofyang.com


  • Prev by Date: Re: Units question
  • Next by Date: Re: Caution on BibTex references with Mathematica
  • Previous by thread: Help with non-linear system of equations (Trigonometric and Polynomial)
  • Next by thread: Re: Help with non-linear system of equations (Trigonometric and Polynomial)