| 查看: 575 | 回复: 2 | ||
[求助]
非线性优化 已有1人参与
|
请各位大神帮忙求解一下以上问题，matlab，mathematica，1stopt都可以，谢谢！ |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有138人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
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
