[go: up one dir, main page]

CN109033015B - Device for executing calculus operation on optical signal - Google Patents

Device for executing calculus operation on optical signal Download PDF

Info

Publication number
CN109033015B
CN109033015B CN201710445240.6A CN201710445240A CN109033015B CN 109033015 B CN109033015 B CN 109033015B CN 201710445240 A CN201710445240 A CN 201710445240A CN 109033015 B CN109033015 B CN 109033015B
Authority
CN
China
Prior art keywords
calculus
signal
arithmetic unit
optical signal
signal source
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
CN201710445240.6A
Other languages
Chinese (zh)
Other versions
CN109033015A (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.)
Nanjing Agricultural University
Original Assignee
Nanjing Agricultural 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 Nanjing Agricultural University filed Critical Nanjing Agricultural University
Priority to CN201710445240.6A priority Critical patent/CN109033015B/en
Publication of CN109033015A publication Critical patent/CN109033015A/en
Application granted granted Critical
Publication of CN109033015B publication Critical patent/CN109033015B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • 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/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Optical Integrated Circuits (AREA)

Abstract

A device for executing calculus operation on optical signals is composed of a signal source 1, a calculus operator 2 and a signal receiver 3. The method is characterized in that: the optical signal sent by the signal source 1 can directly irradiate on the calculus arithmetic unit 2, the calculus arithmetic unit 2 realizes the regulation and control of the electric field intensity amplitude and phase by utilizing the special physical property (graded refractive index) and ingenious structure of the material, and further realizes the first derivative, the second derivative and the integral operation of the input signal profile in the optical signal propagation process, and the signal profile sent by the signal source 1 and the signal profile received by the signal receiver 3 form a mathematical calculus relation.

Description

一种对光信号执行微积分运算的装置A device for performing calculus operations on optical signals

技术领域technical field

本装置利用材料的物理属性和特殊结构在光信号传播过程中实现对输入信号振幅和相位的控制,从而实现对输入光信号轮廓的微积分运算,属于超材料领域。The device utilizes the physical properties and special structure of the material to control the amplitude and phase of the input signal during the propagation of the optical signal, thereby realizing the calculus operation on the profile of the input optical signal, which belongs to the field of metamaterials.

背景技术Background technique

当前的微积分运算要靠微积分电路来实现。输出电压与输入电压成微分关系的电路为微分电路,通常由电容和电阻组成;输出电压与输入电压成积分关系的电路为积分电路,通常由电阻和电容组成。微积分运算电路广泛用于计算机、自动控制和电子仪器中。但电路存在需要外加能源、会发热、速度慢等缺点,为了适应信息化时代高速计算的需要,急需一种新的运算机制。The current calculus operation is realized by the calculus circuit. The circuit whose output voltage is in a differential relationship with the input voltage is a differential circuit, usually composed of capacitors and resistors; the circuit whose output voltage is in an integral relationship with the input voltage is an integral circuit, usually composed of resistors and capacitors. Calculus operation circuits are widely used in computers, automatic control and electronic instruments. However, the circuit has disadvantages such as needing external energy, heating, and slow speed. In order to meet the needs of high-speed computing in the information age, a new computing mechanism is urgently needed.

发明内容Contents of the invention

本发明的目的正是为了解决上述微积分运算电路的不足,利用渐变折射率材料的性质,通过巧妙的结构,实现对输入光信号轮廓的微积分运算。本装置具有原理新颖、结构简单、运算快速、无外加能源、没有中间过程等优点。The purpose of the present invention is to solve the deficiency of the above-mentioned calculus operation circuit, and to realize the calculus operation on the profile of the input optical signal by utilizing the properties of the graded-index material and through an ingenious structure. The device has the advantages of novel principle, simple structure, fast calculation, no external energy, no intermediate process and the like.

本发明是通过以下技术方案实现的:The present invention is achieved through the following technical solutions:

如图1所示,该装置由信号源1、微积分运算器2、信号接收器3组成。其特征在于:信号源1发出的光信号可以直接照射在微积分运算器2上,微积分运算器2在光信号传播过程中实现对输入信号的一阶导数、二阶导数及积分运算。该装置利用材料的物理属性(折射率)和结构尺寸实现对输入光信号振幅和相位的控制,从而实现微积分运算。As shown in Figure 1, the device consists of a signal source 1, a calculus operator 2, and a signal receiver 3. It is characterized in that the optical signal sent by the signal source 1 can directly irradiate the calculus operator 2, and the calculus operator 2 realizes the first-order derivative, second-order derivative and integral operation of the input signal during the propagation of the optical signal. The device uses the physical properties (refractive index) and structure size of the material to control the amplitude and phase of the input optical signal, thereby realizing the calculus operation.

