| 查看: 1315 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
wangdeqiang金虫 (小有名气)
|
[求助]
1stopt拟合中出现的问题 已有1人参与
|
|
|
Title "DIC-K"; Parameters kI,kII,b; constant n=3,g=8.1*10^4,u=0.26,k=2.17; Variable x,y,z; Function z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; Data; 376.78243 -112.14027 9.73371 374.92666 -111.4325 9.61182 373.12866 -110.7178 9.49729 371.38928 -109.99629 9.39441 369.70935 -109.26811 9.3051 368.08966 -108.5334 9.22794 366.53104 -107.79231 9.15739 365.03425 -107.04503 9.08779 363.60006 -106.29172 9.0149 本人想通过1stopt求解kI，kII，但是运行时系统提示参数kI数值范围定义有误! 请问该如何写程序？ |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有170人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有26人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有34人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有6人回复
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 1stopt拟合出参数为0，为什么？
已经有6人回复
- 求助--用1stopt拟合蓝宝石衬底的GaN薄膜红外反射光谱
已经有6人回复
- 求助：以下非线性拟合方程可以用1stopt求解么？用什么方法
已经有16人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 1stopt求最小值
已经有37人回复
- 在用origin拟合的过程中，程序编写这一块实在不会，求指教！
已经有10人回复
- origin中如何自定义Birch–Murnaghan三阶拟合
已经有40人回复
- 生产中出了点问题 求解
已经有11人回复
- 1stopt中使用sharedmodel拟合曲线出现问题，求指示！！
已经有6人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 运算过程有复数出现，求高版本1stOpt
已经有6人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
wangdeqiang
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1025.8
- 散金: 10
- 帖子: 67
- 在线: 25.5小时
- 虫号: 1931929
- 注册: 2012-08-09
- 专业: 金属材料的力学行为
3楼2014-09-10 20:16:08
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wangdeqiang: 金币+20, ★★★很有帮助 2014-09-10 20:11:54
感谢参与，应助指数 +1
wangdeqiang: 金币+20, ★★★很有帮助 2014-09-10 20:11:54
|
公式： z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum((kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; 改为： z=sum(i=1:n)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum(i=1:n)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; 结果： 均方差(RMSE):0.0334465479380339 残差平方和(SSE):0.0100680441207408 相关系数(R): 0.989298144852019 相关系数之平方(R^2): 0.978710819407646 确定系数(DC): 0.978710819407646 卡方系数(Chi-Square): 0.000543327950684942 F统计(F-Statistic): 140.916649515266 参数 最佳估算 ---------- ------------- ki 0.245068724509857 kii -5.08642699323206 b 9.4238457540924 |
2楼2014-09-10 20:00:19
wangdeqiang
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1025.8
- 散金: 10
- 帖子: 67
- 在线: 25.5小时
- 虫号: 1931929
- 注册: 2012-08-09
- 专业: 金属材料的力学行为
4楼2014-09-10 21:14:35
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
公式写成如下： Function z=sum(i=1:n, kI)(kI/2/g*x^(i/2)*(k*sin(i*y/2)+i/2*sin((i/2-2)*y)-(i/2+(-1)^i)*sin(i*y/2))) -sum(i=1:n, KII)(kII/2/g*x^(i/2)*(-k*cos(i*y/2)+i/2*cos((i/2-2)*y)+(i/2-(-1)^n)*cos(i*y/2)))+b; 均方差(RMSE):0.000319000884196416 残差平方和(SSE):9.15854077062857E-7 相关系数(R): 0.999999031699104 相关系数之平方(R^2): 0.999998063399146 确定系数(DC): 0.999998063399146 卡方系数(Chi-Square): 4.91879921856916E-8 F统计(F-Statistic): 0 参数 最佳估算 ---------- ------------- b 9.37438411982016 ki1 240.991105375969 ki2 45.2897205090565 ki3 -0.00759010072910402 kii1 1285.30644313843 kii2 -40.4189608273874 kii3 -4.56724263105769 |
5楼2014-09-10 21:53:48
50