| 查看: 1836 | 回复: 13 | |||
[求助]
使用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 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有281人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大侠能帮忙跑下代码！！！！！1stopt下载安装后使用不了。。。。。
已经有4人回复
- 1stopt自动搜索功能怎么老卡住不动？？？
已经有12人回复
- 1stOpt软件微分方程参数拟合，运行后没有任何反应。
已经有5人回复
- 编程模式下1stopt如何使用积分拟合
已经有7人回复
- 1stopt程序运行
已经有3人回复
- 1stopt文本程序或matlab程序求助
已经有0人回复
- 用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来
已经有4人回复
- 各位前辈，利用1stOpt对非线性方程拟合每次运算结果都不一样怎么办呢？
已经有9人回复
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 求助，有1stopt 5.0的虫友们，帮忙跑一下代码，谢谢了！！
已经有2人回复
- 紧急求助，有1stopt高版本的师兄师姐帮忙做一个 SharedModel 拟合
已经有2人回复
- 哪位朋友有注册版1stOpt软件？急需
已经有5人回复
- 1stopt 安装
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
8楼2016-04-22 07:26:40
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★
海绵宝宝mm: 金币+5, ★★★很有帮助 2016-04-22 19:17:05
海绵宝宝mm: 金币+5, ★★★很有帮助 2016-04-22 19:17:05
|
下面代码供参考： 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 |
13楼2016-04-22 15:53:18
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2016-04-21 16:10:54
3楼2016-04-21 16:59:26
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2016-04-21 17:26:33
5楼2016-04-21 17:31:52
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2016-04-21 20:02:17
|
就是因为解方程比较烦才想要用算法来拟合。。我用这个软件中的pso算法成功算过另一个模型的两个参数，它的计算公式稍微简单一点。。就是不知道这个为什么不行。N只有下限，上限未知。 发自小木虫Android客户端 |
7楼2016-04-21 22:33:54
9楼2016-04-22 09:48:53
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
10楼2016-04-22 09:59:27
10