- 注册时间
- 2021-11-19
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 9011
- 在线时间
- 小时
|
楼主 |
发表于 2024-5-6 09:18:17
|
显示全部楼层
- Clear["Global`*"];
- (*子函数,利用三边计算角的余弦值,角是c边所对的角*)
- cs[a_,b_,c_]:=((a^2+b^2-c^2)/(2*a*b))
- (*子函数,四面体体积公式,a,b,c分别是从一个顶点出发的三条棱,x,y,z分别是对棱*)
- fun[a_,b_,c_,x_,y_,z_]:=Sqrt[Det[{{0,1,1,1,1},{1,0,a^2,b^2,c^2},{1,a^2,0,z^2,y^2},{1,b^2,z^2,0,x^2},{1,c^2,y^2,x^2,0}}]/288]
- (*余弦定理,约束条件1*)
- cond1=cs[PA,6,PB]-Cos[120*Degree]//Together//Numerator
- (*余弦定理,约束条件2*)
- cond2=fun[PA,PB,PC,6,6,6]^2*4
- (*目标函数*)
- f=PB/PC+x1*cond1+x2*cond2
- (*求偏导数*)
- fx=D[f,{{PA,PB,PC,x1,x2}}]//Simplify
- (*解方程组*)
- ans=Solve[{fx==0,PA>0,PB>0,PC>0},{PA,PB,PC,x1,x2}]
- Grid[ans,Alignment->Left](*列表显示*)
- aaa=(f/.ans)//Simplify (*得到目标函数值*)
复制代码
方程组
\[\left\{2 \left(\text{PA} \text{x2} \left(-2 \text{PA}^2+\text{PB}^2+\text{PC}^2+36\right)+(\text{PA}+3) \text{x1}\right),2 \text{PB} \text{x2} \left(\text{PA}^2+\text{PC}^2+36\right)-4 \text{PB}^3 \text{x2}-2 \text{PB} \text{x1}+\frac{1}{\text{PC}},2 \text{PC} \text{x2} \left(\text{PA}^2-2 \text{PC}^2+36\right)+2 \text{PB}^2 \text{PC} \text{x2}-\frac{\text{PB}}{\text{PC}^2},\text{PA}^2+6 \text{PA}-\text{PB}^2+36,-\text{PA}^4+\text{PA}^2 \left(\text{PB}^2+\text{PC}^2+36\right)-\text{PB}^4+\text{PB}^2 \left(\text{PC}^2+36\right)-\text{PC}^4+36 \text{PC}^2-1296\right\}\]
方程组求解结果
\[\begin{array}{lllll}
\text{PA}\to 6 & \text{PB}\to 6 \sqrt{3} & \text{PC}\to 6 & \text{x1}\to -\frac{1}{36 \sqrt{3}} & \text{x2}\to \frac{1}{2592 \sqrt{3}} \\
\text{PA}\to 6 & \text{PB}\to 6 \sqrt{3} & \text{PC}\to 12 & \text{x1}\to \frac{1}{144 \sqrt{3}} & \text{x2}\to -\frac{1}{20736 \sqrt{3}} \\
\end{array}\]
目标函数值
\[\left\{\sqrt{3},\frac{\sqrt{3}}{2}\right\}\] |
|