24小时热门版块排行榜    

查看: 974  |  回复: 3
【奖励】 本帖被评价1次,作者dingd增加金币 0.8

dingd

铁杆木虫 (职业作家)


[资源] 过指定点的微分方程拟合

1stopt 5.0目前可计算过指定点的微分方程拟合,代码格式与风格与一般拟合区别不大,非常容易理解和实现:
已知微分方程:dy/dx=-k*y^n;
k、n为待求参数,x、y数据如下:
x=[0,20,40,60,120,180,300];
y=[10,8,6,5,3,2,1];
要求必须过第二和最后一个点,(20,8)和(300,1)

1stOpt代码
CODE:
Variable x,y;
SubjectTo y[20]-8=0,y[300]-1=0;
ODEFunction y'=-k*y^n;
Data;
x=[0,20,40,60,120,180,300];
y=[10,8,6,5,3,2,1];

结果:
均方差(RMSE): 0.263804091013414
残差平方和(SSE): 0.417555590612481
相关系数(R): 0.997392048218919
相关系数之平方(R^2): 0.99479089785033
决定系数(DC): 0.988012758164235
F统计(F-Statistic): 333.688103969476

参数                  最佳估算
--------------------        -------------
k        0.00534339757163354
n        1.33612562125854

微分方程拟合约束(SubjectTo):
y[20]-8-0: 7.4892980705954E-11
y[300]-1-0: -3.49031470392447E-10

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

文件:数据文件 - 1
No        x        目标y        计算y
1        20        8        8.00000000007489
2        40        6        6.50046220604109
3        60        5        5.35401725694302
4        120        3        3.19568964029055
5        180        2        2.05891109748135
6        300        1        0.999999999650969
回复此楼

» 收录本帖的淘帖专辑推荐

科研工具专辑

» 猜你喜欢

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

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

葡萄酱油茶

捐助贵宾 (初入文坛)


★★★★★ 五星级,优秀推荐

学生党要用一下1stopt拟合微分方程,正版买不起啊,大神能帮泡一下么,万分感谢
2楼2015-05-02 21:44:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葡萄酱油茶

捐助贵宾 (初入文坛)


关于微分方程参数拟合,数据都有
3楼2015-05-02 22:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葡萄酱油茶

捐助贵宾 (初入文坛)


大神求助
4楼2015-05-02 22:52:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dingd 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上项目filecode邪修 +3 西山十月 2026-08-09 3/150 2026-08-09 19:39 by Tide man
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见