24小时热门版块排行榜    

Znn3bq.jpeg
查看: 846  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
信息提示
请填处理意见