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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 使用MATLAB拟合动力学参数问题
192/2返回列表上一页12
查看: 1994  |  回复: 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
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