| 查看: 2259 | 回复: 7 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】求助三元函数条件极值(已解决)
|
|||
|
要求有详细执行程序及结果（推荐用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 ] |
回复此楼» 猜你喜欢
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
物理学I论文润色/翻译怎么收费?
已经有54人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有0人回复
hakuna
木虫 (知名作家)
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
2楼2009-05-27 21:01:30
3楼2009-05-27 21:28:27
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 15090.2
- 散金: 10
- 红花: 3
- 帖子: 2367
- 在线: 149.8小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
4楼2009-05-28 10:01:54
5楼2009-05-28 11:45:49
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 15090.2
- 散金: 10
- 红花: 3
- 帖子: 2367
- 在线: 149.8小时
- 虫号: 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
8楼2009-05-29 10:02:06
20