应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 在用origin拟合的过程中，程序编写这一块实在不会，求指教！
111/2返回列表12下一页
查看: 2471  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cjj1001

新虫 (初入文坛)

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

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


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


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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-17 22:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

cjj1001

新虫 (初入文坛)
第二个公式中的参数d,e,f为固定值，分别等于8.314,303,6.37
赞一下(1人) 回复此楼
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是无法拟合的
赞一下(1人) 回复此楼
泉涸,鱼相与处于陆,相呴以湿,相濡以沫,不如相忘于江湖。
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这块内容我不会写
赞一下(1人) 回复此楼
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 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料类284调剂 +34 想换手机不想解� 2026-04-08 41/2050 2026-04-10 11:51 by wp06
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 求调剂 +3 电气300求调剂不 2026-04-08 3/150 2026-04-10 10:19 by LHGeng
[考研] 342电子信息专硕求调剂 +6 你让我怎么荔枝 2026-04-10 7/350 2026-04-10 10:18 by 314126402
[考研] 070300化学学硕311分求调剂 +12 梁富贵险中求 2026-04-04 14/700 2026-04-10 08:07 by kangsm
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 332，085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 280求调剂 +5 兮兮夜夜 2026-04-09 8/400 2026-04-09 11:15 by 兮兮夜夜
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 复试调剂 +14 呼呼？~+123456 2026-04-05 14/700 2026-04-06 22:50 by chenzhimin
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
信息提示
关闭
请填处理意见
关闭