24小时热门版块排行榜

返回列表

deepseawdd

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 128.1
帖子: 27
在线: 20.5小时
虫号: 1797383
注册: 2012-05-05
专业: 结构工程

[求助] 能帮忙看一下吧，为什么用1stopt5.5试用版，点运行没反应啊已有1人参与

能帮忙看一下吧，为什么用1stopt5.5试用版，点运行没反应啊
ConstStr k=t^b*a*0.473/(2.081-1.608*t);
Parameters a,b;
Variable t,r;
ODEFunction r''=-r'*(1+10*k'/k)/10;
Data;
0.228398899 8.32
0.242900416 13.47
0.255672483 15.64
0.267082228 17.69
0.277390305 18.64
0.286788667 18.96
0.295422884 19.35
0.303406043 20.17
0.310827786 21.69

回复此楼

» 猜你喜欢

基元I理论下三大核心空间现象精准推导与细节解析已经有0人回复
基于基元 I 统一理论的反重力理论推导已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有78人回复
基于基元I统一理论的量子力学本源推导已经有1人回复
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费已经有20人回复
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）已经有0人回复
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）已经有0人回复

1楼 2017-06-30 21:43:30

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
deepseawdd: 金币+5 2017-07-02 20:17:07

先把k的一阶导数求出来代替k':
k'=((((b*(t^(b-1)))*a)*0.473)*(2.081-1.608*t)-(((t^b)*a)*0.473)*(-1.608))/((2.081-1.608*t)^2);

k'由dk代替表示：

CODE:

ConstStr k=t^b*a*0.473/(2.081-1.608*t),

         dk=((((b*(t^(b-1)))*a)*0.473)*(2.081-1.608*t)-(((t^b)*a)*0.473)*(-1.608))/((2.081-1.608*t)^2);

Parameters a,b;

Variable t,r;

ODEFunction r''=-r'*(1+10*dk/k)/10;

Data;

0.228398899        8.32

0.242900416        13.47

0.255672483        15.64

0.267082228        17.69

0.277390305        18.64

0.286788667        18.96

0.295422884        19.35

0.303406043        20.17

0.310827786        21.69

a参数似乎无用了，dk/k把a约掉了，也即a可以为任意值。

均方差(RMSE): 0.474551966680738
残差平方和(SSE): 1.80159655264445
相关系数(R): 0.981564892013995
相关系数之平方(R^2): 0.963469637234446
修正R平方(Adj. R^2): 0.948857492128224
确定系数(DC): 0.961944496067107
F统计(F-Statistic): 68.2782010981496

参数最佳估算
-------------------- -------------
a -67.2023509470802
b 7.18754485098689
r' 初值 406.49264251293