24小时热门版块排行榜

返回列表

changhw

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 210.6
散金: 300
帖子: 69
在线: 96.2小时
虫号: 3921951
注册: 2015-06-09
专业: 等离子体物理

[求助] 求帮忙跑1stopt程序，计算参数k,Q的值已有2人参与

需要跑下拟合计算参数k,Q的值，没有高版本的1stopt，主要就是那两个微分方程和一组变量数据，谢谢

@dingd；

Constant m=0.0142,w0=40,b=4.26,g=9.8;
Parameters k,Q;
Variable t,x,y;
ODEFunction x''=-k/m*(x')^2-Q/m*(w0-b*t)*y';
y''=-k/m*(y')^2+Q/m*(w0-b*t)*x'-g;
data;
t,x,y
0.000000 0.174032 0.076237
0.031250 0.264437 0.217998
0.062500 0.336230 0.336128
0.097917 0.400702 0.473140
0.139583 0.434516 0.605399
0.164583 0.454967 0.671532
0.225000 0.472212 0.806135
0.270833 0.471369 0.874608
0.345833 0.451653 0.940701
0.410417 0.446928 0.940696
0.452083 0.440362 0.898188
0.493750 0.443045 0.872218
0.550000 0.453310 0.806116
0.589583 0.472939 0.747107
0.631250 0.499860 0.671577
0.666667 0.521939 0.605486
0.706250 0.541626 0.541755
0.768750 0.609412 0.409598
0.804167 0.650393 0.343527
0.866667 0.739415 0.213752
0.908333 0.806386 0.147707
0.960417 0.897451 0.043908
1.000000 0.969207 -0.026855
1.045833 1.036121 -0.088178
1.112500 1.136578 -0.187246
1.162500 1.213060 -0.258004
1.206250 1.282394 -0.324047
1.256250 1.356367 -0.383002
1.304167 1.418468 -0.437246
1.397917 1.549904 -0.557533
1.468750 1.645723 -0.663689
1.533333 1.724625 -0.739167
1.610417 1.815632 -0.838244
1.710417 1.944850 -0.970339

回复此楼

» 猜你喜欢

散金币祈福已经有186人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有273人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有114人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有39人回复

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

1楼 2015-06-13 14:16:54

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

changhw

铁虫 (小有名气)

@dingd

2楼2015-06-13 14:18:39

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
changhw: 金币+8, ★★★★★最佳答案 2015-06-16 00:47:34

Root of Mean Square Error (RMSE):0.193317106130155
Sum of Square Error:2.46651923248748
Correlation Coef. (R): 0.972400442408196
R-Square: 0.945562620395654
Determination Coef. (DC): 0.87077858107221
F-Statistic: 88.807625755435

Parameter Best Estimate
-------------------- -------------
k -0.0281819861554189
q -0.000387446301479405
x' Initial Value -0.243025584163639
y' Initial Value 2.17962701790004

赞一下(1人)

回复此楼

3楼2015-06-13 22:50:14

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

changhw

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 210.6
散金: 300
帖子: 69
在线: 96.2小时
虫号: 3921951
注册: 2015-06-09
专业: 等离子体物理

引用回帖:

3楼: Originally posted by dingd at 2015-06-13 22:50:14
Root of Mean Square Error (RMSE):0.193317106130155
Sum of Square Error:2.46651923248748
Correlation Coef. (R): 0.972400442408196
R-Square: 0.945562620395654
Determination Coef. (DC): 0.8707785810 ...

他这个跑起来好像有迭代次数的限制，到后面的直接报是试用版了

赞一下

回复此楼

4楼2015-06-14 12:41:20

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