应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 【求助】求助三元函数条件极值(已解决)
81/1返回列表
查看: 2353  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 式中log指自然对数，^2是2次方，以此类推^3,^4....,

求出y极小值点及对应的x(1),x(2),x(3),求出的y值应该是负数。
通过电子邮件发给我：zsgandyliu@126.com

[ Last edited by woshilsh on 2009-6-13 at 22:24 ]
回复此楼

» 猜你喜欢
1楼 2009-05-27 18:21:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)
★
mingdong(金币+1,VIP+0):谢谢！ 5-29 06:37
引用回帖:
Originally posted by zsgandy at 2009-5-27 18:21:
要求有详细执行程序及结果（推荐用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^ ...

你可以找个遗传算法求函数极值的程序，把里面的函数改一下试试看。
赞一下(10人) 回复此楼
2楼2009-05-27 21:01:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zsgandy

至尊木虫 (职业作家)
请应助者直接给出结果，不需要建议，本人小白！！呵呵
赞一下 回复此楼
3楼2009-05-27 21:28:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

q68

铁杆木虫 (著名写手)
★ ★ ★ ★ ★
zsgandy(金币+5,VIP+0):谢谢 5-28 11:48
我来试一试！
QQ531550249联系一下你的表达式有问题！

[ Last edited by q68 on 2009-5-28 at 10:18 ]
赞一下(20人) 回复此楼
4楼2009-05-28 10:01:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zsgandy

至尊木虫 (职业作家)
引用回帖:
Originally posted by q68 at 2009-5-28 10:01:
我来试一试！
QQ531550249联系一下你的表达式有问题！

[ Last edited by q68 on 2009-5-28 at 10:18 ]

已改正，少了个加号，谢谢！！
赞一下 回复此楼
5楼2009-05-28 11:45:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

q68

铁杆木虫 (著名写手)
★ ★ ★ ★ ★
zsgandy(金币+5,VIP+0):虽然没错误了，但是很感动您的热心！ 5-28 13:29
函数仍然有问题！
回复此楼
6楼2009-05-28 13:15:55
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sldsunny

木虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
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
略加调整自变量，使之在迭代范围内收敛。
赞一下(18人) 回复此楼
7楼2009-05-28 21:00:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zsgandy

至尊木虫 (职业作家)
求助已结束，请版主删帖，谢谢！
赞一下 回复此楼
8楼2009-05-29 10:02:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zsgandy 的主题更新
81/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗？ 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