24小时热门版块排行榜    

查看: 553  |  回复: 2

haibing

铜虫 (正式写手)

[求助] 求1stopt高版本软件帮忙运行程序 已有1人参与

那位大虾帮我用高版本的1stopt运行一下下面两段的程序:微分方程参数拟合:
Variable X,S,P,T;
Parameter r=[0,1],m=[0,1],a=[0,1],b=[0,1],rr=[0,1];
ODEFunction X'=r*X*(1-X/m)-(c*S)/(S+d);
            S'=a*X*(0.008343*P+1.258)-b*(67.39*T/(T+13.53))*a*X*(0.008343*P+1.258)-rr*X;
Data;
0.024        0.0322        108.2        1
0.026        0.0284        196.3        2
0.048        0.0233        168.3        3
0.052        0.0235        281.4        4
0.087        0.0224        144.1        5
0.174        0.018        119.7        6
0.13        0.017        173.6        7
0.164        0.0167        187.1        8
0.1632        0.0166        214.9        9
0.1542        0.0152        121        10
0.1952        0.0142        180.7        11
0.2289        0.0132        181.8        12
0.2585        0.016        201.8        13
0.2888        0.0263        219.9        14
0.3522        0.0198        215.4        15
0.287        0.0145        157.8        16
0.2222        0.0137        200.5        17
0.2323        0.0109        141.2        18
0.2138        0.0113        151.6        19
0.1972        0.0128        225.1        20
0.2105        0.01635        286.6        21
0.1794        0.0151        246.6        22
0.1661        0.01215        136.6        23
0.1822        0.01215        74.4        24
0.1738        0.01185        124.6        25
0.1904        0.0112        271.2        26
0.1612        0.01085        137.9        27
0.1497        0.01045        153        28
0.1422        0.0111        137        29
0.1233        0.0127        185.8        30
0.1026        0.0145        178.6        31
0.1408        0.0146        215.4        32
0.1349        0.01335        157.8        33
0.1368        0.0157        200.5        34
0.1322        0.01        141.2        35
0.12015        0.0103        151.6        36
0.11775        0.0135        225.1        37
0.1083        0.01875        286.6        38
0.1014        0.02045        246.6        39
0.10965        0.0163        136.6        40
0.1035        0.0134        74.4        41
0.09215        0.01595        124.6        42
0.0895        0.0137        271.2        43
0.0924        0.01365        137.9        44
0.09105        0.0173        153        45
0.09765        0.0193        137        46
0.0952        0.0202        185.8        47
0.09495        0.0157        178.6        48
0.09        0.0086        74.4        49
0.09        0.0083        124.6        50
0.091        0.0109        271.2        51
0.093        0.0159        137.9        52
0.093        0.0148        153        53
0.095        0.0143        137        54
0.093        0.0177        185.8        55
0.096        0.0109        178.6        56
=========================================================
第二段程序
Variable X,S,P,T;
Parameter r=[0,1],b=[0,1],rr=[0,1];
ODEFunction X'=r*X*(1-X/0.38)-(c*S)/(S+d);
            S'=0.08*X*(0.008343*P+1.258)-b*(67.39*T/(T+13.53))*0.08*X*(0.008343*P+1.258)-rr*X;
Data;
0.024        0.0322        108.2        1
0.026        0.0284        196.3        2
0.048        0.0233        168.3        3
0.052        0.0235        281.4        4
0.087        0.0224        144.1        5
0.174        0.018        119.7        6
0.13        0.017        173.6        7
0.164        0.0167        187.1        8
0.1632        0.0166        214.9        9
0.1542        0.0152        121        10
0.1952        0.0142        180.7        11
0.2289        0.0132        181.8        12
0.2585        0.016        201.8        13
0.2888        0.0263        219.9        14
0.3522        0.0198        215.4        15
0.287        0.0145        157.8        16
0.2222        0.0137        200.5        17
0.2323        0.0109        141.2        18
0.2138        0.0113        151.6        19
0.1972        0.0128        225.1        20
0.2105        0.01635        286.6        21
0.1794        0.0151        246.6        22
0.1661        0.01215        136.6        23
0.1822        0.01215        74.4        24
0.1738        0.01185        124.6        25
0.1904        0.0112        271.2        26
0.1612        0.01085        137.9        27
0.1497        0.01045        153        28
0.1422        0.0111        137        29
0.1233        0.0127        185.8        30
0.1026        0.0145        178.6        31
0.1408        0.0146        215.4        32
0.1349        0.01335        157.8        33
0.1368        0.0157        200.5        34
0.1322        0.01        141.2        35
0.12015        0.0103        151.6        36
0.11775        0.0135        225.1        37
0.1083        0.01875        286.6        38
0.1014        0.02045        246.6        39
0.10965        0.0163        136.6        40
0.1035        0.0134        74.4        41
0.09215        0.01595        124.6        42
0.0895        0.0137        271.2        43
0.0924        0.01365        137.9        44
0.09105        0.0173        153        45
0.09765        0.0193        137        46
0.0952        0.0202        185.8        47
0.09495        0.0157        178.6        48
0.09        0.0086        74.4        49
0.09        0.0083        124.6        50
0.091        0.0109        271.2        51
0.093        0.0159        137.9        52
0.093        0.0148        153        53
0.095        0.0143        137        54
0.093        0.0177        185.8        55
0.096        0.0109        178.6        56
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

天高任鸟飞,海阔凭鱼跃
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haibing

铜虫 (正式写手)

坛子里没有人有高版本的1stopt吗?求帮忙
天高任鸟飞,海阔凭鱼跃
2楼2014-11-06 09:20:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
微分方程拟合计算比较费时间,下面是第一段代码的参考结果(有时间的话结果还会更好):

