- 注册时间
- 2008-11-26
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 149507
- 在线时间
- 小时
|
发表于 2020-4-24 13:29:50
|
显示全部楼层
本帖最后由 mathematica 于 2020-4-24 13:31 编辑
- Clear["Global`*"];
- (*拉格朗日乘子法*)
- f=1/(1+x^2)+1/(1+y^2)+a*(x+y-2)
- (*求导数*)
- fx=D[f,x]//FullSimplify
- fy=D[f,y]//FullSimplify
- fa=D[f,a]//FullSimplify
- (*求解零点*)
- ans=Solve[{fx,fy,fa}=={0,0,0},{x,y,a}]//FullSimplify//ToRadicals;
- Grid[ans]
- N[%,10]
- (*带入得到函数值*)
- FullSimplify[f/.ans]
- N[%,10]
复制代码
求解零点
\[\begin{array}{ccc}
x\to 1 & y\to 1 & a\to \frac{1}{2} \\
x\to 1-\sqrt{2 \left(\sqrt{2}-1\right)} & y\to \sqrt{2 \left(\sqrt{2}-1\right)}+1 & a\to \frac{1}{4 \sqrt{2}} \\
x\to \sqrt{2 \left(\sqrt{2}-1\right)}+1 & y\to 1-\sqrt{2 \left(\sqrt{2}-1\right)} & a\to \frac{1}{4 \sqrt{2}} \\
x\to 1-i \sqrt{2 \left(\sqrt{2}+1\right)} & y\to 1+i \sqrt{2 \left(\sqrt{2}+1\right)} & a\to -\frac{1}{4 \sqrt{2}} \\
x\to 1+i \sqrt{2 \left(\sqrt{2}+1\right)} & y\to 1-i \sqrt{2 \left(\sqrt{2}+1\right)} & a\to -\frac{1}{4 \sqrt{2}} \\
\end{array}\]
数值化
\[
\begin{array}{ccc}
x\to 1.000000000 & y\to 1.000000000 & a\to 0.5000000000 \\
x\to 0.08982027888 & y\to 1.910179721 & a\to 0.1767766953 \\
x\to 1.910179721 & y\to 0.08982027888 & a\to 0.1767766953 \\
x\to 1.000000000-2.197368227 i & y\to 1.000000000+2.197368227 i & a\to -0.1767766953 \\
x\to 1.000000000+2.197368227 i & y\to 1.000000000-2.197368227 i & a\to -0.1767766953 \\
\end{array}
\]
带入函数得到函数值
\[\left\{1,\frac{1}{2}+\frac{1}{\sqrt{2}},\frac{1}{2}+\frac{1}{\sqrt{2}},\frac{1}{2}-\frac{1}{\sqrt{2}},\frac{1}{2}-\frac{1}{\sqrt{2}}\right\}\]
数值化
\[\{1.000000000,1.207106781,1.207106781,-0.2071067812,-0.2071067812\}\]
为什么我感觉与这儿的非常的像?
已知a+b=1,求ab(a^4+b^4)的最大值
https://bbs.emath.ac.cn/forum.ph ... 174&fromuid=865
至少求出来的拉格朗日零点很像,
两个值相等,然后两个对称的实数值,两个对称的复数值
|
|