CN113850389B - 一种量子线路的构建方法及装置 - Google Patents
一种量子线路的构建方法及装置 Download PDFInfo
- Publication number
- CN113850389B CN113850389B CN202010598078.3A CN202010598078A CN113850389B CN 113850389 B CN113850389 B CN 113850389B CN 202010598078 A CN202010598078 A CN 202010598078A CN 113850389 B CN113850389 B CN 113850389B
- Authority
- CN
- China
- Prior art keywords
- quantum
- matrix
- bit
- logic gate
- qubit
- 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
- 238000010276 construction Methods 0.000 title description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 334
- 239000002096 quantum dot Substances 0.000 claims abstract description 137
- 238000000034 method Methods 0.000 claims abstract description 67
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 40
- 238000005259 measurement Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 15
- 230000002441 reversible effect Effects 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 9
- 230000002829 reductive effect Effects 0.000 abstract description 13
- 238000004088 simulation Methods 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 230000036961 partial effect Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000005610 quantum mechanics Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000021615 conjugation Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
-
- 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
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
- Logic Circuits (AREA)
Abstract
本发明公开了一种量子线路的构建方法及装置,方法包括:获得N*N维矩阵A和N维向量b;获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,辅助量子比特和第一量子比特的初态置为|0>,第二量子比特的初态置为bj为向量b的第j个元素;确定矩阵A对应的酉矩阵U,将矩阵U分解成r个携带受控信息的单量子逻辑门对应的酉矩阵;输出包含r个携带受控信息的单量子逻辑门的子量子线路,根据子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。利用本发明实施例,能够减少量子线路复杂度,提高量子线路的模拟效率,同时减少硬件资源的占用。
Description
技术领域
本发明属于量子计算技术领域,特别是一种量子线路的构建方法及装置。
背景技术
量子计算机是一类遵循量子力学规律进行高速数学和逻辑运算、存储及处理量子信息的物理装置。当某个装置处理和计算的是量子信息,运行的是量子算法时,它就是量子计算机。量子计算机因其具有相对普通计算机更高效的处理数学问题的能力,例如,能将破解RSA密钥的时间从数百年加速到数小时,故成为一种正在研究中的关键技术。
量子计算模拟是一个借助数值计算和计算机科学来仿真遵循量子力学规律的模拟计算,作为一个仿真程序,它依据量子力学的量子比特的基本定律,利用计算机的高速计算能力,刻画量子态的时空演化。
目前,量子算法如HHL算法通常由量子线路实现,量子线路包括量子逻辑门操作和测量操作等等。HHL算法因其在解线性方程组时的指数级加速效果,成为重要的量子算法之一,但在解决实际问题过程中,如何构建实现HHL算法的量子线路,并实现高效运行,是一个非常普遍且困难的问题。现有的解决方法大多处于理论阶段,限制了在实际应用环境中的应用范围。例如,基于GLOA(Group Leaders Optimization Algorithm)进行矩阵分解的HHL算法的量子线路,其包含的量子逻辑门数量和种类较多,线路复杂度较高,导致量子线路的模拟效率较低,且对硬件资源的占用较多,实际应用价值不高。
发明内容
本发明的目的是提供一种量子线路的构建方法及装置,以解决现有技术中的不足,它能够减少HHL算法对应的量子线路的计算量,提高量子线路的模拟效率,同时减少硬件资源的占用。
本申请的一个实施例提供了一种量子线路的构建方法,包括:
获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,所述,所述为阶单位矩阵;
输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
可选的,所述确定所述矩阵A对应的酉矩阵U,包括:
若所述矩阵A为酉矩阵,则将所述矩阵A直接确定为对应的酉矩阵U;
若所述矩阵A为厄米矩阵,确定对应的酉矩阵;其中,所述t为常量;
若所述矩阵A为非厄米矩阵且非酉矩阵,确定对应厄米矩阵,以确定对应的酉矩阵,同时将所述向量b转换为。
可选的,所述将所述酉矩阵分解成个携带受控信息的单量子逻辑门对应的酉矩阵,包括:
确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序;
针对所述排序中的第i个非对角元素,构造特定量子逻辑门的阶酉矩阵,以使矩阵中与该非对角元素同位置的元素置为0,且不改变已置0的非对角元素;
其中,所述特定量子逻辑门包括操作一比特的单量子逻辑门,所述单量子逻辑门携带受其余比特控制的受控信息,; 当所述时,所述单量子逻辑门的酉矩阵由所述量子线路对应的酉矩阵的元素确定;当时,所述单量子逻辑门的酉矩阵由矩阵的元素确定;并且,当所述第i个非对角元素的排序为所在列最后一个时,同时使矩阵的相同列中的对角元素置为1。
可选的,所述确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序,包括:
当n=1时,所述酉矩阵U中对角元素下方的、待置0的非对角元素排序为(2,1);其中,所述(2,1)表示坐标为第2行第1列的非对角元素;
当n>1时,根据(n-1)比特量子线路对应的酉矩阵的第一列排序,确定n比特量子线路对应的所述酉矩阵U中对角元素下方的、待置0的非对角元素的第一列排序;其中,第一列中坐标(N/2+1,1)的非对角元素的排序位于第一列最后一个;
基于所述n比特量子线路对应的第一列排序,分别确定所述n比特量子线路对应的第2列至第N/2列中对角元素下方的、待置0的非对角元素的排序;
根据(n-1)比特量子线路对应的酉矩阵中对角元素下方的、待置0的非对角元素的排序,对应确定所述n比特量子线路对应的第(N/2+1)列至第N列中对角元素下方的、待置0的非对角元素的排序。
可选的,所述携带受控信息的单量子逻辑门的表示形式包括:
,其中,所述Cm表示0、1、*或单量子逻辑门V,所述m表示量子比特位,,并且,有且只有一个Cm表示单量子逻辑门V,所述单量子逻辑门V的酉矩阵由所述酉矩阵U确定;
当Cm为0时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为0态时,执行单量子逻辑门V;
当Cm为1时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为1态时,执行单量子逻辑门V;
当Cm为*时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为任意态时,均执行单量子逻辑门V。
可选的,所述根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路,包括:
根据所述子量子线路,构建相位估计运算对应的第一部分量子线路,以将|b>在所述矩阵A的特征空间上分解为:,并将所述第一量子比特和所述第二量子比特的初态|0>|b>转化为:;其中,所述第一量子比特的数量k取决于相位估计的精度和成功概率,所述第二量子比特的数量为所述n,所述为所述矩阵A的特征向量,所述为的振幅,所述为所述矩阵A的特征值;
构建受控旋转操作对应的第二部分量子线路,以将基态中的值提取到辅助量子比特的量子态振幅上,得到:,其中,所述辅助量子比特的数量为1,所述C为常数,;
构建相位估计逆运算对应的第三部分量子线路,以消除,得到:;
构建针对所述辅助量子比特的量子测量操作,以使测量到所述辅助量子比特的量子态为时,得到:,所述与为振幅归一化的对应关系;
将所述第一部分量子线路、所述第二部分量子线路、所述第三部分量子线路和所述量子测量操作,依次组成HHL算法对应的量子线路。
本申请的又一实施例提供了一种量子线路的构建装置,包括:
第一获得模块,用于获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
第二获得模块,用于获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
分解模块,用于确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,所述,所述为阶单位矩阵;
构建模块,用于输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
可选的,所述分解模块,具体用于:
若所述矩阵A为酉矩阵,则将所述矩阵A直接确定为对应的酉矩阵U;
若所述矩阵A为厄米矩阵,确定对应的酉矩阵;其中,所述t为常量;
若所述矩阵A为非厄米矩阵且非酉矩阵,确定对应厄米矩阵,以确定对应的酉矩阵,同时将所述向量b转换为。
可选的,所述分解模块,包括:
确定单元,用于确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序;
构造单元,用于针对所述排序中的第i个非对角元素,构造特定量子逻辑门的阶酉矩阵,以使矩阵中与该非对角元素同位置的元素置为0,且不改变已置0的非对角元素;
其中,所述特定量子逻辑门包括操作一比特的单量子逻辑门,所述单量子逻辑门携带受其余比特控制的受控信息,; 当所述时,所述单量子逻辑门的酉矩阵由所述量子线路对应的酉矩阵的元素确定;当时,所述单量子逻辑门的酉矩阵由矩阵的元素确定;并且,当所述第i个非对角元素的排序为所在列最后一个时,同时使矩阵的相同列中的对角元素置为1。
可选的,所述确定单元,具体用于:
当n=1时,所述酉矩阵U中对角元素下方的、待置0的非对角元素排序为(2,1);其中,所述(2,1)表示坐标为第2行第1列的非对角元素;
当n>1时,根据(n-1)比特量子线路对应的酉矩阵的第一列排序,确定n比特量子线路对应的所述酉矩阵U中对角元素下方的、待置0的非对角元素的第一列排序;其中,第一列中坐标(N/2+1,1)的非对角元素的排序位于第一列最后一个;
基于所述n比特量子线路对应的第一列排序,分别确定所述n比特量子线路对应的第2列至第N/2列中对角元素下方的、待置0的非对角元素的排序;
根据(n-1)比特量子线路对应的酉矩阵中对角元素下方的、待置0的非对角元素的排序,对应确定所述n比特量子线路对应的第(N/2+1)列至第N列中对角元素下方的、待置0的非对角元素的排序。
可选的,所述携带受控信息的单量子逻辑门的表示形式包括:
,其中,所述Cm表示0、1、*或单量子逻辑门V,所述m表示量子比特位,,并且,有且只有一个Cm表示单量子逻辑门V,所述单量子逻辑门V的酉矩阵由所述酉矩阵U确定;
当Cm为0时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为0态时,执行单量子逻辑门V;
当Cm为1时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为1态时,执行单量子逻辑门V;
当Cm为*时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为任意态时,均执行单量子逻辑门V。
可选的,所述构建模块,具体用于:
根据所述子量子线路,构建相位估计运算对应的第一部分量子线路,以将|b>在所述矩阵A的特征空间上分解为:,并将所述第一量子比特和所述第二量子比特的初态|0>|b>转化为:;其中,所述第一量子比特的数量k取决于相位估计的精度和成功概率,所述第二量子比特的数量为所述n,所述为所述矩阵A的特征向量,所述为的振幅,所述为所述矩阵A的特征值;
构建受控旋转操作对应的第二部分量子线路,以将基态中的值提取到辅助量子比特的量子态振幅上,得到:,其中,所述辅助量子比特的数量为1,所述C为常数,;
构建相位估计逆运算对应的第三部分量子线路,以消除,得到:;
构建针对所述辅助量子比特的量子测量操作,以使测量到所述辅助量子比特的量子态为时,得到:,所述与为振幅归一化的对应关系;
将所述第一部分量子线路、所述第二部分量子线路、所述第三部分量子线路和所述量子测量操作,依次组成HHL算法对应的量子线路。
本申请的又一实施例提供了一种存储介质,所述存储介质中存储有计算机程序,其中,所述计算机程序被设置为运行时执行上述任一项中所述的方法。
本申请的又一实施例提供了一种电子装置,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行上述任一项中所述的方法。
与现有技术相比,本发明提供的一种量子线路的构建方法,首先获得N*N维矩阵A和N维向量b,其中,矩阵A为可逆矩阵,N=2n, n为正整数;获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,辅助量子比特和第一量子比特的初态置为|0>,第二量子比特的初态置为,为向量b的第j个元素;确定矩阵A对应的酉矩阵U,将矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,,为阶单位矩阵;输出包含个携带受控信息的单量子逻辑门的子量子线路,根据子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。由于输出的子量子线路中量子逻辑门数量有限,且消除了酉矩阵形式复杂的多比特量子逻辑门,使得量子逻辑门形式得以简化,从而减少构建的HHL算法的量子线路复杂度,提高量子线路的模拟效率,同时减少硬件资源的占用。
附图说明
图1为本发明实施例提供的一种量子线路的构建方法的计算机终端的硬件结构框图;
图2为本发明实施例提供的一种量子线路的构建方法的流程示意图;
图3为本发明实施例提供的一种特定量子逻辑门位于量子线路中的示意图;
图4为本发明实施例提供的一种相位估计对应的第一部分量子线路示意图;
图5为本发明实施例提供的一种量子线路的构建装置的结构示意图。
具体实施方式
下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本发明实施例首先提供了一种量子线路的构建方法,该方法可以应用于电子设备,如计算机终端,具体如普通电脑、量子计算机等。
下面以运行在计算机终端上为例对其进行详细说明。图1为本发明实施例提供的一种量子线路的构建方法的计算机终端的硬件结构框图。如图1所示,计算机终端可以包括一个或多个(图1中仅示出一个)处理器102(处理器102可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)和用于存储数据的存储器104,可选地,上述计算机终端还可以包括用于通信功能的传输装置106以及输入输出设备108。本领域普通技术人员可以理解,图1所示的结构仅为示意,其并不对上述计算机终端的结构造成限定。例如,计算机终端还可包括比图1中所示更多或者更少的组件,或者具有与图1所示不同的配置。
存储器104可用于存储应用软件的软件程序以及模块,如本申请实施例中的量子线路的构建方法对应的程序指令/模块,处理器102通过运行存储在存储器104内的软件程序以及模块,从而执行各种功能应用以及数据处理,即实现上述的方法。存储器104可包括高速随机存储器,还可包括非易失性存储器,如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中,存储器104可进一步包括相对于处理器102远程设置的存储器,这些远程存储器可以通过网络连接至计算机终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
传输装置106用于经由一个网络接收或者发送数据。上述的网络具体实例可包括计算机终端的通信供应商提供的无线网络。在一个实例中,传输装置106包括一个网络适配器(Network Interface Controller,NIC),其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中,传输装置106可以为射频(Radio Frequency,RF)模块,其用于通过无线方式与互联网进行通讯。
需要说明的是,真正的量子计算机是混合结构的,它包含两大部分:一部分是经典计算机,负责执行经典计算与控制;另一部分是量子设备,负责运行量子程序进而实现量子计算。而量子程序是由量子语言如QRunes语言编写的一串能够在量子计算机上运行的指令序列,实现了对量子逻辑门操作的支持,并最终实现量子计算。具体的说,量子程序就是一系列按照一定时序操作量子逻辑门的指令序列。
在实际应用中,因受限于量子设备硬件的发展,通常需要进行量子计算模拟以验证量子算法、量子应用等等。量子计算模拟即借助普通计算机的资源搭建的虚拟架构(即量子虚拟机)实现特定问题对应的量子程序的模拟运行的过程。通常,需要构建特定问题对应的量子程序。本发明实施例所指量子程序,即是经典语言编写的表征量子比特及其演化的程序,其中与量子计算相关的量子比特、量子逻辑门等等均有相应的经典代码表示。
量子线路作为量子程序的一种体现方式,也称量子逻辑电路,是最常用的通用量子计算模型,表示在抽象概念下对于量子比特进行操作的线路,其组成包括量子比特、线路(时间线),以及各种量子逻辑门,最后常需要通过量子测量操作将结果读取出来。
不同于传统电路是用金属线所连接以传递电压信号或电流信号,在量子线路中,线路可看成是由时间所连接,亦即量子比特的状态随着时间自然演化,在这过程中按照哈密顿运算符的指示,一直到遇上逻辑门而被操作。
一个量子程序整体上对应有一条总的量子线路,本发明所述量子程序即指该条总的量子线路,其中,该总的量子线路中的量子比特总数与量子程序的量子比特总数相同。可以理解为:一个量子程序可以由量子线路、针对量子线路中量子比特的测量操作、保存测量结果的寄存器及控制流节点(跳转指令)组成,一条量子线路可以包含几十上百个甚至千上万个量子逻辑门操作。量子程序的执行过程,就是对所有的量子逻辑门按照一定时序执行的过程。需要说明的是,时序即单个量子逻辑门被执行的时间顺序。
需要说明的是,经典计算中,最基本的单元是比特,而最基本的控制模式是逻辑门,可以通过逻辑门的组合来达到控制电路的目的。类似地,处理量子比特的方式就是量子逻辑门。使用量子逻辑门,能够使量子态发生演化,量子逻辑门是构成量子线路的基础,量子逻辑门包括单比特量子逻辑门,如Hadamard门(H门,阿达马门)、泡利-X门(X门)、泡利-Y门(Y门)、泡利-Z门(Z门)、RX门、RY门、RZ门等等;多比特量子逻辑门,如CNOT门、CR门、iSWAP门、Toffoli门等等。量子逻辑门一般使用酉矩阵表示,而酉矩阵不仅是矩阵形式,也是一种操作和变换。一般量子逻辑门在量子态上的作用是通过酉矩阵左乘以量子态右矢对应的矩阵进行计算的。
参见图2,图2为本发明实施例提供的一种量子线路的构建方法的流程示意图,可以包括如下步骤:
S201,获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
具体的,可以获得用户输入的一个N*N的矩阵A和一个N维向量b。
目前,线性系统是很多科学和工程领域的核心,由于HHL算法在特定条件下相较于经典算法具有指数加速效果,未来能够在数据处理、机器学习、数值计算等场景具有广泛应用。HHL算法解决了一个求解线性方程的问题:输入一个N*N的矩阵A和一个N维向量b,输出n维向量x,满足Ax=b,即x=A-1b。由此,矩阵A需要满足为可逆矩阵,且由于下述需要将向量b的数据加载到量子线路的缘故,向量b的维数N能表示为2的正整数次幂。
S202,获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
具体的,可以获得输入的一组量子比特,具体可以量子比特位表示。例如,一位量子比特位上为0,表示该位量子比特的量子态为态,为1则表示态。
为了后续计算需要,该组量子比特区分为辅助量子比特、第一量子比特、第二量子比特,比特数可由用户根据需求确定,在计算资源充足的情况下,可设置较大数量的量子比特,满足多种计算需求。
并且,各量子比特的初态可由现有的振幅编码方式制备得到。其中,辅助量子比特和第一量子比特的初态置为|0>态,第二量子比特的初态置为。例如,对于4维向量b=[b0,b1,b2,b3],N=4,可得n=2。则,将向量b的数据编码到量子态振幅上,得到:
从而实现:将向量b的数据加载到量子线路中2个第二量子比特的量子态振幅上。
S203,确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,所述,所述为阶单位矩阵;
具体的,若矩阵A为酉矩阵,则将矩阵A直接确定为对应的酉矩阵U;
若矩阵A为厄米矩阵,可以通过哈密顿量模拟实现厄米矩阵到酉矩阵的转换,得到对应的酉矩阵;其中,t为常量,一般取值;
若矩阵A为非厄米矩阵且非酉矩阵,通过扩维将其转换,确定对应厄米矩阵,以确定对应的酉矩阵,同时将向量b转换为,也可写成,得,x即为前述HHL算法的解。
具体的,任意的酉矩阵都可以通过一系列量子逻辑门进行分解,进而转换为量子线路,因为一段量子线路本身就看成一个酉矩阵,当然,该量子线路不能包含任何不可逆操作,比如Measure(测量)或者Reset等操作。将矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵,可以包括:
S2031,确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序;
在一种实现方式中,酉矩阵中对角元素下方的、待置0的非对角元素的排序可以为:第一列按列号排至最后一列,每一列非对角元素按行号从上至下排序,一种两比特量子线路的4阶酉矩阵的排序示例如表1所示。
表1 一种两比特量子线路的酉矩阵的元素排序
其中,00、01、10、11表示行或列对应的二进制表示,二进制位与量子比特位一一对应;(1,1)、(2,2)、(3,3)、(4,4)表示坐标对应的对角元素,诸如(2,1)、(3,1)、(4,1)表示坐标对应的非对角元素,括号后边的数字1、2、3表示对应的排序。
需要强调的是,由于量子逻辑门的矩阵形式均为酉矩阵,即酉矩阵与酉矩阵的转置共轭的乘积是单位阵,且酉矩阵之间的乘积也是酉矩阵,故只需关注矩阵对角线下方的非对角元素的置0操作即可,对角元素上方的同一列非对角元素在该对角元素被置1同时会被置0,这是酉矩阵的特性决定的,不再进行赘述。同理,只关注矩阵对角线上方的非对角元素的置0操作也是可行的。
优选的,为了便于后续的矩阵构造,在另一种实现方式中,酉矩阵中对角元素下方的、待置0的非对角元素的排序可以为:
当n=1时,所述酉矩阵U中对角元素下方的、待置0的非对角元素排序为(2,1);其中,所述(2,1)表示坐标为第2行第1列的非对角元素;
当n>1时,根据(n-1)比特量子线路对应的酉矩阵的第一列排序,确定n比特量子线路对应的所述酉矩阵U中对角元素下方的、待置0的非对角元素的第一列排序;其中,第一列中坐标(N/2+1,1)的非对角元素的排序位于第一列最后一个;
基于所述n比特量子线路对应的第一列排序,分别确定所述n比特量子线路对应的第2列至第N/2列中对角元素下方的、待置0的非对角元素的排序;
根据(n-1)比特量子线路对应的酉矩阵中对角元素下方的、待置0的非对角元素的排序,对应确定所述n比特量子线路对应的第(N/2+1)列至第N列中对角元素下方的、待置0的非对角元素的排序。
示例性的,对于1比特量子线路的2阶酉矩阵,对角元素下方只有一个元素(2,1),故排序有且只有第一列的非对角元素为:(2,1)。
对于2比特量子线路的4阶酉矩阵,第一列排序采用递归的思想,即参考1比特量子线路的第一列排序,且第(N/2+1)行第1列的非对角元素排在该列最后一个,即(2,1)排在第1个、(3,1)排在最后1个,从而确定(4,1)排在第2个,最终可得第1列排序为:(2,1)、(4,1)、(3,1)。
对于3比特量子线路的8阶酉矩阵,第一列排序则参考2比特量子线路的第一列排序,即第1个至第3个为(2,1)、(4,1)、(3,1),且(5,1)位于最后1个,其余(6,1)、(7,1)、(8,1)参考(2,1)、(3,1)、(4,1)的排序为(6,1)、(8,1)、(7,1),最终可得第1列排序为:(2,1)、(4,1)、(3,1)、(6,1)、(8,1)、(7,1)、(5,1)。
以此类推,4比特量子线路的第1列排序参考3比特量子线路的第1列排序,为:(2,1)、(4,1)、(3,1)、(6,1)、(8,1)、(7,1)、(5,1)、(10,1)、(12,1)、(11,1)、(14,1)、(16,1)、(15,1)、(13,1)、(9,1),同理可得更多比特量子线路的第1列排序。
然后,继续以2比特量子线路的4阶酉矩阵为例,确定第2列的排序:
获取第一列中与第2列(3,2)、(4,2)同行的元素(3,1)、(4,1)的排序,即(4,1)、(3,1),对应的行的二进制表示为11和00,分别与第2列对应的二进制表示01进行异或操作:
11 ⊕ 01 = 10 =(3,2)
10 ⊕ 01 = 11 =(4,2)
可得,2比特量子线路的4阶酉矩阵第2列排序为:(3,2)、(4,2)。
确定第3列至第4列的排序:类比1比特量子线路的酉矩阵排序,为:(4,3),最终得到的排序如表2所示。
表2 另一种2比特量子线路的酉矩阵排序
同理以3比特量子线路为例,先确定第2列至第4列的排序:
第1列中与第2列同行的非对角元素为的排序为:(4,1)、(3,1)、(6,1)、(8,1)、(7,1)、(5,1),对应的行的二进制分别与第2列的二进制进行异或,由2比特量子线路可知(3,2)、(4,2)的排序不变,在此可省略异或操作,即:
101 ⊕ 001 = 100 =(5,2)
111 ⊕ 001 = 110 =(7,2)
110 ⊕ 001 = 111 =(8,2)
100 ⊕ 001 = 101 =(6,2)
可得,3比特量子线路的8阶酉矩阵第2列排序为:(3,2)、(4,2)、(5,2)、(7,2)、(8,2)、(6,2);
第1列中与第3列同行的非对角元素为的排序为:(4,1)、(6,1)、(8,1)、(7,1)、(5,1),对应的行的二进制分别与第3列的二进制进行异或,由2比特量子线路可知(4,3)的排序不变,在此可省略异或操作,即:
101 ⊕ 010 = 110 =(8,3)
111 ⊕ 010 = 101 =(6,3)
110 ⊕ 010 = 100 =(5,3)
100 ⊕ 010 = 110 =(7,3)
可得,3比特量子线路的8阶酉矩阵第3列排序为:(4,3)、(8,3)、(6,3)、(5,3)、(7,3);
第1列中与第4列同行的非对角元素为的排序为:(6,1)、(8,1)、(7,1)、(5,1),对应的行的二进制分别与第4列的二进制进行异或操作,即:
101 ⊕ 011 = 110 =(7,4)
111 ⊕ 011 = 100 =(5,4)
110 ⊕ 011 = 101 =(6,4)
100 ⊕ 011 = 111 =(8,4)
可得,3比特量子线路的8阶酉矩阵第4列排序为:(7,4)、(5,4)、(6,4)、(8,4)。
接着,对于3比特量子线路的8阶酉矩阵第5列至第8列排序,则类比2比特量子线路的4阶酉矩阵第1列至第4列排序,可得:
第5列排序:(6,5)、(8,5)、(7,5);
第6列排序:(7,6)、(8,6);
第7列排序:(8,7);
第8列排序:无。
同理可确定更多比特量子线路的酉矩阵第2列至最后一列排序。由上述,一种3比特量子线路的酉矩阵的部分列排序如表3所示。
表3 一种3比特量子线路的酉矩阵的部分列排序
S2032,针对所述排序中的第i个非对角元素,构造特定量子逻辑门的阶酉矩阵,以使矩阵中与该非对角元素同位置的元素置为0,且不改变已置0的非对角元素。并且,当第i个非对角元素的排序为所在列最后一个时,同时使矩阵的相同列中的对角元素置为1。
为了方便区别,携带受控信息的单量子逻辑门也可理解为特定量子逻辑门,因其酉矩阵不再是普通意义上的单量子逻辑门的2阶酉矩阵,而是阶酉矩阵,该特定量子逻辑门的表示形式可以为:
其中,Cm表示0、1、*或单量子逻辑门V, m表示量子比特位,,并且,有且只有一个Cm表示单量子逻辑门V。该单量子逻辑门V为普通意义上的操作一个量子比特的单量子逻辑门,只不过可额外受量子线路中的其余量子比特控制。对于不同项的待置0的非对角元素,构造的特定量子逻辑门可能不同。
当Cm为0时,表示量子线路运行到单量子逻辑门V之前(即下一步即将执行逻辑门V),当判断出该位的量子比特的量子态为0态时,才执行单量子逻辑门V,简称0控;
当Cm为1时,表示量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为1态时,执行单量子逻辑门V,简称1控;
当Cm为*时,表示量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为任意态时,均执行单量子逻辑门V,简称无控。
例如,一种特定量子逻辑门的表示形式为{10*V},表示单量子逻辑门V作用在最低位量子比特上,*表示逻辑门V不受第2位的量子比特控制(无控),0表示逻辑门V受第3位的量子比特控制(0控),1表示逻辑门V受第4位的量子比特控制(1控)。另外可知,量子线路为4比特量子线路,该特定量子逻辑门的酉矩阵为阶酉矩阵。
具体的,i为正整数,取值范围为:。当时,单量子逻辑门V的酉矩阵由量子线路的酉矩阵的元素确定;当时,单量子逻辑门V的酉矩阵由矩阵的元素确定。
示例性的,对于2比特量子线路,在表2基础上,对应增加特定量子逻辑门的表示形式,可如表4所示,,,,,,。
表4 一种2比特量子线路对应的特定量子逻辑门
其中,具体矩阵形式如下:
其中,。可以补充的是,2比特量子线路相关的和的矩阵形式如下:
一种特定量子逻辑门、、、位于量子线路中的示意图可如图3所示,空心点及与V的连线表示0控,实心点及与V的连线表示1控,上方横线表示低位量子比特的时间线,下方横线表示高一位的量子比特时间线。
假设2比特量子线路的酉矩阵如下:
第一步,将同位置的非对角元素置0:
在1比特量子线路情况中,由矩阵乘法可知,V的元素确定与(1,1)项和(2,1)项相关,可看成利用(1,1)项消除(2,1)项,以使矩阵相乘后(2,1)项为0。以类似方式,首先也利用(1,1)项消除(2,1)项,中各元素、、、根据(1,1)项和(2,1)项确定:
第二步,将同位置的非对角元素置0:
以(1,1)项消除(2,1)项类比,左下半部分,利用中的(3,1)项消除(4,1)项,据此确定:
第三步,将同位置的非对角元素置0,同时将该列对角元素置1:
利用中的(1,1)项消除(3,1)项,据此确定:
由于和均为酉矩阵,酉矩阵与其本身的转置共轭的乘积为单位阵,通过该条件,可计算得到:
以此类推,对于第2列,类比消除(4,1)项和(3,1)项的方式来消除(3,2)项和(4,2)项,同时不改变第一列中的元素为0的项。之后,矩阵的形式为二阶单位矩阵与二阶矩阵直和的形式,该二阶矩阵可看成1比特量子线路的情形,采用特定量子逻辑门来处理,这样已置0的前两列就不会受到影响。
对于3比特量子线路的酉矩阵,对于第一列,消除(2,1)、(4,1)、(3,1)项的方式与2比特量子线路的情况一样,只是所使用的特定量子逻辑门不同;对于下半部分,消除(6,1)、(8,1)、(7,1)项可类比为消除(2,1)、(4,1)、(3,1)项,最后利用(1,1)项消除(5,1)项。对于第二列,消除(3,2)、(4,2)项与2比特量子线路的情况一样,类比消除(6,1)、(8,1)、(7,1)、(5,1)项的方式来消除(5,2)、(7,2)、(8,2)、(6,2)项。其余列同理。
更具体的,利用元素为a的项消除元素为b的项,如果a项位置在b项上方,则:
否则,在a项位于b项下方的情况下:
其中,、表示a、b的共轭。
示例性的,对于3比特量子线路,在表3基础上,对应特定量子逻辑门的表示形式如下表5所示:
表5 一种3比特量子线路对应的特定量子逻辑门
第一列,(2,1):;(4,1):;(3,1):;(6,1):;(8,1):;(7,1):;(5,1):;
第二列,(3,2):;(4,2):;(5,2):;(7,2):;(8,2):;(6,2):;
第三列,(4,3):;(8,3):;(6,3):;(5,3):;(7,3):;
第四列,(7,4):;(5,4):;(6,4):;(8,4):;
第五列,(6,5):;(8,5):;(7,5):;
第六列,(7,6):;(8,6):;
第七列,(8,7):;第八列无。
本领域技术人员可以理解的是,待置0的非对角元素的排序及特定量子逻辑门的阶酉矩阵的构造不限以上方式,具体以实现为准。
矩阵构造存在一些基本规律。例如,2比特量子线路,对原酉矩阵的行列根据对应的量子比特数进行二进制编码(前述的二进制表示),即从00至11,一个矩阵作用在4阶酉矩阵的左侧,则只会影响的00和01部分(即前两行前两列),同理可以得到只会影响的10和11部分,只会影响的00和10部分,只会影响的01和11部分。对于和矩阵,其中不包含任何控制,由其矩阵形式可知,左乘原矩阵时会影响原矩阵的所有行和列。
特定量子逻辑门的矩阵表示形式的构造规律可概括如下:
一、首先说明量子线路酉矩阵第一列对应的矩阵构造:
1,一比特量子线路:
线路酉矩阵仅有一个元素(2,1)待置0,构造特定量子逻辑门即可,以使;
2,两比特量子线路:
采用递归思想,参考1比特量子线路,线路酉矩阵除最后一个待置0元素(3,1)外,对应的特定量子逻辑门;
对于酉矩阵上半部分(2,1),最高位量子比特设为不受控,即(2,1):;
对于下半部分(4,1),判断低位量子比特对应的是否为1,若不为1,则(4,1):,否则;判断可得:
(4,1)对应1比特量子线路的(2,1):;
最后一个待置0元素(3,1)直接设为:;
3,三比特量子线路:
对应构造的特定量子逻辑门,线路酉矩阵上半部分参考2比特量子线路,仍将最高位量子比特设为不受控,即,可得:
(2,1)对应2比特量子线路的(2,1):;
(4,1)对应2比特量子线路的(4,1):;
(3,1)对应2比特量子线路的(3,1):;
对于下半部分除最后一个待置0元素(5,1)外,与上半部分按顺序一一对应,判断上半部分的、低2位量子比特对应的、是否均不为1,若均不为1,则,否则;判断可得:
(6,1)对应的中,、与(2,1)对应的、相同,即为*、V,且均不为1,可得:;
同理,(8,1)对应(4,1):;(7,1)对应(3,1):;
最后一个待置0元素(5,1)直接设为:;
以此类推,可实现任意比特量子线路酉矩阵的第一列对应的矩阵构造;
二、量子线路酉矩阵第二列至第N/2列对应的矩阵构造:
1,两比特量子线路,n=2:
第2列,列下标,二进制表示01,二进制低位、高位;根据预设不等式,求得x=1;下半部分按顺序与前一列下半部分相对应,矩阵构造如下:
(3,2):参考(4,1)对应的:若j=n,且中均不为1,则(3,2)对应的中的;若,且中对应的,则(3,2)对应的中的;否则,(3,2)对应的与中对应的保持一致;判断可得:
j=1时,满足与中对应的保持一致的条件,即;
j=2时,满足与中对应的保持一致的条件,即;
可得,(3,2)对应的;
(4,2):其为本列最后一个待置0元素,参考第一列(3,1)对应的:将中的*视为0,进行二进制加1操作,*变为1,得到(3,2)对应的;
2,三比特量子线路,n=3:
第2列,列下标,二进制表示01,;根据,求得x=1,上半部分(3,2)、(4,2)参考两比特量子线路:
(3,2)对应的中,取值与两比特量子线路的(3,2)对应的相同,设为*,即:(3,2)对应的;
(4,2)对应的中,取值与两比特量子线路的(4,2)对应的相同,设为*,即:(3,2)对应的;
下半部分按顺序与第一列下半部分进行对应,矩阵构造如下:
(5,2):参考(6,1)对应的:若j=n,且中均不为1,则(5,2)对应的中的;若,且中对应的,则(5,2)对应的中的;否则,(5,2)对应中的与中对应的保持一致;判断可得:
j=1时,满足与中对应的保持一致的条件,即;
j=2时,满足与中对应的保持一致的条件,即等于*;
j=3时,满足与中对应的保持一致的条件,即;
可得,(5,2)对应的;
同理可得,(7,2)对应的;(8,2)对应的;
(6,2):其为本列最后一个待置0元素,参考第一列(5,1)对应的:将中的*视为0,进行二进制加1操作,00变为01,即**变为*1,得到(6,2)对应的;
同理可得,第3列:
上半部分:(4,3)对应;下半部分:(8,3)对应,(6,3)对应,(5,3)对应;本列最后一个待置0元素(7,3)对应;
第4列不进行赘述;可见,下半部分除每列最后一个待置0元素外,偶数列与前一列(奇数列)的矩阵构造对应相同,奇数列的矩阵则参考第一列确定;
三、量子线路酉矩阵第(N/2+1)列至最后一列对应的矩阵构造:
参考第1列至第N/2列的上半部分,按顺序一一对应,将最高位的*变为1,其余不变,以上述3比特量子线路为例,可得:
第5列:(6,5)对应(2,1),可得;(8,5)对应(4,1),可得;(7,5)对应(3,1),可得;
第6列:(7,6)对应(3,2),可得;(8,6)对应(4,2),可得;
第7列:(8,7)对应(4,3),可得;第8列无;
以此类推,可实现任意比特量子线路酉矩阵的所有列对应的矩阵构造,在此不对其进行赘述。
具体的,,其中,等于:
,若;,若;,若;,若为*。
S204,输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
具体的,由,可得:,、为、的转置共轭,即分解出的个携带受控信息的单量子逻辑门(特定量子逻辑门)是处于转置共轭dagger状态的。
在特定量子逻辑门的矩阵形式确定后,特定量子逻辑门随即确定(例如图3所示的特定量子逻辑门位于量子线路中的示意图),根据从依次至的执行时序,构建并输出分解后的、包括至的子量子线路。相比动辄包括成百上千量子逻辑门且存在较大数量的多比特量子逻辑门的复杂量子线路,子量子线路结构大为精简,运行量子线路时的计算复杂度和资源占用显著下降。
具体的,根据子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路,可以包括:
S2041,根据子量子线路,构建相位估计运算对应的第一部分量子线路,以将|b>在矩阵A的特征空间上分解为:,并将第一量子比特和第二量子比特的初态|0>|b>转化为:;
本领域技术人员可以理解的是,相位估计 (Phase estimation) 是量子傅里叶变换 QFT 的一个重要应用, 它的重要性体现在它是很多量子算法的基础。如图4所示,相位估计对应的第一部分量子线路,可以包括:H门操作模块、操作(受控U算子操作)模块和量子傅里叶逆变换模块,其中,该U算子即为矩阵A对应的酉矩阵U,前述子量子线路即为该U算子对应的量子线路实现,在此不对其进行赘述。经过第一部分量子线路,辅助量子比特(对应图4最上的时间线)的量子态不变、第一量子比特(对应图4中间的时间线)的初态|0>转化为、第二量子比特(对应图4最下的时间线)的初态分解。
其中,第一量子比特的数量k取决于相位估计的精度和成功概率,第二量子比特的数量优选为n(至少为n,也可大于n),为矩阵A的特征向量,为的振幅,为矩阵A的特征值。
实际上,输出的是估计值,可以通过增加第一量子比特的数量,提高相位估计的输出精度。并且,在实际应用中,可以设置辅助量子寄存器、第一量子寄存器、第二量子寄存器,分别存储辅助量子比特、第一量子比特、第二量子比特的量子态。
S2042,构建受控旋转操作对应的第二部分量子线路,以将基态中的值提取到辅助量子比特的量子态振幅上,得到:,其中,辅助量子比特的数量为1,所述C为常数,;
具体的,受控旋转又可称“提取占比”,因为经过相位估计运算之后,第一量子寄存器会存储一系列的特征值(具体存储在基态中),而第二量子寄存器存储的输入态即初态会在矩阵A的特征空间上进行分解,然后通过受控旋转操作,将基态中的值提取到了振幅上,将辅助量子比特的量子态|0>转化为得到,各量子比特的量子态经过第二部分量子线路由转化为。为了减少资源占用,辅助量子比特可设1位,且C为常数,一般取1。
S2043,构建相位估计逆运算对应的第三部分量子线路,以消除,得到:;
本领域技术人员可以理解的是,相位估计逆运算是前述相位估计的还原过程,或称相位估计的转置共轭运算,目的是消除,具体以将量子态转化为如下为准:
S2044,构建针对所述辅助量子比特的量子测量操作,以使测量到所述辅助量子比特的量子态为时,得到:,所述与为振幅归一化的对应关系;
具体的,量子测量操作施加在辅助量子比特上,以在相位估计逆运算后,对辅助量子比特进行测量。测量后,辅助量子比特的状态会塌缩至一个确定的状态,其中,塌缩至|0>的概率是,塌缩至|1>的概率是。当测量辅助量子比特的量子态为、且时,即可得到确定的量子态:,可见是进行振幅归一化的对应结果。在实际应用中,可以根据用户所需的应用场景,对应得到,或者直接将用于后续的场景计算。
S2045,将所述第一部分量子线路、所述第二部分量子线路、所述第三部分量子线路和所述量子测量操作,依次组成HHL算法对应的量子线路。
具体的,根据第一部分量子线路、第二部分量子线路、第三部分量子线路和量子测量操作的执行时序,依次组成一条完整的量子线路,即为HHL算法对应的总量子线路。如前述,第一部分、第二部分、第三部分及总量子线路包含的量子比特总数均相同,可以为1+k+n。
可见,由于输出的子量子线路中量子逻辑门数量有限,且消除了酉矩阵形式复杂的多比特量子逻辑门,使得量子逻辑门形式得以简化,从而减少构建的HHL算法的量子线路复杂度,提高量子线路的模拟效率,同时减少硬件资源的占用。
参见图5,图5为本发明实施例提供的一种量子线路的构建装置的结构示意图,与图2所示的流程相对应,可以包括:
第一获得模块501,用于获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
第二获得模块502,用于获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
分解模块503,用于确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,所述,所述为阶单位矩阵;
构建模块504,用于输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
具体的,所述分解模块,具体用于:
若所述矩阵A为酉矩阵,则将所述矩阵A直接确定为对应的酉矩阵U;
若所述矩阵A为厄米矩阵,确定对应的酉矩阵;其中,所述t为常量;
若所述矩阵A为非厄米矩阵且非酉矩阵,确定对应厄米矩阵,以确定对应的酉矩阵,同时将所述向量b转换为。
具体的,所述分解模块,包括:
确定单元,用于确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序;
构造单元,用于针对所述排序中的第i个非对角元素,构造特定量子逻辑门的阶酉矩阵,以使矩阵中与该非对角元素同位置的元素置为0,且不改变已置0的非对角元素;
其中,所述特定量子逻辑门包括操作一比特的单量子逻辑门,所述单量子逻辑门携带受其余比特控制的受控信息,; 当所述时,所述单量子逻辑门的酉矩阵由所述量子线路对应的酉矩阵的元素确定;当时,所述单量子逻辑门的酉矩阵由矩阵的元素确定;并且,当所述第i个非对角元素的排序为所在列最后一个时,同时使矩阵的相同列中的对角元素置为1。
具体的,所述确定单元,具体用于:
当n=1时,所述酉矩阵U中对角元素下方的、待置0的非对角元素排序为(2,1);其中,所述(2,1)表示坐标为第2行第1列的非对角元素;
当n>1时,根据(n-1)比特量子线路对应的酉矩阵的第一列排序,确定n比特量子线路对应的所述酉矩阵U中对角元素下方的、待置0的非对角元素的第一列排序;其中,第一列中坐标(N/2+1,1)的非对角元素的排序位于第一列最后一个;
基于所述n比特量子线路对应的第一列排序,分别确定所述n比特量子线路对应的第2列至第N/2列中对角元素下方的、待置0的非对角元素的排序;
根据(n-1)比特量子线路对应的酉矩阵中对角元素下方的、待置0的非对角元素的排序,对应确定所述n比特量子线路对应的第(N/2+1)列至第N列中对角元素下方的、待置0的非对角元素的排序。
具体的,所述携带受控信息的单量子逻辑门的表示形式包括:
,其中,所述Cm表示0、1、*或单量子逻辑门V,所述m表示量子比特位,,并且,有且只有一个Cm表示单量子逻辑门V,所述单量子逻辑门V的酉矩阵由所述酉矩阵U确定;
当Cm为0时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为0态时,执行单量子逻辑门V;
当Cm为1时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为1态时,执行单量子逻辑门V;
当Cm为*时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为任意态时,均执行单量子逻辑门V。
具体的,所述构建模块,具体用于:
根据所述子量子线路,构建相位估计运算对应的第一部分量子线路,以将|b>在所述矩阵A的特征空间上分解为:,并将所述第一量子比特和所述第二量子比特的初态|0>|b>转化为:;其中,所述第一量子比特的数量k取决于相位估计的精度和成功概率,所述第二量子比特的数量为所述n,所述为所述矩阵A的特征向量,所述为的振幅,所述为所述矩阵A的特征值;
构建受控旋转操作对应的第二部分量子线路,以将基态中的值提取到所述辅助量子比特的量子态振幅上,得到:,其中,所述辅助量子比特的数量为1,所述C为常数,;
构建相位估计逆运算对应的第三部分量子线路,以消除,得到:;
构建针对所述辅助量子比特的量子测量操作,以使测量到所述辅助量子比特的量子态为时,得到:,所述与为振幅归一化的对应关系;
将所述第一部分量子线路、所述第二部分量子线路、所述第三部分量子线路和所述量子测量操作,依次组成HHL算法对应的量子线路。
可见,由于输出的子量子线路中量子逻辑门数量有限,且消除了酉矩阵形式复杂的多比特量子逻辑门,使得量子逻辑门形式得以简化,从而减少构建的HHL算法的量子线路复杂度,提高量子线路的模拟效率,同时减少硬件资源的占用。
本发明实施例还提供了一种存储介质,所述存储介质中存储有计算机程序,其中,所述计算机程序被设置为运行时执行上述任一项中方法实施例中的步骤。
具体的,在本实施例中,上述存储介质可以被设置为存储用于执行以下步骤的计算机程序:
S1,获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
S2,获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
S3,确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,所述,所述为阶单位矩阵;
S4,输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
具体的,在本实施例中,上述存储介质可以包括但不限于:U盘、只读存储器(Read-Only Memory,简称为ROM)、随机存取存储器(Random Access Memory,简称为RAM)、移动硬盘、磁碟或者光盘等各种可以存储计算机程序的介质。
可见,由于输出的子量子线路中量子逻辑门数量有限,且消除了酉矩阵形式复杂的多比特量子逻辑门,使得量子逻辑门形式得以简化,从而减少构建的HHL算法的量子线路复杂度,提高量子线路的模拟效率,同时减少硬件资源的占用。
本发明实施例还提供了一种电子装置,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行上述任一项中方法实施例中的步骤。
具体的,上述电子装置还可以包括传输设备以及输入输出设备,其中,该传输设备和上述处理器连接,该输入输出设备和上述处理器连接。
具体的,在本实施例中,上述处理器可以被设置为通过计算机程序执行以下步骤:
S1,获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
S2,获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
S3,确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,,所述,所述为阶单位矩阵;
S4,输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
可见,由于输出的子量子线路中量子逻辑门数量有限,且消除了酉矩阵形式复杂的多比特量子逻辑门,使得量子逻辑门形式得以简化,从而减少构建的HHL算法的量子线路复杂度,提高量子线路的模拟效率,同时减少硬件资源的占用。
以上依据图式所示的实施例详细说明了本发明的构造、特征及作用效果,以上所述仅为本发明的较佳实施例,但本发明不以图面所示限定实施范围,凡是依照本发明的构想所作的改变,或修改为等同变化的等效实施例,仍未超出说明书与图示所涵盖的精神时,均应在本发明的保护范围内。
Claims (10)
1.一种量子线路的构建方法,其特征在于,包括:
获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n,所述n为正整数;
获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,且满足矩阵中与第i个非对角元素同位置的元素置为0,且不改变已置0的非对角元素,,所述,所述为阶单位矩阵;
输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
2.根据权利要求1所述的方法,其特征在于,所述确定所述矩阵A对应的酉矩阵U,包括:
若所述矩阵A为酉矩阵,则将所述矩阵A直接确定为对应的酉矩阵U;
若所述矩阵A为厄米矩阵,确定对应的酉矩阵;其中,所述t为常量;
若所述矩阵A为非厄米矩阵且非酉矩阵,确定对应厄米矩阵,以确定对应的酉矩阵,同时将所述向量b转换为。
3.根据权利要求1所述的方法,其特征在于,将所述酉矩阵分解成个携带受控信息的单量子逻辑门对应的酉矩阵,包括:
确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序;
针对所述排序中的第i个非对角元素,构造特定量子逻辑门的阶酉矩阵,以使矩阵中与该非对角元素同位置的元素置为0,且不改变已置0的非对角元素;
其中,所述特定量子逻辑门包括操作一比特的单量子逻辑门,所述单量子逻辑门携带受其余比特控制的受控信息,; 当所述时,所述单量子逻辑门的酉矩阵由所述量子线路对应的酉矩阵的元素确定;当时,所述单量子逻辑门的酉矩阵由矩阵的元素确定;并且,当所述第i个非对角元素的排序为所在列最后一个时,同时使矩阵的相同列中的对角元素置为1。
4.根据权利要求3所述的方法,其特征在于,所述确定所述酉矩阵中对角元素下方的、待置0的非对角元素的排序,包括:
当n=1时,所述酉矩阵U中对角元素下方的、待置0的非对角元素排序为(2,1);其中,所述(2,1)表示坐标为第2行第1列的非对角元素;
当n>1时,根据(n-1)比特量子线路对应的酉矩阵的第一列排序,确定n比特量子线路对应的所述酉矩阵U中对角元素下方的、待置0的非对角元素的第一列排序;其中,第一列中坐标(N/2+1,1)的非对角元素的排序位于第一列最后一个;
基于所述n比特量子线路对应的第一列排序,分别确定所述n比特量子线路对应的第2列至第N/2列中对角元素下方的、待置0的非对角元素的排序;
根据(n-1)比特量子线路对应的酉矩阵中对角元素下方的、待置0的非对角元素的排序,对应确定所述n比特量子线路对应的第(N/2+1)列至第N列中对角元素下方的、待置0的非对角元素的排序。
5.根据权利要求1所述的方法,其特征在于,所述携带受控信息的单量子逻辑门的表示形式包括:
,其中,所述Cm表示0、1、*或单量子逻辑门V,所述m表示量子比特位,,并且,有且只有一个Cm表示单量子逻辑门V,所述单量子逻辑门V的酉矩阵由所述酉矩阵U确定;
当Cm为0时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为0态时,执行单量子逻辑门V;
当Cm为1时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为1态时,执行单量子逻辑门V;
当Cm为*时,表示所述量子线路运行到单量子逻辑门V之前,当该位的量子比特的量子态为任意态时,均执行单量子逻辑门V。
6.根据权利要求1所述的方法,其特征在于,所述根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路,包括:
根据所述子量子线路,构建相位估计运算对应的第一部分量子线路,以将|b>在所述矩阵A的特征空间上分解为:,并将所述第一量子比特和所述第二量子比特的初态|0>|b>转化为:;其中,所述第一量子比特的数量k取决于相位估计的精度和成功概率,所述第二量子比特的数量为所述n,所述为所述矩阵A的特征向量,所述为的振幅,所述为所述矩阵A的特征值;
构建受控旋转操作对应的第二部分量子线路,以将基态中的值提取到辅助量子比特的量子态振幅上,得到:,其中,所述辅助量子比特的数量为1,所述C为常数,;
构建相位估计逆运算对应的第三部分量子线路,以消除,得到:;
构建针对所述辅助量子比特的量子测量操作,以使测量到所述辅助量子比特的量子态为时,得到:,所述与为振幅归一化的对应关系;
将所述第一部分量子线路、所述第二部分量子线路、所述第三部分量子线路和所述量子测量操作,依次组成HHL算法对应的量子线路。
7.一种量子线路的构建装置,其特征在于,包括:
第一获得模块,用于获得N*N维矩阵A和N维向量b,其中,所述矩阵A为可逆矩阵,所述N=2n, 所述n为正整数;
第二获得模块,用于获得包括辅助量子比特、第一量子比特、第二量子比特的若干量子比特,其中,所述辅助量子比特和所述第一量子比特的初态置为|0>,所述第二量子比特的初态置为,所述为所述向量b的第j个元素;
分解模块,用于确定所述矩阵A对应的酉矩阵U,将所述矩阵U分解成个携带受控信息的单量子逻辑门对应的酉矩阵;其中,满足,所述为第i个携带受控信息的单量子逻辑门对应的酉矩阵,且满足矩阵中与第i个非对角元素同位置的元素置为0,且不改变已置0的非对角元素,,所述,所述为阶单位矩阵;
构建模块,用于输出包含所述个携带受控信息的单量子逻辑门的子量子线路,根据所述子量子线路和各量子比特及其初态,构建HHL算法对应的量子线路。
8.根据权利要求7所述的装置,其特征在于,所述分解模块,具体用于:
若所述矩阵A为酉矩阵,则将所述矩阵A直接确定为对应的酉矩阵U;
若所述矩阵A为厄米矩阵,确定对应的酉矩阵;其中,所述t为常量;
若所述矩阵A为非厄米矩阵且非酉矩阵,确定对应厄米矩阵,以确定对应的酉矩阵,同时将所述向量b转换为。
9.一种存储介质,其特征在于,所述存储介质中存储有计算机程序,其中,所述计算机程序被设置为运行时执行所述权利要求1至6任一项中所述的方法。
10.一种电子装置,包括存储器和处理器,其特征在于,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行所述权利要求1至6任一项中所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010598078.3A CN113850389B (zh) | 2020-06-28 | 2020-06-28 | 一种量子线路的构建方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010598078.3A CN113850389B (zh) | 2020-06-28 | 2020-06-28 | 一种量子线路的构建方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113850389A CN113850389A (zh) | 2021-12-28 |
CN113850389B true CN113850389B (zh) | 2023-08-11 |
Family
ID=78972167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010598078.3A Active CN113850389B (zh) | 2020-06-28 | 2020-06-28 | 一种量子线路的构建方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113850389B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114219048A (zh) * | 2022-02-21 | 2022-03-22 | 合肥本源量子计算科技有限责任公司 | 基于量子计算的谱聚类方法、装置、电子设备及存储介质 |
WO2023169345A1 (zh) * | 2022-03-11 | 2023-09-14 | 本源量子计算科技(合肥)股份有限公司 | 数据模拟任务的处理方法、装置、电子设备及存储介质 |
CN117010506A (zh) * | 2022-04-29 | 2023-11-07 | 腾讯科技(深圳)有限公司 | 量子态制备电路的生成方法、装置、量子芯片和电子设备 |
CN114819166B (zh) * | 2022-05-27 | 2023-05-30 | 北京大学 | 一种量子系统的演化方法及装置 |
CN117852659A (zh) * | 2022-09-30 | 2024-04-09 | 本源量子计算科技(合肥)股份有限公司 | 量子线路生成方法、装置、存储介质及电子装置 |
CN116681138B (zh) * | 2023-05-29 | 2024-06-14 | 本源量子计算科技(合肥)股份有限公司 | 一种数据大小比较任务的处理方法、装置、存储介质及电子装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005332242A (ja) * | 2004-05-20 | 2005-12-02 | Nippon Telegr & Teleph Corp <Ntt> | ユニタリ行列分解方法、ユニタリ行列分解装置、ユニタリ行列分解プログラム及び記録媒体 |
CN101118608A (zh) * | 2007-08-23 | 2008-02-06 | 清华大学 | 任意量子比特门的分解方法 |
CA3133427A1 (en) * | 2017-10-18 | 2019-04-25 | Google Llc | Simulation of quantum circuits |
CN110826719A (zh) * | 2019-10-14 | 2020-02-21 | 合肥本源量子计算科技有限责任公司 | 一种量子程序的处理方法、装置、存储介质和电子装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060123363A1 (en) * | 2004-12-07 | 2006-06-08 | Williams Colin P | Method and apparatus for automated design of quantum circuits |
US10242321B2 (en) * | 2014-04-09 | 2019-03-26 | Microsoft Technology Licensing, Llc | Efficient synthesis of repeat-until-success circuits in clifford + T basis |
-
2020
- 2020-06-28 CN CN202010598078.3A patent/CN113850389B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005332242A (ja) * | 2004-05-20 | 2005-12-02 | Nippon Telegr & Teleph Corp <Ntt> | ユニタリ行列分解方法、ユニタリ行列分解装置、ユニタリ行列分解プログラム及び記録媒体 |
CN101118608A (zh) * | 2007-08-23 | 2008-02-06 | 清华大学 | 任意量子比特门的分解方法 |
CA3133427A1 (en) * | 2017-10-18 | 2019-04-25 | Google Llc | Simulation of quantum circuits |
CN110826719A (zh) * | 2019-10-14 | 2020-02-21 | 合肥本源量子计算科技有限责任公司 | 一种量子程序的处理方法、装置、存储介质和电子装置 |
Non-Patent Citations (1)
Title |
---|
MAHBOOBEH HOUSHMAND etal..Quantum Circuit Synthesis Targeting to Improve One-Way Quantum Computation Pattern Cost Metrics.《ACM Journal on Emerging Technologies in Computing Systems》.2017,第第13卷卷(第第4期期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113850389A (zh) | 2021-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113850389B (zh) | 一种量子线路的构建方法及装置 | |
CN111563599B (zh) | 一种量子线路的分解方法、装置、存储介质及电子装置 | |
CN112114776B (zh) | 一种量子乘法运算方法、装置、电子装置及存储介质 | |
CN112073221B (zh) | 一种实现网络节点排序的方法及装置 | |
CN114418105B (zh) | 一种基于量子线路处理量子应用问题的方法及装置 | |
CN112068798B (zh) | 一种实现网络节点重要性排序的方法及装置 | |
CN114764549B (zh) | 基于矩阵乘积态的量子线路模拟计算方法、装置 | |
CN113222150B (zh) | 一种量子态的变换方法及装置 | |
CN114819163B (zh) | 量子生成对抗网络的训练方法、装置、介质及电子装置 | |
CN113222153A (zh) | 一种量子态的模拟方法、装置、存储介质和电子装置 | |
CN113222151B (zh) | 一种量子态的变换方法及装置 | |
CN116431807B (zh) | 一种文本分类方法、装置、存储介质及电子装置 | |
CN114418104B (zh) | 一种量子应用问题的处理方法及装置 | |
CN115346080B (zh) | 基于量子计算的图像处理方法及相关设备 | |
CN114372539B (zh) | 基于机器学习框架的分类方法及相关设备 | |
CN115809707B (zh) | 量子比较运算方法、装置、电子装置及基础算术组件 | |
CN114819167B (zh) | 一种稀疏线性系统的稀疏近似逆量子预处理方法及装置 | |
CN114881238A (zh) | 量子鉴别器的构造方法、装置、介质及电子装置 | |
CN114764618B (zh) | 一种针对线性系统的量子预处理方法及装置 | |
CN114692879B (zh) | 一种基于稀疏线性系统的量子预处理方法及装置 | |
CN115713122B (zh) | 一种量子数据与经典数据的大小关系的确定方法及装置 | |
CN115879557B (zh) | 基于量子线路的数据大小比较方法、装置及量子计算机 | |
CN116432721B (zh) | 数据处理方法、机器学习框架及相关设备 | |
CN116499466B (zh) | 一种智能体导航方法、装置、存储介质及电子装置 | |
CN116090571A (zh) | 基于广义极小残量的量子线性求解方法、装置及介质 |
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 | ||
CB02 | Change of applicant information |
Address after: 230088 6th floor, E2 building, phase II, innovation industrial park, 2800 innovation Avenue, Hefei high tech Zone, Hefei City, Anhui Province Applicant after: Benyuan Quantum Computing Technology (Hefei) Co.,Ltd. Address before: 230088 6th floor, E2 building, phase II, innovation industrial park, 2800 innovation Avenue, Hefei high tech Zone, Hefei City, Anhui Province Applicant before: ORIGIN QUANTUM COMPUTING COMPANY, LIMITED, HEFEI |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |