- 注册时间
- 2008-11-26
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 149507
- 在线时间
- 小时
|
楼主 |
发表于 2021-4-9 08:57:25
|
显示全部楼层
- Clear["Global`*"];(*mathematica11.2,win7(64bit)Clear all variables*)
- (*海伦公式*)
- heron[a_,b_,c_]:=Module[{p=(a+b+c)/2},Sqrt[p*(p-a)*(p-b)*(p-c)]]
- (*利用第一个正弦定理与余弦定理,简化第一个约束条件*)
- aaa=TrigExpand[Sin[A-B]-(3-4*Cos[A])*Sin[B]]/.{Sin[A]->a,Sin[B]->b,Cos[A]->(b^2+c^2-a^2)/(2*b*c),Cos[B]->(a^2+c^2-b^2)/(2*a*c)};
- bbb=Together[aaa]//Numerator
- (*得到目标函数*)
- ff=c-x1*(bbb)-x2*(heron[a,b,c]^2-5);
- (*求解所有可能的极值点,实数范围内*)
- ccc=Solve[D[ff,{{a,b,c,x1,x2}}]==0,{a,b,c,x1,x2},Reals]//FullSimplify//ToRadicals;
- Grid[ccc,Alignment->Left]
- ddd=(ff/.ccc)//FullSimplify;
- Grid[{ddd},Alignment->Left]
- Maximize[{c,bbb==0&&heron[a,b,c]==Sqrt[5]&&a>=0&&b>=0&&c>=0},{a,b,c}]
- Minimize[{c,bbb==0&&heron[a,b,c]==Sqrt[5]&&a>=0&&b>=0&&c>=0},{a,b,c}]//FullSimplify
复制代码
我只会拉格朗日乘子法求解
求解结果
\[-a^2+b^2-3 b c+2 c^2\]
\[\begin{array}{lllll}
a\to -2 \sqrt{\frac{7}{5}} & b\to -3 \sqrt{\frac{2}{5}} & c\to -\sqrt{10} & \text{x1}\to -\frac{1}{2 \sqrt{10}} & \text{x2}\to -\frac{1}{2 \sqrt{10}} \\
a\to -2 \sqrt{\frac{7}{5}} & b\to 3 \sqrt{\frac{2}{5}} & c\to \sqrt{10} & \text{x1}\to \frac{1}{2 \sqrt{10}} & \text{x2}\to \frac{1}{2 \sqrt{10}} \\
a\to 2 \sqrt{\frac{7}{5}} & b\to -3 \sqrt{\frac{2}{5}} & c\to -\sqrt{10} & \text{x1}\to -\frac{1}{2 \sqrt{10}} & \text{x2}\to -\frac{1}{2 \sqrt{10}} \\
a\to 2 \sqrt{\frac{7}{5}} & b\to 3 \sqrt{\frac{2}{5}} & c\to \sqrt{10} & \text{x1}\to \frac{1}{2 \sqrt{10}} & \text{x2}\to \frac{1}{2 \sqrt{10}} \\
\end{array}\]
\[\begin{array}{llll}
-\sqrt{10} & \sqrt{10} & -\sqrt{10} & \sqrt{10} \\
\end{array}\]
\[\left\{\sqrt{10},\left\{a\to 2 \sqrt{\frac{7}{5}},b\to 3 \sqrt{\frac{2}{5}},c\to \sqrt{10}\right\}\right\}\]
|
|