本发明的有益效果是:该装置结构简单,信号处理速度快,无外加能源、适用范围广,在光子计算机、照片及图像存储、信号传递等方面有着广泛的应用。The beneficial effects of the invention are: the device has simple structure, fast signal processing speed, no external energy source, wide application range, and is widely used in photonic computer, photo and image storage, signal transmission and the like.

附图说明Description of drawings

图1是本发明一种对光信号执行微积分运算装置的结构示意图。FIG. 1 is a schematic structural diagram of a device for performing calculus operations on optical signals according to the present invention.

具体实施方式Detailed ways

如图1所示,该装置由信号源1、微积分运算器2、信号接收器3组成。信号源1发出的光信号可以直接照射在微积分运算器2上,微积分运算器2的材料属性和结构尺寸满足一定的条件,不改变其它参数,只调节微积分运算器2中材料II的折射率就可以实现对输入信号的一阶导数、二阶导数及积分运算,信号接收器3接收到的信号轮廓与信号源1发出的信号轮廓成数学微积分关系。As shown in Figure 1, the device consists of a signal source 1, a calculus operator 2, and a signal receiver 3. The optical signal sent by the signal source 1 can be directly irradiated on the calculus calculator 2. The material properties and structural dimensions of the calculus calculator 2 meet certain conditions, and only the material II in the calculus calculator 2 is adjusted without changing other parameters. The refractive index can realize the first-order derivative, second-order derivative and integral operation of the input signal, and the signal profile received by the signal receiver 3 and the signal profile sent by the signal source 1 form a mathematical calculus relationship.

理论依据Theoretical basis

如果光沿x轴正向在某一介质中传播,该介质的折射率在y方向按的规律变化,这种材料称为渐变折射率材料,其中η1为y=0处的折射率,η21=[π/(2Lg)]2,Lg为介质中传播的距离。当光波在具有渐变折射率的材料中传播时,渐变折射率材料将会在特征长度Lg上对电场进行傅里叶变换,利用该性质,可以设计具有微积分运算功能的超材料。如图1所示,微积分运算器2由三层材料(I,II和III)组成,设入射的光信号波长为λ,且沿x轴正方向传播,则三种材料x方向的尺寸分别为Lg=11.619λ、Δ=λ/3.012和Lg,y方向的尺寸均为W=9.876λ。坐标原点位于左侧中间。材料I和材料III的折射率满足ηI(y)=-ηIII(y)=η(y),分别执行傅里叶变换和逆傅里叶变换,材料II执行的变换是微积分运算器2的核心变换,该变换命名为G(y)。当入射光信号的电场强度大小E(y)输入到材料I的左端时,材料I右端输出的将是E(y)的傅里叶变换,即其中代表傅里叶变换,傅里叶变量ky与y属于同一区域,所以ky与y成正比。因为材料II执行的变换为G(y),所以材料II的输出函数为同时因为材料III执行逆傅里叶变换,所以最终出射端的函数应为若对输入的电场强度实现导数运算,则有If light travels in a medium along the positive direction of the x-axis, the refractive index of the medium in the y-direction is This kind of material is called graded index material, where η 1 is the refractive index at y=0, η 21 =[π/(2Lg)] 2 , and Lg is the propagation distance in the medium. When a light wave propagates in a material with a graded-refractive index, the graded-refractive-index material will perform a Fourier transform on the electric field at the characteristic length Lg. Using this property, metamaterials with calculus functions can be designed. As shown in Figure 1, the calculus operator 2 is composed of three layers of materials (I, II and III), and if the wavelength of the incident optical signal is λ, and it propagates along the positive direction of the x-axis, the dimensions of the three materials in the x-direction are respectively For Lg=11.619λ, Δ=λ/3.012 and Lg, the dimension in the y direction is W=9.876λ. The coordinate origin is in the left middle. The refractive index of material I and material III satisfy η I (y) = -η III (y) = η (y), perform Fourier transform and inverse Fourier transform respectively, and the transformation performed by material II is a calculus operator 2's core transformation, which is named G(y). When the electric field intensity E(y) of the incident light signal is input to the left end of material I, the output from the right end of material I will be the Fourier transform of E(y), namely in Represents the Fourier transform, the Fourier variable ky and y belong to the same area, so ky is proportional to y. Because the transformation performed by Material II is G(y), the output function of Material II is At the same time, because material III performs inverse Fourier transform, the function of the final exit end should be If the derivative operation is implemented for the input electric field strength, then there is

