应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 使用1stopt运行没有报错，但是出不了结果
52/1返回列表
查看: 1795  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

海绵宝宝mm

新虫 (初入文坛)

[求助] 使用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

使用1stopt运行没有报错，但是出不了结果
JM模型参数计算公式.png
回复此楼

» 猜你喜欢
1楼 2016-04-21 14:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
海绵宝宝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
赞一下(1人) 回复此楼
13楼2016-04-22 15:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
看不大明白，想求什么啊？拟合问题吗？哪个是自变量、哪个是因变量、哪个又是待求未知参数？注意1stOpt不区分大小写，因此n和N其中之一最好用不同符号。
赞一下 回复此楼
2楼2016-04-21 16:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

海绵宝宝mm

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by dingd at 2016-04-21 16:10:54
看不大明白，想求什么啊？拟合问题吗？哪个是自变量、哪个是因变量、哪个又是待求未知参数？注意1stOpt不区分大小写，因此n和N其中之一最好用不同符号。

想求参数p和N，t(1:n)也就是下面的136个数据是变量。。。刚我把N换成其他字母了，还是不行

发自小木虫Android客户端
赞一下 回复此楼
3楼2016-04-21 16:59:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

晕啊，看不懂原题，目标函数是什么？有明白的给解释下。
赞一下 回复此楼
4楼2016-04-21 17:26:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +61 AI8RaGaPaSCR 2026-08-07 64/3200 2026-08-14 05:45 by yNk0cTSDcreS
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
关闭
请填处理意见
关闭