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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助，计算方程最大值
61/1返回列表
查看: 1145  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zh327112793

木虫 (小有名气)

[求助] 求助，计算方程最大值

求助，计算方程最大值，请给出结果及方法
0.00038596*x3*x3*x3-1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113
其中x1的范围是0.56-0.8，x2的范围是0.001-0.031，x3的范围是600-1100
如发错板块，请版主处理
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 方程式.doc
    • 2013-10-22 14:17:44, 19.5 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-10-22 14:18:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fegg7502: 金币+1, 应助指数+1, 鼓励交流 2013-10-22 15:12:42
zh327112793: 金币+45, ★★★很有帮助 2013-10-22 15:26:00
1stopt软件计算代码和结果：
CODE:
Parameters x1[0.56,0.8],x2[0.001,0.031],x3[600,1100];
MaxFunction    -1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113;

结果：
目标函数值(最大): 41244.1831170177
x1: 0.8
x2: 0.0158529785212449
x3: 600
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2013-10-22 15:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zh327112793: 金币+50, ★★★★★最佳答案 2013-10-22 15:27:33
csgt0: 金币+1, 谢谢 2013-10-23 16:52:55
迭代数: 31
计算用时(时:分:秒:毫秒): 00:00:00:359
计算中止原因: 达到收敛判定标准
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
函数表达式: 0.00038596*x3*x3*x3-1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696
          *x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113
目标函数值(最大): 124611.543116985
x1: 0.8
x2: 0.0158530344445098
x3: 600

====== 计算结束 ======
赞一下 回复此楼
3楼2013-10-22 15:21:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zh327112793

木虫 (小有名气)
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-10-22 15:03:01
1stopt软件计算代码和结果：
Parameters x1,x2,x3;
MaxFunction    -1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3 ...

把x1 x2 x3代入方程后 函数值等于124832.8627不等于41244.183？
赞一下 回复此楼
4楼2013-10-22 15:22:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
zh327112793: 金币+5, ★★★很有帮助 2013-10-22 15:28:35
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-10-22 15:03:01
1stopt软件计算代码和结果：
Parameters x1,x2,x3;
MaxFunction    -1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3 ...

CODE:
function max_Fmincon
format short

clear all
clc

x0 = [0.6; 0.01;800];

[x,fval] = fmincon(@ObjFunc,x0,[],[],[],[],[],[],@NlinCons)
% ------------------------------------------------------------------
ff=-ObjFunc(x)

function f = ObjFunc(x)
x1=x(1);
x2=x(2);
x3=x(3);
f = -(0.00038596*x3*x3*x3-1793444.01*x1*x1*x1-...
    0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+...
    333.10356*x3+877480.696*x2-2158343.84*x1+...
    7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113);
% ------------------------------------------------------------------

function [c,ceq] = NlinCons(x)
c(1) = x(1)-0.8;
c(2)= 0.56-x(1);
c(3) = x(2)-0.031;
c(4)= 0.001-x(2);
c(5) = x(3)-1100;
c(6)= 600-x(3);
ceq=[];

结果：
x =

    0.8000
    0.0159
  600.0000

x都是一样的，不过函数的最大值：1.2461e+005
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2013-10-22 15:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
不好意思，目标函数很长，2楼的1stopt代码漏掉一些，改正后如下：
CODE:
Parameters x1[0.56,0.8],x2[0.001,0.031],x3[600,1100];
MaxFunction    0.00038596*x3*x3*x3-1793444.01*x1*x1*x1-0.9269*x3*x3-10308400*x2*x2+3746080.65*x1*x1+333.10356*x3+877480.696*x2-2158343.84*x1+7.1376*x1*x3+391.7883*x2*x3-982145*x1*x2+416476.113;

结果：
目标函数值(最大): 124611.543117018
x1: 0.8
x2: 0.0158529780617937
x3: 600

这样的话，MATLAB、1stopt的计算结果是一致的。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2013-10-22 15:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zh327112793 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 求调剂 +16 MAX怅惘 2026-04-09 18/900 2026-04-10 20:24 by 小鬼喔～～
[考研] 一志愿京区985，085401电子信息，本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 求调剂 +9 翩翩一书生 2026-04-09 9/450 2026-04-10 11:43 by wp06
[考研] 机械专硕270求调剂，接受跨专业 +10 老师看看我吧aba 2026-04-09 11/550 2026-04-10 10:41 by cal0306
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 本科郑州大学，一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[考研] 322求调剂，08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 求调剂，262机械专硕 +6 嗯yyl 2026-04-08 6/300 2026-04-09 12:01 by zhouyuwinner
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 复试调剂 +14 呼呼？~+123456 2026-04-05 14/700 2026-04-06 22:50 by chenzhimin
[考研] 085405软件工程301分求调剂，专硕可跨专业，四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 346分的生物与医药08600求调剂 +6 常雨阳上岸 2026-04-05 7/350 2026-04-06 12:36 by lys0704
[考研] 一志愿211生物学280分 求调剂 +5 李rien 2026-04-05 5/250 2026-04-06 10:30 by zhyzzh
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
信息提示
关闭
请填处理意见
关闭