24小时热门版块排行榜    

查看: 939  |  回复: 6

changhw

铁虫 (小有名气)

[求助] 求帮忙跑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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

changhw

铁虫 (小有名气)

2楼2015-06-13 14:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
3楼2015-06-13 22:50:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

changhw

铁虫 (小有名气)

引用回帖:
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的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
4楼: Originally posted by changhw at 2015-06-14 12:41:20
他这个跑起来好像有迭代次数的限制,到后面的直接报是试用版了...

迭代次数是自己设置的,计算停止要么是达到了收敛标准,要么是达到了最大迭代次数,这里的最大迭代次数就是自己设置的。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
5楼2015-06-15 11:47:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

changhw

铁虫 (小有名气)

引用回帖:
5楼: Originally posted by 月只蓝 at 2015-06-15 11:47:42
迭代次数是自己设置的,计算停止要么是达到了收敛标准,要么是达到了最大迭代次数,这里的最大迭代次数就是自己设置的。...

哦哦,那可能是我设置不对,它迭代了四十几次就报 试用版 了
6楼2015-06-15 12:23:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by changhw at 2015-06-15 12:23:00
哦哦,那可能是我设置不对,它迭代了四十几次就报 试用版 了...

你查看一下计算停止的原因是什么,如果四十几次应该是收敛了。试用版的软件,计算结束就是提示是试用版,不代表计算结果不对。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2015-06-15 12:55:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 changhw 的主题更新
信息提示
请填处理意见