Question

求解非线性方程组，用MATLAB求解得到无解析解，用mathematica跑了一上午也没个结果。

matlab程序如下：
syms a1 a2 a3 a4 a5 b31 b32 b41 b42 b43 b51 b52 b53 b54 c1 c2 c3 c4 c5 c6 c7 c8 c9 ;
[a1,a2,a3,a4,a5,b31,b32,b41,b42,b43,b51,b52,b53,b54,c1,c2,c3,c4,c5,c6,c7,c8,c9]=solve('a2-b31-b32','a3-b41-b42-b43','c2*a1+c3*a2+c4*a3-1/2','c2*a1^2+c3*a2^2+c4*a3^2-1/3','c2*a1^3+c3*a2^3+c4*a3^3-1/4','c4*b42*a1^2+c4*b43*a2^2-1/12','c3*b32*a1*a2+c4*b42*a1*a3+c4*b43*a2*a3-1/8','c3*b32*a1+c4*b42*a1+c4*b43*a2-1/6','c1+c2+c3+c4-1','b51+b52+b53+b54-a4','c5+c6+c7+c8+c9-1','c6*a1+c7*a2+c8*a3+c9*a4-1/2','c6*a1^2+c7*a2^2+c8*a3^2+c9*a4^2-1/3','c6*a1^3+c7*a2^3+c8*a3^3+c9*a4^3-1/4','c6*a1^4+c7*a2^4+c8*a3^4+c9*a4^4-1/5','c7*b32*a1+c7*b42*a1+c8*b43*a2+c9*b52*a1+c9*b53*a2+c9*b54*a3-1/6','c7*b32*a1^2+c7*b42*a1^2+c8*b43*a2^2+c9*b52*a1^2+c9*b53*a2^2+c9*b54*a3^2-1/12','c7*b32*a1^3+c7*b42*a1^3+c8*b43*a2^3+c9*b52*a1^3+c9*b53*a2^3+c9*b54*a3^3-1/20','c8*b43*b32*a1+c9*b53*b32*a1+c9*b54*b42*a1+c9*b54*b43*a2-1/24','c7*b32*a1*a2^2+c8*b42*a1*a3^2+c8*b43*a2*a3^2+c9*b52*a1*a4^2+c9*b53*a2*a4^2+c9*b54*a3*a4^2-1/10','c7*b32*a1^2*a2+c8*b42*a1^2*a3+c8*b43*a2^2*a3+c9*b52*a1^2*a4+c9*b53*a2^2*a4+c9*b54*a3^2*a4-1/15','c8*b43*b32*(a1*a2+a1*a3)+c9*b53*b32*(a1*a2+a1*a4)+c9*b54*b42*(a1*a3+a1*a4)+c9*b54*b43*(a2*a3+a2*a4)-7/120','c7*b32^2*a1^2+c8*(b42*a1+b43*a2)^2+c9*(b52*a1+b53*a2+b54*a3)^2-1/20')
请各位大神帮忙！谢谢啦~~
方程做了一点修改：又添加了几个方程：

c9*b54*b43*b32*a1==1/120,
c4*b43*b32*a1==1/24,
c8*b43*b32*a1^2+c9*b53*b32*a1^2+c9*b54*b42*a1^2==1/60

Answer 1

不懂matlab语法的mathematica使用者路过。你用matlab求解析解……算了，这个不是重点。但是感觉上解不出来是因为你的方程阶数太高了，话说你电脑还真撑得住啊，内存满大嘛。尽可能缩小未知数的取值范围再用Reduce试下，再不行就试下FindRoot，再再不行就考虑修正你的模型吧……

追问

搞得不太明白，能具体点吗？万分感谢！

追答

……这还要怎么具体啊，我的意思是，比如，你的各个未知量，是否全是实数？是否有较小的取值范围？有的话，把这些不等式也加进去，然后用Reduce解，就有可能解的出来，FindRoot也是思路之一，当然那要求你对你的方程的解的性质有更多的了解：方程只有一组解吗？有多组解的话，它们是在哪些位置附近？要是这样还不行，就要考虑你对原始问题的分析是不是有纰漏，是不是有什么重要的可用条件没有在你的方程里体现出来？方程是否能列的更为简洁？诸如此类。总之，十分依赖你对原始问题的分析就是了。

Answer 2

解方程组用LinGo吧