24小时热门版块排行榜    

查看: 780  |  回复: 2

CloudYuan

铁虫 (初入文坛)

[求助] 1stOpt求解168个未知数的非线性方程组,求可行解,或者改进建议 已有1人参与

模型是一个有168个未知数的方程组,求教高手有无可行解,或者方程组存在什么问题。具体请见附件,模型结构如下:
@@@@@@@@@
440<=a*(X1^2)*(Y1^1)*Z1+b*(X1^2)*(Y1^2)*Z1+c*(X2^1)*(Y2^4)*Z2+...<=450;
1.待求的Zn和Xn,Yn一一对应,X、Y、Z一共56对,计168个;
2.不同的Zn与Xn、Yn的多次方(0-5次方、且二者之和不大于5)的乘积之和在一个区间内,预计最长的式子包含有56个“+”;
3.a, b, c...为常数。
4.理论上应有(以下模型约束似乎无法添加入1.5版中):
x2>x1;

x4>x3;
x5>x4;
x6>x5;
x7>x6;

x9>x8;
x10>x9;
x11>x10;
x12>x11;
x13>x12;

x15>x14;
x16>x15;
x17>x16;
x18>x17;
x19>x18;
x20>x19;

x22>x21;
x23>x22;
x24>x23;
x25>x24;
x26>x25;
x27>x26;
x28>x27;

x30>x29;
x31>x30;
x32>x31;
x33>x32;
x34>x33;

x36>x35;
x37>x36;
x38>x37;
x39>x38;
x40>x39;

x42>x41;
x43>x42;
x44>x43;
x45>x44;
x46>x45;

x48>x47;
x49>x48;
x50>x49;
x51>x50;

x53>x52;
x54>x53;
x55>x54;
x56>x55;
@@@@@@@@@@@@@
用v1.5尝试过很多次,结果都很怪,不知道高版本是否有可行解,或者高手有什么建议?

v1.5中使用的解法是
差分进化法DE
子模式3-Vi=X3+F(Xbest+X1-X2-X3)
种群数20
变异率0.95
(有时采用局部搜索,步长为2)
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 1stOpt.txt
  • 2014-12-17 23:55:39, 118.6 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
CloudYuan: 金币+10, 有帮助, 感谢! 2015-01-07 18:31:04
x1..x56, y1..y56, z1..z56, mz277..mz534,一共327个未知数?
2楼2014-12-18 10:05:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CloudYuan

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2014-12-18 10:05:20
x1..x56, y1..y56, z1..z56, mz277..mz534,一共327个未知数?

mz277至mz534是我用来表示区间范围的,初衷是在function中表达:
  x*y*z+x^2*y*z...<1500;
  x*y*z+x^2*y*z...>0;

但我不知道在v1.5版怎么书写才能让上面的不等式正确运行,所以就写成:
  parameter中,mz277(0,1500);
  function中,
  x*y*z+x^2*y*z...=mz277;
3楼2014-12-29 18:56:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 CloudYuan 的主题更新
信息提示
请填处理意见