| 查看: 308 | 回复: 2 | ||
[求助]
非线性优化 已有1人参与
|
请各位大神帮忙求解一下以上问题，matlab，mathematica，1stopt都可以，谢谢！ |
回复此楼» 猜你喜欢
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有271人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有1人回复
-
投稿chemical physical letters不送审？
已经有2人回复
onesupeng
金虫 (职业作家)
- 计算强帖: 13
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2336.2
- 散金: 9224
- 红花: 92
- 帖子: 4583
- 在线: 1303.8小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
2楼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]); |
3楼2013-12-23 10:40:27