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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助，计算方程最大值
61/1返回列表
查看: 1146  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 广东省 085601 329分求调剂 +13 Eddieddd 2026-04-10 13/650 2026-04-10 23:27 by 314126402
[考研] 本科西工大 324求调剂 +4 wysyjs25 2026-04-10 4/200 2026-04-10 20:00 by 来看流星雨10
[考研] 计算机类求调剂，22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 269电子信息求调剂，可转专业 +9 独酌wl 2026-04-06 9/450 2026-04-09 20:55 by laoshidan
[考研] 考研调剂-材料类-284 +28 想换手机不想解� 2026-04-08 28/1400 2026-04-09 20:08 by 倒数321?
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 336求调剂，一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