上式两边求傅里叶变换,有Finding the Fourier transform on both sides of the above formula, we have

根据傅里叶变换的微分性质上式可化为According to the differential properties of the Fourier transform The above formula can be transformed into

由方程(2)和(3)可求得材料II的变换函数G(y)=(iky)n∝(-iy)n(根据调试结果,取ky∝-y)。因材料沿y方向的尺寸为W,坐标原点又位于W中心,所以y方向的最大值为y0=W/2,即归一化的函数G(y)∝(-iy/y0)n。若要实现对光信号的微积分运算,还需求出材料II折射率的表达式,可以利用电磁波的一些性质,具体如下:The transformation function G(y)=(ik y ) n ∝(-iy) n of material II can be obtained from equations (2) and (3) (according to the debugging results, take k y ∝-y). Since the size of the material along the y direction is W, and the coordinate origin is located at the center of W, the maximum value in the y direction is y 0 =W/2, which is the normalized function G(y)∝(-iy/y 0 ) n . In order to realize the calculus operation of the optical signal, the expression of the refractive index of the material II is also required, and some properties of electromagnetic waves can be used, as follows:

当平面光波沿x轴正向传播时,x处的电场强度大小E(x)满足如下亥姆霍兹方程When a plane light wave propagates forward along the x-axis, the electric field intensity E(x) at x satisfies the following Helmholtz equation

其中,λ为入射波波长,εr,μr分别为介质的相对介电常数和相对磁导率。因为相对介电常数和相对磁导率与折射率之间满足关系εrμr=(η-iκ)2,其中η和κ分别代表折射率的实部和虚部,所以k也可以表示为k=(2π/λ)(η-iκ)。解方程(4),得in, λ is the wavelength of the incident wave, ε r and μ r are the relative permittivity and relative permeability of the medium, respectively. Because the relative permittivity and relative permeability satisfy the relationship ε r μ r = (η-iκ) 2 with the refractive index, where η and κ represent the real and imaginary parts of the refractive index respectively, so k can also be expressed as k=(2π/λ)(η-iκ). Solving equation (4), we get

E(x)=E0eikx      (5)E(x)=E 0 e ikx (5)

其中E0为x=0处的电场强度大小,E(x)为x处的电场强度大小。因为材料II左端的电场强度为根据方程(5),材料II中经过Δ距离的传播后,电场强度为因为则有Where E 0 is the magnitude of the electric field intensity at x=0, and E(x) is the magnitude of the electric field strength at x. Because the electric field strength at the left end of material II is According to Equation (5), after the propagation of Δ distance in material II, the electric field intensity is because then there is

Right now

为了使输出电场正比于输入电场的一阶导数,取G(y)=(-iy/y0),根据方程(7),材料II的折射率应满足In order to make the output electric field proportional to the first derivative of the input electric field, take G(y)=(-iy/y 0 ), according to equation (7), the refractive index of material II should satisfy

同理,若对输入函数进行二阶导数运算,则取G(y)=(-iy/y0)2,此时材料II的折射率应满足Similarly, if the second-order derivative operation is performed on the input function, G(y)=(-iy/y 0 ) 2 is taken, and the refractive index of material II should satisfy

Right now

其中n=1代表一阶导数运算,n=2代表二阶导数运算。Wherein, n=1 represents a first-order derivative operation, and n=2 represents a second-order derivative operation.

对于积分运算,(1)式应修改为For the integral operation, formula (1) should be modified as

上式两边求傅里叶变换,有Finding the Fourier transform on both sides of the above formula, we have

根据傅里叶变换的积分性质由方程(12)可得,执行积分运算的材料II的函数G(y)=(iky)∝(-iy)-1,设归一化常数为d,并取d=λ/4,则对于积分运算可取G(y)=(-iy/d)-1,为了避免|y|<d区域的透射系数大于1,将|y|=d设置为转折点,并假设|y|<d区域折射率的绝对值为常数,因此材料II可以进行积分运算的折射率表达式为According to the integral property of the Fourier transform From Equation (12), it can be obtained that the function G(y)=(ik y )∝(-iy) -1 of the material II performing the integral operation, setting the normalization constant as d, and taking d=λ/4, then For the integral operation, G(y)=(-iy/d) -1 can be taken. In order to avoid the transmission coefficient of |y|<d area being greater than 1, set |y|=d as the turning point, and assume that |y|<d area The absolute value of the refractive index is constant, so the refractive index expression of material II that can be integrated is

其中sign(.)为符号函数。Where sign(.) is a sign function.

