CN112906899A - 一种基于量子计算的最小均方误差检测方法 - Google Patents
一种基于量子计算的最小均方误差检测方法 Download PDFInfo
- Publication number
- CN112906899A CN112906899A CN202110318825.8A CN202110318825A CN112906899A CN 112906899 A CN112906899 A CN 112906899A CN 202110318825 A CN202110318825 A CN 202110318825A CN 112906899 A CN112906899 A CN 112906899A
- Authority
- CN
- China
- Prior art keywords
- quantum
- formula
- arrival direction
- equation
- state
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title description 3
- 239000011159 matrix material Substances 0.000 claims abstract description 59
- 239000013598 vector Substances 0.000 claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims description 9
- 238000002360 preparation method Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000005610 quantum mechanics Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Databases & Information Systems (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Radio Transmission System (AREA)
Abstract
本发明公开了一种基于量子计算的混合大规模MIMO到达方向估计方法,包括以下步骤:S1:根据量子奇异值估计构造量子态形式的空间协方差矩阵;S2:根据变分原理实现密度矩阵特征分解算法并对步骤S1中的协方差矩阵进行特征分解得到特征向量;S3:根据步骤S2得到的特正向量实现量子标记操作并实现到达方向搜索。本发明针对混合大规模MIMO通信系统中传统的到达方向估计算法,采用量子奇异值估计算法、量子可获取数据结构以及变分量子算法,降低了传统到达方向估计问题的复杂度,从而得到更稳定的性能。
Description
技术领域
本发明涉及通信领域,特别是涉及一种基于量子计算的混合大规模MIMO到达方向估计方法。
背景技术
1982年,Feynman描述了量子计算的巨大潜力,并且建议在量子力学原理的基础上构造量子计算机,以此来挖掘量子计算的潜力。进一步的,1994年Shor提出素数因子分解和离散对数问题的多项式时间量子算法;1995年Grover提出了在没有结构的搜索空间上进行搜索的量子算法。这些量子计算算法都展示了量子计算的特殊潜力,对传统的算法提供了加速。另一方面,随着技术和需求的增加,通信系统发展迅速,特别是大规模多输入多输出系统(Massive MIMO:Massive Multiple-Input Multiple-Output)系统作为移动通信领域的关键技术之一,现如今越来越受到青睐。由于大规模MIMO系统的研究存在大量基站天线来联合服务于多个用户的现象,因此这类研究涉及了很多高计算复杂度的问题,并且其中一部分至今也没有较好的处理方式。例如信号达到方向估计在混合大规模MIMO系统中有较多的应用,但为了避免复杂度为O(poly(MQ))的大规模矩阵特征分解和协方差矩阵重构,研究人员只好牺牲部分检测性能来降低计算复杂度。目前,量子计算还没有被应用到混合大规模MIMO到达方向估计方法的先例。
发明内容
技术问题:本发明的目的是提供一种基于量子计算的混合大规模MIMO到达方向估计方法,将传统混合大规模MIMO到达方向估计方法大规模矩阵特征分解,协方差重构以及到达方向查找问题的复杂度降低至O(poly(L log(MQ))),从而更加适用于大规模多用户系统应用场景中并且得到更快的运算速度。
技术方案:为达到此目的,本发明的一种基于量子计算的混合大规模MIMO到达方向估计方包括以下步骤:
S1:根据量子奇异值估计构造量子态形式的空间协方差矩阵;
S2:根据变分原理实现密度矩阵特征分解算法并对步骤S1中的协方差矩阵进行特征分解得到特征向量;
S3:根据步骤S2得到的特征向量实现量子标记操作并实现到达方向搜索。
其中:
所述能量矢量P根据式(2)构造:
式(2)中,UP为构造能量矢量的量子过程且能够通过复杂度O(poly(log Q))实现,i为索引值且表示向量的i个元素;矩阵A根据式(3),(4)得到:
式(3),(4)中,UM为构造公式(3)量子态的量子过程且能够以复杂度O(poly(logM))实现,UN为构造公式(4)量子态的量子过程且能够以复杂度O(poly(log Q))实现,j为索引值且表示矩阵A的第j个列。
所述步骤S1中构造量子态形式的空间协方差矩阵,其特定的空间协方差矩阵量子态为:
其中,αi=<ui|P>,γ为A的秩,|ui>、|vi>和σi分别为A的左、右奇异向量和奇异值。
所述空间协方差矩阵的量子态通过以下方法获得:
式(6)中,αi=<ui|P>;
S1.3:增加一个初始态为|0>的寄存器并执行由奇异值寄存器控制的受控旋转操作,通过式(8)得到量子态:
S1.4:执行逆奇异值估计操作并对第一个寄存器进行测量,得到式(9)的量子态:
S1.5:执行一个酉操作实现变换T:|ui>→|vi>,得到式(10)的量子态:
所述步骤S2具体包括以下步骤:
S2.1:根基步骤S1中生成的量子态,构造空间协方差矩阵相关的密度矩阵
式(11)中,该密度矩阵和空间协方差矩阵具有相同的特征向量且特征值存在平方归一化关系,|ui>和σi分别为A的左奇异向量和奇异值;
S2.2:利用冯诺依曼定理和变分量子算法设计成本函数,得到式(12)的函数:
S2.3:重新表述式(12)的成本函数,得到式(13):
S2.4:在量子计算机上计算式(13)成本函数,在传统计算机上通过梯度算法或者免梯度算法对成本函数进行优化,并将新的参数θ反馈到量子计算机上重新计算成本函数,重复执行上述传统-量子计算机上运算,直到参数的更新对成本函数的变化小于预设的阈值,得到参数θ*;
所述步骤S3具体包括以下步骤:
S3.1:制备由搜索空间方向向量构成的量子叠加态|φS>,如下式(14):
S3.2:对存储特征态的寄存器执行酉操作VH(θ*),得到式(15)量子态:
S3.4:对式(15)执行标记映射,得到式(18)的量子态:
式(18)中,方向向量在信号子空间和噪声子空间的投影可通过标记寄存器值为1或者0来有效地区分;
S3.5:测量标记寄存器并且对第一个计数寄存器进行抽样可获得式(19)的量子态:
式(19)中,Us为由步骤S2特征分解得到的特征向量构成的矩阵,由于概率越高,样本越容易被抽到,抽到次数最多的样本即为要预测的最优到达方向。
有益效果:本发明公开了一种基于量子计算的混合大规模MIMO到达方向估计方法,针对大规模通信系统中传统的混合大规模MIMO到达方向估计问题,采用量子方法构造空间协方差矩阵、执行密度矩阵特征分解和达到方向搜索,降低了传统混合大规模MIMO到达方向估计方法大规模矩阵特征分解、协方差矩阵重构和到达方向搜索的复杂度,从而得到更准确的方向估计,能更好的适用于如大规模MIMO系统在内的应用场景。
附图说明
图1为本发明具体实施方式中的方法流程图;
图2为本发明具体实施方式中步骤S1的流程图;
图3为本发明具体实施方式中步骤S2和S3的流程图。
具体实施方式
下面结合具体实施方式对本发明的技术方案作进一步的介绍。
本具体实施方式公开了一种基于量子计算的混合大规模MIMO到达方向估计方法,如图1所示,包括以下步骤:
S1:根据量子奇异值估计构造量子态形式的空间协方差矩阵;
S2:根据变分原理实现密度矩阵特征分解算法并对步骤S1中的协方差矩阵进行特征分解得到特征向量;
S3:根据步骤S2得到的特正向量实现量子标记操作并实现到达方向搜索。
步骤S1中,能量矢量P根据式(3)制备:
式(2)中,UP为构造能量矢量的量子过程且能够通过复杂度O(poly(log Q))实现,i为索引值且表示向量的第i个元素。矩阵A根据式(3),(4)得到:
式(3),(4)中,UM为构造公式(3)量子态的量子过程且能够以复杂度O(poly(logM))实现,UN构造公式(4)量子态的量子过程且能够以复杂度O(poly(log Q))实现,j为索引值且表示矩阵A的第j个列。
步骤S1中,制备得到的特定的空间协方差矩阵量子态为:
其中,αi=<ui|P>和γ为A的秩,|ui>、|vi>和σi分别为A的左、右奇异向量和奇异值。
空间协方差矩阵的量子态表示通过以下方法获得:
式(6)中,αi=<ui|P>;
S1.3:执行由奇异值寄存器控制的受控旋转操作,通过式(8)得到量子态:
S1.4:对第一个寄存器使用算子|0><0|进行测量,得到式(9)的量子态:
S1.5:执行一个酉操作实现变换T:|ui>→|vi>,得到式(10)的量子态:
步骤S2和S3如图3所示。步骤S2中具体包括以下步骤:
S2.1:根基步骤S1中生成的量子态,构造一个与空间协方差矩阵相关的密度矩阵
式(11)中,该密度矩阵和空间协方差矩阵具有相同的特征向量且特征值存在平方归一化关系,|ui>和σi分别为A的左奇异向量和奇异值;
S2.2:利用冯诺依曼定理和变分量子算法设计成本函数,得到式(12)的函数:
S2.3:重新表述式(12)的成本函数,得到式(13):
S2.4:在量子计算机上计算式(13)成本函数,在传统计算机上通过梯度算法或者免梯度算法对成本函数进行优化,并将新的参数θ反馈到量子计算机上重新计算成本函数,重复执行上述传统-量子计算机上运算,直到参数的更新对成本函数的变化小于预设的阈值,得到参数θ*。
步骤S3具体包括以下步骤:
S3.1:制备由搜索空间方向向量构成的量子叠加态|φS>,如下式(14):
S3.2:对存储特征态的寄存器执行酉操作VH(θ*),得到式(15)量子态:
S3.4:对(15)执行标记映射,得到式(18)的量子态:
式(18)中,方向向量在信号子空间和噪声子空间的投影可通过标记寄存器值为1或者0来有效地区分;
S3.5:测量标记寄存器并且对第一个计数寄存器进行抽样可获得式(19):
式(19)中,Us为由步骤S2特征分解得到的特征向量构成的矩阵。由于概率越高,样本越容易被抽到,抽到次数最多的样本即为要预测的最优到达方向。
方法的复杂度如下:第一,量子态|P>能够以复杂度O[poly(log Q)]有效地制备,酉操作UM和UN也能够通过适当的存储结构以及量子线路有效地构造,其复杂度可分别近似地估计为O[poly(log M)]和O[poly(log Q)],空间协方差矩阵的量子态形式能够以复杂度其中,κA和||A||F是矩阵A的条件数和F范数;第二,量子计算机可通过O[1/ε2]次测量、抽样计算成本函数。不是一般性,假设传统计算机上参数更新阶段的迭代步数为TI,那么总共需要TI/ε2个密度矩阵备份实现变分量子态特征分解算法,复杂度可近似为第三,以复杂度O(L)可实现L个受控操作CU(λi),i=1,...,L和进行O(polyL)次抽样可最终确定L个到达方向;最后,将上述步骤合并以总的复杂实现基于量子计算的混合大规模MIMO到达方向估计方法。
Claims (7)
1.一种基于量子计算的混合大规模MIMO到达方向估计方法,其特征在于:该估计方法包括以下步骤:
S1:根据量子奇异值估计构造量子态形式的空间协方差矩阵;
S2:根据变分原理实现密度矩阵特征分解算法并对步骤S1中的协方差矩阵进行特征分解得到特征向量;
S3:根据步骤S2得到的特征向量实现量子标记操作并实现到达方向搜索。
6.根据权利要求1所述的基于量子计算的混合大规模MIMO到达方向估计方法,其特征在于:所述步骤S2具体包括以下步骤:
S2.1:根基步骤S1中生成的量子态,构造空间协方差矩阵相关的密度矩阵
式(11)中,该密度矩阵和空间协方差矩阵具有相同的特征向量且特征值存在平方归一化关系,|ui>和σi分别为A的左奇异向量和奇异值;
S2.2:利用冯诺依曼定理和变分量子算法设计成本函数,得到式(12)的函数:
S2.3:重新表述式(12)的成本函数,得到式(13):
S2.4:在量子计算机上计算式(13)成本函数,在传统计算机上通过梯度算法或者免梯度算法对成本函数进行优化,并将新的参数θ反馈到量子计算机上重新计算成本函数,重复执行上述传统-量子计算机上运算,直到参数的更新对成本函数的变化小于预设的阈值,得到参数θ*;
7.根据权利要求1所述的基于量子计算的混合大规模MIMO到达方向估计方法,其特征在于:所述步骤S3具体包括以下步骤:
S3.1:制备由搜索空间方向向量构成的量子叠加态|φS>,如下式(14):
S3.2:对存储特征态的寄存器执行酉操作VH(θ*),得到式(15)量子态:
S3.4:对式(15)执行标记映射,得到式(18)的量子态:
式(18)中,方向向量在信号子空间和噪声子空间的投影可通过标记寄存器值为1或者0来有效地区分;
S3.5:测量标记寄存器并且对第一个计数寄存器进行抽样可获得式(19)的量子态:
式(19)中,Us为由步骤S2特征分解得到的特征向量构成的矩阵,由于概率越高,样本越容易被抽到,抽到次数最多的样本即为要预测的最优到达方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110318825.8A CN112906899B (zh) | 2021-03-25 | 2021-03-25 | 基于量子计算的混合大规模mimo到达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110318825.8A CN112906899B (zh) | 2021-03-25 | 2021-03-25 | 基于量子计算的混合大规模mimo到达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112906899A true CN112906899A (zh) | 2021-06-04 |
CN112906899B CN112906899B (zh) | 2024-04-12 |
Family
ID=76106530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110318825.8A Active CN112906899B (zh) | 2021-03-25 | 2021-03-25 | 基于量子计算的混合大规模mimo到达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112906899B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113434646A (zh) * | 2021-06-08 | 2021-09-24 | 天津大学 | 基于量子测量与自注意力机制的问答任务匹配模型及方法 |
CN115022130A (zh) * | 2022-05-16 | 2022-09-06 | 海南大学 | 基于量子力学去噪的毫米波大规模mimo信道估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109767007B (zh) * | 2018-12-10 | 2023-04-18 | 东南大学 | 一种基于量子计算的最小均方误差检测方法 |
CN112183756B (zh) * | 2020-09-30 | 2024-04-23 | 西北大学 | 一种可选择的精确量子主成分分析方法及应用 |
-
2021
- 2021-03-25 CN CN202110318825.8A patent/CN112906899B/zh active Active
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113434646A (zh) * | 2021-06-08 | 2021-09-24 | 天津大学 | 基于量子测量与自注意力机制的问答任务匹配模型及方法 |
CN115022130A (zh) * | 2022-05-16 | 2022-09-06 | 海南大学 | 基于量子力学去噪的毫米波大规模mimo信道估计方法 |
CN115022130B (zh) * | 2022-05-16 | 2024-02-06 | 海南大学 | 基于量子力学去噪的毫米波大规模mimo信道估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112906899B (zh) | 2024-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109767007B (zh) | 一种基于量子计算的最小均方误差检测方法 | |
CN106980106B (zh) | 阵元互耦下的稀疏doa估计方法 | |
De Lathauwer et al. | Independent component analysis and (simultaneous) third-order tensor diagonalization | |
CN112906899A (zh) | 一种基于量子计算的最小均方误差检测方法 | |
CN110109050B (zh) | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 | |
CN112731275B (zh) | 一种基于零化插值的互质阵部分极化信号参数估计方法 | |
CN112731278B (zh) | 一种部分极化信号的角度与极化参数欠定联合估计方法 | |
Liu et al. | Understanding and comparing scalable Gaussian process regression for big data | |
CN104809475A (zh) | 基于增量线性判别分析的多类标场景分类方法 | |
Cheng et al. | Towards overfitting avoidance: Tuning-free tensor-aided multi-user channel estimation for 3D massive MIMO communications | |
CN112149045A (zh) | 一种适用于大规模数据的降维、关联分析方法 | |
CN110174658A (zh) | 基于秩一降维模型和矩阵补全的波达方向估计方法 | |
Ji et al. | FastVGBS: A fast version of the volume-gradient-based band selection method for hyperspectral imagery | |
Meng et al. | Quantum algorithm for MUSIC-based DOA estimation in hybrid MIMO systems | |
Lee et al. | Speed-up of the Matrix Computation on the Ridge Regression. | |
Cheng et al. | Tensor subspace Tracking via Kronecker structured projections (TeTraKron) for time-varying multidimensional harmonic retrieval | |
Liu et al. | A novel architecture to eliminate bottlenecks in a parallel tiled QRD algorithm for future MIMO systems | |
CN107645460A (zh) | 实值平行因子分解的多径参数估算方法 | |
Yu et al. | Dvsai: Diverse view-shared anchors based incomplete multi-view clustering | |
CN114201453A (zh) | 基于云计算的大规模高光谱数据张量环分解随机抽样方法 | |
Bartelmaos et al. | Fast principal component extraction using givens rotations | |
CN107801149B (zh) | 实值平行因子分解的多径参数估算方法 | |
CN113052319B (zh) | 一种量子计算的克莱姆-斯密特正交化方法 | |
Vandecappelle et al. | Updating the multilinear utv decomposition | |
Wang et al. | Compressed Sensing Based Data Acquisition Via Intelligent Vehicles |
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 |