Origin使用自定义函数拟合曲线函数
(2019年2月19日注:这篇文章原先发在自己github那边的博客,时间是2016年10月28日)
最近应该是六叔的物化理论作业要交了吧,很多人问我六叔的作业里面有两道题要怎么进行图像函数的拟合。综合起来的问题主要有两个:
1. 知道图像的准确拟合方程,但是不知道怎么拟合。(这个是本文的重点)
2. 不知道图像的准确拟合方程,也不知道怎么拟合,这个我可以稍微提供一个拟合的方向。
先从不知道图像的准确拟合方程开始吧,当你拿到一组数据时,你肯定是先看横纵坐标写的是什么,在19题的第二问里面,要求蒸发焓和热力学温度的关系,学过物化的都知道,肯定是克劳修斯-克拉贝龙方程才能出现这两个变量的关系。因此这里我们就要利用这个方程对我们的数据进行变形。$$ln(\frac{p}{p^{\Theta}}) = \frac{\Delta_{p}H^{m}}{R}(\frac{1}{373}-\frac{1}{T})$$
这就很明显了嘛,把我们得到的蒸气压除以标准状态下的气压,它的对数和$-1/T$成线性关系,接着就不用我说了吧。数学建模的其中一个思路是根据已有的经验公式进行改进和修正啊,因为建模也要符合基本法啊,不能随便就上个算符吧。
18题我留着后面说,先把19题的第一问的完整拟合过程在这里写一遍。
1. 根据已有的拟合形式将所需要的数据处理好,19题的公式是:$$lg(p) = A-\frac{B}{T-C}$$要的自变量是热力学温度T,因变量是$log_{10}(p),先整理好,这个在Excel里面是很容易的事情。
2. 打开Origin,我的是2015,不同版本操作大同小异,打开Analysis→ Fitting→ Nonlinear Curve Fit→ Open dialog。
3. 在Category里面选择User Defined,可见Origin的拟合函数也是类似于Matlab的M文件一样可以自己编写。
4. 当然是选择New一个函数咯,不是New一个对象。
5. 从这张图开始,注意看左边的Hints,这个比网上的所谓的教程不知道高到哪里去了。从Hints中可以看到,要输入函数名,后缀名是FDF,然后是函数模型,Explicit是指明确的,Function Type是表达式,也就是所谓的y=f(x)的形式,可以看到还能选Equlations,左边的Hints告诉我们,这种类型是处理只有一个自变量的。
6. 点击Next,这个界面是输入自变量,因变量,参数的名字,感天谢地我的数理统计是英文版的,这里直接阅读无压力,可以看到,多个参数时用逗号分隔。
7. 这一个界面是整个过程最关键的部分,自定义函数的拟合过程中最重要的就是初值的选取,这里我们先按照要求把函数的形式输入。注意到左边的Hints,Fixed是说如果初值选取的不好的时候也不能变参数的值,显然不符合我们的意思,所以我们这里都不勾选。Fixed用在已经确定了某些参数的范围的时候。接着是初值的尝试,看到了下面的Quick Check了吗,这里是要我们选取不同的ABC作为初值,使得我们输入一个T的时候,跑出来的结果接近我们测试出来的数据值。
8. 经过漫长的测试,我测试出了一组值,只要接近就可以,不需要完全吻合,然后点击Next。
9. 这个界面也是非常的重要,刚才是测试一个点,现在就是需要将一个点变成一组数值,使得我们的函数尽可能匹配多的测试值,我们点击红色框框中的按钮,打开代码界面。
10. 看到这个界面,一股浓浓的C++代码风格袭来……,确实长的很像C++,但是又混有Matlab的风格,根据Matlab的思路,很容易就想到这个Vector存储的就是Matlab的矩阵中的值,至于这些double类型的,明显就是数字啦,这里提示我们,Code to be executed to initialized parameters,代码将在初始化参数以后执行。我们要把我们的表达式输入进去,使得vector能存储尽可能多而接近的测试值,根据变量存储的名字,很容易就知道输入的是什么东西。
11. 点击左上方的Compile,也就是编译,果然是C++风格的代码……,看到编译成功,如果失败的话请仔细检查是不是少打了什么东西,不过以逗号作为每行结尾也是反人类……
12. 然后关闭代码窗口,看到这里多来一串这样的符号,然后点击Next,一路Next直到Finish。
13. Finish以后回到这个界面,然后点击红色方框中的按钮,这个是指Fit untiled Converged,也就是根据Chi-Square Test,直到满足条件以后才能函数收敛,点击一次就可以了。可以看到提示,Fit may be Converged。
14. 点击Fit,我们就能看到一条很漂亮的拟合曲线了,R-Square还是0.99997哦,相信这个拟合应该满足要求了,接下来小伙伴就知道怎么做了吧。
好了开始研究18题,18题里面原来是没有公式的,但是不少小伙伴拿到答案,答案里面是一种很奇怪的拟合方式,并且还有迭代什么的,其实本质上也不难,也是利用刚才的方法,只不过这里就多了好几步,其中一个是在第5步输入参数名时要把常量名字输入,这里比如说是sigma,然后在尝试初值的时候要将Constant界面的sigma的值输入为0.07275,才能开始尝试初值。
按照上述的步骤开始做自定义函数拟合,得到的界面如下。这里的Iteration Algorithm就是答案里所写的方法,不过这个Fit did not Converge是什么情况啊,我试了很多次都不行,包括答案给的初值以及自己试的初值,后来硬着头皮直接Fit,发现也能够Fit的出来,只是。
第一个是根据上述的Fit出来的结果,不觉得这个0.996瞬间low了很多了吗,而且它还是未收敛的,换句话说,肯定存在一个形式比它更好的式子满足要求。
我们回头看看式子,你看到了ln(1+x),你看到这个难道真的不会想到用泰勒展开吗,用一下不就啥问题都解决了吗,而且多项式拟合在Origin里是有快捷键的啊!
根据Analysis→ Fitting→ Polynormal Fit,随手试了2,3,4三个数,至少这3个数拟合出来的R-Square都比那个自定义函数的高,其中3又是最高的那一个,因此我们只要选3次多项式作为拟合值就好了啊,当然咯,能够用泰勒展开式的前提,是余项要趋于0,因此,泰勒展开式也不能随便乱用。
Origin使用自定义函数拟合曲线函数的更多相关文章
- 我的MYSQL学习心得(十) 自定义存储过程和函数
我的MYSQL学习心得(十) 自定义存储过程和函数 我的MYSQL学习心得(一) 简单语法 我的MYSQL学习心得(二) 数据类型宽度 我的MYSQL学习心得(三) 查看字段长度 我的MYSQL学习心 ...
- PHP 错误与异常 笔记与总结(8)自定义错误处理函数 set_error_handler()
通过 Set_error_handler() 函数设置用户自定义的错误处理函数. 步骤: ① 创建错误处理函数 ② 设置不同级别调用函数 ③ Set_error_handler() 函数制定接管错误处 ...
- 在String()构造器不存在的情况下自定义一个MyString()函数,实现如下内建String()方法和属性:
在String()构造器不存在的情况下自定义一个MyString()函数,实现如下内建String()方法和属性: var s = new MyString("hello"); s ...
- Qt自定义sleep延时函数(巧妙的使用时间差,但这样似乎CPU满格,而不是沉睡)
Qt不像VC++的win32/MFC编程那样,提供了现成的sleep函数可供调用.Qt把sleep函数封装在QThread类中.子线程可以调用sleep函数.但是如果用户想在主线程实现延时功能,该怎么 ...
- C程序中引用自定义的C函数模块
原文:C程序中引用自定义的C函数模块 我们知道,刚开始接触C语言编程,一般都是在一个.c或者.cpp(以下只说.c)的文件中编写代码,其中一定会有一个入口函数, 也就是main()函数,你可以将程序代 ...
- 自定义一个EL函数
自定义一个EL函数 一般就是一下几个步骤,顺便提供一个工作常用的 案例: 1.编写一个java类,并编写一个静态方法(必需是静态方法),如下所示: public class DateTag { pri ...
- 自定义JS Map 函数
// 自定义JS Map 函数 function Map() { var map = function (key, value) {//键值对 this.key = key; this.value = ...
- GRDB自定义的纯函数
GRDB自定义的纯函数 在GRDB中,用户可以自定义SQlite函数.这样,在SQL语句中,可以直接调用这些函数.但是在定义的时候,用户需要指定函数的pure属性,表示该函数是否为纯函数.纯函数是 ...
- sql server 系统常用函数:聚合函数 数学函数 字符串函数 日期和时间函数和自定义函数
一.系统函数 1.聚合函数 聚合函数常用于GROUP BY子句,在SQL Server 2008提供的所有聚合函数中,除了COUNT函数以外,聚合函数都会忽略空值AVG.COUNT.COUNT_BIG ...
随机推荐
- Linux安全审计-基础篇
安全审计这块我能想到的有两种方案可以解决,一种是在Linux中配置实现,一种是使用Python开发堡垒机实现,我先实现了第一种比较简单的:后面会开发堡垒机: 一.首先我们需要在/etc/profi ...
- Linux 常用命令:文本查看篇
前言 Linux常用命令中,除了cat还有很多其他用于文本查看的命令.本文将简单介绍一下这些文本查看的命令. 全文本显示--cat cat可能是常用的一个文本查看命令了,使用方法也很简单: cat f ...
- Android Button.getWidth()为0的问题
View在onCreate的时候,没有渲染组件,所以获取到的宽度和高度为0, 需要添加一个观察者,在layout渲染后再去取宽高.代码如下: private Button btn_icon; @Ove ...
- BCZM : 1.8
问题: 所有的员工均在1楼进电梯的时候,选择所要到达的楼层.然后计算出停靠的楼层i,当到达楼层i的时候,电梯停止.所有人走出电梯,步行到所在的楼层中.求所有人爬的楼层数目和的最小值. 解法一 ...
- Java中的常量池
JVM中有: Class文件常量池.运行时常量池.全局字符串常量池.基本类型包装类对象 常量池 Class文件常量池: class文件是一组以字节为单位的二进制数据流,在java代码的编译期间,编写的 ...
- CIE XYZ
了解CIE XYZ的来龙去脉,看维基之前,先读这两篇文章: https://medium.com/hipster-color-science/a-beginners-guide-to-colorime ...
- vue框架介绍
vue框架介绍 一.vue 概念 vue 是一种开发用户界面的渐进式开发框架.渐进式指的是:你可以将vue作为一部分嵌入到web应用中,带来丰富的交互体验 二.vue特点及常见开发中的高级功能 1.解 ...
- 修改sql server表字段的字符串
网站标题被注入黑链接,使用sql脚本update修改字段内的字符串截取UPDATE [qds0460132_db].[dbo].[Blood_News] SET [Blood_Name] = SU ...
- CSS表格属性
border-collapse:表格边框线合并,取值:collapse.
- MT Call来电流程分析????