以上已经公开了本发明设计方案,但并非用以限制本发明,凡采取同等思路和方法获得的其它微积分运算装置,均在本发明的保护范围内。The design scheme of the present invention has been disclosed above, but it is not intended to limit the present invention. All other calculus computing devices obtained by adopting equivalent ideas and methods are within the scope of protection of the present invention.

Claims (1)

1. The device for executing the calculus operation on the optical signal consists of a signal source (1), a calculus operator (2) and a signal receiver (3); the light signal emitted by the signal source (1) is directly irradiated on the calculus arithmetic unit (2), the calculus arithmetic unit (2) is formed by tightly connecting three layers of materials I, II and III, the dimensions of the three layers of materials in the x direction are lg= 11.619 λ, delta=λ/3.012 and lg= 11.619 λ respectively, the dimensions of the three layers of materials in the y direction are w= 9.876 λ respectively, and the structural parameters are based on the wavelength λ of the incident light signal; the material I and the material III are graded index materials, and the refractive index of the material II is different according to the different operation functions; the signal source (1), the calculus arithmetic unit (2) and the signal receiver (3) are directly connected by optical signals, and the whole device is in an air environment.
CN201710445240.6A 2017-06-09 2017-06-09 Device for executing calculus operation on optical signal Active CN109033015B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710445240.6A CN109033015B (en) 2017-06-09 2017-06-09 Device for executing calculus operation on optical signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710445240.6A CN109033015B (en) 2017-06-09 2017-06-09 Device for executing calculus operation on optical signal

Publications (2)

Publication Number Publication Date
CN109033015A CN109033015A (en) 2018-12-18
CN109033015B true CN109033015B (en) 2023-05-09

Family

ID=64630123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710445240.6A Active CN109033015B (en) 2017-06-09 2017-06-09 Device for executing calculus operation on optical signal

Country Status (1)

Country Link
CN (1) CN109033015B (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2797331B1 (en) * 1999-08-03 2002-07-19 Cit Alcatel OPTICAL DIFFERENTIATOR
CN104849878A (en) * 2015-06-03 2015-08-19 东南大学 Silicon nitride waveguide calorescence switch array chip based on Mach-Zahnder structure and production method thereof

Also Published As

Publication number Publication date
CN109033015A (en) 2018-12-18

Similar Documents

Publication Publication Date Title
Zangeneh-Nejad et al. Performing mathematical operations using high-index acoustic metamaterials
Gómez-Díaz et al. Graphene-based plasmonic switches at near infrared frequencies
Ma et al. An omnidirectional retroreflector based on the transmutation of dielectric singularities
Ju et al. Controllable and enhanced nanojet effects excited by surface plasmon polariton
Rejaei et al. Scattering of surface plasmons on graphene by a discontinuity in surface conductivity
Maksimov et al. Coupled mode theory for acoustic resonators
Li et al. An equivalent realization of coherent perfect absorption under single beam illumination
Lin et al. Transverse-electric Brewster effect enabled by nonmagnetic two-dimensional materials
Yu et al. Plasmonic interaction between nanospheres
Ren et al. High-order unified symplectic FDTD scheme for the metamaterials
Wu et al. The modern high frequency methods for solving electromagnetic scattering problems
Huang et al. Broadband radar absorbing characteristic based on periodic hollow truncated cone structure
Asgari et al. Tunable midinfrared wavelength selective structures based on resonator with antisymmetric parallel graphene pair
CN109033015B (en) Device for executing calculus operation on optical signal
Liu et al. Nonideal ultrathin mantle cloak for electrically large conducting cylinders
Jang et al. Perfectly matched layer formulation of the INBC-FDTD algorithm for electromagnetic analysis of thin film materials
Xu et al. Logarithm conformal mapping brings the cloaking effect
Kamandi et al. Integral equation analysis of multilayered waveguide bends using complex images Green's function technique
Karagounis et al. Full-wave simulation of optical waveguides via truncation in the method of moments using PML absorbing boundary conditions
Emadi et al. Use of an epsilon-near-zero region comprised of a graphene strip-silica stack for designing cloaking and reflection devices
Lu et al. Verification and improvement of equivalent refractive index models for evaluating the shielding effectiveness of high-transmittance double-layer metallic meshes
Belkhir et al. Implementation of the FDTD method in cylindrical coordinates for dispersive materials: Modal study of C-shaped nano-waveguides
Kim et al. Three-dimensional efficient dispersive alternating-direction-implicit finite-difference time-domain algorithm using a quadratic complex rational function
Tu et al. Single optical element to generate a meter-scale thz diffraction-free beam
Giloan et al. Optical imaging system based on transformation-optics lenses

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