| 查看: 2146 | 回复: 20 | ||||
[求助]
请问，有谁用matlab拟合过J-C材料本构方程吗，求帮助
|
||||
我做了不同温度下和不同速率下的拉伸试验，想拟合如图中参数B C n m ，用之前的方法拟合太差，文献里有用matlab拟合的，貌似还要编程，因为没有接触过matlab，说以请教大家  QQ截图20160324145330.png |
回复此楼» 本帖已获得的红花（最新10朵）
木虫小生» 本帖@通知
@拐点
» 猜你喜欢
香港科技大学全奖招收博士/博士后，光子、纳米材料、液晶显示方向
已经有15人回复
-
法国博士招聘
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有169人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有6人回复
-
论文接收，散金币了！
已经有149人回复
-
vasp算居里温度
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 67793.9
- 散金: 1938
- 红花: 442
- 沙发: 4
- 帖子: 4373
- 在线: 3284.6小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
|
Parameter m[0.001,0.01]; Function y=(425+b*x^c)*(1+d*ln(m/0.001)); Data; 0 432 0.013 442 0.014 452 0.031 527 0.049 582 0.058 600 0.08 635 0.102 661 0.116 672 0.13 682 以上代码正确。 此外，公式y=(425+b*x^c)*(1+d*ln(m/0.001)) 是过拟合的，式中，d*ln(m/0.001)完全可以合并为一个参数，存在过拟合的情况下，d和m会有无穷多组结果。 关于过拟合，参见该链接中第7条：http://muchong.com/bbs/viewthread.php?tid=9176191 |
17楼2016-04-17 18:22:03
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 67793.9
- 散金: 1938
- 红花: 442
- 沙发: 4
- 帖子: 4373
- 在线: 3284.6小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
拐点: 金币+20, ★★★★★最佳答案, 谢谢，虽然我回复的太晚了！抱歉，之前弄出来了，所以就没上木虫了。抱歉。 2016-06-03 14:24:37
拐点: 金币+20, ★★★★★最佳答案, 谢谢，虽然我回复的太晚了！抱歉，之前弄出来了，所以就没上木虫了。抱歉。 2016-06-03 14:24:37
|
抱歉，之前时间仓促，17楼代码不正确。为了实现变参数的拟合，代码修正如下： VarConstant m=[0.001,0.01]; Function y=(425+b*x^c)*(1+d*ln(m/0.001)); Data; 0 432 0.013 442 0.014 452 0.031 527 0.049 582 0.058 600 0.08 635 0.102 661 0.116 672 0.13 682 Data; 0 435 0.01 442 0.014 466 0.026 522 0.04 571 0.055 607 0.066 627 0.08 648 0.094 664 0.109 678 但因为原问题公式过拟合问题仍然存在，解释见17楼，所以计算结果仍不唯一： 相关系数(R): 0.97706015366352 相关系数之平方(R^2): 0.954646543876981 决定系数(DC): 0.954646543876981 卡方系数(Chi-Square): 7.19594508693612 F统计(F-Statistic): 378.882653246445 参数 最佳估算 ---------- ------------- b 1187.73192654974 c 0.667411209648257 d 0.000499331780954299 m 2.7014340335103E-24 参数 最佳估算 ---------- ------------- b 1187.73192701695 c 0.667411210301741 d 2.66058047772146 m 0.000991150872922775 |
18楼2016-04-18 10:08:46
9楼2016-03-26 22:31:50
| 文件在这里，前面那个作废 @月只蓝 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : 数据加说明.xlsx
2016-04-15 08:39:11, 11.74 K
15楼2016-04-15 08:39:40
2楼2016-03-24 15:04:09
3楼2016-03-24 15:38:47
4楼2016-03-24 15:39:20
5楼2016-03-24 15:57:12
6楼2016-03-24 16:29:20
7楼2016-03-25 15:37:18
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 67793.9
- 散金: 1938
- 红花: 442
- 沙发: 4
- 帖子: 4373
- 在线: 3284.6小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
8楼2016-03-25 15:57:03
你好，我看了，之前用origin拟合，，但是我好几组在不同条件下的数据，拟合到一个公式里，不知道怎么弄了，最近看了点matlab以前也没有接触过这方面，也是没有头绪不知道改怎么弄了，
10楼2016-04-14 10:38:19
10