eqn1 = 2 b^2 + b^4 + 4 a^3 (3 + b) + 6 a^2 (-1 + 2 b + b^2 - 4 c) +
    4 a (-b^2 + b^3 - 6 b c + 3 c (1 + c)) ==
   4 b^3 + 6 c^2 + 4 b (-3 + c) c^2;(*用此等式解出c的实数解，其中0<a<c<a+b<1*)
exp1 = ((-2 + 4 a + 4 b) (-a + c)^3 -
     b (a - c) (3 a^3 + 2 b - 4 b^2 + b^3 + a^2 (6 b - c) - c^3 +
        a (-4 b + 4 b^2 - c^2)))/b;
(*然后将c代入上述表达式中得到表达式【1】*)
exp2 = -(3 + 3 a^4 + 6 a^3 (-2 + b) - 8 b + 6 b^2 - b^3 + 2 a^2 (9 - 11 b + 2 b^2) +
     a (-12 + 24 b - 12 b^2 + b^3));
exp3 = -a (3 a^3 + 6 a^2 b + 4 a (-1 + b) b + b (2 - 4 b + b^2));
(*表达式【1】和上面两个表达式相等，得到两个方程，求a,b的近似值(小数点后保留4位)使得表达式【1】值<0*)

cond = {exp1 < 0, 0 < a < c < a + b < 1}

NMinimize[{( exp1 - exp3)^2 + (Subtract @@ eqn1)^2 + (exp1 - exp2)^2, cond}, {a, b, c},
  Method -> "DifferentialEvolution"] // AbsoluteTiming

(* {1.871953, {3.71924*10^-17, {a -> 0.223149, b -> 0.553702, c -> 0.5}}} *)

2楼: Originally posted by zho123855 at 2016-01-10 19:54:17
这个应该没太大难度吧，把代码贴出来看看

4楼: Originally posted by zho123855 at 2016-01-10 20:49:48
懂你的意思了。
你需要首先求解一个符号方程得到C的解析表达是吧
这个就有难度了，一方面计算量确实很大，另一方面，C一般不会只有一个解，你这个思路可能不大行得通，换个思路？

7楼: Originally posted by zho123855 at 2016-01-11 11:58:41
早上过来用MATLAB尝试了一下，五分钟搞定，解出来a=0.2232,b=0.5537,c=0.5
不过这个方程组确实是多解的，需要一些小技巧来排除不需要的解。

7楼: Originally posted by zho123855 at 2016-01-11 11:58:41
早上过来用MATLAB尝试了一下，五分钟搞定，解出来a=0.2232,b=0.5537,c=0.5
不过这个方程组确实是多解的，需要一些小技巧来排除不需要的解。

9楼: Originally posted by zho123855 at 2016-01-11 12:13:40
对我而言，绝大多数时候会优先考虑MATLAB，Mathematica确实是优秀的数学软件，尤其是在符号推导方面，但做数值计算，尤其是大规模数值计算，目前看不出相对于MATLAB的优势。
具体到你这个问题，既然你知道C的取值 ...