| 查看: 1835 | 回复: 10 | ||
| 【悬赏金币】回答本帖问题，作者hitsz.civil将赠送您 100 个金币 | ||
hitsz.civil金虫 (正式写手)
|
[求助]
求助1stOpt，Bouc-wen模型识别参数，非常感谢，已有2人参与
|
|
|
各位大神，请帮我运算下下列表达式(可能有错)和数据，非常感谢， parameters a,b,c,d,e,n; variable t,x,y,f; odefunction diff(z,t)=a*y-b*y*(abs(z))^(n)-c*abs(y)*z*(abs(z))^(n-1); function f=d*x+e*z; data; 数据太多，请见附件，谢谢 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : 数据.xlsx
2019-06-12 11:14:41, 108.53 K
» 猜你喜欢
氧化硅薄膜玻璃 求购
已经有3人回复
-
vasp极化强度
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有96人回复
-
vasp计算的一些问题
已经有0人回复
-
国科温州研究院 项目聘用有意“二战”的优秀本科毕业生
已经有31人回复
-
SCAPS-1D学习交流
已经有0人回复
-
深圳大学跟哈工大联合博士
已经有5人回复
独孤神宇
版主 (知名作家)
- 应助: 488 (硕士)
- 贵宾: 0.008
- 金币: 29532.8
- 散金: 802
- 红花: 120
- 沙发: 1
- 帖子: 5337
- 在线: 850.2小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
你这两个方程没有什么耦合关系，可以分开拟合，先拟合第一个，然后拟合第二个。这样参数都不超过4个，试用版就可以运行了。 http://muchong.com/t-8050245-1 QuickReg = 10; Parameter a,b,c,n; Variable t,y,z; ODEFunction z'=a*y-b*y*(abs(z))^(n)-c*abs(y)*z*(abs(z))^(n-1); // f=d*x+e*z; Data; //t,y,z; |
2楼2019-06-12 16:47:38
hitsz.civil
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 334.9
- 散金: 295
- 帖子: 349
- 在线: 438.5小时
- 虫号: 1267111
- 注册: 2011-04-16
- 专业: 结构工程
3楼2019-06-12 17:05:02
hitsz.civil
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 334.9
- 散金: 295
- 帖子: 349
- 在线: 438.5小时
- 虫号: 1267111
- 注册: 2011-04-16
- 专业: 结构工程
4楼2019-06-12 17:11:43
独孤神宇
版主 (知名作家)
- 应助: 488 (硕士)
- 贵宾: 0.008
- 金币: 29532.8
- 散金: 802
- 红花: 120
- 沙发: 1
- 帖子: 5337
- 在线: 850.2小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
5楼2019-06-12 18:20:55
hitsz.civil
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 334.9
- 散金: 295
- 帖子: 349
- 在线: 438.5小时
- 虫号: 1267111
- 注册: 2011-04-16
- 专业: 结构工程
6楼2019-06-13 09:43:19
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2019-06-13 14:34:01
感谢参与，应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2019-06-13 14:34:01
|
这种微分方程与代数方程混合的拟合问题比较复杂，加之数据多，计算量会很大的。均匀取了十分之一的数据，用4阶龙格库塔算法替代缺省的RKF45算法（RKF45精度高但更费时间）试了下，结果参考如下。要想获得更好结果，其他人无法帮，只能自己慢慢算。 ODE Algorithms: Fourth Order Runge-Kutta Method - RK4 Algorithms: Universal Global Optimization(UGO1) Stop Reason of Computation: Convergence tolerance reached Elapsed Time (Hr:Min:Sec:Msec): 00:04:53:648 Root of Mean Square Error (RMSE): 1.13858352091496 Sum of Squared Residual: 217.790568928649 Correlation Coef. (R): 0.993685383888104 R-Square: 0.98741064215285 Adjusted R-Square: 0.987258043875914 Determination Coef. (DC): 0.987399390830626 F-Statistic: 2089.35690605282 Parameter Best Estimate -------------------- ------------- a 4.43611964934065 b 0.171990099795075 c -0.192529680045407 d 0.136012768179973 e -0.484114290416918 n 0.680155516905957 z Initial Value 27.5072080867539  c177.jpg |
» 本帖已获得的红花（最新10朵）
hitsz.civil7楼2019-06-13 10:23:41
hitsz.civil
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 334.9
- 散金: 295
- 帖子: 349
- 在线: 438.5小时
- 虫号: 1267111
- 注册: 2011-04-16
- 专业: 结构工程
8楼2019-06-13 11:57:14
hitsz.civil
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 334.9
- 散金: 295
- 帖子: 349
- 在线: 438.5小时
- 虫号: 1267111
- 注册: 2011-04-16
- 专业: 结构工程
9楼2019-06-14 06:45:19
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
10楼2019-06-14 10:29:26
