- 注册时间
- 2021-11-19
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 8641
- 在线时间
- 小时
|
发表于 2023-3-20 14:16:12
|
显示全部楼层
本帖最后由 nyy 于 2023-3-20 14:41 编辑
https://bbs.emath.ac.cn/forum.ph ... 0&fromuid=14149
引用这儿的结果,化简问题
\[f=\sqrt{a^2-6 a+5}+\sqrt{b^2-6 b+5}+\sqrt{c^2-6 c+5}+t \left(a^2+b^2+c^2+2 a b c-1\right)\]
求偏导数,解方程组,求解结果如下:
\[
\begin{array}{llll}
a\to -1.000000000 & b\to 3.000000000 & c\to 3.000000000 & t\to 0.07216878365 \\
a\to 0.5000000000 & b\to 0.5000000000 & c\to 0.5000000000 & t\to 1.111111111 \\
a\to 3.000000000 & b\to -1.000000000 & c\to 3.000000000 & t\to 0.07216878365 \\
a\to 3.000000000 & b\to 3.000000000 & c\to -1.000000000 & t\to 0.07216878365 \\
a\to 0.2674923595+0.2608869774 i & b\to 0.9930197051-0.2791410926 i & c\to 0.2674923595+0.2608869774 i & t\to 0.8201469710-0.4808912947 i \\
a\to 0.9930197051-0.2791410926 i & b\to 0.2674923595+0.2608869774 i & c\to 0.2674923595+0.2608869774 i & t\to 0.8201469710-0.4808912947 i \\
a\to 0.2674923595+0.2608869774 i & b\to 0.2674923595+0.2608869774 i & c\to 0.9930197051-0.2791410926 i & t\to 0.8201469710-0.4808912947 i \\
a\to 0.2674923595-0.2608869774 i & b\to 0.9930197051+0.2791410926 i & c\to 0.2674923595-0.2608869774 i & t\to 0.8201469710+0.4808912947 i \\
a\to 0.9930197051+0.2791410926 i & b\to 0.2674923595-0.2608869774 i & c\to 0.2674923595-0.2608869774 i & t\to 0.8201469710+0.4808912947 i \\
a\to 0.2674923595-0.2608869774 i & b\to 0.2674923595-0.2608869774 i & c\to 0.9930197051+0.2791410926 i & t\to 0.8201469710+0.4808912947 i \\
\end{array}
\]
很显然取第二行的解。
代入前四个解,得到函数值
\[
\left\{2 \sqrt{3}+4 i,\frac{9}{2},2 \sqrt{3}+4 i,2 \sqrt{3}+4 i\right\}
\]
全部代码如下:
- Clear["Global`*"];
- cond=Det[{{-1,a,b},{a,-1,c},{b,c,-1}}](*约束条件*)
- (*拉格朗日乘子法建立目标函数*)
- f=Sqrt[a^2-6a+5]+Sqrt[b^2-6b+5]+Sqrt[c^2-6c+5]+t*cond
- (*求偏导数,解方程组*)
- ans=Solve[GroebnerBasis@D[f,{{a,b,c,t}}]==0,{a,b,c,t}]
- aaa=N[ans,10](*结果数值化*)
- Grid[aaa,Alignment->Left](*列表显示*)
- bbb=f/.ans[[1;;4]](*只代入前几个解,因为后面是复数解*)
复制代码
|
|