| 查看: 2372 | 回复: 7 | |||
| 当前主题已经存档。 | |||
zsgandy至尊木虫 (职业作家)
|
[交流]
【求助】求助三元函数条件极值(已解决)
|
||
|
要求有详细执行程序及结果（推荐用MATLAB），题目如下： y=x(1)*(-271.21+211.206579*x(3)-38.5844296*x(3)*log(x(3))+0.018531982*x(3)^2-0.000005764227*x(3)^3+74092*x(3)^(-1)+7934*10^(-20)*x(3)^7)+(1-x(1)-x(2))*(-5677.958+146.176046*x(3)-32.4913959*x(3)*log(x(3))+0.00154613*x(3)^2)+x(2)*(9496.31-9.809114*x(3)-8.2590486*x(3)*log(x(3))-0.016814429*x(3)^2+2.623132*10^(-6)*x(3)^3-1081244*x(3)^(-1))+8.314*x(3)*(x(1)*log(x(1))+x(2)*log(x(2))+(1-x(1)-x(2))*log(1-x(1)-x(2))) +x(1)*(1-x(1)-x(2))*(47933.6-10.71995*x(3)+(14407.33-6.65287*x(3))*(2*x(1)+x(2)-1)+(4742.6-0.72034*x(3))*(2*x(1)+x(2)-1)^2)+x(1)*x(2)*(16330-4.983*x(3)+(4112-1.151*x(3))*(x(1)-x(2))+(1765-0.574*x(3))*(x(1)-x(2))^2) +x(2)*(1-x(1)-x(2))*(6200-0.418*x(3)+(1-x(1)-2*x(2))*(790-1.914*x(3))) +x(1)*x(2)*(1-x(1)-x(2))*((-100000+57.7*x(3))*x(1)-3286*(1-x(1)-x(2))-364*x(2)) 其中， 0 求出y极小值点及对应的x(1),x(2),x(3),求出的y值应该是负数。 通过电子邮件发给我：zsgandyliu@126.com [ Last edited by woshilsh on 2009-6-13 at 22:24 ] |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有75人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
hakuna
木虫 (知名作家)
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
2楼2009-05-27 21:01:30
zsgandy
至尊木虫 (职业作家)
- 应助: 1 (幼儿园)
- 金币: 19782.8
- 散金: 70
- 红花: 5
- 帖子: 3682
- 在线: 331.3小时
- 虫号: 427932
- 注册: 2007-07-30
- 性别: GG
- 专业: 金属材料的制备科学与跨学
3楼2009-05-27 21:28:27
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 16098.2
- 散金: 10
- 红花: 3
- 帖子: 2529
- 在线: 149.9小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
4楼2009-05-28 10:01:54
zsgandy
至尊木虫 (职业作家)
- 应助: 1 (幼儿园)
- 金币: 19782.8
- 散金: 70
- 红花: 5
- 帖子: 3682
- 在线: 331.3小时
- 虫号: 427932
- 注册: 2007-07-30
- 性别: GG
- 专业: 金属材料的制备科学与跨学
5楼2009-05-28 11:45:49
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 16098.2
- 散金: 10
- 红花: 3
- 帖子: 2529
- 在线: 149.9小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
6楼2009-05-28 13:15:55
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
zsgandy(金币+90,VIP+0):太谢谢了！！！！！！！ 5-28 23:32
zsgandy(金币+90,VIP+0):太谢谢了！！！！！！！ 5-28 23:32
|
-38.5844296*x(3)*log(x(3))+0.018531982*x(3)^2-0.000005764227*x(3)^3+74092*x(3)^(-1)+7934*10^(-20)*x(3)^7)+(1-x(1)-x(2))*(-5677.958+146.176046*x(3)-32.4913959*x(3)*log(x(3))+0.00154613*x(3)^2)+x(2)*(9496.31-9.809114*x(3)-8.2590486*x(3)*log(x(3))-0.016814429*x(3)^2+2.623132*10^(-6)*x(3)^3-1081244*x(3)^(-1))+8.314*x(3)*(x(1)*log(x(1))+x(2)*log(x(2))+(1-x(1)-x(2))*log(1-x(1)-x(2)))+x(1)*(1-x(1)-x(2))*(47933.6-10.71995*x(3)+(14407.33-6.65287*x(3))*(2*x(1)+x(2)-1)+(4742.6-0.72034*x(3))*(2*x(1)+x(2)-1)^2)+x(1)*x(2)*(16330-4.983*x(3)+(4112-1.151*x(3))*(x(1)-x(2))+(1765-0.574*x(3))*(x(1)-x(2))^2)+x(2)*(1-x(1)-x(2))*(6200-0.418*x(3)+(1-x(1)-2*x(2))*(790-1.914*x(3)))+x(1)*x(2)*(1-x(1)-x(2))*((-100000+57.7*x(3))*x(1)-3286*(1-x(1)-x(2))-364*x(2)) ','x') LB=[0.004 0.7355 456.15]; UB=[1 1 753.15]; x0=[0.003 0 400] [x,fval,exitflag]=fmincon(Y,x0,[],[],[],[],LB,UB) Y = Inline function: Y(x) = x(1)*(-271.21+211.206579*x(3)-38.5844296*x(3)*log(x(3))+0.018531982*x(3)^2-0.000005764227*x(3)^3+74092*x(3)^(-1)+7934*10^(-20)*x(3)^7)+(1-x(1)-x(2))*(-5677.958+146.176046*x(3)-32.4913959*x(3)*log(x(3))+0.00154613*x(3)^2)+x(2)*(9496.31-9.809114*x(3)-8.2590486*x(3)*log(x(3))-0.016814429*x(3)^2+2.623132*10^(-6)*x(3)^3-1081244*x(3)^(-1))+8.314*x(3)*(x(1)*log(x(1))+x(2)*log(x(2))+(1-x(1)-x(2))*log(1-x(1)-x(2)))+x(1)*(1-x(1)-x(2))*(47933.6-10.71995*x(3)+(14407.33-6.65287*x(3))*(2*x(1)+x(2)-1)+(4742.6-0.72034*x(3))*(2*x(1)+x(2)-1)^2)+x(1)*x(2)*(16330-4.983*x(3)+(4112-1.151*x(3))*(x(1)-x(2))+(1765-0.574*x(3))*(x(1)-x(2))^2)+x(2)*(1-x(1)-x(2))*(6200-0.418*x(3)+(1-x(1)-2*x(2))*(790-1.914*x(3)))+x(1)*x(2)*(1-x(1)-x(2))*((-100000+57.7*x(3))*x(1)-3286*(1-x(1)-x(2))-364*x(2)) x0 = 0.0030 0 400.0000 Warning: Large-scale (trust region) method does not currently solve this type of problem, switching to medium-scale (line search). > In fmincon at 260 Optimization terminated: first-order optimality measure less than options.TolFun and maximum constraint violation is less than options.TolCon. Active inequalities (to within options.TolCon = 1e-006): lower upper ineqlin ineqnonlin 1 3 2 x = 0.0040 0.7355 753.1500 fval = -5.3350e+004 exitflag = 1 略加调整自变量，使之在迭代范围内收敛。 |
7楼2009-05-28 21:00:08
zsgandy
至尊木虫 (职业作家)
- 应助: 1 (幼儿园)
- 金币: 19782.8
- 散金: 70
- 红花: 5
- 帖子: 3682
- 在线: 331.3小时
- 虫号: 427932
- 注册: 2007-07-30
- 性别: GG
- 专业: 金属材料的制备科学与跨学
8楼2009-05-29 10:02:06
