CN116224604A - 一种实现矩阵求逆的光学结构和方法 - Google Patents

一种实现矩阵求逆的光学结构和方法 Download PDF

Info

Publication number
CN116224604A
CN116224604A CN202211675054.9A CN202211675054A CN116224604A CN 116224604 A CN116224604 A CN 116224604A CN 202211675054 A CN202211675054 A CN 202211675054A CN 116224604 A CN116224604 A CN 116224604A
Authority
CN
China
Prior art keywords
optical
matrix
matrix inversion
path structure
phase grating
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.)
Pending
Application number
CN202211675054.9A
Other languages
English (en)
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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202211675054.9A priority Critical patent/CN116224604A/zh
Publication of CN116224604A publication Critical patent/CN116224604A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/09Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
    • G02B27/0938Using specific optical elements
    • G02B27/0944Diffractive optical elements, e.g. gratings, holograms
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/09Beam shaping, e.g. changing the cross-sectional area, not otherwise provided for
    • G02B27/0938Using specific optical elements
    • G02B27/0994Fibers, light pipes
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B6/00Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings
    • G02B6/02Optical fibres with cladding with or without a coating
    • G02B6/028Optical fibres with cladding with or without a coating with core or cladding having graded refractive index
    • G02B6/0288Multimode fibre, e.g. graded index core for compensating modal dispersion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Dispersion Chemistry (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Optical Couplings Of Light Guides (AREA)

Abstract

本发明公开了一种实现矩阵求逆的光学结构和方法,矩阵求逆的光路结构由n个相互关联且完全相同的单元模块、1块透镜、1个位相光栅组成,每个所述单元模块均由n‑1个光栅衰减器、1个多模光纤阵列组成。本发明采用上述的一种实现矩阵求逆的光学结构和方法,具有新颖、光学元件少,可拓展性强的优点,只需一个位相光栅可以实现2×2到N×N的逆矩阵计算,在数值分析、光计算,光通信网络和深度学习神经网络等方面有着广泛的应用前景。

Description

一种实现矩阵求逆的光学结构和方法
技术领域
本发明涉及光计算技术领域,尤其是涉及一种实现矩阵求逆的光学结构和方法。
背景技术
随着全球化以及科技的高速发展,需要处理的数据量在急剧增加,相应的数据处理模型和算法也在不断增加,带来的结果就是对算力和功耗的要求不断提高。而目前冯诺依曼架构的电子计算机存在传输瓶颈、功耗增加以及算力瓶颈等问题,越来越难以满足大数据时代对算力与功耗的需求,因此提高运算速度同时降低运算功耗是目前面临的紧要问题。面向未来的计算系统,光计算具有天然优势:光子具有光速传播、抗电磁干扰、任意叠加等特性,因此光学计算具有天然的并行计算特性,因而运算速度极快,且十分适合做并行运算。
目前针对光学矩阵计算的研究主要有多平面衍射光转换,片上微环谐振器阵列和片上马赫曾德尔干涉仪网络等实现方式。英国南安普敦大学的Kan Wu,Cesare Soci等人提出了一种在光纤网络上实现3×3矩阵求逆的方案,采用了光纤,9个1×2光分路器,9个光纤调制器,当矩阵维度扩充到N×N时,需要N^2个1×2光分路器,因此此方案虽然能实现矩阵求逆,但所需元件多,结构复杂,耦合损耗高,系统庞大,并且光路的可拓展性差,在矩阵维度增加时,需要增加大量光学元件,光路结构复杂,调试困难。本发明提出的矩阵求逆结构,基于位相光栅,避免大量使用二元分光元件,对于任意N维维度都仅用一块位相光栅和一块透镜即可实现多维矩阵求逆运算,使得系统更加简单,减少了所需的光学元件,可拓展性强,使光计算架构朝着微型化、集成化的方向发展。
发明内容
本发明的目的在于提供一种实现矩阵求逆的光学结构和方法,具有新颖、光学元件少、可拓展性强、可以实现n*n矩阵求逆功能的特点。
为实现上述目的,本发明提供了一种实现矩阵求逆的光路结构,矩阵求逆的光路结构由n个相互关联且完全相同的单元模块、1块透镜、1个位相光栅组成,每个所述单元模块均由n-1个光栅衰减器、1个多模光纤阵列组成。
一种矩阵求逆的光路结构实现方法,具体步骤如下:
S1、从多模光纤阵列输出的n路光信号经过1*n均匀分光的位相光栅,得到n组平行光束沿不同方向输出;
S2、由位相光栅输出的n组平行光经过透镜汇聚后耦合进入多模光纤阵列中;
S3、由多模光纤阵列中的光信号经过各自的光衰减器衰减指定倍率,分别到不同模块成为输入信号之一;
S4、依次改变输入光信号I的光强值,即得n*n矩阵的逆矩阵;
S5、整个n×n矩阵求逆光路结构采用n个与步骤S1-S4完全相同的单元模块,每一个所述单元模块的工作方式相同。
优选的,步骤S4中,光强值输入时,依次令Ij
Figure BDA0004016859890000021
(j=1,2……n),其他值为零,其中,n*k为位相光栅总分光效率;每次得到每个模块的Xi信号值即为矩阵M的逆矩阵的j列。
优选的,步骤S5中,单元模块之间采用固定光纤通道形式连接。
优选的,步骤S5中,每一个单元模块光信号的输出为所有单元模块光信号的输入的其中一个。
因此,本发明采用上述一种实现矩阵求逆的光学结构和方法,具有新颖、光学元件少、可拓展性强、可以实现n*n矩阵求逆功能的特点。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1是已有方案中单个节点的结构图;
图2是已有方案中节点间的相互连接关系图;
图3是为本方案全光方式实现3×3矩阵求逆光路的光路图;
图4为本方案全光方式实现n×n矩阵求逆光路的一个单元模块的结构图;
图5为本方案全光方式实现n×n矩阵求逆光路总体结构图;
图6入射光经过位相光栅S后分光的过程;
图7不同方向入射平行光经过透镜L2耦合进入多模光纤阵列过程;
图8为本方案全光方式实现n×n矩阵求逆光路的一个单元模块的光路图;
图9为本方案全光方式实现n×n矩阵求逆光路总体光路图。
附图标记
L、透镜;S、位相光栅;Dij、可调光纤衰减器;F、多模光纤列阵。
具体实施方式
以下通过附图和实施例对本发明的技术方案作进一步说明。
除非另外定义,本发明使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。
全光方式实现矩阵求逆的光路,采用自由空间光和光纤混合互连结构,选用少量的光学元件和模块化结构,以3×3矩阵求逆光路为例,在南安普敦大学的Kan Wu等人提出的方案中,系统分为3个相互连接的节点,每个节点的结构如图1所示,通过光衰减器D对光强调节,通过3个2×2光耦合器实现分光和和光功能。不同节点间的连接方式如图2所示。采用这种方案组建的光路需要用到大量的2×2光耦合器,系统复杂,不易扩展。本发明在此基础上进行了改进,采用一块位相光栅和一块透镜实现了同样的功能,如图3所示。
光路结构由上下三个相同的光路模块组成,模块之间通过光纤进行连接。在上方的模块中,光从光纤阵列F出射,经过1×3分光的位相光栅分光后,经过透镜L后汇聚进入光纤阵列F中,随后经过光衰减器。其中位相光栅总分光效率为3×k,I1、I2、I3为信号输入口,X1、X2、X3为信号输出口,光线经过光栅分光,透镜汇聚后的光强分别为X1、X2、X3
因此我们可以得到如下关系:
k*(X2*D2,1+X3*D3,1+I1)=X1 (1)
k*(X1*D1,2+X3*D3,2+I2)=X2 (2)
k*(X1*D1,3+X2*D2,3+I3)=X3 (3)
连列(1)(2)(3)式可得矩阵方程:
Figure BDA0004016859890000041
即:
k(DX+I)=X (5)
X=k*(E-k*D)-1I (6)
令E-k*D=M,
Figure BDA0004016859890000042
I2=0,I3=0可以得到X为矩阵M的逆矩阵第一列,令I1=0,
Figure BDA0004016859890000043
I1=0可以得到X为矩阵M的逆矩阵第二列,再令I1=0,I2=0,/>
Figure BDA0004016859890000044
可以得到X为矩阵M的逆矩阵第三列。
如图5所示,从3*3矩阵拓展到n*n矩阵,光路结构由n个完全相同的光路模块组成,模块之间采用光纤连接的方式。单个模块的结构如图4所示。整个光路仅由一块1×n均匀分光的位相光栅S,一个透镜L,和可调光纤衰减器Dij(j=1,2…n)组成,i表示第i个模块。输入的n束光信号阵列入射经过1*n均匀分光位相光栅S后,每一束光信号都会分为n束相同光强,沿不同出射方向的平行光,因此我们得到了n组平行光。透镜L2将每组平行光汇聚后进入多模光纤阵列F中的不同光纤,随后经过光纤衰减器Di对光强按照要求进行衰减,分别成为其他n-1个单元模块的输入信号之一,最后一路输出作监测作用。
图6为位相光栅S对入射光进行分光的光路图,位相光栅是一种可将一束入射的单色平行光分为一维或二维等间距(或非等间距)的光束列阵,具有不等间距、周期重复结构的傅里叶型分束器件。此处使用的是1*n均匀分光的位相光栅,会将入射光束均匀地分成n束平行光沿不同方向出射。可以通过使用光刻的方法在位相光栅分光器表面周期性地刻蚀出深浅不一的台阶状结构,当平行入射的激光束经过位相光栅时,通过光栅厚度的周期性变化对入射的单色平面波进行空间位相调制,经透镜变换,在后焦面成像,形成周期性分布的光点。
图7是透镜L,将位相光栅S出射的n束平行光汇聚到焦平面上,然后耦合进入多模光纤阵列。为多模光纤选择耦合镜片需要透镜的数值孔径(NA)和光纤的数值孔径接近,使光源的焦点大小和光纤的纤芯大小匹配,并且使入射圆锥角不超过光纤数值孔径的反正弦,来达到较高的耦合效率。
可调光纤衰减器是通过光信号的吸收、反射、扩散、散射、偏转、衍射、色散等来降低光功率。本发明中可调光纤衰减器设置数值由需要求逆的矩阵决定。
n×n矩阵求逆的光路结构一个单元模块的光信号处理过程,如图8所示,输入光信号入射到光学系统中,光束入射到位相光栅S上,位相光栅S会将入射平行光束均分为n束沿不同方向出射的平行光,如图6所示。平行光经过透镜L汇聚到焦平面上。因为透镜一个数值孔径(NA)和光纤的数值孔径较为接近,使光源的焦点大小和光纤的纤芯大小匹配,并且入射圆锥角不超过光纤数值孔径的反正弦,光能较好地耦合进入多模光纤中,可参考图7。
在光纤中传输的光信号,经过可调光纤衰减器衰减指定倍率后,分别成为n个单元模块输入信号之一,如此循环往复,系统趋于稳定。
整体的矩阵求逆光路图如图9所示:
Xi表示第i个模块中光信号在光纤阵列L2中的光强,因此在第1个模块中有如下关系:
k*(X2*D2,1+…+Xn*Dn,1+I1)=X1 (7)
k等于光栅的总分光效率除以n。
同理,可以得到在第i个模块中有如下的关系:
k*(X1*D1,i+…+Xi-1*Di-l,i+Xi+1*Di+1,i+…+Xn*Dn,i+Ii)=Xi (8)
由(8)式可得矩阵方程:
Figure BDA0004016859890000063
k(D*X+I)=X (10)
化简可得:
X=k*(E-k*D)-1I (11)
其中E为单位矩阵。
进一步令E-k*D=M,并带入(11)式可得:
Figure BDA0004016859890000064
由此,首先令
Figure BDA0004016859890000061
Ii=0(i≠1)可以得到/>
Figure BDA0004016859890000062
为M逆矩阵的第一列。依此类推令/>
Figure BDA0004016859890000071
其余值为0,可以以此得到矩阵M的逆矩阵的第j列。
综上所述,本发明采用上述一种实现矩阵求逆的光学结构和方法,具有新颖、光学元件少、可拓展性强、可以实现n*n矩阵求逆功能的特点。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。

Claims (5)

1.一种实现矩阵求逆的光路结构,其特征在于:矩阵求逆的光路结构由n个相互关联且完全相同的单元模块、1块透镜、1个位相光栅组成,每个单元模块均由n-1个光栅衰减器、1个多模光纤阵列组成。
2.一种矩阵求逆的光路结构实现方法,其特征在于,具体步骤如下:
S1、从多模光纤阵列输出的n路光信号经过1*n均匀分光的位相光栅,得到n组平行光束沿不同方向输出;
S2、由位相光栅输出的n组平行光经过透镜汇聚后耦合进入多模光纤阵列中;
S3、由多模光纤阵列中的光信号经过各自的光衰减器衰减指定倍率,分别到不同模块成为输入信号之一;
S4、依次改变输入光信号I的光强值,即得n*n矩阵的逆矩阵;
S5、整个n×n矩阵求逆光路结构采用n个与步骤S1-S4完全相同的单元模块,每一个单元模块的工作方式相同。
3.根据权利要求2所述的一种矩阵求逆的光路结构实现方法,其特征在于:步骤S4中,光强值输入时,依次令Ij
Figure FDA0004016859880000011
其他值为零,其中,n*k为位相光栅总分光效率;每次得到每个模块的Xi信号值即为矩阵M的逆矩阵的j列。
4.根据权利要求2所述的一种矩阵求逆的光路结构实现方法,其特征在于:步骤S5中,单元模块之间采用固定光纤通道形式连接。
5.根据权利要求2所述的一种矩阵求逆的光路结构实现方法,其特征在于:步骤S5中,每一个单元模块光信号的输出为所有单元模块光信号的输入的其中一个。
CN202211675054.9A 2022-12-26 2022-12-26 一种实现矩阵求逆的光学结构和方法 Pending CN116224604A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211675054.9A CN116224604A (zh) 2022-12-26 2022-12-26 一种实现矩阵求逆的光学结构和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211675054.9A CN116224604A (zh) 2022-12-26 2022-12-26 一种实现矩阵求逆的光学结构和方法

Publications (1)

Publication Number Publication Date
CN116224604A true CN116224604A (zh) 2023-06-06

Family

ID=86577688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211675054.9A Pending CN116224604A (zh) 2022-12-26 2022-12-26 一种实现矩阵求逆的光学结构和方法

Country Status (1)

Country Link
CN (1) CN116224604A (zh)

Similar Documents

Publication Publication Date Title
Liu et al. Holylight: A nanophotonic accelerator for deep learning in data centers
Bai et al. Photonic multiplexing techniques for neuromorphic computing
WO2021227767A1 (zh) 在光波导和/或光芯片上实现的全光衍射神经网络及系统
Huang et al. Programmable matrix operation with reconfigurable time-wavelength plane manipulation and dispersed time delay
Soref Reconfigurable integrated optoelectronics
CN110908428B (zh) 一种高效实现大规模矩阵运算的并行光计算系统
Miller Analyzing and generating multimode optical fields using self-configuring networks
Shiflett et al. Albireo: Energy-efficient acceleration of convolutional neural networks via silicon photonics
US12001947B2 (en) Ultra-wide data band optical processor
Peserico et al. Integrated photonic tensor processing unit for a matrix multiply: a review
Shao et al. Generalized robust training scheme using genetic algorithm for optical neural networks with imprecise components
Tanomura et al. Scalable and robust photonic integrated unitary converter based on multiplane light conversion
CN117077750B (zh) 一种并行光子卷积运算芯片及系统
CN113935493A (zh) 一种基于集成光学的可编程高维量子计算芯片结构
SeyedinNavadeh et al. Determining the optimal communication channels of arbitrary optical systems using integrated photonic processors
Guo et al. StarLight: a photonic neural network accelerator featuring a hybrid mode-wavelength division multiplexing and photonic nonvolatile memory
US20240078419A1 (en) Optical neuron unit and network of the same
CN116224604A (zh) 一种实现矩阵求逆的光学结构和方法
Gur et al. Optical implementation of fuzzy-logic controllers. Part I
Pinkston et al. Design of an optical reconfigurable shared-bus–hypercube interconnect
WO2022060908A1 (en) Frequency multiplexed photonic neural networks
WO2022136146A1 (en) Optical computing and reconfiguring with spatiotemporal nonlinearities in waveguides
CN116112088B (zh) 多输入多输出型光子器件的自适应调制方法
Nakajima Parallel neuromorphic computing on space and wavelength division multiplexed photonic processor
Morris et al. Reconfigurable optical interconnects by a combined computer-generated hologram and spatial light modulator method

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