24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

123-smiling

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 168.5
帖子: 14
在线: 5.7小时
虫号: 2022771
注册: 2012-09-23
专业: 工业工程与管理

[求助] 非线性优化已有1人参与

请各位大神帮忙求解一下以上问题，matlab，mathematica，1stopt都可以，谢谢！

回复此楼

» 猜你喜欢

新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有117人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复

1楼 2013-12-23 06:27:55

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

123-smiling

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 168.5
帖子: 14
在线: 5.7小时
虫号: 2022771
注册: 2012-09-23
专业: 工业工程与管理

引用回帖:

2楼: Originally posted by onesupeng at 2013-12-23 08:36:53
很多软件可以做这个问题，但是你要给代码这就为难大家了。因为你这个没有本质的困难，软件找来研究一下就能做到了

是麻烦大家了，我是新手，不是很会这些软件，我用mathematica和1stopt解决不了，可能递归函数里面出了问题，这是我的程序，麻烦大家纠正。
1.
Mathematica程序：

Maximize[{w*S[12]*(1 - P[12]) + (1-w)*S[16]*(1 - P[16]),
  S[1] == 10 && R == 100 && k == 0.6 &&
S[2] == S[1]*E^-(S[1]/(R*(r[1]))) &&
S[3] == S[2]*E^-(S[2]/(R*(r[2]))) &&
S[4] == S[3]*E^-(S[3]/(R*(r[3]))) &&
S[5] == S[4]*E^-(S[4]/(R*(r[4]))) &&
S[6] == S[5]*E^-(S[5]/(R*(r[5]))) &&
S[7] == S[6]*E^-(S[6]/(R*(r[6]))) &&
S[8] == S[7]*E^-(S[7]/(R*(r[7]))) &&
S[9] == S[8]*E^-(S[8]/(R*(r[8]))) &&
S[10] == S[9]*E^-(S[9]/(R*(r[9]))) &&
S[11] == S[10]*E^-(S[10]/(R*(r[10]))) &&
S[12] == S[11]*E^-(S[11]/(R*(r[11]))) &&
S[13] == k*S[12]*E^-(S[12]/(R*(r[12]))) &&
S[14] == S[13]*E^-(S[13]/(R*(r[13]))) &&
S[15] == S[14]*E^-(S[14]/(R*(r[14]))) &&
S[16] == S[15]*E^-(S[15]/(R*(r[15]))) &&
P[1] == 1 - E^-(S[1]/(R*(r[1]))) &&
P[2] == 1 - E^-(S[2]/(R*(r[2]))) &&
P[3] == 1 - E^-(S[3]/(R*(r[3]))) &&
P[4] == 1 - E^-(S[4]/(R*(r[4]))) &&
P[5] == 1 - E^-(S[5]/(R*(r[5]))) &&
P[6] == 1 - E^-(S[6]/(R*(r[6]))) &&
P[7] == 1 - E^-(S[7]/(R*(r[7]))) &&
P[8] == 1 - E^-(S[8]/(R*(r[8]))) &&
P[9] == 1 - E^-(S[9]/(R*(r[9]))) &&
P[10] == 1 - E^-(S[10]/(R*(r[10]))) &&
P[11] == 1 - E^-(S[11]/(R*(r[11]))) &&
P[12] == 1 - E^-(S[12]/(R*(r[12]))) &&
P[13] == 1 - E^-(S[13]/(R*(r[13]))) &&
P[14] == 1 - E^-(S[14]/(R*(r[14]))) &&
P[15] == 1 - E^-(S[15]/(R*(r[15]))) &&
P[16] == 1 - E^-(S[16]/(R*(r[16]))) &&
r[1] + r[2] + r[3] + r[4] + r[5] + r[6] + r[7] + r[8] + r[9] +
   r[10] + r[11] + r[12] + r[13] + r[14] + r[15] + r[16] == 1 &&
0 <= r[1] <= 1 && 0 <= r[2] <= 1 && 0 <= r[3] <= 1 &&
0 <= r[4] <= 1 && 0 <= r[5] <= 1 && 0 <= r[6] <= 1 &&
0 <= r[7] <= 1 && 0 <= r[8] <= 1 && 0 <= r[9] <= 1 &&
0 <= r[10] <= 1 && 0 <= r[11] <= 1 && 0 <= r[12] <= 1 &&
0 <= r[13] <= 1 && 0 <= r[14] <= 1 && 0 <= r[15] <= 1 &&
0 <= r[16] <= 1}, {r[1], r[2], r[3], r[4], r[5], r[6], r[7], r[8],
  r[9], r[10], r[11], r[12], r[13], r[14], r[15], r[16]}]

2.1stopt程序
Constant n=16,RR=4000,k=0.6;
ConstStr S[1]=200,for i := 1 to n do P=1-exp(-S/RR/r);
Parameter r(1:n)[0,1], S(2:16)[0,],w;
MaxFunction w*S12*(1-P12)+(1-w)*S16*(1-P16);
         Sum(i=1:n)(r)=1;
         S[2]=S[1]*(1-P[1]);
         S[3]=S[2]*(1-P[2]);
         S[4]=S[3]*(1-P[3]);
         S[5]=S[4]*(1-P[4]);
         S[6]=S[5]*(1-P[5]);
         S[7]=S[6]*(1-P[6]);
         S[8]=S[7]*(1-P[7]);
         S[9]=S[8]*(1-P[8]);
         S[10]=S[9]*(1-P[9]);
         S[11]=S[10]*(1-P[10]);
         S[12]=S[11]*(1-P[11]);
         S[13]=S[12]*k*(1-P[12]);
         S[14]=S[13]*(1-P[13]);
         S[15]=S[14]*(1-P[14]);
         S[16]=S[15]*(1-P[15]);