24小时热门版块排行榜    

查看: 2416  |  回复: 10

cjj1001

新虫 (初入文坛)

[求助] 在用origin拟合的过程中,程序编写这一块实在不会,求指教! 已有4人参与

文献的原始公式是第一个,然后我改写成第二个形式的进行自定义函数拟合,第三个图为我的原始数据
在用origin拟合过程中公式程序编写实在不会了,看小木虫上好多帖子都说用1stOpt比origin更简单,论坛里哪位大侠可以教我拟合一下,两种方法都行,因为后续还有好多组数据需要拟合,所以必须学会了,小女子先谢过各位!
在用origin拟合的过程中,程序编写这一块实在不会,求指教!
QQ图片20131217215153.jpg


在用origin拟合的过程中,程序编写这一块实在不会,求指教!-1
QQ图片20131217215220.jpg


在用origin拟合的过程中,程序编写这一块实在不会,求指教!-2
QQ图片20131217215913.jpg
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

cjj1001

新虫 (初入文坛)

第二个公式中的参数d,e,f为固定值,分别等于8.314,303,6.37
2楼2013-12-17 22:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

polypro

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2013-12-19 07:31:50
d,e,f为常数,定义上就行,a,b,c为参数,第二个公式在origin写成
a*exp(-1.0*(d*e*log(f/x)/b)^c);
如上说,x=0是无法拟合的
泉涸,鱼相与处于陆,相呴以湿,相濡以沫,不如相忘于江湖。
6楼2013-12-18 11:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjj1001

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by polypro at 2013-12-18 11:01:05
d,e,f为常数,定义上就行,a,b,c为参数,第二个公式在origin写成
a*exp(-1.0*(d*e*log(f/x)/b)^c);
如上说,x=0是无法拟合的

我已经把零点去掉了,origin公式输入这一步我会的,就是parameter initialization这块内容我不会写
7楼2013-12-18 14:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2013-12-19 07:29:52
妹纸,首先X列有一个零蛋,你放分母上啥意思呢?其次,方程的非线性程度,你不给出较接近的initial guess,ORIGIN也算不准的。ORIGIN做非线性拟合用的是Levenberg-Marquardt算法,这类优化算法都需要给出合适的初值。给出了合适的初值,应该能拟合出来。但本猪不保证能求出全局最优解(海塞矩阵不一定正定或半正定)。这是数学模型的问题,不能归咎于拟合软件。要么你可以尝试Monte-Carlo型算法,就是比较慢。
3楼2013-12-17 23:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2013-12-19 07:30:02
内容已删除
4楼2013-12-18 09:11:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjj1001

新虫 (初入文坛)

我把数据和公式的文本上传了,刚接触拟合这块内容,有好多错误的地方,请大家多批评指正,同时也希望能给我一些帮助,谢谢大家

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 数据.xls
  • 2013-12-18 09:46:10, 18.5 K
  • 附件 2 : 公式.docx
  • 2013-12-18 09:46:50, 12.29 K
5楼2013-12-18 09:47:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


fegg7502: 金币+1, 应助指数+1, 3ks 2013-12-19 07:31:59
1stOpt求解:
CODE:
Constant [d,e,f]=[8.314,303,6.37];
Function  y=a*exp(-1.0*(d*e*ln(f/x)/b)^c);
Data;
x        y
0.1718        243
0.4707        297
0.8448        330
1.2799        374
1.5762        386

均方差(RMSE):3.40899033271159
残差平方和(SSE):58.1060754426055
相关系数(R): 0.997869302680017
相关系数之平方(R^2): 0.995743145231104
确定系数(DC): 0.995743144656219
卡方系数(Chi-Square): 0.0823308572481187
F统计(F-Statistic): 234.91519425508

参数        最佳估算
----------        -------------
a        686.017513775682
b        8594.0482022177
c        0.631733566013568
在用origin拟合的过程中,程序编写这一块实在不会,求指教!-3
f3.jpg

8楼2013-12-18 17:30:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjj1001

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by dingd at 2013-12-18 17:30:11
1stOpt求解:

Constant =;
Function  y=a*exp(-1.0*(d*e*ln(f/x)/b)^c);
Data;
x        y
0.1718        243
0.4707        297
0.8448        330
1.2799        374
1.5762        386

均方差(RMSE):3.40899033271159
残差平方和(SSE):58.10 ...

太感谢了,这个软件能发我一份吗,我想自己学艺下,邮箱jiejingcao@163.com
9楼2013-12-18 17:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiw123

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 3ks 2013-12-19 07:32:07
1stopt拟合一条线很好,但是一系列的多条线就不行了

[ 发自小木虫客户端 ]
10楼2013-12-18 18:05:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cjj1001 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 一志愿武理材料工程348求调剂 +3  ̄^ ̄゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
信息提示
请填处理意见