24小时热门版块排行榜    

查看: 995  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂323材料与化工 +3 1124361 2026-03-24 3/150 2026-03-24 20:24 by peike
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +15 吃吃吃才有意义 2026-03-19 16/800 2026-03-24 19:51 by 了了了了。。
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
信息提示
请填处理意见