应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 使用MATLAB拟合动力学参数问题
北京石油化工学院2026年研究生招生接收调剂公告
192/2返回列表上一页12
查看: 2121  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +15 GENJIOW 2026-04-07 16/800 2026-04-08 06:23 by lijunpoly
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 复试调剂 +14 呼呼？~+123456 2026-04-05 14/700 2026-04-06 22:50 by chenzhimin
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 22408 总分320，一篇论文二作，两个国三，求调剂 +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 326分求调剂 +3 于是乎呢 2026-04-01 5/250 2026-04-03 14:23 by 于是乎呢
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 348求调剂 +6 吴彦祖24k 2026-04-02 6/300 2026-04-02 14:07 by 给你你注意休息
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 324求调剂 +5 想上学求调 2026-04-01 6/300 2026-04-02 10:16 by sanrepian
信息提示
关闭
请填处理意见
关闭