[交流] 用1stopt进行非线性常微分方程的参数拟合问题 已有3人参与

本人苦逼小硕一枚，研究传染病传播模型，将模型进行数学变化，得到一个非线性常微分方程组，这个方程组有1个自变量t，7个因变量，3个待估计参数，我下载了1stopt 5.0试用版，但提示试用版无法进行估计，请大神们帮我看看代码是否正确并帮助运行。 Parameter k1,k2,k3 ; Variable t,y2,y3,y4,y5,y7; ODEFunction y1'=y7-k2*y2*y1/y3-0.000019*y1; y2'=y4+y5-0.0349*y2;                        y4'=k1*k2*y2*y1/y3-(0.000019+k1)*y4; y5'=k1*k3*y2*y6/y3-(0.000019+k1)*y5; y6'=0.0348*y2-0.000019*y6-k3*y2*y6/y3; data; 0,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,280,290,300,310,320,330,340,350,360,370,380,390,400,410,420,430,440,450,460,470,480,490,500,510,520,530,540,550,560,570,580,590,600,610,620,630,640,650,660,670,680,690,700,710,720,730,740,750,760,770,780,790,800,810,820,830,840,850,860,870,880,890,900,910,920,930,940,950,960,970,980,990,1000,1010,1020,1030,1040,1050,1060,1070,1080; 37224,38270,38468,38165,39722,40084,39877,40129,40438,39374,40797,40896,40498,41401,41658,41326,41729,41728,41453,41789,41897,41882,43244,43330,42963,43291,43064,42526,43102,42549,41965,41953,41633,41209,40973,40097,39397,40023,40587,40979,41619,41913,42004,42853,43492,43772,44376,44846,44964,45195,45436,45311,45516,45714,45563,45965,46187,46621,47246,47491,46961,46891,46940,46299,45576,45301,44858,43833,43892,43393,42267,42453,41233,40223,40989,41107,41113,41465,42015,41849,42836,42556,42062,43105,43076,43120,43653,43903,43732,44209,44177,43954,44837,44901,45035,45346,46150,45995,46334,46188,45690,46429,45848,45345,45490,44878,44462,44105,44004; 65700000,65707123,65714246,65721369,65728492,65735615,65742738,65749861,65756984,65764107,65771230,65778353,65785476,65792599,65799722,65806845,65813968,65821091,65828214,65835337,65842460,65849583,65856706,65863829,65870952,65878075,65885198,65892321,65899444,65906567,65913690,65920813,65927936,65935059,65942182,65949305,65956428,65965984,65977733,65989482,66001231,66012980,66024729,66036478,66048227,66059976,66071725,66083474,66095223,66106972,66118721,66130470,66142219,66153968,66165717,66177466,66189215,66200964,66212713,66224462,66236211,66247960,66259709,66271458,66283207,66294956,66306705,66318454,66330203,66341952,66353701,66365450,66377199,66388948,66402825,66417072,66431319,66445566,66459813,66474060,66488307,66502554,66516801,66531048,66545295,66559542,66573789,66588036,66602283,66616530,66630777,66645024,66659271,66673518,66687765,66702012,66716259,66730506,66744753,66759000,66773247,66787494,66801741,66815988,66830235,66844482,66858729,66872976,66887223; 2558,1480,855,1843,1784,1849,1990,1595,1002,2294,1634,959,2137,1558,1012,1706,1415,1308,1737,1317,1254,1782,1476,1161,1547,1345,1168,1647,1433,1161,1500,1219,1075,1565,1213,1018,2278,1785,1300,1987,1864,1509,2256,1924,1455,1920,1758,1491,1912,1425,1421,1620,1548,1340,1747,1504,1265,1567,1451,1339,1290,1458,1300,1224,1566,1285,912,1587,1061,828,1537,1102,719,2290,1476,1315,1743,1391,1560,2483,1636,1174,2184,1640,1382,2090,1429,1386,1731,1278,1180,1826,1437,1168,1716,1351,1218,1679,1319,1304,1889,1480,1188,1570,1255,1071,1636,1168,1047; 128,61,26,84,101,108,77,63,44,121,53,36,115,63,34,72,56,45,71,53,57,75,48,40,56,46,58,66,41,50,52,35,33,45,48,34,117,79,51,89,82,68,100,98,77,82,90,65,97,63,52,74,61,52,56,71,53,64,62,40,51,80,51,52,60,46,34,74,22,25,72,55,38,116,62,59,66,52,79,100,60,43,80,62,54,95,50,49,69,48,41,71,40,34,50,50,33,58,40,43,71,44,27,52,35,24,52,31,35; 24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24030,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24540,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555,24555; 1.jpg

回复此楼