股票直线拟合怎么做?

阿伦尼斯(Arrhenius Mode)应该是所有经验加速模型中被最常采用的了,大部分资料中给出的加速因子计算式一般为:EXP[Ea/K*(1/T0-1/T1)],式中有个未知量E......

股票直线拟合怎么做

接下来具体说说

顶点式线性平均内插拟合算法

最近用Matlab做了各种不同模型的三维数据的仿真实验 ,不少需要对三维数据进行拟合,matlab 自带一个拟合工具箱(cftool),确实强大,使用工具箱做数据拟合省不少事,对数据集足够的原始模型而言,cftool能很准确的拟合出我们想要的结果。当然我们不讲Matlab ,本篇通篇跟Matlab毫无关系,~~话讲多了,下面开始。


一、数据原型

首先导入一张完整的数据表格,我们可以看到这是一个y=9 x=15 ,大小是9x15的矩阵,矩阵每个元素对应一个数值z,第二张图是我们在靠前张图里面的数据对应位置挑出来的数据。大小是3x5的有效矩阵。

这样我们就在原始数据虚拟出了一个残缺度高达90%左右的数据,那我们后面的内插拟合就是针对这个残缺数据而做相应的计算。

股票直线拟合怎么做?
股票直线拟合怎么做?

二、插值拟合

1、概念导入

给你两个数,(129 192)让你在这两个数中间线性**一个值,那我们**的方式其实有很多,例如在二维平面线性 B2 =(A1+A2)/2, 或三维空间线性,我们用到如下。

A: 192 129

B :192 157 129

股票直线拟合怎么做?

2、一阶拟合

经过上述的概念导入,我们接下对我们残缺数据进行一阶插值,参考下面公式

股票直线拟合怎么做?

黑色框为有效的顶点数据,我们用算法遍历任意两个顶点中,只要有一个是有效数据(黑框的,这点很重要,很重要,很重要,重要的事情说三遍)插一个数据(蓝色框数据),如行相邻的129与192两个顶点内插157,129与267两个相邻顶点内插186。结果如下图所示。

靠前阶内插完了吗? 不,再看上面的数据。插完靠前轮后新增相邻顶点还需要内插,由于我的数据矩阵不大,第二轮插完后我的数据已经找不到相邻空白点了。如果数据矩阵非常大,在算法实现上,这个节点我们的代码可以写成递归形式,直到插完才开始退栈(效果如下)。

>3、二阶拟合

我们的二阶拟合完全是根据一阶拟合的结果来调整的,我的二阶设计是通过四个有效顶点(黑色框框为有效顶点,重要的事情说一遍)内插一个绿框数据。AD对角,BC对角,参考如下公式

注:插完之后不执行一阶拟合,一阶是需要黑框数据支持,任意两个蓝色、绿色数据不需要。

为了好看些,我把一阶的内容拿掉,如果一阶拟合出来的结果不是我下面的模型,需要进行模型转换。具体怎么弄,自己想想看,或给我留言吧(核心私聊)。

4、三阶拟合

三阶插值跟二阶插值基本是一样的,唯一不一样的地方:任意四个顶点(不限定是否黑框),参考上面公式**。得到第二个图。

5、删除辅助点

最后告诉大家,上面除了黑色框框的为有效数据,其他颜色框框的为辅助数据,类似于解几何题作的辅助线。现在把他删掉吧。参考删除之后的结果如下。

6、与数据原型比较偏差

平均偏差7%最大偏差22%,相对来说拟合度挺好的,别忘了我只有10%的数据,要拟合90%的数据。

7、结果输出

8、总结

95%时间在想办法,5%的时间把办法写成代码。看完本篇,是否曲面拟合可以轻松搞定了?

阿伦尼斯线性模型拟合

阿伦尼斯(Arrhenius Mode)应该是所有经验加速模型中被最常采用的了,大部分资料中给出的加速因子计算式一般为:EXP[Ea/K*(1/T0-1/T1)],式中有个未知量Ea(激活能/活化能),标准中一般要求实测该值(也有部分试验人员会引用一些大厂提供的参考值,有文献研究认为这些参考值可能会成为一个误导,可信度存疑)。

那该怎么来计算这个Ea值呢?其实这里就用到了阿伦尼斯一个线性表达式:LN(L)=A+B/T,式中A和B(Ea/K)是未知常数,L和T分别表示寿命和开氏温度( T在500K以内,Ea值相对稳定,电子产品一般不会超过此限值 )。由此可见,只要有一些失效数据即可求解出拟合出寿命和温度的关系表达式,具体操作步骤如下:

①将待计算Ea的样品分成3组或者4组(每组推荐不少于3个样品,对于部分价值较高的样品1个也行,一些标准文献上给的就是1个)

②确定合适的试验温度点及梯度(一般会在极限温度以下,温度点数等于组数,梯度尽可能快暴露缺陷为宜)

以上就是股票直线拟合怎么做?的详细内容,希望通过阅读小编的文章之后能够有所收获!

版权:本文由用户自行上传,观点仅代表作者本人,本站仅供存储服务。如有侵权,请联系管理员删除,了解详情>>

发布
问题