24小时热门版块排行榜    

Znn3bq.jpeg
查看: 832  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +3 qut2026 2026-04-09 5/250 2026-04-10 22:52 by qut2026
[考研] 生物学调剂 可调剂到生物与医药 +7 李政莹 2026-04-06 8/400 2026-04-10 21:15 by zhouxiaoyu
[考研] 284求调剂 +11 archer.. 2026-04-10 12/600 2026-04-10 18:57 by HPUCZ
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 求调剂 +7 璃茉一定上岸 2026-04-10 7/350 2026-04-10 15:25 by jiajinhpu
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 267求调剂 +5 再忙也要吃饭啊 2026-04-09 5/250 2026-04-09 18:47 by stone_128
[考研] 求调剂希望还是希望在山河四省附近 +3 快乐的小白鸽 2026-04-05 3/150 2026-04-09 17:36 by wp06
[考研] 085600材料与化工,求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 求调剂 +15 熊二想上岸 2026-04-06 15/750 2026-04-08 04:53 by 无际的草原
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 一志愿北京化工085600 310分求调剂 +20 0856材料与化工3 2026-04-04 22/1100 2026-04-07 15:14 by 上岸快快
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
信息提示
请填处理意见