| 查看: 509 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
非线性优化 已有1人参与
|
|||
请各位大神帮忙求解一下以上问题，matlab，mathematica，1stopt都可以，谢谢！ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有60人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
|
是麻烦大家了，我是新手，不是很会这些软件，我用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
onesupeng
金虫 (职业作家)
- 计算强帖: 13
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2116.7
- 散金: 9437
- 红花: 92
- 帖子: 4585
- 在线: 1305小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
2楼2013-12-23 08:36:53
30