| 查看: 24203 | 回复: 72 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
drk123新虫 (著名写手)
|
[交流]
利用Excel数据前处理，结合1stopt数据自定义拟合初值，美美origin数据作图分析！已有49人参与
|
||||
前言：本人也菜鸟，接粗数据处理模拟之后，自学了小木虫里面的软件教程，看到有些人的数据不会处理，想必自己也会遇到，所以在大家求助帖子里面都说把数据弄上来，让我来处理或者拟合，所以现在弄这个帖子来一个总说，好下一次谁不会数据拟合处理作图分析的，可以引用本帖进行学习，本帖也不是很高水平，希望高手们多多包涵，觉得好的虫友多多献上几朵小红花！ 因为Excel进行数据倒置，数据的一些加减乘除等运算成新数据，数值转换批量处理等比较快，而且可以直接在原始数据中处理，所以我觉得Excel的数据处理是非常必要学会的，然后1stopt拟合自定义数据基本上能够获得自己需要的参数值，但是软件作图本身还是有点欠缺，但是参数值可以作为origin拟合的初值，origin本身作图高手，但是拟合自定义函数时候没有初始值就显得没那么好用，所以要结合1stopt的自定义拟合的参数值作为初始值，那样就可以完美结合。 首先，当然你要有Excel软件，1stopt软件，和origin软件，我把教程例子的所有数据上传到附件，大家可以下载进行学习哈，顺便弄个origin的自定义拟合的视频教程在附件有兴趣的可以下载看看。可以参考我这帖子下载：http://muchong.com/bbs/viewthread.php?tid=9900656&_from=share&u=2470864&h=38a3eda0 接着当然实例解说一下一个友友求助过的帖子： 图片中的公式拟合我的数据，其中t是自变量，Y是因变量；A,B,C,a,b,c是参数，但是有些约束条件， 0<A,B,C<1; 0.1<a<1, 0.01<b<1, 0<c<0.01  最好能用origin或者excel进行拟合。 虫友给出的数据如下，大家需要可以见附件下载：  好接下来我们看Excel的处理，可以看到数据本身第一行是x,第二三四五是y1y2y3y4数据，我们需要将他们倒置过来 全选复制  在新的数据表格中右键--选择性粘贴  选择倒置，从图里面我们可以看到，其实选择性粘贴里面有很多功能，大家别忽略了这些作用大大的功能  好了数据顺利变成我们的XYYYY的形式了：  在其他的数据转换中，记得别忘了表格右下角的鼠标变成黑色实心的十字光标的批量处理的大用处哈，好，数据的前处理结束。 接下来我们看看1stopt的数据自定义拟合初值，右键软件以管理员身份打开，  根据条件和函数我们可以输入成如图所示，并将xy1的数据复制进去。  我想说明一下就是这些代码的书写，这个parameter是可以自己定义的参数，参数可以给定范围（形式如图所示），也可以不给定范围。然后函数方面除本例外，大家可以输入任意自己需要的函数，但是一定要对应variable的变量；数据的输入可以参考图中所示，大家拟合其他自己需要的自定义公式都可以参考这图，是一个万能模板，有什么错误的书写，按提示修修改改就行啦。我们在这个软件上不需要用到那么高深，刚刚号够用就行了。 然后点击上面的三角符号开始拟合（旁边的闪电按钮是自动全函数拟合，如果不知道自己的函数拟合类型，可以用自动拟合按钮，参考软件给的拟合函数和参数值），接着软件跳到拟合结果界面  可以看到软件给出来了参考值和相关系数和决定系数，但是拟合是一种迭代收敛的过程，随着迭代的数量和次数不同，收敛程度不同，参数值就会发生变化，有时候迭代出来也有可能超出了我们所定义条件范围，我们可以按F9进行重新迭代，可以带到我们所需要的参数范围值。这里不能着急，其实得到上面的结果基本上就能够满足虫友所提出的要求，但是我们还是需要进一步用origin作图美化。同理我们将Y2Y3Y4的数据按同样的方法进行，得出我们需要的ABCa1b1c1参数值。 根据条件和函数我们可以输入成如图所示，并将xy2的数据复制进去，并进行拟合得出结果。   根据条件和函数我们可以输入成如图所示，并将xy3的数据复制进去，并进行拟合得出结果。   根据条件和函数我们可以输入成如图所示，并将xy4的数据复制进去，并进行拟合得出结果。   好，到这里，我们进入origin作图并拟合。保留刚刚1stopt拟合的ABCa1b1c1参数值作为origin的参数值。 复制所有数据粘贴到origin中，然后所有数据做散点线图看看   在这里我也要强调一下输入的设置，如上图，我们有时候要注意论文期刊对图表的分辨度的要求，关于其他的图片如何调节分辨率可以使用Photoshop，可以参考我写的这帖子http://muchong.com/bbs/viewthread.php?tid=9925119&_from=share&u=2470864&h=e3f9faaa 好啦，我们看到了整体的趋势了，那么我们开始分别作图拟合，一次性拟合的话增加难度，建议一组一组来，熟练之后就会很快的。 选择xy1数据做散点图，然后我们先设置函数，点击origin的函数拟合管理器如下图所示  然后按照下图输入需要的条件，最后第六步中需要进行验证所输入的函数是否能让origin识别   好接下来我们需要完整输入我们的条件，不过这里也可以点击验证我们的条件   设置好了以后点击save,然后OK  进行拟合，如下图所示  接下来选择我们刚刚设置好的newcategory  点击parameters界面，然后将1stopt的参数值作为origin的初始值输入，点击一次迭代，再点击fit。   需要注意的是：点击一次迭代或者fit以后，要看看其他参数值是否满足我们的参数条件，如果不满足，重新拟合，我们需要固定某些origin必然超出范围值的参数，比如本例中，a1和c1值是会超出我们限定的范围，那么我们在输入1stopt的初始参数值的时候就点击fixed的复选框让它将值固定不变，拟合迭代其他参数，这里基本需要点击一次迭代按钮就足够了，因为从1stopt的参数数据已经迭代收敛过了。 同理我们可以用同样的方法得出x-y2,x-y3,x-y4的拟合数据图    好，基本弄完了，看看整体效果，对比可以知道数据Y4的拟合效果比其他的数据要好。  哈哈，最后推荐一下我其他两篇帖子，希望都对你们有用，记得给几朵小红花，哈哈 与origin有关的 简单利用PPT和origin放映数据曲线作图 http://muchong.com/bbs/viewthread.php?tid=9772838&_from=share&u=2470864&h=c64d4512 论文参考文献格式导入技巧--中国知网、百度学术、谷歌学术、web of science、endnote http://muchong.com/bbs/viewthread.php?tid=9579680&_from=share&u=2470864&h=4daaed44 再推荐一发自写的草根教程，别见外 红外图谱(XRD图谱)从数据到origin作图(叠堆光滑标峰)简单教程 http://muchong.com/bbs/viewthread.php?tid=9931830&_from=share&u=2470864&h=cb7abff6   彩虹色线性填充紫外可见光谱图http://muchong.com/bbs/viewthread.php?tid=9990955&_from=share&u=2470864&h=ba2c8ef4 [ Last edited by drk123 on 2016-2-23 at 18:07 ] 彩虹色线性填充紫外可见光谱图 http://muchong.com/bbs/viewthread.php?tid=9990955&_from=share&u=2470864&h=ba2c8ef4 [ Last edited by drk123 on 2016-2-23 at 18:08 ] 荧光三维作图XYZ三维等高线图origin作图 http://muchong.com/bbs/viewthread.php?tid=10363741&_from=share&u=2470864&h=f39d7356 [ Last edited by drk123 on 2016-5-31 at 11:39 ] 是啦，如果大家想作带误差棒拟合效果图，下载附件中的教程吧。  [ Last edited by drk123 on 2016-6-2 at 19:13 ]origin线性拟合及添加置信区间预测区间范围图 http://muchong.com/bbs/viewthread.php?tid=10467152&_from=share&u=2470864&h=0d701fea [ Last edited by drk123 on 2016-7-7 at 22:33 ] 哈哈，楼主分享科研的公众号平台来了，微信公众号：olddrivergo。欢迎关注。  [ Last edited by drk123 on 2016-8-19 at 18:08 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 数据前处理.xls
- 附件 2 : origin工程文件.opj
- 附件 3 : 34.拟合(一)_线性与多项式拟合.wmv
- 附件 4 : 35.拟合(二)_拟合报表分析.wmv
- 附件 5 : 36.拟合(三)_非线性拟合.wmv
- 附件 6 : 37.拟合(四)_自定义函数拟合.wmv
- 附件 7 : 38.拟合(五)_峰值拟合,指数拟合与反曲拟合.wmv
- 附件 8 : 39.拟合(六)_比对.wmv
- 附件 9 : Y1数据1stopt拟合.mff
- 附件 10 : Y2数据1stopt拟合.mff
- 附件 11 : Y3数据1stopt拟合.mff
- 附件 12 : Y4数据1stopt拟合.mff
- 附件 13 : 原始数据.xlsx
- 附件 14 : 在origin中添加误差棒的方法与误差棒的计算方法.txt
2016-01-29 12:34:29, 10.5 K
2016-01-29 12:34:39, 1.39 M
2016-01-29 12:38:48, 39.9 M
2016-01-29 12:41:31, 46.04 M
2016-01-29 12:44:20, 48.18 M
2016-01-29 12:46:36, 38.83 M
2016-01-29 12:48:36, 33.88 M
2016-01-29 12:51:09, 43.87 M
2016-01-29 13:28:20, 282 bytes
2016-01-29 13:28:20, 282 bytes
2016-01-29 13:28:20, 282 bytes
2016-01-29 13:28:21, 282 bytes
2016-01-29 13:31:14, 6.75 K
2016-06-02 19:13:18, 378 bytes
» 本帖已获得的红花（最新10朵）
optimas91» 猜你喜欢
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
-
RSC ADV状态问题
已经有4人回复
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
1stopt代码可以一次列出计算： Parameters A1[0,1],B1[0,1],C1[0,1],a11[0.1,1],b11[0.01,1],c11[0,0.01], A2[0,1],B2[0,1],C2[0,1],a22[0.1,1],b22[0.01,1],c22[0,0.01], A3[0,1],B3[0,1],C3[0,1],a33[0.1,1],b33[0.01,1],c33[0,0.01], A4[0,1],B4[0,1],C4[0,1],a44[0.1,1],b44[0.01,1],c44[0,0.01]; Variable t,y1,y2,y3,y4; Function Y1=A1*exp(-a11*t)+B1*exp(-b11*t)+C1*exp(-c11*t); Y2=A2*exp(-a22*t)+B2*exp(-b22*t)+C2*exp(-c22*t); Y3=A3*exp(-a33*t)+B3*exp(-b33*t)+C3*exp(-c33*t); Y4=A4*exp(-a44*t)+B4*exp(-b44*t)+C4*exp(-c44*t); data; //t,y1,y2,y3,y4 1 0.642766547 0.658512643 0.670991891 0.800273001 2 0.464945956 0.5518936 0.571903848 0.738213852 4 0.428643115 0.513411445 0.550921127 0.712506579 6 0.426348447 0.508989708 0.548478241 0.70907828 12 0.41089034 0.45090639 0.53753663 0.705228105 24 0.395117563 0.431672713 0.535204853 0.703349159 48 0.383968102 0.422623233 0.53369585 0.700954519 96 0.381161451 0.420238554 0.531561544 0.699314008 120 0.379980721 0.417965224 0.530789817 0.698852127 240 0.378106999 0.416326383 0.530141909 0.698309918 |
11楼2016-02-01 09:42:19
xiaoshm
至尊木虫 (文坛精英)
Eugénie Grandet
- ACI: 1
- 应助: 44 (小学生)
- 贵宾: 0.355
- 金币: 45876.3
- 红花: 173
- 沙发: 1
- 帖子: 12292
- 在线: 2760.5小时
- 虫号: 654237
- 注册: 2008-11-14
2楼2016-01-30 17:04:11
txtzhimeng
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2.5
- 帖子: 14
- 在线: 10.6小时
- 虫号: 3222437
- 注册: 2014-05-21
- 专业: 气象观测原理、方法及数据
3楼2016-01-30 22:36:16
txtzhimeng
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2.5
- 帖子: 14
- 在线: 10.6小时
- 虫号: 3222437
- 注册: 2014-05-21
- 专业: 气象观测原理、方法及数据
4楼2016-01-30 22:48:07