24小时热门版块排行榜    

查看: 2344  |  回复: 18

allengjn

新虫 (正式写手)

引用回帖:
10楼: Originally posted by 月只蓝 at 2014-11-06 21:39:29
代码参见9楼给出的链接中第4楼 dingd 的回复。...

这个一阶导可以用x',如果是二阶导,比如我标题里的式子,那怎么写?是x后边跟两个单引号,还是x后边一个双引号?
写成x''=(Q-C*x'-K*x)/M    这种形式吗?
还有,使用1stopt拟合,按那个键开始拟合?是”执行“、”快速公式拟合搜索“、”全执行“、”热执行“这些按键的哪个?需要先保存在运行吗?
11楼2014-11-06 22:08:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
allengjn: 金币+10, ★★★★★最佳答案 2014-11-06 22:27:33
1stOpt代码:
CODE:
Constant M=0.3986,Q=100;
Variable t,x,x',x'';
ODEFunction x''= (Q-C*x'-K*x)/M;
Data;
//时间        位移x        速度dx        加速度ddx
0        0        0        26663.2
0.000500005        0.000137226        0.382478        281.033
0.00100001        0.000360148        0.505711        229.476
0.00150001        0.000640082        0.612797        212.787
0.00200002        0.000965359        0.664756        -9.53666
0.00250002        0.00130017        0.676495        -53.0477
0.00300003        0.0016186        0.587299        -131.129
....

均方差(RMSE):57.5493822281378
残差平方和(SSE):1987158.83690418
相关系数(R): 0.291877273942345
相关系数之平方(R^2): 0.0851923430440146
确定系数(DC): -2.20653576106198
F统计(F-Statistic): 25.27741280379

参数                  最佳估算
--------------------        -------------
c        49.6465860963742
k        70046293.1591487

其实如果不考虑加速度ddx的话,效果似乎更好:
CODE:
Constant M=0.3986,Q=100;
Variable t,x,x';
ODEFunction x''= (Q-C*x'-K*x)/M;
Data;
//时间        位移x        速度dx
0        0        0
0.000500005        0.000137226        0.382478
0.00100001        0.000360148        0.505711
0.00150001        0.000640082        0.612797
0.00200002        0.000965359        0.664756
0.00250002        0.00130017        0.676495
0.00300003        0.0016186        0.587299
0.00350003        0.00191261        0.610209
.....

均方差(RMSE):0.134482150106052
残差平方和(SSE):7.23417947885864
相关系数(R): 0.905464687918125
相关系数之平方(R^2): 0.819866301066668
确定系数(DC): -0.231632439076251
F统计(F-Statistic): 90.6285375843916

参数                  最佳估算
--------------------        -------------
c        24.0628939603168
k        11361.0933692284
使用MATLAB拟合动力学参数问题
c122.jpg


使用MATLAB拟合动力学参数问题-1
c123.jpg

12楼2014-11-06 22:21:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allengjn

新虫 (正式写手)

引用回帖:
12楼: Originally posted by dingd at 2014-11-06 22:21:38
1stOpt代码:

Constant M=0.3986,Q=100;
Variable t,x,x',x'';
ODEFunction x''= (Q-C*x'-K*x)/M;
Data;
//时间        位移x        速度dx        加速度ddx
0        0        0        26663.2
0.000500005        0.000137226        0.382478        281.033
0.0010 ...

哦,非常感谢,还有最后一个小问题,将代码输入到1stopt之后怎么运行才是对它的函数拟合,我试了一下午总是不灵啊,您的1stopt是哪个版本的,是版本的问题吗?我输入的代码,之后点”执行“,就全都是0;点”快速公式拟合搜索“,就是拟合出来一个sin函数。这个怎么搞?
Parameter K,C;
Variable t, x;
ODEFunction x''= (100-C*x'-K-x)/0.39;
Data;
0        0
0.000500005        0.000137226
0.00100001        0.000360148
0.00150001        0.000640082
0.00200002        0.000965359
0.00250002        0.00130017
0.00300003        0.0016186
0.00350003        0.00191261
0.00400004        0.00220283
0.00450004        0.0024554
0.00500005        0.00266415
0.00550005        0.00288092
0.00600005        0.00314881
0.00650006        0.00344612
13楼2014-11-06 22:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allengjn

新虫 (正式写手)

送红花一朵
悬赏分给了楼上的那位了,不好意思就只能送您一朵花了,非常感谢,希望以后还能得到您的帮助
14楼2014-11-06 22:29:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
13楼: Originally posted by allengjn at 2014-11-06 22:27:15
哦,非常感谢,还有最后一个小问题,将代码输入到1stopt之后怎么运行才是对它的函数拟合,我试了一下午总是不灵啊,您的1stopt是哪个版本的,是版本的问题吗?我输入的代码,之后点”执行“,就全都是0;点”快速公 ...

如果用的是1.5破解版就别试了,没法运行的。
15楼2014-11-06 23:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allengjn

新虫 (正式写手)

内容已删除
16楼2014-11-06 23:04:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allengjn

新虫 (正式写手)

内容已删除
17楼2014-11-06 23:34:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dikeway

新虫 (小有名气)

引用回帖:
17楼: Originally posted by allengjn at 2014-11-06 23:34:01
您的意思是不是基本上网上能下到的,都不能拟合?您的是花钱买的吗?...

都告诉你了,低版本的不能拟合微分方程,所以你输入代码结果全是0。大神是高版本的,据悉相当之贵,1stopt6.0加密狗版1.5万+
18楼2014-11-07 09:04:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allengjn

新虫 (正式写手)

版主,我看了你回复的http://muchong.com/bbs/viewthread.php?tid=6425538&authorid=1122189这个帖子,我现在有些混乱了,我的二阶微分方程如何在KineticsEqs方程中写出啊,因为我的拟合数据有三组,位移、速度、加速度,这三组是相关的,如何编写这个程序,希望版主能够指导一下
19楼2014-11-07 10:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 allengjn 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 范进中举一文的中心思想 +5 炎黄贵胄 2026-08-22 6/300 2026-08-24 06:34 by sincosx
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见