24小时热门版块排行榜    

查看: 1392  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

522142149

新虫 (正式写手)

[求助] 求助一个origin 8.0 自定义非线性拟合的具体操作

小弟最近看了一篇文献讲到用origin做非线性拟合,但不知道具体怎么操作,希望哪位大神现身相助,帮忙写出origin的具体操作,并给出Kg值,感激不尽。这个程序的算法是这样的,我需要拟合的参数是Kg,自变量是x,因变量是y, 我自己试过,但是不知道怎么定义a,b,c,d,H,step这些东西,下面是具体的算法流程,

a=0.515*Kg;
b=0.515+Kg+0.515*Kg*50+0.515*Kg*x-0.515*Kg*50;
c=1+0.515*50+Kg*x-(0.515+Kg)*50;
d=-50;       
for(H=50,step=1;abs(step)>1E-15;H=H-step){
step=(a*H*H*H+b*H*H+c*H+d)/(3*a*H*H+2*b*H+c)};
y=50*(0.004142+0.02*0.515*H)/(1+0.515*H)

x,y的参考数据为:

x              y
0                0.40692
8.16                0.46934
16.29        0.52351
24.39        0.57062
32.48        0.60668
40.52        0.64045
48.54        0.67391
56.54        0.70164
65.52        0.71935
72.46        0.73048
80.39        0.75476
88.28        0.73918
96.15        0.7396
104        0.73894
111.82        0.73845
119.62        0.73769

[ 来自小组 界面化学 ]
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
522142149: 金币+10, ★★★很有帮助, 谢谢,我已经用会用origin解决这个问题了,我拟合的值为37,跟文献差不多,现在我想用其他软件拟合下能发我一份1stopt的软件吗?在网上找的1.5pro版本运行起来有问题,另外顺便求matlab做非线性拟合的教程,视频文件都好,多谢啦 2013-08-06 11:39:00
关系式是否有误啊,如:
b=0.515+Kg+0.515*Kg*50+0.515*Kg*x-0.515*Kg*50;

里面既有“0.515*Kg*50”,也有“-0.515*Kg*50”,何不简化为:

b=0.515+Kg+0.515*Kg*x;

用1stOpt算了个结果,仅一个参数,效果不太好,参考下:

均方差(RMSE):0.18549434787848
残差平方和(SSE):0.550530449517803
相关系数(R): 0.954945715891251
相关系数之平方(R^2): 0.911921320299054
确定系数(DC): -2.05214149870227

参数                  最佳估算
--------------------        -------------
kg        0.0814361786948646

====== 结果输出 ======

文件: 数据文件-1
No        目标 y        计算 y
1        0.40692        0.858385919521254
2        0.46934        0.834627451860141
3        0.52351        0.812864072583636
4        0.57062        0.79312215539788
5        0.60668        0.775164214328853
6        0.64045        0.758842124357923
7        0.67391        0.743872372783508
8        0.70164        0.730071238766839
9        0.71935        0.71574555836118
10        0.73048        0.705415121301159
11        0.75476        0.694309183742899
12        0.73918        0.683919717474956
13        0.7396        0.674143733588132
14        0.73894        0.664918341827109
15        0.73845        0.65620062096851
16        0.73769        0.647932335500046
5楼2013-08-05 08:51:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
522142149: 金币+1, 有帮助, 谢谢,那个教程我已经看过,另外我还看了一些视频教程,那些都是一些很直接的拟合,但是我这个方程有点特殊,就是有个H,这个H 是随x变化的,我不知道怎么处理这个H,把它当成自变量? 2013-08-03 21:13:23
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-08-03 14:54:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

522142149

新虫 (正式写手)

教程上的拟合只有公式,这里还有一些C语言在里面,像是那个for 语句,所以很困扰这个到底怎么在origin上运行
3楼2013-08-03 21:16:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

522142149

新虫 (正式写手)

千万别沉了,急着用呢。。。。。哪位大神快来帮忙吧
4楼2013-08-04 05:16:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 06:40 by gy1nBQXYQJqL
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 06:21 by gy1nBQXYQJqL
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 05:13 by gy1nBQXYQJqL
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +4 gouxfjh 2026-08-28 4/200 2026-08-29 04:46 by alongwaytogo
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见