24小时热门版块排行榜

返回列表

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

海绵宝宝mm

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 13
帖子: 8
在线: 1.7小时
虫号: 4606271
注册: 2016-04-21
专业: 计算机软件

[求助] 使用1stopt运行没有报错，但是出不了结果已有1人参与

刚接触1stopt软件，运行下面代码时没有报错但是出不了结果，不知道是啥问题，请各位大神给看看，感激不尽
具体公式见附件
代码如下：
Title "JM模型参数估计";
Constant n=136;
Parameters p[0,1];N[136,];
Variable t(1:n);
Function  p=n/(N*Sum(i=1:n)(t)-Sum(i=1:n)((i-1)*t));
      Sum(i=1:n)(1/(N-(i-1)))=n/(N-(1/Sum(i=1:n)(t))*Sum(i=1:n)((i-1)*t));
data;
1 3
2 30
3 113
4 81
5 115
6 9
7 2
8 91
9 112
10 15
11 138
12 50
13 77
14    24
15 108
16 88
17 670
18 120
19 26
20 114
21 325
22 55
23 242
24 68
25 422
26 180
27 10
28 1146
29 600
30 15
31 36
32 4
33 0
34 8
35 227
36 65
37 176
38 58
39 457
40 300
41 97
42 263
43 452
44 255
45 197
46 193
47 6
48 79
49 816
50 1351
51 148
52 21
53 233
54 134
55 357
56 193
57 236
58 31
59 369
60 748
61 0
62 232
63 330
64 365
65 1222
66 543
67 10
68 16
69 529
70 379
71 44
72 129
73 810
74 290
75 300
76 529
77 281
78 160
79 828
80 1011
81 445
82 296
83 1755
84 1064
85 1783
86 860
87 983
88 707
89 33
90 868
91 724
92 2323
93 2930
94 1461
95 843
96 12
97 261
98 1800
99 865
100 1435
101 30
102 143
103 108
104 0
105 3110
106 1247
107 943
108    700
109 875
110 245
111 729
112 1897
113 447
114 386
115 446
116 122
117 990
118    948
119 1082
120 22
121 75
122 482
123 5509
124 100
125 10
126 1071
127 371
128 790
129 6150
130 3321
131 1045
132 648
133 5485
134 1160
135 1864
136 4116

JM模型参数计算公式.png

回复此楼

» 猜你喜欢

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

1楼 2016-04-21 14:00:35

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★
海绵宝宝mm: 金币+5, ★★★很有帮助 2016-04-22 19:17:05

下面代码供参考：

CODE:

Constant n=136;

Constant t(n)=[3,30,113,81,115,9,2,91,112,15,138,50,77,24,108,88,670,120,26,114,325,55,242,68,422,180,10,1146,600,15,36,4,0,8,227,65,176,58,457,300,97,263,452,255,197,193,6,79,816,1351,148,21,233,134,357,193,236,31,369,748,0,232,330,365,1222,543,10,16,529,379,44,129,810,290,300,529,281,160,828,1011,445,296,1755,1064,1783,860,983,707,33,868,724,2323,2930,1461,843,12,261,1800,865,1435,30,143,108,0,3110,1247,943,700,875,245,729,1897,447,386,446,122,990,948,1082,22,75,482,5509,100,10,1071,371,790,6150,3321,1045,648,5485,1160,1864,4116];

Parameter M;

PassParameter p=n/(M*Sum(i=1:n,t)(t)-Sum(i=1:n,t)((i-1)*t));

Function Sum(i=1:n,t)(1/(M-(i-1)))=n/(M-(1/Sum(i=1:n,t)(t))*Sum(i=1:n,t)((i-1)*t));

N用M代替。

如果M没有范围约束的话，多解：
1：
目标函数值(最小): 0
m: -1.22575693225421E17

传递参数(PassParameter):
p: -1.25112030157401E-20

2：
目标函数值(最小): 0
m: -1.65191601797123E16

传递参数(PassParameter):
p: -9.28357958912305E-20

3：
目标函数值(最小): 0
m: -3.41327139364123E16

传递参数(PassParameter):
p: -4.49296058202489E-20

如果M范围为[0,1]，结果稳定唯一：
目标函数值(最小): 4.93038065763132E-32
m: 0.220715983682821

传递参数(PassParameter):
p: -1.56768295994417E-5