24小时热门版块排行榜

返回列表

tircool

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6
帖子: 2
在线: 2.7小时
虫号: 4010122
注册: 2015-08-06

[求助] 1stOpt5.0代码问题，一运行就未响应，新手求指教已有1人参与

我的代码如下：

Title "优化180908";

Constant n = 2, al = 0.6, ak = 0.4, a0 = 0, rho2 = 0.8376, sigxi2 = 0.1643;

Parameter p1(1:n) = [0.001, 5], p2(1:n) = [0.001, 5];

Variable A, B, C, D, E, F;
DataSet;
A, B, C, D, E, F =
0.247392986224052 -0.0250052158625110 2.87032657566461 -0.0355787864255006 -0.0351595036550584 5.38103159249501
0.605786280504227 0.00230152237385910 2.47760642515506 0.432865574172738 -0.111377081083603 8.05827308706437
EndDataSet;

ConstStr l = ((sigxi2^2)/2 + ln(al) + a0 + rho2*A - B + ln(exp(p1)-1) + ak*C)/(1-al);
ConstStr y = a0 + al*l + ak*C + D + E;
ConstStr PQ = sum(i = 1:n)(exp(y + p1));
ConstStr Q = PQ/exp(p2);
ConstStr p11 = (-1/F)*y + p2 + (1/F)*ln(Q);

MinFunction sum(i = 1:n)((p2 - p11)*y);

新手编程虚心求教！

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有226人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

1楼 2018-09-10 20:52:19

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+1 2018-09-10 23:25:03
独孤神宇: 鼓励交流 2018-10-23 08:42:21
tircool(独孤神宇代发): 金币+5, 最佳答案 2018-10-23 08:43:42

1stOpt不区分大小写，改下如下：

CODE:

Title "优化180908";

Constant n = 2, al = 0.6, ak = 0.4, a0 = 0, rho2 = 0.8376, sigxi2 = 0.1643;

Parameter p1(1:n) = [0.001, 5], p2(1:n) = [0.001, 5];

DataSet;

AA, BB, CC, DD, EE, FF =

0.247392986224052        -0.0250052158625110        2.87032657566461        -0.0355787864255006        -0.0351595036550584        5.38103159249501

0.605786280504227        0.00230152237385910        2.47760642515506        0.432865574172738        -0.111377081083603        8.05827308706437

EndDataSet;

ConstStr l =  ((sigxi2^2)/2 + ln(al) + a0 + rho2*AA - BB + ln(exp(p1)-1) + ak*CC)/(1-al);

ConstStr y = a0 + al*l + ak*CC + DD + EE;

ConstStr PQ = sum(i = 1:n)(exp(y + p1));

ConstStr Q = PQ/exp(p2);

ConstStr pp = (-1/FF)*y +  p2 + (1/FF)*ln(Q);

MinFunction sum(i = 1:n,p1,p2,AA,BB,CC,DD,EE,FF)((p2 - pp)*y);

Objective Function (Min.): -17.7397811486867
p11: 5
p12: 5
p21: 0.00100000000000166
p22: 0.001

赞一下

回复此楼

2楼2018-09-10 21:48:26

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

tircool

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6
帖子: 2
在线: 2.7小时
虫号: 4010122
注册: 2015-08-06

引用回帖:

2楼: Originally posted by dingd at 2018-09-10 21:48:26
1stOpt不区分大小写，改下如下：

Title "优化180908";
Constant n = 2, al = 0.6, ak = 0.4, a0 = 0, rho2 = 0.8376, sigxi2 = 0.1643;
Parameter p1(1:n) = , p2(1:n) = ;
DataSet;
AA ...

我想问一下为什么我带入了您的代码之后，仍然未响应呢？是电脑版本问题吗？还是因为我安装的是5.0不正式？