| 查看: 2126 | 回复: 6 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
dlionel银虫 (正式写手)
|
[求助]
1stopt 非线性约束拟合 已有1人参与
|
|||
请教大家一个关于多非线性约束下拟合的问题， 自变量： theta 应变量：V 参数： c11, c12, c13, c33, c44 常数： ro 1stopt (v1.5) 的代码如下： Parameter c11=53.872,c12=10,c13=15,c33=47.424,c44=19.968; Constant ro=2080; Function y=sqrt((c44+c11*sin(x)^2+c33*cos(x)^2+sqrt(((c11-c44)^2*sin(x)^2+(c44-c33)^2*cos(x)^2)^2+(c13+c44)^2*sin(2*x)^2))/(2*ro)); c11>c12; c12>0; c44>0; (c11+c12)*c33>2*c13^2; Data; 0 22.8 0.382227106186758 22.2 0.911061869541040 29.2 1.22696646415201 25.0 1.54461638801498 22.8 1.57079632679490 25.9 拟合结果： 迭代数: 16 计算用时(时:分:秒:微秒): 00:00:00:63 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 23.7738301499779 残差平方和(SSE): 3391.17 相关系数(R): 0 相关系数之平方(R^2): 0 决定系数(DC): -95.2443593018306 卡方系数(Chi-Square): 130.843543616615 F统计(F-Statistic): 0 参数 最佳估算 ---------- ------------- c11 53.872 c12 10 c13 15 c33 47.424 ====== 结果输出 ===== No 实测值y 计算值y 1 22.8 1.0000000 2 22.2 1.0000000 3 29.2 1.0000000 4 25 1.0000000 5 22.8 1.0000000 6 25.9 1.0000000 参数拟合结果明显是错误的，而且没有参数c44的估算，计算值y都是1.00。 请问是多个约束条件命令不对吗？ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有177人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 非线性拟合
已经有16人回复
- 请教一个非线性方程非线性项转化为多项式的方法
已经有31人回复
- 关于有约束非线性函数的优化问题
已经有8人回复
- 怎么做matlab多元非线性拟合
已经有8人回复
- 不可线性化的非线性曲线方程拟合
已经有5人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 求助1stopt含约束非线性拟合
已经有5人回复
- ]紧急求助复介电常数的模型拟合-Matlab or 1stOpt
已经有9人回复
- 在进行多元非线性拟合的时候，用matlab拟合的值残差都大于0
已经有14人回复
- 非线性拟合得到的模型可靠性怎么验证？
已经有8人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 用1stopt或者origin非线性拟合磁导率
已经有5人回复
- 非线性数据拟合的数学原理！！！！！？？？？
已经有30人回复
- 求指点 langmuir线性拟合 与 非线性拟合 不一致
已经有12人回复
- 自定义非线性拟合中两个自变量，如何拟合？
已经有4人回复
dsp410327
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8816.6
- 帖子: 401
- 在线: 83.6小时
- 虫号: 1281931
- 注册: 2011-04-30
- 性别: GG
- 专业: 遥感机理与方法
7楼2018-09-24 16:19:32
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
¾ù·½²î(RMSE): 1.64771329728459 ²Ð²îƽ·½ºÍ(SSE): 16.2897546602908 Ïà¹ØϵÊý(R): 0.735295008053769 Ïà¹ØϵÊý֮ƽ·½(R^2): 0.540658748868792 ¾ö¶¨ÏµÊý(DC): 0.537682569595835 ¿¨·½ÏµÊý(Chi-Square): 0.329591175354763 Fͳ¼Æ(F-Statistic): 4.70812275220068 ²ÎÊý ×î¼Ñ¹ÀËã ---------- ------------- c11 1570.27237574273 c44 2.01763943470314E-9 c33 -1401.85923757196 c13 2126665.49354796 ====== ½á¹ûÊä³ö ===== No ʵ²âÖµy ¼ÆËãÖµy 1 22.8 21.7271513 2 22.2 24.5656304 3 29.2 27.1722123 4 25 25.7695871 5 22.8 24.3644782 6 25.9 24.3537906 建议你看看该软件的指导手册：http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank |
2楼2015-01-06 11:19:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★
dlionel: 金币+5, ★★★很有帮助 2015-01-06 13:52:09
dlionel: 金币+5, ★★★很有帮助 2015-01-06 13:52:09
|
修正乱码： 均方差(RMSE): 1.64771329728459 残差平方和(SSE): 16.2897546602908 相关系数(R): 0.735295008053769 相关系数之平方(R^2): 0.540658748868792 决定系数(DC): 0.537682569595835 卡方系数(Chi-Square): 0.329591175354763 F统计(F-Statistic): 4.70812275220068 参数 最佳估算 ---------- ------------- c11 1570.27237574273 c44 2.01763943470314E-9 c33 -1401.85923757196 c13 2126665.49354796 ====== 结果输出 ===== No 实测值y 计算值y 1 22.8 21.7271513 2 22.2 24.5656304 3 29.2 27.1722123 4 25 25.7695871 5 22.8 24.3644782 6 25.9 24.3537906 |
3楼2015-01-06 11:20:25
dlionel
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3424.6
- 散金: 109
- 红花: 4
- 帖子: 842
- 在线: 156.5小时
- 虫号: 761805
- 注册: 2009-05-02
- 专业: 固体力学
4楼2015-01-06 13:59:39
5