技术领域
[0001] 本发明涉及一种量子拉普拉斯特征映射方法。
相关背景技术
[0002] 机器学习与数据分析在降维，预测和分类中扮演着越来越重要的角色。在许多例子中原有的数据在高维的特征空间，如一张有n平方像素的图片(每个像素作为一个特征)。所以为了分析这些高维度的特征数据，我们需要将自然结构视为低维流形嵌入高维空间的数据降维。
[0003] 为了将高维的数据降维，无论我们选择哪种方式，我们都需要考虑所需要的时间。正如我们所知的，一个设计良好的量子算法可以极大的改善我们的经典算法。劳埃德等人提出量子版本的PCA，可以指数级的提高算法速度。cong等人泛化了HHL算法，使之可以应用于量子判别式分析。尽管如此，这里依旧没有非线性的量子版本的降维方法。
具体实施方式
[0038] 下面进一步详细描述本发明的技术方案：
[0039] 对于经典拉普拉斯特征映射算法：拉普拉斯特征映射算法为假设在高维空间中的数据有对应的低维结构。利用数据的位置信息建立一个图G，顶点V是数据，边E是不同领域数据的相似性。
[0040] 为了使数据降维，我们要最小化目标函数J(u)，通过下面的方程：
[0041]
[0042] 其中，yi是数据点xi的低维表现，wij对应xi与xj的权重，L代表图G的拉普拉斯矩阵。
[0043] 而对于优化min(2YTLY)的问题可以转化为一个广义特征值问题：
[0044] Lv＝λDv
[0045] 其中D是一个对角矩阵，Dii＝∑jW(i,j)，Wij对应与xi与xj的权重，λ表示特征值。
[0046] 最后特征向量v的最小非零特征值可以构造出数据的低维表示Y。
[0047] 而在本实施例中的量子拉普拉斯特征映射方法，将拉普拉斯矩阵视为数据集的协方差矩阵，则可以简便地得到一个密度矩阵，即L＝I·IT，其中I是图G＝(V,E)的关联矩阵。关联矩阵存储每个节点及其连接边之间的关系。如果一个有向边j从点i出发，则Iij＝1，如果在点i结束，则Iij＝-1，否则Iij＝0(边j与点i之间没有连接)。
[0048] 最终的特征向量问题可以转换为：
[0049] D-1I·ITv＝λv
[0050] 我们的任务就是将拉普拉斯矩阵转换为密度矩阵，并解出上面的方程。
[0051] 解决问题之前，需要将上述中提到的I与D转化为可以在量子随机存储器QRAM中输入 的形式，并假设I的列为ai，D的列为di。于是我们可以通过下面的一些过程来得到关联矩阵I＝∑i|ai||ai＞＜i|，对角矩阵D＝∑i|di||di＞＜i|的量子形式。
[0052] (1)访问QRAM，通过QRAM构造关联矩阵I和对角矩阵D的量子态：
[0053] O(|i＞|0＞|0＞)→|i＞|di＞||di|＞
[0054] O(|i＞|0＞|0＞)→|i＞|ai＞||ai|＞
[0055] 其中O()表示Operation，存取操作，即分别表示对QRAM进行从0到di以及对QRAM进行从0到ai的获取操作。
[0056] 并且构造|ψ1＞和|ψ2＞的状态：
[0057]
[0058]
[0059] 式中，I是|ψ1＞的密度矩阵，D是|ψ2＞的密度矩阵，m表示列的数量，每一个|＞均表示一个寄存器； 和 为归一化的系数，归一化的原因为为了满足量子力学的基本假设。
[0060] (2)得到对应量子态的密度矩阵，由于现在第二个寄存器中的密度矩阵在两种情况下与I和D都是成正比的：
[0061]
[0062]
[0063] 归一化的系数也是成正比的原因。
[0064] (3)由于D和I都是半正定的厄米算符(Hermitian matrices)，令：
[0065] u＝I-1/2v
[0066] 于是，原问题转换为：
[0067] I1/2D-1/2D-1/2I1/2u＝λu
[0068] 令L＝I1/2D-1/2D-1/2I1/2，并使用文献[1]中提出的共轭链式乘法(Hermitian Chain Product)，转化L：
[0069]
[0070] 其中：
[0071] f2(X)＝X-1/2
[0072] f1(X)＝X1/2。
[0073] 文献[1]为Cong I,Duan L.Quantum Discriminant Analysis for Dimensionality Reduction and Classification[J].arXiv preprint arXiv:1510.00113,2015。具体地，这里文献[1]就是提出了将D-1I·ITv＝λv转化成I1/2D-1/2D-1/2I1/
2
u＝λu。
[0074] (4)最后的任务只剩下解出(3)中得到的最后方程组，方法为：使用文献[2]中的矩阵运算技术(matrix exponentiation technique，就是一种方法计算e-iX，这里X是矩阵)，并应用量子相位估计得到|φ＞的状态：
[0075]
[0076] 式中，φ表示最终的结果态，最后求的结果都保存在这个态里面；通过对|φ＞的取样，得到特征向量u；根据u＝I-1/2v进一步得到特征向量v；根据Lv＝λDv得到L。
[0077] 文献[2]为Harrow A W,Hassidim A,Lloyd S.Quantum algorithm for linear systems of equations[J].Physical review letters,2009,103(15):150502。
