CN106919797B - 一种量子拉普拉斯特征映射方法 - Google Patents
一种量子拉普拉斯特征映射方法 Download PDFInfo
- Publication number
- CN106919797B CN106919797B CN201710122846.6A CN201710122846A CN106919797B CN 106919797 B CN106919797 B CN 106919797B CN 201710122846 A CN201710122846 A CN 201710122846A CN 106919797 B CN106919797 B CN 106919797B
- Authority
- CN
- China
- Prior art keywords
- matrix
- quantum
- data
- laplace
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种量子拉普拉斯特征映射方法,在现有的拉普拉斯特征映射算法之上,将拉普拉斯矩阵视为数据集的协方差矩阵,可以简便地得到一个密度矩阵,同时将现有的特征向量问题进行相应的转换,采用量子的方式进行计算。本发明提出了一种量子版本的拉普拉斯特征映射方法‑QLE(Quantum Laplacian Eigenmaps),应用了共轭链以及矩阵运算以解决非线性的降维问题。比起经典的拉普拉斯特征映射需要的多项式时间,本发明可以提供指数级的加速。
Description
技术领域
本发明涉及一种量子拉普拉斯特征映射方法。
背景技术
机器学习与数据分析在降维,预测和分类中扮演着越来越重要的角色。在许多例子中原有的数据在高维的特征空间,如一张有n平方像素的图片(每个像素作为一个特征)。所以为了分析这些高维度的特征数据,我们需要将自然结构视为低维流形嵌入高维空间的数据降维。
为了将高维的数据降维,无论我们选择哪种方式,我们都需要考虑所需要的时间。正如我们所知的,一个设计良好的量子算法可以极大的改善我们的经典算法。劳埃德等人提出量子版本的PCA,可以指数级的提高算法速度。cong等人泛化了HHL算法,使之可以应用于量子判别式分析。尽管如此,这里依旧没有非线性的量子版本的降维方法。
发明内容
本发明的目的在于克服现有技术的不足,提供一种量子拉普拉斯特征映射方法,应用了共轭链以及矩阵运算,解决非线性的降维问题,指数级地加快原有的拉普拉斯特征算法。
本发明的目的是通过以下技术方案来实现的:一种量子拉普拉斯特征映射方法,包括以下步骤:
S1:利用数据的位置信息建立一个图G,顶点V是数据,边E是不同领域数据的相似性;为了使数据降维,需要最小化目标函数J(u):
式中,yi是数据点xi的低维表现,wij对应xi与xj的权重,L代表图G的拉普拉斯矩阵;
S2:将目标函数min(2YTLY)求解转化为广义特征值求解:
Lv=λDv
式中,D是一个对角矩阵,Dii=∑jW(i,j),特征向量v的最小非零特征值构造出数据的低维表示Y,λ表示特征值;
S3:将拉普拉斯矩阵L视为数据集的协方差矩阵,得到一个密度矩阵,即L=I·IT;其中I是图G=(V,E)的关联矩阵;所述的关联矩阵I存储每个节点及其连接边之间的关系,如果一个有向边j从点i出发,则Iij=1,如果在点i结束,则Iij=-1,否则Iij=0;
S4:将步骤S2中的广义特征值求解转换为:
D-1I·ITv=λv;
S5:将关联矩阵I和对角矩阵D转化为可以在量子随机存储器QRAM中输入的形式,ai为关联矩阵I的列,di为对角矩阵D的列;
S6:访问QRAM以得到关联矩阵I和对角矩阵D的量子态:
O(|i>|0>|0>)→|i>|di>||di|>
O(|i>|0>|0>)→|i>|ai>||ai|>
S7:通过QRAM构造|ψ1>和|ψ2>的状态:
式中,I是|ψ1>的密度矩阵,D是|ψ2>的密度矩阵,m表示为列的数量;
S8:得到对应量子态的密度矩阵,由于|ψ1>的第二个寄存器中的密度矩阵|ai>与I、以及|ψ2>的第二个寄存器中的密度矩阵|di>和D都是成正比的,因此:
S9:由于D和I都是半正定的厄米算符,令:
u=I-1/2v
于是,步骤S4中的求解转换为:
I1/2D-1/2D-1/2I1/2u=λu
令L=I1/2D-1/2D-1/2I1/2,采用共轭链式乘法,转化L:
其中:
S10:采矩阵运算技术,并应用量子相位估计得到|φ>的状态:
式中,φ表示最终的结果态;通过对|φ>的取样,得到特征向量u,并进一步得到特征向量v与L。
本发明的有益效果是:本发明提出了一种量子版本的拉普拉斯特征映射方法-QLE(Quantum Laplacian Eigenmaps),应用了共轭链以及矩阵运算以解决非线性的降维问题。比起经典的拉普拉斯特征映射需要的多项式时间,本发明可以提供指数级的加速。QLE可以在O(poly(log(mn))/ε3)进行非线性降维,其中m是数据点的数量,n是边的数量,ε为容许误差。
具体实施方式
下面进一步详细描述本发明的技术方案:
对于经典拉普拉斯特征映射算法:拉普拉斯特征映射算法为假设在高维空间中的数据有对应的低维结构。利用数据的位置信息建立一个图G,顶点V是数据,边E是不同领域数据的相似性。
为了使数据降维,我们要最小化目标函数J(u),通过下面的方程:
其中,yi是数据点xi的低维表现,wij对应xi与xj的权重,L代表图G的拉普拉斯矩阵。
而对于优化min(2YTLY)的问题可以转化为一个广义特征值问题:
Lv=λDv
其中D是一个对角矩阵,Dii=∑jW(i,j),Wij对应与xi与xj的权重,λ表示特征值。
最后特征向量v的最小非零特征值可以构造出数据的低维表示Y。
而在本实施例中的量子拉普拉斯特征映射方法,将拉普拉斯矩阵视为数据集的协方差矩阵,则可以简便地得到一个密度矩阵,即L=I·IT,其中I是图G=(V,E)的关联矩阵。关联矩阵存储每个节点及其连接边之间的关系。如果一个有向边j从点i出发,则Iij=1,如果在点i结束,则Iij=-1,否则Iij=0(边j与点i之间没有连接)。
最终的特征向量问题可以转换为:
D-1I·ITv=λv
我们的任务就是将拉普拉斯矩阵转换为密度矩阵,并解出上面的方程。
解决问题之前,需要将上述中提到的I与D转化为可以在量子随机存储器QRAM中输入 的形式,并假设I的列为ai,D的列为di。于是我们可以通过下面的一些过程来得到关联矩阵I=∑i|ai||ai><i|,对角矩阵D=∑i|di||di><i|的量子形式。
(1)访问QRAM,通过QRAM构造关联矩阵I和对角矩阵D的量子态:
O(|i>|0>|0>)→|i>|di>||di|>
O(|i>|0>|0>)→|i>|ai>||ai|>
其中O()表示Operation,存取操作,即分别表示对QRAM进行从0到di以及对QRAM进行从0到ai的获取操作。
并且构造|ψ1>和|ψ2>的状态:
(2)得到对应量子态的密度矩阵,由于现在第二个寄存器中的密度矩阵在两种情况下与I和D都是成正比的:
归一化的系数也是成正比的原因。
(3)由于D和I都是半正定的厄米算符(Hermitian matrices),令:
u=I-1/2v
于是,原问题转换为:
I1/2D-1/2D-1/2I1/2u=λu
令L=I1/2D-1/2D-1/2I1/2,并使用文献[1]中提出的共轭链式乘法(Hermitian ChainProduct),转化L:
其中:
f2(X)=X-1/2
f1(X)=X1/2。
文献[1]为Cong I,Duan L.Quantum Discriminant Analysis forDimensionality Reduction and Classification[J].arXiv preprint arXiv:1510.00113,2015。具体地,这里文献[1]就是提出了将D-1I·ITv=λv转化成I1/2D-1/2D-1/2I1/ 2u=λu。
(4)最后的任务只剩下解出(3)中得到的最后方程组,方法为:使用文献[2]中的矩阵运算技术(matrix exponentiation technique,就是一种方法计算e-iX,这里X是矩阵),并应用量子相位估计得到|φ>的状态:
式中,φ表示最终的结果态,最后求的结果都保存在这个态里面;通过对|φ>的取样,得到特征向量u;根据u=I-1/2v进一步得到特征向量v;根据Lv=λDv得到L。
文献[2]为Harrow A W,Hassidim A,Lloyd S.Quantum algorithm for linearsystems of equations[J].Physical review letters,2009,103(15):150502。
Claims (1)
1.一种量子拉普拉斯特征映射方法,用于对具有n平方像素的图片数据进行降维,每个像素作为一个特征数据,其特征在于:包括以下步骤:
S1:利用数据的位置信息建立一个图G,顶点V是数据,边E是不同数据的相似性;为了使数据降维,需要最小化目标函数J(Y):
式中,yi是数据点xi的低维表现,wij对应xi与xj的权重,L代表图G的拉普拉斯矩阵;
S2:将目标函数min(2YTLY)求解转化为广义特征值求解:
Lv=λDv
式中,D是一个对角矩阵,Dii=∑jW(i,j),特征向量v的最小非零特征值构造出数据的低维表示Y,λ表示特征值;
S3:将拉普拉斯矩阵L视为数据集的协方差矩阵,得到一个关联矩阵I,即L=I·IT;其中I是图G=(V,E)的关联矩阵;所述的关联矩阵I存储每个节点及其连接边之间的关系,如果一个有向边从点i出发,则Iij=1,如果在点i结束,则Iij=-1,否则Iij=0;
S4:将步骤S2中的广义特征值求解转换为:
D-1I·ITv=λv;
S5:将关联矩阵I和对角矩阵D转化为可以在量子随机存储器QRAM中输入的形式,ai为关联矩阵I的列,di为对角矩阵D的列;
S6:访问QRAM以得到关联矩阵I和对角矩阵D的量子态:
O(|i>|0>|0>)→|i>|di>||di|>
O(|i>|0>|0>)→|i>|ai>||ai|>
S7:通过QRAM构造|ψ1>和|ψ2>的状态:
式中,I是|ψ1>的密度矩阵,D是|ψ2>的密度矩阵,m表示为列的数量;
S8:得到对应量子态的密度矩阵,由于|ψ1>的第二个寄存器中的密度矩阵|ai>与I、以及|ψ2>的第二个寄存器中的密度矩阵|di>和D都是成正比的,因此:
S9:由于D和I都是半正定的厄米算符,令:
u=I-1/2v
于是,步骤S4中的求解转换为:
I1/2D-1/2D-1/2I1/2u=λu
其中:
S10:采矩阵运算技术,并应用量子相位估计得到|φ>的状态:
式中,φ表示最终的结果态;通过对|φ>的取样,得到特征向量u,并进一步得到特征向量v与L。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710122846.6A CN106919797B (zh) | 2017-03-03 | 2017-03-03 | 一种量子拉普拉斯特征映射方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710122846.6A CN106919797B (zh) | 2017-03-03 | 2017-03-03 | 一种量子拉普拉斯特征映射方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106919797A CN106919797A (zh) | 2017-07-04 |
CN106919797B true CN106919797B (zh) | 2020-04-07 |
Family
ID=59460877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710122846.6A Active CN106919797B (zh) | 2017-03-03 | 2017-03-03 | 一种量子拉普拉斯特征映射方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106919797B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109740758B (zh) * | 2019-01-09 | 2023-04-07 | 电子科技大学 | 一种基于量子计算的核方法 |
CN112068798B (zh) * | 2020-08-14 | 2023-11-03 | 本源量子计算科技(合肥)股份有限公司 | 一种实现网络节点重要性排序的方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372719A (zh) * | 2016-08-31 | 2017-02-01 | 浙江莱达信息技术有限公司 | 一种医学x光影像处理控制参量的智能优选方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1937002B1 (en) * | 2006-12-21 | 2017-11-01 | Rohde & Schwarz GmbH & Co. KG | Method and device for estimating the image quality of compressed images and/or video sequences |
CN102867171B (zh) * | 2012-08-23 | 2014-11-19 | 山东师范大学 | 一种基于标记传播与邻域嵌入保持的人脸表情识别方法 |
EP3311140B8 (en) * | 2015-06-19 | 2021-04-21 | The Government of the United States of America as represented by the Secretary of the Navy | Determining extracellular analyte concentration with nanoplasmonic sensors |
-
2017
- 2017-03-03 CN CN201710122846.6A patent/CN106919797B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372719A (zh) * | 2016-08-31 | 2017-02-01 | 浙江莱达信息技术有限公司 | 一种医学x光影像处理控制参量的智能优选方法 |
Non-Patent Citations (3)
Title |
---|
Quantum algorithm for linear systems of equations;Harrow A W,Hassidim A,Lloyd S;《Physical review letters》;20091231;第103卷(第15期);全文 * |
Quantum support vector machine for big data classification;Rebentrost P, Mohseni M,Lloyd S;《Physical Review Letters》;20141231;第113卷(第13期);全文 * |
遥感影像的张量表达与流形学习方法研究;张乐飞;《中国博士学位论文全文数据库 信息科技辑》;20140615(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106919797A (zh) | 2017-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cichocki et al. | Tensor networks for dimensionality reduction and large-scale optimization: Part 2 applications and future perspectives | |
Cao et al. | The sparse matrix transform for covariance estimation and analysis of high dimensional signals | |
Chakraborty et al. | An image denoising technique using quantum wavelet transform | |
CN107180013B (zh) | 一种量子d(4)小波变换实现量子线路设计的方法 | |
CN108520495B (zh) | 基于聚类流形先验的高光谱图像超分辨重建方法 | |
Luttinen et al. | Transformations in variational Bayesian factor analysis to speed up learning | |
Ma et al. | Sparse representation and position prior based face hallucination upon classified over-complete dictionaries | |
Maggu et al. | Kernel transform learning | |
CN106919797B (zh) | 一种量子拉普拉斯特征映射方法 | |
Maggu et al. | Greedy deep transform learning | |
Freris et al. | Compressed sensing of streaming data | |
Sterck et al. | An adaptive algebraic multigrid algorithm for low-rank canonical tensor decomposition | |
Song et al. | Robust exponential squared loss-based variable selection for high-dimensional single-index varying-coefficient model | |
Lu et al. | Linear programming support vector regression with wavelet kernel: A new approach to nonlinear dynamical systems identification | |
Shimodaira | Cross-validation of matching correlation analysis by resampling matching weights | |
Qi et al. | Two dimensional synthesis sparse model | |
Wang et al. | Hyperspectral unmixing via plug-and-play priors | |
Chaurra-Gutierrez et al. | QIST: One-dimensional quantum integer wavelet S-transform | |
Zhang et al. | Low-rank tensor completion via combined tucker and tensor train for color image recovery | |
Song et al. | A parameterized fuzzy adaptive K-SVD approach for the multi-classes study of pursuit algorithms | |
Menon et al. | Novel EM based ML Kalman estimation framework for superresolution of stochastic three-states microtubule signal | |
Montalto et al. | A linear approach for sparse coding by a two-layer neural network | |
Gkillas et al. | Efficient coupled dictionary learning and sparse coding for noisy piecewise-smooth signals: Application to hyperspectral imaging | |
Yang et al. | Principal whitened gradient for information geometry | |
Khader et al. | A model-guided deep convolutional sparse coding network for hyperspectral and multispectral image fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |