版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dlionel

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 3424.6
散金: 109
红花: 4
帖子: 842
在线: 156.5小时
虫号: 761805
注册: 2009-05-02
专业: 固体力学

[求助] 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。
请问是多个约束条件命令不对吗？

回复此楼

» 猜你喜欢

求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有51人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

Good times, hard times, but never bad times.

1楼2015-01-06 10:31:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dlionel

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 3424.6
散金: 109
红花: 4
帖子: 842
在线: 156.5小时
虫号: 761805
注册: 2009-05-02
专业: 固体力学

引用回帖:

3楼: Originally posted by 月只蓝 at 2015-01-06 11:20:25
修正乱码：
Constant ro=2080;
Parameters c11,c12,c44;
Variable x,y;
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)*c33>2*c13^2 怎么没有在你的代码里给出啊？应该怎么加上去呢？

根据你的回复，我又修改了一下原问题：c12是不需要拟合的，拟合函数形式不变，
自变量：theta
应变量：V
拟合参数：c11, c13, c33, c44
已知常数：c12,ro
约束条件：c11>c12;
            c33>0;
            c44>0;
            (c11+c12)*c33>2*c13^2;

赞一下

回复此楼

Good times, hard times, but never bad times.

4楼2015-01-06 13:59:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dlionel

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 3424.6
散金: 109
红花: 4
帖子: 842
在线: 156.5小时
虫号: 761805
注册: 2009-05-02
专业: 固体力学

引用回帖:

4楼: Originally posted by dlionel at 2015-01-06 13:59:39
多谢版主，手册确实很有帮助。

请问约束 (c11+c12)*c33>2*c13^2 怎么没有在你的代码里给出啊？应该怎么加上去呢？

根据你的回复，我又修改了一下原问题：c12是不需要拟合的，拟合函数形式不变，
自变量： ...

且已知常数
ro = 2080
c12 = 850000

参数initial guess
c11 = 1384531.2
c13 = 8000000
c33 = 1081267.2
c44 = 191692.8

赞一下

回复此楼

Good times, hard times, but never bad times.

5楼2015-01-06 14:08:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 dlionel 的主题更新

返回列表