一种相对完整的卷积计算

前几天写了篇不完整的去卷积方法,那里所谓的不完整主要是指未进行相应的消除速度计算(该值从文献中引用),本文进行一些补充,主要的编程软件依然是python,数据同上篇文章。说明下:毛毛很懒,没有心思把细节做的尽善尽美,所以这里面的代码可能会报错,我只是展示个思路而已。

1. 消除速度计算

关于消除速度计算我需要说一句,我这里采用的是一种懒人算法,即通过插值计算出消除速度(该计算也仅仅针对非房室模型)。

首先绘制半对数图,后面直线的斜率为消除速度,因此只需要差分就可以求出来。

代码为:先用插值将时间间隔换为一致,均为一小时间隔(取平稳序列的均值),计算出ke = -0.21,与文献报道-1/3.2 = -0.21一致。

吸收速度

吸收速度求解的难度略微大于消除速度,需要解高阶方程。此处求出吸收速度为2.9.根据速度公式,我们进行一次重新拟合PK曲线,结果:

图1:为未处理的原图,图二未乘以放大系数之后的PK曲线。存在该问题的主要原因是

这个值我当作1处理,您可以将该值看作是缩放因子。

5. 卷积模拟

该卷积属于弱智卷积,在python numpy数据库中又np.convolve 函数可以直接调用,其卷积结果为:

红色为numpy卷积后的数据,绿色为实测值,存在差异的原因:是在计算过程中,我并未对时间序列进行处理。此处介绍一种新的自己编写的卷积函数。具体代码为

模拟的PK 曲线为:

红色为模拟曲线,蓝色为实测曲线

  • 发表于:
  • 原文链接https://kuaibao.qq.com/s/20180730G0H2HY00?refer=cp_1026
  • 腾讯「云+社区」是腾讯内容开放平台帐号(企鹅号)传播渠道之一,根据《腾讯内容开放平台服务协议》转载发布内容。

扫码关注云+社区

领取腾讯云代金券