Aquí hay una manera. Supongo que el -0.8 es precisamente -4/5 aquí.
Primero, use identidades de suma de ángulos para expandir cosas como [math] \ cos (5 \ theta_3) [/ math] en polinomios en [math] \ sin (\ theta_3) [/ math] y [math] \ cos (\ theta_3) [/ math]. Deberías obtener:
[matemáticas] -2 \ cos \ left (\ theta _1 \ right) +2 \ cos \ left (\ theta _2 \ right) -2 \ cos \ left (\ theta_3 \ right) +1 = \ dfrac {-4} {5} [/ matemáticas]
[matemáticas] -2 \ cos ^ 5 \ left (\ theta _1 \ right) +2 \ cos ^ 5 \ left (\ theta _2 \ right) -2 \ cos ^ 5 \ left (\ theta _3 \ right) +20 \ sin ^ 2 \ left (\ theta _1 \ right) \ cos ^ 3 \ left (\ theta_1 \ right) -20 \ sin ^ 2 \ left (\ theta _2 \ right) \ cos ^ 3 \ left (\ theta _2 \ right) +20 \ sin ^ 2 \ left (\ theta _3 \ right) \ cos ^ 3 \ left (\ theta _3 \ right) -10 \ sin ^ 4 \ left (\ theta _1 \ right) \ cos \ left (\ theta _1 \ right) +10 \ sin ^ 4 \ left (\ theta_2 \ right) \ cos \ left (\ theta _2 \ right) -10 \ sin ^ 4 \ left (\ theta _3 \ right) \ cos \ left (\ theta _3 \ right) + 1 = 0 [/ math]
[matemáticas] -2 \ cos ^ 7 \ left (\ theta _1 \ right) +2 \ cos ^ 7 \ left (\ theta _2 \ right) -2 \ cos ^ 7 \ left (\ theta _3 \ right) +42 \ sin ^ 2 \ left (\ theta _1 \ right) \ cos ^ 5 \ left (\ theta_1 \ right) -42 \ sin ^ 2 \ left (\ theta _2 \ right) \ cos ^ 5 \ left (\ theta _2 \ right) +42 \ sin ^ 2 \ left (\ theta _3 \ right) \ cos ^ 5 \ left (\ theta _3 \ right) -70 \ sin ^ 4 \ left (\ theta _1 \ right) \ cos ^ 3 \ left (\ theta _1 \ right) +70 \ sin ^ 4 \ left (\ theta_2 \ right) \ cos ^ 3 \ left (\ theta _2 \ right) -70 \ sin ^ 4 \ left (\ theta _3 \ right ) \ cos ^ 3 \ left (\ theta _3 \ right) +14 \ sin ^ 6 \ left (\ theta _1 \ right) \ cos \ left (\ theta_1 \ right) -14 \ sin ^ 6 \ left (\ theta _2 \ right) \ cos \ left (\ theta _2 \ right) +14 \ sin ^ 6 \ left (\ theta _3 \ right) \ cos \ left (\ theta _3 \ right) + 1 = 0 [/ math]
Ahora observe que, dado que 5 y 7 son impares, [math] \ sin (\ theta_i) [/ math] siempre aparece con un exponente par. Esto significa que podemos sustituir [math] x_i = \ cos (\ theta_i) [/ math] y obtener un polinomio en [math] x_i [/ math] ‘s, porque entonces [math] \ sin ^ 2 (\ theta_i) = 1 – x_i ^ 2 [/ matemáticas]. Cuando haces esto, obtienes:
[matemáticas] -2 x_1 + 2 x_2-2 x_3 + 1 = \ dfrac {-4} {5} [/ matemáticas]
[matemáticas] -32 x_1 ^ 5 + 40 x_1 ^ 3-10 x_1 + 32 x_2 ^ 5-32 x_3 ^ 5-40 x_2 ^ 3 + 40 x_3 ^ 3 + 10 x_2-10 x_3 + 1 = 0 [/ matemáticas]
[matemáticas] -128 x_1 ^ 7 + 224 x_1 ^ 5-112 x_1 ^ 3 + 14 x_1 + 128 x_2 ^ 7-128 x_3 ^ 7-224 x_2 ^ 5 + 224 x_3 ^ 5 + 112 x_2 ^ 3-112 x_3 ^ 3-14 x_2 + 14 x_3 + 1 = 0 [/ matemáticas]
Ahora, resolver un sistema de ecuaciones polinómicas no es exactamente sencillo, pero se puede hacer a través de métodos básicos de Groebner. Calculé una base de Groebner aquí, y da una descripción concisa de las soluciones: [matemáticas] x_3 [/ matemáticas] es una de las raíces del polinomio irreducible de noveno grado
en general. ]
Habiendo elegido [math] x_3 [/ math], [math] x_2 [/ math] es entonces una raíz de la cuadrática
[Matemáticas] 157415574206075150 x_2 ^ 2 + \ left (141674016785467635-157415574206075150 x_3 \ right) x_2 + (-9122024017968750000 x_3 ^ 8 + 16419643232343750000 x_3 ^ 7 + 8916240158083125000 x_3 ^ 6-20229093108665437500 x_3 ^ 5-2040509184748059375 x_3 ^ 4 + 6629981099774130000 x_3 ^ 3-73395994436560600 x_3 ^ 2-457457311406778480 x_3 + 11558978705814424) [/ matemáticas]
y finalmente
[matemáticas] x_1 = x_2 – x_3 + \ dfrac {9} {10} [/ matemáticas]
Las nueve raíces del polinomio que dan x_3 son reales y se encuentran entre -1 y 1. Además, al conectar cada una de ellas en la definición cuadrática x_2 se obtienen soluciones reales entre -1 y 1.
Por supuesto, no hay una manera visiblemente más simple de describir nuestras soluciones exactamente que “el coseno inverso de un cierto número algebraico de noveno grado”. Sin embargo, si aproximamos estas soluciones, obtenemos
[matemáticas] \ begin {array} {ccc} \ theta _1 \ approx 0.759441 & \ theta _2 \ approx 2.88862 & \ theta _3 \ approx 2.48716 \\ \ theta _1 \ approx 0.252972 y \ theta _2 \ approx 2.38215 y \ theta _3 \ aprox. 2.48716 \\ \ theta _1 \ aprox 1.4003 & \ theta _2 \ aprox 2.9857 & \ theta _3 \ aprox 1.83127 \\ \ theta _1 \ aprox 0.155894 & \ theta _2 \ aprox 1.74129 & \ theta _3 \ aprox 1.83127 \\\ theta _1 \ aprox 1.25446 & \ theta _2 \ aprox 2.25617 & \ theta _3 \ aprox 1.61486 \\ \ theta _1 \ aprox 0.885427 & \ theta _2 \ aprox 1.88714 & \ theta _3 \ aprox 1.61486 \\ \ theta _1 \ aprox 1.83127 & \ theta _2 \ approx 2.9857 & \ theta _3 \ approx 1.4003 \\ \ theta _1 \ approx 0.155894 & \ theta _2 \ approx 1.31032 & \ theta _3 \ approx 1.4003 \\ \ theta _1 \ approx 1.61486 y \ theta _2 \ approx 2.25617 & \ theta _3 \ aprox 1.25446 \\ \ theta _1 \ aprox 0.885427 & \ theta _2 \ aprox 1.52673 y \ theta _3 \ aprox 1.25446 \\ \ theta _1 \ aprox 1.61486 & \ theta _2 \ aprox 1.88714 & \ theta _3 \ aprox 0.885427 \\ \ theta _1 \ aprox 1.25446 y \ theta _2 \ aprox 1.52673 y \ theta _3 \ aprox 0.88 5427 \\ \ theta _1 \ approx 2.48716 & \ theta _2 \ approx 2.88862 & \ theta _3 \ aprox 0.759441 \\ \ theta _1 \ aprox 0.252972 & \ theta _2 \ aprox 0.654432 & \ theta _3 \ aprox 0.759441 \\ \ theta _1 \ approx 2.48716 & \ theta _2 \ approx 2.38215 & \ theta _3 \ approx 0.252972 \\ \ theta _1 \ approx 0.759441 & \ theta _2 \ approx 0.654432 & \ theta _3 \ approx 0.252972 \\ \ theta _1 \ approx 1.83127 y \ theta _2 \ aprox 1.74129 & \ theta _3 \ aprox 0.155894 \\ \ theta _1 \ aprox 1.4003 & \ theta _2 \ aprox 1.31032 & \ theta _3 \ aprox 0.155894 \\ \ end {array} [/ math]
Hice estos cálculos en Mathematica usando TrigExpand
, TrigExpand
, GroebnerBasis
, CountRoots
, las instalaciones básicas de sustitución y programación funcional del lenguaje, TeXForm
para salida y NSolve
para las soluciones aproximadas. En una situación en la que el sistema de ecuaciones era más complicado, habría utilizado Mathematica para hacer los cálculos iniciales y luego poner el cálculo de la base de Groebner en singular
.