11楼: Originally posted by xzczd at 2016-01-11 20:04:50
答题前姑且问一句，你的方程只有一组解是运气好还是必如此？如果必然只有一组解，那这个问题会简单很多。

11楼: Originally posted by xzczd at 2016-01-11 20:04:50
答题前姑且问一句，你的方程只有一组解是运气好还是必如此？如果必然只有一组解，那这个问题会简单很多。

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}}} *)

14楼: Originally posted by xzczd at 2016-01-11 20:55:23
必然只有一组解的话，那么有一个很简单的方法：

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 ...

15楼: Originally posted by 青沙山 at 2016-01-11 21:13:59
太棒了！  万分感谢啊！

代码还在学习模仿中，请问”Subtract @@ eqn1“这个表示什么意思呢？还有”差分进化算法“——好高大上啊！...

16楼: Originally posted by xzczd at 2016-01-11 21:29:27
你直接在帮助里搜 Subtract 和 @@ 就明白了。...

17楼: Originally posted by 青沙山 at 2016-01-11 22:11:03
懂了！谢谢！

刚才在nb文件中可能输入了”?@“，楞是没看明白！...

18楼: Originally posted by xzczd at 2016-01-11 22:17:24
用？查函数效率太低，要养成按把光标停在函数上然后按F1的习惯。...