一种相对完整的卷积计算

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

1. 消除速度计算

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

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

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

吸收速度

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

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

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

5. 卷积模拟

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

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

模拟的PK 曲线为：

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