| 查看: 744 | 回复: 3 | ||||
洋洋sunnyy铁虫 (初入文坛)
|
[求助]
多条曲线进行拟合，有两个参数 已有1人参与
|
|
http://muchong.com/bbs/viewthread.php?tid=6793259&fpage=1&target=blank 这是有人之前发的帖子，我的麻烦和这个差不多，希望大神能以这个为例用1stopt告知怎么编写的代码？ 1、横坐标S表示比例，分别为1/4.8,1/3.2,1/2.4，为作图方便所以乘了4.8。 2、各条曲线表示模型比例与弯矩的关系，五条曲线分别为不同转交下的模型比例-弯矩示意图。 3、因此图中有两个未知参数：模型比例与转角。 现在需要拟合一条曲线，使图中五条曲线与其都比较接近且方差最小，因此所拟合的曲线方程中会有两个参数（模型比例和转角） 达到的目的就是根据上述实验曲线中的点拟合出一个方程，当别人给你一个模型的比例跟转角时，直接带入就可以计算出弯矩的大小。 请各位大神救救我。 横坐标 0.02rad 0.04rad 0.06rad 0.08rad 0.1rad 1 0.2466 0.405 0.537 0.77 1.29 1.5 0.417 0.6827 1.0 1.412 1.64 2 0.517 1.07 1.62 2.107 2.388  1.jpg |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有50人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 靠方程来描述复杂的真实世界，科学家是不是太天真了？
已经有12人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 多条曲线进行拟合，有两个参数
已经有9人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- 关于房室模型、权重的选择
已经有0人回复
- 怎样写好论文---一个大学教授、审稿专家的写作经验
已经有43人回复
- 怎么写好论文——一位大学教授、审稿专家写的经验文章【转载】
已经有1人回复
- 【sobereva个人文集】乱弹DFT-D（从sobereva的博客中转载过来的）
已经有27人回复
- 400金币急求一篇文献翻译中译英
已经有9人回复
- 生物学软件介绍及其下载地址【学生物的会很有用】
已经有24人回复
- 一个大学教授、审稿专家的写作经验
已经有33人回复
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
洋洋sunnyy: 金币+5, ★★★★★最佳答案, 谢谢 2016-11-19 13:21:02
感谢参与，应助指数 +1
洋洋sunnyy: 金币+5, ★★★★★最佳答案, 谢谢 2016-11-19 13:21:02
|
Title ""; Variable x,y,z; Function //z=(p1*x+p2)*y+p3*x^2+p4*x+p5; //z=(p1+p2*x+p3*y+p4*y^2+p5*y^3)/(1+p6*x+p7*x^2+p8*x^3+p9*y+p10*y^2); z=(a+b*ln(x)+c*ln(y)+d*(ln(y))^2+e*(ln(y))^3)/(1+f*ln(x)+g*(ln(x))^2+h*ln(y)+i*(ln(y))^2); data; 1.0 0.02 0.2466 1.5 0.02 0.4170 2.0 0.02 0.5170 1.0 0.04 0.4050 1.5 0.04 0.6827 2.0 0.04 1.0700 1.0 0.06 0.5370 1.5 0.06 1.0000 2.0 0.06 1.6200 1.0 0.08 0.7700 1.5 0.08 1.4120 2.0 0.08 2.1070 1.0 0.10 1.2900 1.5 0.10 1.6400 2.0 0.10 2.3880 迭代数: 82 计算用时(时:分:秒:微秒): 00:00:13:718 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 0.0264649583888136 残差平方和(SSE): 0.0105059103378245 相关系数(R): 0.999122427680878 相关系数之平方(R^2): 0.99824562549493 决定系数(DC): 0.998245065023016 卡方系数(Chi-Square): 0.00663031725872216 F统计(F-Statistic): 7397.04840325372 参数 最佳估算 ---------- ------------- a 2.52863283294756 b 0.195280522431928 c 2.41501095860775 d 0.759609301305623 e 0.0768577202998456 f 0.156561812928556 g -0.151373317264167 h 0.796925855227041 i 0.165879175786186 ====== 结果输出 ===== No 实测值z 计算值z 1 0.2466 0.2485435 2 0.417 0.3999500 3 0.517 0.5253854 4 0.405 0.4048510 5 0.6827 0.7357092 6 1.07 1.0433878 7 0.537 0.4964134 8 1 1.0445746 9 1.62 1.5985498 10 0.77 0.8019357 11 1.412 1.3763012 12 2.107 2.1159353 13 1.29 1.2802362 14 1.64 1.6385647 15 2.388 2.3955566 ====== 输出结果 ===== |
2楼2016-11-14 09:15:47
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
3楼2016-11-14 09:19:39
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
4楼2016-11-14 09:20:58