均方差(RMSE):0.017640728970135
残差平方和(SSE):0.0342314850457537
相关系数(R): 0.602461362176299
相关系数之平方(R^2): 0.362959692915321
确定系数(DC): -3.36028539491201
F统计(F-Statistic): 69.8177377311799

参数                  最佳估算
--------------------        -------------
c        0.00106428654192994
d        -0.0349011012047928
r        3.60972530286868E-10
m        9.33157757082492E-10
a        0.00332753093676446
b        0.0278157876334486
rr        0.0014446672737485

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 x        计算 x        目标 s        计算 s
1        2        0.026        0.036819251843059        0.0284        0.0323928918894434
2        3        0.048        0.0504793431206788        0.0233        0.0326048751860358
3        4        0.052        0.0655576197566643        0.0235        0.0329407245701933
4        5        0.087        0.0825052174306297        0.0224        0.0331619435163049
5        6        0.174        0.100795442428006        0.018        0.033346624951942
6        7        0.13        0.120787975446028        0.017        0.0335716045574381
7        8        0.164        0.143535218303242        0.0167        0.0337889386727208
8        9        0.1632        0.170114630181354        0.0166        0.0339994979630296
9        10        0.1542        0.19858412512844        0.0152        0.0340509079602557
10        11        0.1952        0.225043607501893        0.0142        0.0340911826461687
11        12        0.2289        0.247476562408349        0.0132        0.0340519379474219
12        13        0.2585        0.262343553713506        0.016        0.033930360088174
13        14        0.2888        0.268735755210552        0.0263        0.0337193453012847
14        15        0.3522        0.267836891537942        0.0198        0.0334143145239424
15        16        0.287        0.261864960294848        0.0145        0.0330303388838698
16        17        0.2222        0.25297900847369        0.0137        0.0325839453327286
17        18        0.2323        0.242727942646882        0.0109        0.0321090741112015
18        19        0.2138        0.232090832828099        0.0113        0.0316021888326738
19        20        0.1972        0.221521024683804        0.0128        0.0310278125687688
20        21        0.2105        0.211226796686975        0.01635        0.0303814698534273
21        22        0.1794        0.201384721932077        0.0151        0.029737451164964
22        23        0.1661        0.192168163806734        0.01215        0.0291740022818871
23        24        0.1822        0.183664080011589        0.01215        0.0286729929308619
24        25        0.1738        0.175829539643054        0.01185        0.0281343466862331
25        26        0.1904        0.168535998576585        0.0112        0.0274510459221761
26        27        0.1612        0.161724726074676        0.01085        0.0268766834354867
27        28        0.1497        0.155410096276067        0.01045        0.02630701093047
28        29        0.1422        0.149543711382441        0.0111        0.0257545687657379
29        30        0.1233        0.144079761827116        0.0127        0.0251582642253836
30        31        0.1026        0.138970741825484        0.0145        0.0245667266219661
31        32        0.1408        0.134183341186696        0.0146        0.0239425583475144
32        33        0.1349        0.12969106566469        0.01335        0.0233783820402628
33        34        0.1368        0.12547546471084        0.0157        0.0227831318167776
34        35        0.1322        0.121511828246697        0.01        0.0222496977090231
35        36        0.12015        0.1177872212671        0.0103        0.021719627760067
36        37        0.11775        0.114272485486526        0.0135        0.0211174136925277
37        38        0.1083        0.110936188795761        0.01875        0.020447330994152
38        39        0.1014        0.107762174781625        0.02045        0.0198169763772046
39        40        0.10965        0.104752119949238        0.0163        0.0193123707600452
40        41        0.1035        0.101909009228637        0.0134        0.0188928965660205
41        42        0.09215        0.0992215745590295        0.01595        0.0184365715356579
42        43        0.0895        0.0966637264522639        0.0137        0.017831130642767
43        44        0.0924        0.0942207813135252        0.01365        0.017346513857017
44        45        0.09105        0.0918949943467024        0.0173        0.0168708559987182
45        46        0.09765        0.0896772250290444        0.0193        0.0164149611221576
46        47        0.0952        0.0875586327235844        0.0202        0.0159191627761857
47        48        0.09495        0.0855293652404816        0.0157        0.0154300542654557
48        49        0.09        0.0835880867476915        0.0086        0.015048237961479
49        50        0.09        0.0817331351387703        0.0083        0.0146369814626851
50        51        0.091        0.0799502666950306        0.0109        0.0140859026699785
51        52        0.093        0.0782316922833754        0.0159        0.0136476948874128
52        53        0.093        0.0765801833523736        0.0148        0.013217663923481
53        54        0.095        0.0749912461179193        0.0143        0.0128058267597263
54        55        0.093        0.0734603614707211        0.0177        0.0123566752817365
55        56        0.096        0.0719821261000118        0.0109        0.0119135640050082
3楼2014-11-06 10:55:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 haibing 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +3 上岸6666@ 2026-03-20 3/150 2026-03-21 16:49 by barlinike
[考研] 336求调剂 +4 rmc8866 2026-03-21 4/200 2026-03-21 16:48 by macy2011
[考研] 求调剂 +4 十三加油 2026-03-21 4/200 2026-03-21 16:20 by barlinike
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 15/750 2026-03-21 10:16 by 脱颖而出
[考研] 313求调剂 +3 肆叁贰壹22 2026-03-19 3/150 2026-03-21 08:01 by JourneyLucky
[考研] 08工科 320总分 求调剂 +6 梨花珞晚风 2026-03-17 6/300 2026-03-21 03:40 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +4 一瓶苯甲酸 2026-03-14 4/200 2026-03-20 20:36 by fen_rao
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
信息提示
请填处理意见