CN109033015B - Device for executing calculus operation on optical signal - Google Patents
Device for executing calculus operation on optical signal Download PDFInfo
- 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
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 16
- 239000000463 material Substances 0.000 claims abstract description 36
- 230000005684 electric field Effects 0.000 abstract description 11
- 238000000034 method Methods 0.000 abstract description 5
- 230000000704 physical effect Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 description 4
- 239000003990 capacitor Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009699 differential effect Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
Images
Classifications
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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/11—Complex 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
Description
技术领域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
本发明的有益效果是:该装置结构简单,信号处理速度快,无外加能源、适用范围广,在光子计算机、照片及图像存储、信号传递等方面有着广泛的应用。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
理论依据Theoretical basis
如果光沿x轴正向在某一介质中传播,该介质的折射率在y方向按的规律变化,这种材料称为渐变折射率材料,其中η1为y=0处的折射率,η2/η1=[π/(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, η 2 /η 1 =[π/(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)
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)
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 |
-
2017
- 2017-06-09 CN CN201710445240.6A patent/CN109033015B/en active Active
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 |