在Mathematica中求解36个方程和24个变量的系统

时间:2015-01-23 20:04:53

标签: wolfram-mathematica mathematica-8 mathematica-frontend

我想解决下面的方程组,但是当我使用Solve或NRoots时,我的机器似乎没有响应。有什么想法可以提供帮助吗?你可以在你的机器上运行代码并给我结果。

Solve[{a122 * a223 + a123*a323 + a124*a423 == 
   a222*a123 + a322*a133 + a422*a134, 
  a122*a224 + a123*a324 + a124*a424 == 
   a124*a222 + a322*a134 + a422*a144, 
  a123*a224 + a133*a324 + a134*a424 == 
   a223*a124 + a323*a134 + a423*a144, 
  a122*a233 + a123*a333 + a124*a433 == 
   a223*a123 + a323*a133 + a423*a134, 
  a122*a234 + a123*a334 + a124*a434 == 
   a223*a124 + a323*a134 + a423*a144, 
  a123*a234 + a133*a334 + a134*a434 == 
   a233*a124 + a333*a134 + a433*a144, 
  a122*a234 + a123*a334 + a124*a434 == 
   a224*a123 + a324*a133 + a424*a134, 
  a122*a244 + a123*a344 + a124*a444 == 
   a224*a124 + a324*a134 + a424*a144, 
  a123*a244 + a133*a344 + a134*a444 == 
   a234*a124 + a334*a134 + a434*a144, 
  a123 + a223*a323 + a224*a423 == a322*a233 + a422*a234, 
  a124 + a223*a324 + a224*a424 == a322*a234 + a422*a224, 
  a233*a324 + a234*a424 == a323*a234 + a423*a244, 
  a133 + a222*a233 + a223*a333 + a224*a433 == 
   a223*a223 + a323*a233 + a423*a234, 
  a134 + a222*a234 + a223*a334 + a224*a434 == 
   a223*a224 + a323*a234 + a423*a244, 
  a223*a234 + a233*a334 + a234*a434 == 
   a233*a224 + a333*a234 + a433*a244, 
  a134 + a222*a234 + a223*a334 + a224*a434 == 
   a224*a223 + a324*a233 + a424*a234, 
  a144 + a222*a244 + a223*a344 + a224*a444 == 
   a224*a224 + a324*a234 + a424*a244, 
  a223*a244 + a233*a344 + a234*a444 == 
   a234*a224 + a334*a234 + a434*a244, 
  a322*a223 + a323*a323 + a324*a423 == 
   a122 + a222*a323 + a322*a333 + a422*a334, 
  a322*a224 + a323*a324 + a324*a424 == 
   a222*a324 + a322*a334 + a422*a344, 
  a124 + a323*a224 + a333*a324 + a334*a424 == 
   a223*a324 + a323*a334 + a423*a334, 
  a123 + a223*a323 + a423*a334 == a322*a233 + a324*a433, 
  a322*a234 + a324*a434 == a223*a324 + a423*a344, 
  a134 + a323*a234 + a334*a434 == a233*a324 + a433*a344, 
  a322*a234 + a323*a334 + a324*a434 == 
   a124 + a224*a323 + a324*a333 + a424*a334, 
  a322*a244 + a323*a344 + a324*a444 == 
   a224*a324 + a324*a334 + a424 * a344, 
  a144 + a323*a244 + a333*a344 + a334*a444 == 
   a234*a324 + a334*a334 + a434*a344, 
  a422*a234 + a423*a334 == a224*a423 + a324*a433, 
  a422*a244 + a423*a344 == a224*a424 + a324*a434, 
  a423*a244 + a433*a344 == a134 + a234*a424 + a334*a434, 
  a422*a223 + a423*a323 + a424*a423 == 
   a122 + a222*a423 + a322*a433 + a422*a434, 
  a422*a224 + a423*a324 + a424*a424 == 
   a122 + a222*a424 + a323*a434 + a422*a444, 
  a423*a224 + a433*a324 + a434*a424 == 
   a123 + a223*a424 + a323*a434 + a423*a444, 
  a422*a233 + a423*a333 + a424*a433 == 
   a222*a423 + a323*a433 + a423*a434, 
  a422*a234 + a423*a334 + a424*a434 == 
   a123 + a222*a424 + a323*a434 + a423*a444, 
  a423*a234 + a433*a334 + a434*a434 == 
   a133 + a233*a424 + a333*a433 + a433*a444}, {a122, a123, a124, a133,
   a134, a144, a222, a223, a224, a233, a234, a244, a322, a323, a324, 
  a333, a334, a344, a422, a423, a424, a433, a434, a444}]

1 个答案:

答案 0 :(得分:0)

校正。最后我向我解释了我的错误。谢谢你坚持不懈。有时候我很难认清事物。所以现在我正确地定位了Norm。 (是的还有其他方法可以做到这一点)

NMinimize[
 Norm[a122*a223+a123*a323+a124*a423-(a222*a123+a322*a133+a422*a134)] +
 Norm[a122*a224+a123*a324+a124*a424-(a124*a222+a322*a134+a422*a144)] +
  ... +
 Norm[a423*a234+a433*a334+a434*a434-(a133+a233*a424+a333*a433+a433*a444)]],
 {a122, a123, a124, a133, a134, a144, a222, a223, a224, a233, a234, a244,
  a322, a323, a324, a333, a334, a344, a422, a423, a424, a433, a434, a444}]

返回

{1.11167*10^-6, {a122 -> -7.49406*10^-8, a123 -> -4.98174*10^-7, a124 -> 2.50543*10^-7,
a133 -> 2.0337*10^-7, a134 -> -1.16196*10^-6, a144 -> 0.0219087, a222 -> 0.174636,
a223 -> 0.0743701, a224 -> 0.138138, a233 -> 0.0398568, a234 -> -0.0165303,
a244 -> 0.0390774, a322 -> 0.150348, a323 -> 0.0805803, a324 -> -0.0334195,
a333 -> 0.0239558, a334 -> 0.159051, a344 -> -0.182692, a422 -> 4.5244*10^-6,
a423 -> 1.40132*10^-7, a424 -> 9.09125*10^-7, a433 -> 2.59763*10^-6,
a434 -> 4.41731*10^-6, a444 -> 0.0324925}}

在几分之一秒内。

系统似乎至少有一个非平凡的(并非全是零)解决方案。