[求助] 使用MATLAB拟合动力学参数问题已有2人参与

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 数据.xlsx

2014-11-06 11:07:36, 18.31 K

1楼2014-11-06 11:07:57

allengjn

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3772.5
红花: 1
帖子: 487
在线: 28.3小时
虫号: 3367570
注册: 2014-08-15
专业: 机械动力学

引用回帖:

2楼: Originally posted by dingd at 2014-11-06 11:27:29
1stOpt试试：

均方差(RMSE):58.1680578466045
残差平方和(SSE):2030113.77218755
相关系数(R): 0.269739459742299
相关系数之平方(R^2): 0.0727593761420674
确定系数(DC): -1.78401536046027
F统计(F-Stat ...

数据是没有问题的，同时从仿真软件里提取出来的
1stOpt这个是matlab的计算程序吗？我上网查好像是一个软件？
您这里的c和k是怎么求出来的，希望您能详细的帮我解释一下，非常感谢了

3楼2014-11-06 12:26:40

allengjn

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3772.5
红花: 1
帖子: 487
在线: 28.3小时
虫号: 3367570
注册: 2014-08-15
专业: 机械动力学

引用回帖:

4楼: Originally posted by dikeway at 2014-11-06 20:00:58
1stopt跟matlab不同，是个软件，C和K是通过拟合寻参的...

请问，您会1stopt这个软件吗？楼上他是怎么拟合的，这个过程您能贴出来吗？

5楼2014-11-06 20:07:42

allengjn

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3772.5
红花: 1
帖子: 487
在线: 28.3小时
虫号: 3367570
注册: 2014-08-15
专业: 机械动力学

引用回帖:

6楼: Originally posted by dikeway at 2014-11-06 20:23:13
楼上通过高版本的1stopt拟合的，有相应代码。网上有低版本的，但没有拟合微分方程功能
代码大致如下
Parameter c,k;
Variable t,x;
ODEFunction .......
Data ..........

哦，好吧，我看能不能等到大神回来，问问他怎么得到的，看能不能把代码贴上来

7楼2014-11-06 20:55:02

allengjn

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3772.5
红花: 1
帖子: 487
在线: 28.3小时
虫号: 3367570
注册: 2014-08-15
专业: 机械动力学

引用回帖:

大神，求回复啊？代码能贴出来吗？我一整天都在搞这个东西，等回复等了一天了

8楼2014-11-06 21:08:59

allengjn

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3772.5
红花: 1
帖子: 487
在线: 28.3小时
虫号: 3367570
注册: 2014-08-15
专业: 机械动力学

引用回帖:

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

allengjn

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3772.5
红花: 1
帖子: 487
在线: 28.3小时
虫号: 3367570
注册: 2014-08-15
专业: 机械动力学

引用回帖:

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