CN113190998B - 一种波茨模型的高维复用计算方法及装置 - Google Patents

一种波茨模型的高维复用计算方法及装置 Download PDF

Info

Publication number
CN113190998B
CN113190998B CN202110477851.5A CN202110477851A CN113190998B CN 113190998 B CN113190998 B CN 113190998B CN 202110477851 A CN202110477851 A CN 202110477851A CN 113190998 B CN113190998 B CN 113190998B
Authority
CN
China
Prior art keywords
node
model
mode
independent
optical
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
Application number
CN202110477851.5A
Other languages
English (en)
Other versions
CN113190998A (zh
Inventor
刘洁
李振华
余思远
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202110477851.5A priority Critical patent/CN113190998B/zh
Publication of CN113190998A publication Critical patent/CN113190998A/zh
Application granted granted Critical
Publication of CN113190998B publication Critical patent/CN113190998B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10—Noise analysis or noise optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Communication System (AREA)

Abstract

本发明涉及波茨模型计算技术领域，公开了一种波茨模型的高维复用计算方法及装置，通过形成光电环路的信息传递，以光场的强度编码节点的状态，实现了波茨模型的简易、高效求解。并且，本发明通过利用了光场的两个正交的偏振态X和Y以及两个正交相位0(I)和π(Q)进行计算通道的复用，实现了XI、XQ、YI、YQ四路独立运行的二值化节点伊辛机，且能够实现可取相位状态数目任意的节点的波茨机，使得算法复杂度得到降低；并且，由于四个独立光场分量形成的四个通道，可根据模型类型选择不同的计算模式，具有较高的灵活性；另外，由于基于正交偏振和相位的光场可以相互叠加而不产生串扰，所以本发明的装置更加紧凑且更加易于实现。

Description

一种波茨模型的高维复用计算方法及装置
技术领域
本发明涉及波茨模型计算技术领域，特别是涉及一种波茨模型的高维复用计算方法及装置。
背景技术
波茨模型(Potts model)是节点取值为二元相位的伊辛模型(Ising model)的拓展，描述了一个由N个具有m种可取相位状态
Figure GDA0003692495050000011
且以一定强度Jij相互耦合的节点组成的数学模型。该模型的哈密顿量可以表示为：
Figure GDA0003692495050000012
在此，σi和σj、φi和φj分别表示第i和第j个节点状态的单位矢量和相位。
而求解该哈密顿量的最小值及其对应的节点状态，是一种典型的NP-难/完全(NP-hard/complete)问题。在传统电子计算机上，这类问题的求解，所需的计算时间或硬件资源随着节点数N的增加而指数增长，所以面对大规模波茨模型时，传统计算机的算力遇到了瓶颈。波茨机作为一种新型的求解器，通过模拟谐振腔受激振荡并稳定在能量最低的模式的过程，可以高效地求解波茨模型的最低哈密顿量及其对应的节点状态。以二值化节点的伊辛机为例，每一个节点的初始值设置为0，在系统的噪声作用、相位敏感增益作用与节点之间相互耦合作用的三重作用下，会趋向于向+1或者-1演化。任一时刻中，所有节点的符号组合(即模式)决定了哈密顿量的大小，也对应着该模式下的系统损耗。而每个节点的的演化方向，将使系统的哈密顿量趋于最低，即：系统将趋于演化到并最终稳定在损耗达到最小的模式上。
近年来，随着伊辛机的实现与成熟，波茨机也受到关注。基于非线性光学和光学参量振荡的波茨机实现了对波茨模型求解的演示，但是存在算法复杂度增加、节点数目规模受限等缺点。
中国发明专利申请CN111373344A(公开日为2020年07月03日)公开了一种波茨模型的计算装置，具备：伊辛模型的计算装置；计算结果保持判定部，保持在所述伊辛模型的计算装置中设定某个耦合系数时得到的伊辛模型的自旋的值，并且判定计算结束；以及耦合系数设定部，针对所述伊辛模型的计算装置，设定基于所述保持的自旋的值生成的耦合系数，所述波茨模型的计算装置的特征在于，所述耦合系数设定部根据在所述伊辛模型的计算装置中作为与第m次设定的耦合系数对应的计算结果而获得的伊辛模型的自旋的值，针对所述伊辛模型计算装置重新生成并设定第m+1次设定的耦合系数，在多值自旋的可取的值为Si＝0、1、2……M－1(M为自然数)且M≤2Ms的情况下，在所述耦合系数设定部中的耦合系数的设定次数到达与Ms(Ms为自然数)对应的次数为止的情况下，所述计算结果判定部判定为计算结束，基于第m次设定的耦合系数计算出将作为计算结果而获得的自旋的值σim代入下述算式而获得的值Si，由此利用伊辛模型来计算映射到波茨模型的问题。该专利仅仅通过光学传递信息，存在算法复杂、节点数目规模受限的问题。
发明内容
本发明的目的是提供一种简易、高效的波茨模型的高维复用计算方法及装置。
为了实现上述目的，本发明提供了一种波茨模型的高维复用计算方法，包括如下步骤：
S1、获取待求解的模型及其模型参数，根据待求解的模型类型选择计算模式；
S2、将激光光源发出的连续激光分为两路，一路作为承载节点状态的载波，另一路作为本振参考光以供四路独立的复用光场节点状态的相干检测；
S3、将步骤S2作为承载节点状态的载波的激光分成XI、XQ、YI、YQ四个独立光场分量，并把四个独立光场分量对应的通道的节点状态信息中的模拟电域信息
Figure GDA0003692495050000021
分别加载到上述四个独立光场分量的强度
Figure GDA0003692495050000022
上，其中，
Figure GDA0003692495050000023
表示第k个循环中的第i节点的模拟电域信息，
Figure GDA0003692495050000024
表示第k个循环中的第i节点的独立光场分量的强度；
S4、通过光路复用系统传播XI、XQ、YI、YQ四个独立光场分量；
S5、接收经过传播的XI、XQ、YI、YQ四个独立光场分量，并将四个独立光场分量与步骤S2作为本振参考光的激光进行对比、干涉，分离成四路独立的光场信号，并转换成包含节点状态信息的四路模拟电信号；
S6、将步骤S5得到的包含节点状态信息的四路模拟电信号转换为数字电域信息，根据步骤S1获得的模型参数和选定的计算模式进行计算并更新节点状态信息，将更新后的节点状态信息转换为模拟电信号反馈到步骤S3中用于加载到四个独立光场分量上，循环执行步骤S3、S4、S5和S6；步骤S1获得的模型参数包括增益系数α、耦合强度系数β和耦合常数Jij，在步骤S6中，根据接收到的节点状态信息、获得的模型参数和选定的计算模式进行计算，包括：
S6.1、节点状态信息的采集与重整:采集每个节点的状态信息，然后根据选择的计算模式重整节点的形式；
S6.2、节点增益计算：根据模型参数α，计算增益项αxi，其中，α为模型的增益系数，xi为第i节点的状态值；
S6.3、节点耦合计算，根据模型参数β、Jij，计算耦合项β∑j≠iJijxj，其中，β为模型的耦合强度系数，Jij为模型耦合常数，xj为第j节点的状态值；
S6.4、节点状态信息反馈信号生成：根据
f0i＝αxi+β∑j≠iJijxj，
求得f0i，f0i为节点状态信息反馈信号，再根据各计算模式下取值空间的要求，将f0i投影到其最接近的矢量方向上，更新节点状态信息反馈信号，生成更新的节点状态信息反馈信号
Figure GDA0003692495050000031
作为优选方案，在步骤S1中，选择计算模式包括选择求解器模式和循环迭代次数，其中，求解器模式包括：
(1)伊辛机模式：该模式下，节点的状态空间是一维的；
(2)波茨机模式：该模式下，节点的状态空间是二维的，可由两路通道的信号重整为复数；
(3)三维节点状态空间模式：该模式下，节点的状态空间是三维的，可由三路通道的信号重整化为一个三元数；
(4)四维节点状态空间模式：该模式下，节点的状态空间是四维的，可由四路通道的信号重整化为一个四元数。
作为优选方案，在步骤S3中，四个独立光场分量通过双偏振IQ电光调制器加载，对于双偏振IQ电光调制器的每一个独立调制通道，其直流偏置点都应在线性区的-π/4线性区工作点，使得对第k个循环中的第i节点在反馈信号
Figure GDA0003692495050000041
的调制后为：
Figure GDA0003692495050000042
其中，外加的常数项
Figure GDA0003692495050000043
使节点状态取值为[-0.5,0.5]，包含正值或者负值，再取
Figure GDA0003692495050000044
符号以实现节点的二值化。
作为优选方案，在步骤S6.1中，在伊辛机模式下，每个伊辛机通道的节点都不需要变换形式；在波茨机模式下，同一偏振态X或Y下的I、Q两路的节点重整合并为一路：xX＝xXI+j·xXQ和xY＝xYI+j·xYQ，其中j是虚数单位；在三维节点状态空间模式和四维节点状态空间模式下，多个通道的数据重整为一路。
作为优选方案，对于波茨机模式、三维节点状态空间模式和四维节点状态空间模式，节点状态信息反馈信号为：
Figure GDA0003692495050000045
其中单位矢量
Figure GDA0003692495050000046
的相位角为
Figure GDA0003692495050000047
且与f0i的相位角最接近。
作为优选方案，在步骤S4中，在XI、XQ、YI、YQ四个独立光场分量的传播过程中，在每个独立光场分量上引入独立的随机噪声。
本发明还提供一种波茨模型的高维复用计算装置，包括：
模型与任务获取模块，用于获取待求解的模型及模型参数，并根据待求解的模型类型选择计算模式；
激光光源，用于发出的连续激光，且将激光分为两路，一路作为承载节点状态的载波，输入到电-光节点更新模块进行处理；另一路作为本振参考光输入到光-电节点测量模块以供四路独立的复用光场节点状态的相干检测；
电-光节点更新模块，用于接收激光光源发出的作为承载节点状态的载波的激光，并将该激光分成XI、XQ、YI、YQ四个独立光场分量，并把节点信息处理模块反馈的四个通道的节点状态信息中的模拟电域信息
Figure GDA0003692495050000048
分别加载到上述四个独立光场分量的强度
Figure GDA0003692495050000051
上，其中，
Figure GDA0003692495050000052
表示第k个循环中的第i节点的模拟电域信息，
Figure GDA0003692495050000053
表示第k个循环中的第i节点的独立光场分量的强度；
光路复用系统，用于连接电-光节点更新模块与光-电节点测量模块，传播XI、XQ、YI、YQ四个独立光场分量；
光-电节点测量模块，用于接收光路复用系统传播的XI、XQ、YI、YQ四个独立光场分量以及激光光源发出的作为本振参考光的激光，且将XI、XQ、YI、YQ四个独立光场分量和本振参考光进行干涉、对比，并分离成四路独立的光场信号，将光场信号转换成包含节点状态信息的四路模拟电信号；
节点信息处理模块，用于接收光-电节点测量模块的包含节点状态信息的四路模拟电信号并转换为数字电域信息且暂时储存，接收模型与任务获取模块的计算模式和模型参数，对节点信息进行计算、更新和暂时储存，将计算更新后的节点信息转换成模拟电信号反馈给电-光节点更新模块；所述节点信息处理模块包括数学运算单元，数学运算单元用于根据接收到的模型与任务获取模块的计算模式和模型参数，对节点信息进行计算、更新，数学运算单元包括：
节点信息采集与重整子单元，用于采集每个节点的状态值，且根据所指定的计算模式，重整节点的形式；
节点增益计算子单元，用于根据模型参数α，计算增益项αxi，其中，α为模型的增益系数，xi为第i节点的状态值；
节点耦合计算子单元，用于模型参数β、Jij，计算耦合项β∑j≠iJijxj，其中，β为模型的耦合强度系数，Jij为模型耦合常数，nj为第j节点的状态值；
节点状态更新反馈信号生成子单元，用于获取增益项αxi和耦合项β∑j≠iJijxj，并计算f0i＝αxi+β∑j≠iJijxj得到节点状态信息反馈信号f0i，再根据各计算模式下取值空间的要求，将f0i投影到其最接近的矢量方向上，更新节点状态信息反馈信号，生成更新的节点状态信息反馈信号
Figure GDA0003692495050000054
作为优选方案，所述电-光节点更新模块包括双偏振IQ电光调制器。
与现有技术相比，本发明的有益效果在于：
本发明通过激光传播节点状态信息，接收后将光场信号转换成模拟电信号，在计算时转换成数据信号，再将计算更新后的反馈信号转换成模拟电信号加载到光场分量上，形成光电环路的信息传递，以光场的强度编码节点的状态，实现了波茨模型的简易、高效求解，且解决了现有的计算装置复杂、维持装置稳定成本高的问题。并且，本发明通过利用了光场的两个正交的偏振态X和Y以及两个正交相位0(I)和π(Q)进行计算通道的复用，实现了XI、XQ、YI、YQ四路独立运行的二值化节点伊辛机，且能够实现可取相位状态数目任意的节点的波茨机，使得算法复杂度得到降低；并且，由于四个独立光场分量形成的四个通道，可根据模型类型选择不同的计算模式，具有较高的灵活性；另外，由于基于正交偏振和相位的光场可以相互叠加而不产生串扰，所以本发明的装置更加紧凑且更加易于实现。
附图说明
图1是本发明实施例的波茨模型的高维复用计算方法的流程图。
图2是本发明实施例的波茨模型的高维复用计算方法步骤6的模型计算更新的流程图。
图3是本发明实施例的波茨模型的高维复用计算装置的原理框图。
图4是本发明实施例的利用的四路光场独立分量叠加示意图。
图5是本发明实施例的利用的X偏振方向光场分量示意图。
图6是本发明实施例的利用的Y偏振方向光场分量示意图。
图7是本发明实施例的利用的I和Q光场分量叠加示意图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
如图1所示，本实施例的波茨模型的高维复用计算方法，包括如下步骤：
S1、获取待求解的模型及其模型参数，根据待求解的模型类型选择计算模式。本实施例获取的模型参数包括增益系数α、耦合强度系数β和耦合常数Jij。选择计算模式包括选择求解器模式和循环迭代次数，其中，求解器模式包括：
(1)伊辛机模式：该模式下，节点的状态空间是一维的；
(2)波茨机模式：该模式下，节点的状态空间是二维的，可由两路通道的信号重整为复数；
(3)三维节点状态空间模式：该模式下，节点的状态空间是三维的，可由三路通道的信号重整化为一个三元数；
(4)四维节点状态空间模式：该模式下，节点的状态空间是四维的，可由四路通道的信号重整化为一个四元数。
S2、将激光光源发出的连续激光分为两路，一路作为承载节点状态的载波，另一路作为本振参考光以供四路独立的复用光场节点状态的相干检测。
S3、将步骤S2作为承载节点状态的载波的激光分成XI、XQ、YI、YQ四个独立光场分量，并把四个独立光场分量对应的通道的节点状态信息中的模拟电域信息
Figure GDA0003692495050000071
分别加载到上述四个独立光场分量的强度
Figure GDA0003692495050000072
上，其中，
Figure GDA0003692495050000073
表示第k个循环中的第i节点的模拟电域信息，
Figure GDA0003692495050000074
表示第k个循环中的第i节点的独立光场分量的强度。
在本实施例中，四个独立光场分量通过双偏振IQ电光调制器加载，对于双偏振IQ电光调制器的每一个独立调制通道，其直流偏置点都应在线性区的-π/4工作点，使得对第k个循环中的第i节点在反馈信号
Figure GDA0003692495050000075
的调制后为：
Figure GDA0003692495050000076
其中，外加的常数项
Figure GDA0003692495050000077
使节点状态取值为[-0.5,0.5]，包含正值或者负值，再取
Figure GDA0003692495050000078
符号以实现节点的二值化。
四个独立光场分量对应的通道的节点状态信息中的模拟电域信息
Figure GDA0003692495050000079
是通过对模型进行计算和更新后反馈得到的，在第一次循环迭代时，所有节点的初值都设置为0。
S4、通过光路复用系统传播XI、XQ、YI、YQ四个独立光场分量。光路复用系统采用现有技术的系统，一般就是对偏振不敏感的光纤，或者其他波导器件、自由空间光等，需要强调的是，这一段光路，可以支持光的高维度传输，而对具体实现形式没有特殊要求。通过光路复用系统传播，可以支持XI、XQ、YI、YQ四个独立光场分量共享同一个物理通道，并且同步且无串扰地从电-光节点更新模块到光-电节点测量模块传输节点状态信息。另外，在步骤S4中，在XI、XQ、YI、YQ四个独立光场分量的传播过程中，在每个独立光场分量上引入独立的随机噪声。
如图4所示，一列平面电磁波(光场)最多可以由4个独立分量叠加，且完全在同一物理空间上传播，这四个分量分别为XI、XQ、YI、YQ。通过现有成熟的方法，这四个分量可以容易地被区分开并被测量。如图5和图6所示，分别为X偏振方向光场分量和Y偏振方向光场分量，这两个分量上的电磁波相互正交。如图7所示，每个偏振方向的分量中又分别含有I和Q两个分量，这两个分量之间存在的相位差，在复平面上相互正交，所以复平面上的任意一点都可以分解成正交的I、Q两个分量，如图2中矢量A可分解为AI、AQ，矢量B可分解为BI、BQ。
因此，可将光场分成四个分量，形成四个通路，并且这些通路可以单独或一起使用，对光场进行高维复用，高维复用是指充分利用光的多个维度。本实施例的通过利用光场在偏振、相位上的独立分量，构建了多个计算模式：对于伊辛机模式，XI、XQ、YI、YQ四个正交的分量可以分别用于加载四路不同的节点，故可以独立求解4个伊辛模型；对于波茨机模式，正如图6中矢量A或矢量B，利用I、Q两个分量可以构建可取相位数大于2个的节点，故利用X和Y偏振分量上的I、Q分量，可以加载两路不同的节点，故可以独立求解2个波茨模型；对于三维节点状态空间模式和四维节点状态空间模式这样的高维求解器模式，与可以选用四路独立分量中的3路或4路，构建可取值的参数空间维度更高的节点，以求解更高阶的模型。
S5、接收经过传播的XI、XQ、YI、YQ四个独立光场分量，并将四个独立光场分量与步骤S2作为本振参考光的激光进行对比、干涉，分离成四路独立的光场信号，并转换成包含节点状态信息的四路模拟电信号。本实施例通过相干光信号接收机接收四个独立光场分量。通过光电探测器阵列接收已经分离开的四路独立的光场信号并转换成四路模拟电信号。
S6、将步骤S5得到的包含节点状态信息的四路模拟电信号转换为数字电域信息，根据步骤S1获得的模型参数和选定的计算模式进行计算并更新节点状态信息，将更新后的节点状态信息转换为模拟电信号反馈到步骤S3中用于加载到四个独立光场分量上，循环执行步骤S3、S4、S5和S6，直至循环至步骤S1选择的循环迭代次数。如图2所示，其中，根据接收到的节点状态信息、获得的模型参数和选定的计算模式进行计算，包括：
S6.1、节点状态信息的采集与重整:采集每个节点的状态信息，然后根据选择的计算模式重整节点的形式，如：伊辛机模式下，每个伊辛机通道的节点都不需要变换形式；波茨机模式下，同一偏振态X或Y下的I、Q两路的节点重整合并为一路：xX＝xXI+j·xXQ和xY＝xYI+j·xYQ，其中j是虚数单位；在三维节点状态空间模式和四维节点状态空间模式的高维求解器模式下，多个通道的数据重整为一路，如：xX＝xXI+j·xXQ+k·xYI+l·xYQ，其中j,k,l为四元数的三个虚数单位；
S6.2、节点增益计算：根据模型参数α，计算增益项αxi，其中，α为模型的增益系数，xi为第i节点的状态值；
S6.3、节点耦合计算，根据模型参数β、Jij，计算耦合项β∑j≠iJijxj，其中，β为模型的耦合强度系数，Jij为模型耦合常数，xj为第j节点的状态值；
S6.4、节点状态信息反馈信号生成：根据
f0i＝αxi+β∑j≠iJijxj，
求得f0i，f0i为节点状态信息反馈信号，再根据各计算模式下取值空间的要求，将f0i投影到其最接近的矢量方向上，更新节点状态信息反馈信号，生成更新的节点状态信息反馈信号
Figure GDA0003692495050000091
对于波茨机模式、三维节点状态空间模式和四维节点状态空间模式，节点状态信息反馈信号为：
Figure GDA0003692495050000092
其中单位矢量
Figure GDA0003692495050000093
的相位角为
Figure GDA0003692495050000094
且与f0i的相位角最接近。
如图3所示，本实施例还提供一种波茨模型的高维复用计算装置，包括：
模型与任务获取模块，用于获取待求解的模型及模型参数，并根据待求解的模型类型选择计算模式，本实施例的模型与任务获取模块还用于接收并存储来自节点信息处理模块的节点状态信息。计算模式有四个选项：(1)伊辛机模式：该模式下，节点的状态空间是一维的；(2)波茨机模式：该模式下，节点的状态空间是二维的，可由两路通道的信号重整为复数；(3)三维节点状态空间模式：该模式下，节点的状态空间是三维的，可由三路通道的信号重整化为一个三元数；(4)四维节点状态空间模式：该模式下，节点的状态空间是四维的，可由四路通道的信号重整化为一个四元数。
激光光源，用于发出的连续激光，且将激光分为两路，一路作为承载节点状态的载波，输入到电-光节点更新模块进行处理；另一路作为本振参考光输入到光-电节点测量模块以供四路独立的复用光场节点状态的相干检测。
电-光节点更新模块，用于接收激光光源发出的作为承载节点状态的载波的激光，并将该激光分成XI、XQ、YI、YQ四个独立光场分量，并把节点信息处理模块反馈的四个通道的节点状态信息中的模拟电域信息
Figure GDA0003692495050000101
分别加载到上述四个独立光场分量的强度
Figure GDA0003692495050000102
上，其中，
Figure GDA0003692495050000103
表示第k个循环中的第i节点的模拟电域信息，
Figure GDA0003692495050000104
表示第k个循环中的第i节点的独立光场分量的强度；在进行第一次加载时，所有节点的初值都设置为0。本实施例的电-光节点更新模块主要由一个双偏振IQ电光调制器组成，使双偏振IQ电光调制器的每一个独立调制通道，其直流偏置点都应在-π/4(线性区)工作点，使得对第k个循环中的第i节点在反馈信号
Figure GDA0003692495050000105
的调制后为：
Figure GDA0003692495050000106
在此，外加的常数项
Figure GDA0003692495050000107
使节点状态取值为[-1,1]，包含正值或者负值，再取
Figure GDA0003692495050000108
符号以实现节点的二值化。
光路复用系统，用于连接电-光节点更新模块与光-电节点测量模块，传播XI、XQ、YI、YQ四个独立光场分量；本实施例的光路复用系统采用对偏振不敏感的光纤，或者其他波导器件、自由空间光等。需要强调的是，这一段光路，可以支持光的高维度传输。通过光路复用系统传播，可以支持XI、XQ、YI、YQ四个独立光场分量共享同一个物理通道，并且同步且无串扰地从电-光节点更新模块到光-电节点测量模块传输节点状态信息。另外，光路复用系统还用于在传播时，在每个独立分量上引入独立的随机噪声。
光-电节点测量模块，用于接收光路复用系统传播的XI、XQ、YI、YQ四个独立光场分量以及激光光源发出的作为本振参考光的激光，且将XI、XQ、YI、YQ四个独立光场分量和本振参考光进行干涉、对比，并分离成四路独立的光场信号，将光场信号转换成包含节点状态信息的四路模拟电信号；本实施例的光-电节点测量模块包括相干光信号接收机和光电探测器阵列，相干光信号接收机用于接收四个独立光场分量以及本振参考光，光电探测器阵列用于将分离的四路独立的光场信号转换成四路模拟电信号。
节点信息处理模块，用于接收光-电节点测量模块的包含节点状态信息的四路模拟电信号并转换为数字电域信息且暂时储存，接收模型与任务获取模块的计算模式和模型参数，对节点信息进行计算、更新和暂时储存，将计算更新后的节点信息转换成模拟电信号反馈给电-光节点更新模块。本实施例的节点信息处理模块包括：模拟-数字转换器，用于将接收到的光-电节点测量模块的包含节点状态信息的四路模拟电信号并转换为数字电域信息；数字-模拟转换器，用于将计算更新后的节点信息转换成模拟电信号并反馈给电-光节点更新模块。
节点信息处理模块包括数学运算单元，数学运算单元用于根据接收到的模型与任务获取模块的计算模式和模型参数，对节点信息进行计算、更新，数学运算单元包括：
节点信息采集与重整子单元，用于采集每个节点的状态值，且根据所指定的计算模式，重整节点的形式；
节点增益计算子单元，用于根据模型参数α，计算增益项αxi，其中，α为模型的增益系数，xi为第i节点的状态值；
节点耦合计算子单元，用于模型参数β、Jij，计算耦合项β∑j≠iJijxj，其中，β为模型的耦合强度系数，Iij为模型耦合常数，xj为第j节点的状态值；
节点状态更新反馈信号生成子单元，用于获取增益项αxi和耦合项β∑j≠iJijxj，并计算f0i＝αxi+β∑j≠iJijxj得到节点状态信息反馈信号f0i，再根据各计算模式下取值空间的要求，将f0i投影到其最接近的矢量方向上，更新节点状态信息反馈信号，生成更新的节点状态信息反馈信号
Figure GDA0003692495050000111
在本实施例中，若模型与任务获取模块选择的计算模式为波茨机模式、三维节点状态空间模式或四维节点状态空间模式，更新的节点状态反馈信号应为：
Figure GDA0003692495050000112
其中单位矢量
Figure GDA0003692495050000113
的相位角为
Figure GDA0003692495050000114
且与f0i的相位角最接近。
综上，本发明实施例提供一种波茨模型的高维复用计算方法及装置，通过激光传播节点状态信息，接收后将光场信号转换成模拟电信号，在计算时转换成数据信号，再将计算更新后的反馈信号转换成模拟电信号加载到光场分量上，形成光电环路的信息传递，以光场的强度编码节点的状态，实现了波茨模型的简易、高效求解，且解决了现有的计算装置复杂、维持装置稳定成本高的问题。并且，本发明通过利用了光场的两个正交的偏振态X和Y以及两个正交相位0(I)和π(Q)进行计算通道的复用，实现了XI、XQ、YI、YQ四路独立运行的二值化节点伊辛机，且能够实现可取相位状态数目任意的节点的波茨机，使得算法复杂度得到降低；并且，由于四个独立光场分量形成的四个通道，可根据模型类型选择不同的计算模式，具有较高的灵活性；另外，由于基于正交偏振和相位的光场可以相互叠加而不产生串扰，所以本发明的装置更加紧凑且更加易于实现。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和替换，这些改进和替换也应视为本发明的保护范围。

Claims (8)

1.一种波茨模型的高维复用计算方法，其特征在于，包括如下步骤：
S1、获取待求解的模型及其模型参数，根据待求解的模型类型选择计算模式；
S2、将激光光源发出的连续激光分为两路，一路作为承载节点状态的载波，另一路作为本振参考光以供四路独立的复用光场节点状态的相干检测；
S3、将步骤S2作为承载节点状态的载波的激光分成XI、XQ、YI、YQ四个独立光场分量，并把四个独立光场分量对应的通道的节点状态信息中的模拟电域信息
Figure FDA0003692495040000011
分别加载到上述四个独立光场分量的强度
Figure FDA0003692495040000012
上，其中，
Figure FDA0003692495040000013
表示第k个循环中的第i节点的模拟电域信息，
Figure FDA0003692495040000014
表示第k个循环中的第i节点的独立光场分量的强度；
S4、通过光路复用系统传播XI、XQ、YI、YQ四个独立光场分量；
S5、接收经过传播的XI、XQ、YI、YQ四个独立光场分量，并将四个独立光场分量与步骤S2作为本振参考光的激光进行对比、干涉，分离成四路独立的光场信号，并转换成包含节点状态信息的四路模拟电信号；
S6、将步骤S5得到的包含节点状态信息的四路模拟电信号转换为数字电域信息，根据步骤S1获得的模型参数和选定的计算模式进行计算并更新节点状态信息，将更新后的节点状态信息转换为模拟电信号反馈到步骤S3中用于加载到四个独立光场分量上，循环执行步骤S3、S4、S5和S6；步骤S1获得的模型参数包括增益系数α、耦合强度系数β和耦合常数Jij，在步骤S6中，根据接收到的节点状态信息、获得的模型参数和选定的计算模式进行计算，包括：
S6.1、节点状态信息的采集与重整:采集每个节点的状态信息，然后根据选择的计算模式重整节点的形式；
S6.2、节点增益计算：根据模型参数α，计算增益项αxi，其中，α为模型的增益系数，xi为第i节点的状态值；
S6.3、节点耦合计算，根据模型参数β、Jij，计算耦合项β∑j≠iJijxj，其中，β为模型的耦合强度系数，Jij为模型耦合常数，xj为第j节点的状态值；
S6.4、节点状态信息反馈信号生成：根据
f0i＝αxi+β∑j≠iJijxj，
求得f0i，f0i为节点状态信息反馈信号，再根据各计算模式下取值空间的要求，将f0i投影到其最接近的矢量方向上，更新节点状态信息反馈信号，生成更新的节点状态信息反馈信号
Figure FDA0003692495040000021
2.根据权利要求1所述的波茨模型的高维复用计算方法，其特征在于，在步骤S1中，选择计算模式包括选择求解器模式和循环迭代次数，其中，求解器模式包括：
(1)伊辛机模式：该模式下，节点的状态空间是一维的；
(2)波茨机模式：该模式下，节点的状态空间是二维的，可由两路通道的信号重整为复数；
(3)三维节点状态空间模式：该模式下，节点的状态空间是三维的，可由三路通道的信号重整化为一个三元数；
(4)四维节点状态空间模式：该模式下，节点的状态空间是四维的，可由四路通道的信号重整化为一个四元数。
3.根据权利要求1所述的波茨模型的高维复用计算方法，其特征在于，在步骤S3中，四个独立光场分量通过双偏振IQ电光调制器加载，对于双偏振IQ电光调制器的每一个独立调制通道，其直流偏置点都应在线性区的-π/4线性区工作点，使得对第k个循环中的第i节点在反馈信号
Figure FDA0003692495040000022
的调制后为：
Figure FDA0003692495040000023
其中，外加的常数项
Figure FDA0003692495040000024
使节点状态取值为[-0.5,0.5]，包含正值或者负值，再取
Figure FDA0003692495040000025
符号以实现节点的二值化。
4.根据权利要求2所述的波茨模型的高维复用计算方法，其特征在于，在步骤S6.1中，在伊辛机模式下，每个伊辛机通道的节点都不需要变换形式；在波茨机模式下，同一偏振态X或Y下的I、Q两路的节点重整合并为一路：xX＝xXI+j·xXQ和xY＝xYI+j·xYQ，其中j是虚数单位；在三维节点状态空间模式和四维节点状态空间模式下，多个通道的数据重整为一路。
5.根据权利要求2所述的波茨模型的高维复用计算方法，其特征在于，对于波茨机模式、三维节点状态空间模式和四维节点状态空间模式，节点状态信息反馈信号为：
Figure FDA0003692495040000031
其中单位矢量
Figure FDA0003692495040000032
的相位角为
Figure FDA0003692495040000033
且与f0i的相位角最接近。
6.根据权利要求1所述的波茨模型的高维复用计算方法，其特征在于，在步骤S4中，在XI、XQ、YI、YQ四个独立光场分量的传播过程中，在每个独立光场分量上引入独立的随机噪声。
7.一种波茨模型的高维复用计算装置，其特征在于，包括：
模型与任务获取模块，用于获取待求解的模型及模型参数，并根据待求解的模型类型选择计算模式；
激光光源，用于发出的连续激光，且将激光分为两路，一路作为承载节点状态的载波，输入到电-光节点更新模块进行处理；另一路作为本振参考光输入到光-电节点测量模块以供四路独立的复用光场节点状态的相干检测；
电-光节点更新模块，用于接收激光光源发出的作为承载节点状态的载波的激光，并将该激光分成XI、XQ、YI、YQ四个独立光场分量，并把节点信息处理模块反馈的四个通道的节点状态信息中的模拟电域信息
Figure FDA0003692495040000034
分别加载到上述四个独立光场分量的强度
Figure FDA0003692495040000035
上，其中，
Figure FDA0003692495040000036
表示第k个循环中的第i节点的模拟电域信息，
Figure FDA0003692495040000037
表示第k个循环中的第i节点的独立光场分量的强度；
光路复用系统，用于连接电-光节点更新模块与光-电节点测量模块，传播XI、XQ、YI、YQ四个独立光场分量；
光-电节点测量模块，用于接收光路复用系统传播的XI、XQ、YI、YQ四个独立光场分量以及激光光源发出的作为本振参考光的激光，且将XI、XQ、YI、YQ四个独立光场分量和本振参考光进行干涉、对比，并分离成四路独立的光场信号，将光场信号转换成包含节点状态信息的四路模拟电信号；
节点信息处理模块，用于接收光-电节点测量模块的包含节点状态信息的四路模拟电信号并转换为数字电域信息且暂时储存，接收模型与任务获取模块的计算模式和模型参数，对节点信息进行计算、更新和暂时储存，将计算更新后的节点信息转换成模拟电信号反馈给电-光节点更新模块；所述节点信息处理模块包括数学运算单元，数学运算单元用于根据接收到的模型与任务获取模块的计算模式和模型参数，对节点信息进行计算、更新，数学运算单元包括：
节点信息采集与重整子单元，用于采集每个节点的状态值，且根据所指定的计算模式，重整节点的形式；
节点增益计算子单元，用于根据模型参数α，计算增益项αxi，其中，α为模型的增益系数，xi为第i节点的状态值；
节点耦合计算子单元，用于模型参数β、Jij，计算耦合项β∑j≠iJijxj，其中，β为模型的耦合强度系数，Jij为模型耦合常数，xj为第j节点的状态值；
节点状态更新反馈信号生成子单元，用于获取增益项axi和耦合项β∑j≠iJijxj，并计算f0i＝αxi+β∑j≠iJijxj得到节点状态信息反馈信号f0i，再根据各计算模式下取值空间的要求，将f0i投影到其最接近的矢量方向上，更新节点状态信息反馈信号，生成更新的节点状态信息反馈信号
Figure FDA0003692495040000041
8.根据权利要求7所述的波茨模型的高维复用计算装置，其特征在于，所述电-光节点更新模块包括双偏振IQ电光调制器。
CN202110477851.5A 2021-04-29 2021-04-29 一种波茨模型的高维复用计算方法及装置 Active CN113190998B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110477851.5A CN113190998B (zh) 2021-04-29 2021-04-29 一种波茨模型的高维复用计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110477851.5A CN113190998B (zh) 2021-04-29 2021-04-29 一种波茨模型的高维复用计算方法及装置

Publications (2)

Publication Number Publication Date
CN113190998A CN113190998A (zh) 2021-07-30
CN113190998B true CN113190998B (zh) 2022-11-18

Family

ID=76982885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110477851.5A Active CN113190998B (zh) 2021-04-29 2021-04-29 一种波茨模型的高维复用计算方法及装置

Country Status (1)

Country Link
CN (1) CN113190998B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111373344A (zh) * 2017-10-19 2020-07-03 日本电信电话株式会社 波茨模型的计算装置
CN111999914A (zh) * 2020-08-18 2020-11-27 华中科技大学 一种集成全维度高速光场调控的方法和装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006033967A2 (en) * 2004-09-16 2006-03-30 Hewlett-Packart Development Company, L.P. A method and apparatus for modeling systems
US11017309B2 (en) * 2017-07-11 2021-05-25 Massachusetts Institute Of Technology Optical Ising machines and optical convolutional neural networks
JP6684259B2 (ja) * 2017-10-19 2020-04-22 日本電信電話株式会社 ポッツモデルの計算装置
CN110445550A (zh) * 2019-08-08 2019-11-12 深圳市深光谷科技有限公司 一种相干光直接检测方法、装置及光通信系统
CN110516378B (zh) * 2019-08-30 2022-08-09 桂林电子科技大学 一种锡基二元共晶相分离的有限元模拟方法
CN112486244A (zh) * 2019-09-11 2021-03-12 华为技术有限公司 一种光计算设备以及光信号处理方法
CN111006675B (zh) * 2019-12-27 2022-10-18 西安理工大学 基于高精度重力模型的车载激光惯导系统自标定方法
CN111277336B (zh) * 2020-01-19 2023-06-27 山东大学 一种实现数据中心内部互连的高速光交换模块及其运行方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111373344A (zh) * 2017-10-19 2020-07-03 日本电信电话株式会社 波茨模型的计算装置
CN111999914A (zh) * 2020-08-18 2020-11-27 华中科技大学 一种集成全维度高速光场调控的方法和装置

Also Published As

Publication number Publication date
CN113190998A (zh) 2021-07-30

Similar Documents

Publication Publication Date Title
CN103678258B (zh) 一种提高硅基光学矩阵处理器数据分辨率的方法
TWI825452B (zh) 光電計算系統
Rafayelyan et al. Large-scale optical reservoir computing for spatiotemporal chaotic systems prediction
Hu et al. Progress in quantum teleportation
CN109784486B (zh) 一种光学神经网络处理器及其训练方法
US11398871B2 (en) Systems and methods for analog computing using a linear photonic processor
TW202201165A (zh) 計算系統、計算裝置及計算系統的操作方法
CN114819132B (zh) 一种基于时间-波长交织的光子二维卷积加速方法及系统
TW202147060A (zh) 光電計算系統
van Loock Quantum communication with continuous variables
US20240078421A1 (en) Two-dimensional photonic convolutional acceleration system and device for convolutional neural network
CN109784485A (zh) 一种光学神经网络处理器及其计算方法
CN112823359A (zh) 光电计算系统
CN111753977A (zh) 光学神经网络卷积层芯片、卷积计算方法和电子设备
Calvanese Strinati et al. All-optical scalable spatial coherent Ising machine
CN113722667A (zh) 基于伊辛机的数据处理方法、装置及伊辛机
CN113190998B (zh) 一种波茨模型的高维复用计算方法及装置
Bendjaballah Introduction to photon communication
Hong et al. Generation of two-photon hybrid-entangled W state with photonic qubit and time-bin via cross-Kerr nonlinearities
JP3999752B2 (ja) エンタングル光子対発生装置及びその方法
WO2023174072A1 (zh) 数据处理方法及系统
Zhang et al. Experimental observation of quantum nonlocality in general networks with different topologies
Cohen Synchronization and prediction of chaotic dynamics on networks of optoelectronic oscillators
CN109117117B (zh) 一种自适应的源无关量子随机数发生器
CN115118342B (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
GR01 Patent grant
GR01 Patent grant