[发明专利]一种基于广义N体问题的高斯过程回归计算方法

摘要

本发明公开了一种基于广义N体问题的高斯过程回归计算方法，包括基于双kd树的数据集划分法、双kd树的遍历剪枝法、基于高阶分治的核函数矩阵求解法和Cholesky分解算法，主要步骤包括运用kd树对输入数据集进行空间上的划分，再运用双kd树遍历剪枝法求得简化欧氏距离并存储；根据简化欧式距离及相关索引，采用平方指数核函数，求得训练数据集与测试数据集的核函数矩阵K*以及训练数据集自身的核函数矩阵K；运用Cholesky分解算法快速求得核函数矩阵K的逆矩阵K‑1；最后，由核函数矩阵K的逆矩阵K‑1、核函数矩阵K*以及输入的目标函数值，可以求得预测结果向量。通过该方法提高了高斯过程回归的预测效率和处理大数据的能力，推动高斯过程回归在大数据分析中的广泛应用。

主权项

一种基于广义N体问题的高斯过程回归计算方法，其特征在于，该方法是利用基于双kd树的数据集划分法、双kd树的遍历剪枝法、基于高阶分治的核函数矩阵求解法和Cholesky分解算法协同实现高斯过程回归计算；高斯过程回归的计算公式是：f*‾=K*T(K+σn2I)-1y]]>其中K*T是训练数据集与测试数据集之间的核函数矩阵，K是训练数据集自身的核函数矩阵，σn2是噪音，I代表单位矩阵，y是输入的训练目标函数值向量，是预测目标值向量；具体包括下述步骤：运用基于双kd树的数据集划分法对输入数据集进行空间上的划分，再运用双kd树的遍历剪枝法求得简化欧氏距离并存储；根据简化欧式距离及相关索引，采用平方指数核函数，求得训练数据集与测试数据集的核函数矩阵K*以及训练数据集自身的核函数矩阵K；运用Cholesky分解算法快速求得核函数矩阵K的逆矩阵K‑1；最后，由训练数据集自身的核函数矩阵K的逆矩阵K‑1、训练数据集与测试数据集之间的核函数矩阵K*以及输入的目标函数值，求得预测结果向量；所述基于双kd树的数据集划分法，用于空间划分数据，在高斯过程回归的计算中，对于训练数据集R和测试数据集Q，分别按照相同的规则构建R和Q各自的kd树TR、TQ；通过对此两个kd树同时进行遍历、计算及剪枝，最终求得核函数矩阵；所述双kd树的遍历剪枝法，用于双kd树的遍历、计算及剪枝；从两个kd树各自的根节点开始，分别计算两个kd树对应节点的距离；设一个阈值∈，当计算得到的节点间距离大于∈的值时，根据核函数的计算性质，则认为此两节点间的距离为无穷大，对此两节点进行剪枝处理，不再递归这两节点及其所有的子节点；对于距离值小于∈的节点，继续递归计算直到递归完毕所有节点并且每个叶子节点都至少被计算过一次；所述基于高阶分治的核函数矩阵求解法，用于求解核函数矩阵；根据上述双kd树的计算结果，求得核函数矩阵的值；核函数是定义数据点间相似性或距离的函数，因此根据两数据点间的距离x‑x′计算相应的核函数值；所述Cholesky分解算法，用于核函数矩阵求逆；即应用Cholesky分解的核函数矩阵快速求逆，对于矩阵求逆的过程，运用Cholesky分解算法加快求逆的效率，加速高斯过程回归的预测效率。