24小时热门版块排行榜    

CyRhmU.jpeg
查看: 755  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见