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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助，计算方程最大值
61/1返回列表
查看: 1043  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