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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助，计算方程最大值
61/1返回列表
查看: 1134  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 280求调剂 +7 咕噜晓晓 2026-03-18 8/400 2026-03-21 01:27 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 316求调剂 +5 梁茜雯 2026-03-19 5/250 2026-03-20 21:26 by 你好你好666
[考研] 求调剂 +5 Mqqqqqq 2026-03-19 5/250 2026-03-20 20:46 by zhukairuo
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
关闭
请填处理意见
关闭