24小时热门版块排行榜    

查看: 824  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by 余麟余
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
信息提示
请填处理意见