- 注册时间
- 2021-11-19
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 8641
- 在线时间
- 小时
|
楼主 |
发表于 2024-4-15 11:08:34
|
显示全部楼层
- Clear["Global`*"];
- (*计算余弦值子函数,利用三边计算余弦值*)
- cs[a_,b_,c_]:=((a^2+b^2-c^2)/(2*a*b))
- aaa=Solve[cs[3,25/8,BC]==1/2&&BC>0,{BC}] (*用余弦定理求解BC的长度*)
- (*子函数,四面体体积公式,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]
- {AB,BC,AC}={3,(BC/.aaa[[1]]),25/8};(*线段长度赋值*)
- f=7*AP+5*BP+8*CP+t*(fun[AP,BP,CP,BC,AC,AB]^2-0)//Simplify (*定义目录函数,约束条件:四面体的体积等于零*)
- fx=D[f,{{AP,BP,CP,t}}]//Simplify (*对四个变量求解偏导数*)
- ans=Solve[{fx==0},{AP,BP,CP,t}] (*偏导数等于零,求解方程组*)
- out=f/.ans[[2]] (*得到目标函数的极值*)
复制代码
我只会拉格朗日乘子法,
求解结果:
BC长度
\[\left\{\left\{\text{BC}\to \frac{\sqrt{601}}{8}\right\}\right\}\]
目标函数是
\[-\frac{t \left(38464 \text{AP}^4-8 \text{AP}^2 \left(5200 \text{BP}^2+4416 \text{CP}^2+45075\right)+40000 \text{BP}^4-600 \text{BP}^2 \left(64 \text{CP}^2+575\right)+9 \left(4096 \text{CP}^4-41600 \text{CP}^2+375625\right)\right)}{589824}+7 \text{AP}+5 \text{BP}+8 \text{CP}\]
四个偏导数是
\[\left\{-\frac{601 \text{AP}^3 t}{2304}+\frac{\text{AP} t \left(5200 \text{BP}^2+4416 \text{CP}^2+45075\right)}{36864}+7,\frac{25 \text{BP} t \left(208 \text{AP}^2+192 \text{CP}^2+1725\right)}{36864}-\frac{625 \text{BP}^3 t}{2304}+5,\frac{1}{768} \text{CP} t \left(92 \text{AP}^2+100 \text{BP}^2+975\right)-\frac{\text{CP}^3 t}{4}+8,\frac{-38464 \text{AP}^4+8 \text{AP}^2 \left(5200 \text{BP}^2+4416 \text{CP}^2+45075\right)-40000 \text{BP}^4+600 \text{BP}^2 \left(64 \text{CP}^2+575\right)-9 \left(4096 \text{CP}^4-41600 \text{CP}^2+375625\right)}{589824}\right\}\]
求解结果是
\[\left\{\left\{\text{AP}\to -\frac{15}{8},\text{BP}\to -\frac{117}{56},\text{CP}\to -\frac{10}{7},t\to \frac{702464}{219375}\right\},\left\{\text{AP}\to \frac{15}{8},\text{BP}\to \frac{117}{56},\text{CP}\to \frac{10}{7},t\to -\frac{702464}{219375}\right\}\right\}\]
最小值是35
|
|