24小时热门版块排行榜    

查看: 949  |  回复: 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
回复此楼

» 猜你喜欢

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

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

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by lipenggg at 2016-11-14 09:15:47
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
QQ图片20161114091805.png


多条曲线进行拟合,有两个参数-2
QQ图片20161114091831.png

一团和气
3楼2016-11-14 09:19:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by lipenggg at 2016-11-14 09:19:39
图形

QQ图片20161114091805.png

QQ图片20161114091831.png
...

图形
多条曲线进行拟合,有两个参数-3
QQ图片20161114092004.png

一团和气
4楼2016-11-14 09:20:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] FileCode能看出啥? +6 要乐观耀哥 2026-08-10 12/600 2026-08-10 18:37 by wmfsnow
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] filecode +7 documentary 2026-08-10 7/350 2026-08-10 18:32 by GZUZM
[基金申请] 为什么网上很多人说本周 12号出结果 +3 瞬息宇宙 2026-08-10 3/150 2026-08-10 18:27 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见