Modulator photoelectric joint simulation method and device, computer equipment and storage medium

Abstract

The invention relates to a modulator photoelectric joint simulation method and device, computer equipment and a storage medium. The method comprises the following steps: acquiring an input radio frequency signal of a modulator; inputting the radio frequency signal into a radio frequency time domain coupling mode equation to obtain forward propagation radio frequency wave amplitude and reverse propagation radio frequency wave amplitude in the modulator; acquiring radio frequency voltage distribution in the modulator according to the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude; obtaining a first light wave amplitude and a second light wave amplitude in the modulator according to the radio frequency voltage distribution and a light wave time domain coupling mode equation; according to the first light wave amplitude and the second light wave amplitude, the light signal output by the modulator is obtained, and a simulated eye pattern is obtained according to the light signal. According to the invention, electrical simulation and optical simulation in the modulator can be realized at the same time.

Description

Technical Field

[0001] The present invention relates to the field of simulation, and in particular to a modulator optoelectronic joint simulation method, device, computer equipment and storage medium. Background Art

[0002] In the field of optical communication and optical information processing, modulators are key components for encoding signals onto light waves. Modulators can be made of a variety of materials, and simulation is crucial for the design and manufacture of modulators. The performance of modulators can be pre-evaluated based on eye diagrams before the modulator is manufactured, providing guidance for design improvements of modulators.

[0003] The modulator itself involves technical fields such as optics and circuits, so the simulation of the modulator requires electrical simulation and optical simulation, which is also called optoelectronic joint simulation.

[0004] In the prior art, simulation of modulators mainly relies on two methods, namely, segmented circuit method and transmission line circuit method, but these methods have significant defects.

[0005] In the segmented circuit method, optical devices are converted into equivalent circuit devices, and optical waveguides and electrodes are segmented into a set of equivalent circuit device arrays to model the propagation of RF waves and optical waves, and then simulated in software such as Verilog-A or SPICE to achieve optoelectronic joint simulation. However, due to the significant physical differences between current-voltage (IV) signals and wave signals, this method is difficult to handle the characteristics of wave signals such as multi-dimensional and bidirectional propagation. In the transmission line circuit method, electrodes are processed using the equivalent circuit model of the transmission line. For the dispersion of the electrode, additional equivalent circuit devices need to be introduced in the transmission line, and convolution calculations are performed in each step of the calculation. Each convolution requires access to all stored data, which is extremely time-consuming. For silicon optical modulators, it is also necessary to introduce the nonlinear IV equation of the PN junction, so the computational overhead and complexity are increased. Therefore, this method is very inefficient in handling the nonlinear IV relationship and dispersion of the electrode.

[0006] In addition, both of the above methods are material-dependent, as a specific circuit model needs to be established for each material, which lacks universality.

[0007] It can be seen that the optoelectronic simulation method in the prior art has low efficiency and high dependence on materials. Summary of the invention

[0008] In order to solve the above technical problem or at least partially solve the above technical problem, the present invention provides a modulator optoelectronic joint simulation, device, computer equipment and storage medium.

[0009] In a first aspect, the present invention provides a modulator optoelectronic co-simulation method, wherein the modulator receives a radio frequency signal and outputs an optical signal, and the method comprises:

[0010] Acquire a radio frequency signal input to the modulator;

[0011] Inputting the radio frequency signal into the radio frequency time domain coupled mode equation to obtain the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave in the modulator;

[0012] Acquiring a radio frequency voltage distribution in the modulator according to the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave;

[0013] Obtaining a first light wave amplitude and a second light wave amplitude in the modulator according to the RF voltage distribution and the light wave time domain coupling mode equation;

[0014] The optical signal output by the modulator is acquired according to the first optical wave amplitude and the second optical wave amplitude, and a simulated eye diagram is obtained according to the optical signal.

[0015] Optionally, before inputting the RF signal into the RF time-domain coupled mode equation to obtain the amplitude of the forward propagating RF wave and the amplitude of the reverse propagating RF wave in the modulator, the method further includes:

[0016] Get the initial RF equation,

[0017] Determine the nonlinear coefficients and boundary conditions of the initial RF equations,

[0018] Using the initial RF equation for determining the nonlinear coefficient as the RF time-domain coupled mode equation;

[0019] Before obtaining the first lightwave amplitude and the second lightwave amplitude in the modulator according to the RF voltage distribution and the lightwave time-domain coupling mode equation, the method further includes:

[0020] Get the initial light wave equation,

[0021] determining nonlinear coefficients and boundary conditions in the initial light wave equation,

[0022] The initial light wave equation for determining the nonlinear coefficient is used as the light wave time domain coupled mode equation.

[0023] Optionally, the initial radio frequency equation is:

[0024]

[0025]

[0026]

[0027]

[0028]

[0029]

[0030]

[0031] V rf，F (z, t) = A rf，F (z,t)exp(iβ rf z-iω rf t)+cc

[0032] V rf，R (z, t) = A rf，R (z,t)exp(-iβ rf z-iω rf t)+cc

[0033] Among them, α rf is the radio frequency wave attenuation constant, A rf，F (z, t) is the amplitude of the forward propagating RF wave, v g，rf is the RF wave group velocity, A rf，R (z, t) is the amplitude of the reverse propagating RF wave, V dc is the electrostatic field voltage, V rf，F (z, t) is the forward propagation RF wave voltage, V rf，R (z, t) is the reverse propagation RF wave voltage, β rf is the phase constant of the radio frequency wave, z is the spatial coordinate, t is the time, ω rf is the angular frequency of the radio frequency wave, cc represents the complex conjugate of the previous expression, is the second-order nonlinearity of the forward RF wave, is the third-order nonlinearity of the forward RF wave, is the fourth-order nonlinearity of the forward RF wave, is the second-order nonlinearity of the reflected RF wave, is the third-order nonlinearity of the reflected RF wave, is the fourth-order nonlinearity of the reflected RF wave, Λ dr，12 ,Λ r，4 ,Λ dr，22 ,Λ dr，22 ,Λ dr，14 and Λ dr，32 is the nonlinear coefficient.

[0034] Optionally, the initial light wave equation is:

[0035]

[0036]

[0037] Among them, α opt is the light wave attenuation constant, v g，opt is the light wave group velocity, A opt (z, t) is the light wave amplitude, is the second-order nonlinearity of the light wave, is the third-order nonlinearity of the light wave, V dc is the electrostatic field voltage, V rf，F is the forward propagating RF wave voltage, V rf，R is the reverse propagation RF wave voltage, ε0 is the vacuum dielectric constant, and h is The power index, j is The power index, k is The power index of is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, Λ dp，12 ,Λ rp，12 ,Λ dp，22 ,Λ rdp，112 and Λ rp，22 is the nonlinear coefficient.

[0038] Optionally, determining the nonlinear coefficients and boundary conditions of the initial radio frequency equation, and determining the nonlinear coefficients and boundary conditions in the initial light wave equation, are performed in the following manner:

[0039]

[0040]

[0041] Λ dp，12 =Λ rp，12

[0042] Λ dp，22 =Λ rp，22 =Λ drp，112

[0043] Optionally, the boundary conditions are as follows:

[0044] A rf，R (L, t) = Γ L A rf，F (L, t)exp(2iβ rf L)

[0045]

[0046]

[0047]

[0048] Λ drp，hjk is the nonlinear coefficient, is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, h is The power index, j is The power index, k is The power index of h, j and k are positive integers and h+j+k-1≥2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δn eff，rf is the effective refractive index of radio frequency waves, k 0，rf is the vacuum wave vector of the radio frequency wave, A rf，R (L, t) is A at the coordinate (L, t) rf，R , A rf，F (L, t) is A at the coordinate (L, t) rf，F , β rf is the RF wave phase constant, L is the length of the modulator traveling wave region, ω rf is the angular frequency of the RF wave, V in (t) is the RF voltage signal input to the modulator, Z c is the impedance of the equivalent RF waveguide of the electrode, R g is the internal resistance of the equivalent RF waveguide of the electrode, R L is the load resistance of the equivalent RF waveguide of the electrode, Λ dp，12 ,Λ rp，12 ,Λ dp，22 ,Λ rdp，112 ,Λ rp，22 ,Λ drp，hjk ,Λ dr，12 ,Λ dr，22 and Λ dr，32 is the nonlinear coefficient.

[0049] Optionally, the RF signal is input into a RF time-domain coupled mode equation to obtain the amplitude of the forward propagating RF wave and the amplitude of the reverse propagating RF wave in the modulator in the following manner:

[0050]

[0051]

[0052] Δz rf =Δt / β 1，rf

[0053] Δz opt=Δt / β 1，opt

[0054]

[0055]

[0056] Δt is the time interval after the time t is discretized, α rf is the RF wave attenuation constant, V dc is the electrostatic field voltage, β 1，rf and β 1，opt is the first derivative of the propagation constant, v g，rf is the group velocity of RF, v g，opt is the group velocity of the light wave, Λ dr，12 ,Λ r，4 ,Λ dr，22 ,Λ dr，32 and Λ dr，14 is the nonlinear coefficient.

[0057] Optionally, the first lightwave amplitude and the second lightwave amplitude in the modulator are obtained according to the RF voltage distribution and the lightwave time-domain coupling mode equation in the following manner:

[0058] A opt，1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt，1 (z, t)

[0059] A opt，2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt，2 (z, t)

[0060] Δφ1(z, t)=-α opt Δz opt +2iΛ dp，12 (V rf，F (z,t)+V rf，R (z,t)+V dc )Δz opt +3iΛ dp，22 (V rf，F (z,t)+V rf，R (z,t)+V dc ) 2 Δz opt

[0061] Δφ2(z, t)=-α opt Δz opt +2iΛ dp，12 (-V rf，F (z, t)-V rf，R(z,t)+V dc )Δz opt

[0062] Among them, A opt，1 (z, t) is the amplitude of the first light wave, A opt，2 (z, t) is the amplitude of the second light wave, Δz opt represents the discretization interval of the spatial coordinate z, α opt is the light wave attenuation constant, V rf，F (z, t) is the forward propagation RF wave voltage, V rf，R (z, t) is the reverse propagation RF wave voltage, V dc is the electrostatic field voltage, Λ dp，12 ,Λ dp，22 and Λ dp，12 is the nonlinear coefficient.

[0063] In a second aspect, a modulator optoelectronic co-simulation device is provided, wherein the modulator receives a radio frequency signal and outputs an optical signal, and the device comprises:

[0064] A signal acquisition unit, used to acquire the input radio frequency signal of the modulator;

[0065] A radio frequency simulation unit, used for inputting the radio frequency signal into a radio frequency time domain coupled mode equation to obtain the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave in the modulator;

[0066] The RF simulation unit is further used to obtain the RF voltage distribution in the modulator according to the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude;

[0067] A light wave simulation unit, used for obtaining the first light wave amplitude and the second light wave amplitude in the modulator according to the RF voltage distribution and the light wave time domain coupling mode equation;

[0068] The eye diagram generating unit is used to obtain the optical signal output by the modulator according to the first optical wave amplitude and the second optical wave amplitude, and obtain a simulated eye diagram according to the optical signal.

[0069] According to a third aspect, a computer device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements any of the above methods when executing the computer program.

[0070] In a fourth aspect, a computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the method as described in any one of the above items is implemented.

[0071] The present invention provides a modulator optoelectronic joint simulation, device, computer equipment and storage medium, wherein the modulator receives a radio frequency signal and outputs an optical signal, and the method comprises: obtaining the input radio frequency signal of the modulator; inputting the radio frequency signal into the radio frequency time domain coupling mode equation to obtain the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude in the modulator; obtaining the radio frequency voltage distribution in the modulator according to the forward propagation radio frequency wave amplitude and the reverse propagation radio frequency wave amplitude; obtaining the first light wave amplitude and the second light wave amplitude in the modulator according to the radio frequency voltage distribution and the optical wave time domain coupling mode equation; obtaining the optical signal output by the modulator according to the first light wave amplitude and the second light wave amplitude, and obtaining a simulated eye diagram according to the optical signal. The method of the embodiment of the present invention realizes electrical simulation and optical simulation at the same time in the design process of the modulator, thereby realizing efficient and accurate design and simulation of the modulator, and promoting the development of optical communication and optical computing systems. In addition, the method of the embodiment of the present invention does not need to rely on specific materials or establish specific circuit models, and has wider applicability and higher efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0072] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0073] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0074] Figure 1 The figure shows an application environment diagram of the modulator optoelectronic joint simulation method according to an embodiment of the present invention;

[0075] Figure 2 FIG. 1 is a flow chart of a modulator optoelectronic co-simulation method according to an embodiment of the present invention;

[0076] Figure 3 FIG. 1 is a schematic diagram of a radio frequency circuit model according to an embodiment of the present invention;

[0077] Figure 4 The calculation A of the embodiment of the present invention is shown. rf，F and A rf，R The numerical calculation process of the RF wave time domain coupled mode equation;

[0078] Figure 5 The calculation A of the embodiment of the present invention is shown. opt The numerical calculation process of the lightwave time-domain coupled mode equation;

[0079] Figure 6 The figure shows the comparison between the eye diagram obtained by the simulation method of the embodiment of the present invention and the eye diagram measured by the experiment.

[0080] Figure 7 It is a structural block diagram of a modulator optoelectronic joint simulation device according to an embodiment of the present invention;

[0081] Figure 8 1 is a diagram showing the internal structure of a computer device in an embodiment of the present invention. DETAILED DESCRIPTION

[0082] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0083] Figure 1 FIG. 1 is an application environment diagram of a modulator optoelectronic joint simulation method in an embodiment. Figure 1 The modulator optoelectronic co-simulation method is applied to a modulator optoelectronic co-simulation system. The modulator optoelectronic co-simulation method includes a terminal 110 and / or a server 120. The terminal 110 and the server 120 are connected via a network. The terminal 110 may be a desktop terminal or a mobile terminal, and the mobile terminal may be at least one of a mobile phone, a tablet computer, a laptop computer, etc. The server 120 may be implemented as an independent server or a server cluster consisting of multiple servers.

[0084] like Figure 2 As shown, in one embodiment, a modulator optoelectronic co-simulation method is provided. This embodiment mainly applies this method to the above Figure 1 The server 120 in FIG. 1 is used as an example. The modulator receives a radio frequency signal and outputs an optical signal. The method of the embodiment of the present invention is used to simulate the output optical signal of the modulator.

[0085] Reference Figure 2 , the modulator optoelectronic joint simulation method includes:

[0086] Step 210, obtaining an input radio frequency signal of the modulator;

[0087] Step 220, inputting the RF signal into the RF time domain coupled mode equation to obtain the amplitude of the forward propagating RF wave and the amplitude of the reverse propagating RF wave in the modulator;

[0088] Step 230, obtaining a radio frequency voltage distribution in the modulator according to the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave;

[0089] Step 240, obtaining a first lightwave amplitude and a second lightwave amplitude in the modulator according to the RF voltage distribution and the lightwave time-domain coupling mode equation;

[0090] Step 250: Acquire the optical signal output by the modulator according to the first optical wave amplitude and the second optical wave amplitude, and obtain a simulated eye diagram according to the optical signal.

[0091] Time-domain coupled-mode equations (TCME for short), in the embodiment of the present invention, the time-domain coupled-mode equations include radio frequency time-domain coupled-mode equations and light wave time-domain coupled-mode equations.

[0092] The modulator usually includes two parts: a radio frequency electrode and an optical waveguide. In an embodiment of the present invention, the optical waveguide is designed as a Mach-Zehnder interferometer (MZI for short). Due to the phase difference between the two arms of the MZI, the amplitude of the optical signal can be modulated according to the interference of light. In order to generate the phase difference, the refractive index of the optical waveguide of the two arms needs to change. Specifically, the silicon optical modulator is to apply a voltage between the G and S electrodes, and the voltage on the PN junction changes the distribution of electrons and holes, and the refractive index of the optical waveguide changes according to the plasma dispersion effect. These PN junctions are arranged in a lateral push-pull configuration, thereby halving the capacitance contributed by the PN junction. The PN junction is formed by doping N-type and P-type regions in the optical waveguide. The N++ and P++ doped regions realize the electrical connection between the electrode and the PN junction. Similarly, the lithium niobate modulator is to apply a voltage between the G and S electrodes, and the refractive index of the optical waveguide changes according to the Pockels effect of the lithium niobate material. In an embodiment of the present invention, the time domain coupled mode equation (TCME) is derived from Maxwell's equations to describe the nonlinear effects of radio frequency waves and optical waves. Secondly, the nonlinear coefficients in TCME are obtained by solving and fitting the eigenmodes of RF and optical waves. Finally, the differential form of TCME on a two-dimensional time-space grid is introduced to realize the numerical calculation of TCME.

[0093] In the method of the embodiment of the present invention, the input RF signal of the modulator is obtained; the RF signal is input into the RF time-domain coupling mode equation to obtain the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude in the modulator; the RF voltage distribution in the modulator is obtained according to the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; the first light wave amplitude and the second light wave amplitude in the modulator are obtained according to the RF voltage distribution and the light wave time-domain coupling mode equation; the optical signal output by the modulator is obtained according to the first light wave amplitude and the second light wave amplitude, and a simulated eye diagram is obtained according to the optical signal. The method of the embodiment of the present invention realizes electrical simulation and optical simulation at the same time in the design process of the modulator, so that efficient and accurate design and simulation of the modulator can be realized, and the development of optical communication and optical computing systems can be promoted. In addition, the method of the embodiment of the present invention does not need to rely on specific materials or establish specific circuit models, and has wider applicability and higher efficiency.

[0094] In the embodiment of the present invention, before step 220, before inputting the RF signal into the RF time domain coupled mode equation to obtain the amplitude of the forward propagating RF wave and the amplitude of the reverse propagating RF wave in the modulator, the method further includes:

[0095] Get the initial RF equation,

[0096] Determine the nonlinear coefficients and boundary conditions of the initial RF equations,

[0097] Using the initial RF equation for determining the nonlinear coefficient as the RF time-domain coupled mode equation;

[0098] The method of the embodiment of the present invention, before step 240, before obtaining the first lightwave amplitude and the second lightwave amplitude in the modulator according to the RF voltage distribution and the lightwave time-domain coupling mode equation, the method further includes:

[0099] Get the initial light wave equation,

[0100] determining nonlinear coefficients and boundary conditions in the initial light wave equation,

[0101] The initial light wave equation for determining the nonlinear coefficient is used as the light wave time domain coupled mode equation.

[0102] In the embodiment of the present invention, the initial radio frequency equation is:

[0103]

[0104]

[0105]

[0106]

[0107]

[0108]

[0109]

[0110]

[0111] V rf，F (z, t) = A rf，F (z,t)exp(iβ rf z-iω rf t)+cc

[0112] V rf，R (z, t) = A rf，R (z,t)exp(-iβ rf z-iω rf t)+cc

[0113] Among them, α rf is the radio frequency wave attenuation constant, A rf，F (z, t) is the amplitude of the forward propagating RF wave, v g，rf is the RF wave group velocity, A rf，R (z, t) is the amplitude of the reverse propagating RF wave, V dc is the electrostatic field voltage, V rf，F (z, t) is the forward propagation RF wave voltage, V rf，R (z, t) is the reverse propagation RF wave voltage, β rf is the phase constant of the radio frequency wave, z is the spatial coordinate, t is the time, ω rf is the angular frequency of the radio frequency wave, cc represents the complex conjugate of the previous expression, is the second-order nonlinearity of the forward RF wave, is the third-order nonlinearity of the forward RF wave, is the fourth-order nonlinearity of the forward RF wave, is the second-order nonlinearity of the reflected RF wave, is the third-order nonlinearity of the reflected RF wave, is the fourth-order nonlinearity of the reflected RF wave, Λ dr，12 ,Λ r，4 ,Λ dr，22 ,Λ dr，22 ,Λ dr，14 and Λ dr，32 is the nonlinear coefficient.

[0114] The electromagnetic field in the modulator can be approximately considered to consist only of the static voltage Vdc The static electric field distribution E caused by dc , light wave amplitude A opt , Mode field distribution E opt , the amplitude of the forward propagating RF wave A rf，F And the mode field distribution E rf , the amplitude of the reverse propagating RF wave A rf，R The RF mode field distribution of the reverse propagating RF wave is the same as the RF mode field distribution of the forward propagating RF wave. In the embodiment of the present invention, the time domain coupled mode equation describing the nonlinear effect in the modulator can be derived from the Maxwell equations, and the time domain coupled mode equation without confirming the nonlinear parameters can be the initial RF equation.

[0115] In the embodiment of the present invention, the initial light wave equation is:

[0116]

[0117]

[0118]

[0119] Among them, α opt is the light wave attenuation constant, v g，opt is the light wave group velocity, A opt (z, t) is the light wave amplitude, is the second-order nonlinearity of the light wave, is the third-order nonlinearity of the light wave, V dc is the electrostatic field voltage, V rf，F is the forward propagating RF wave voltage, V rf，R is the reverse propagation RF wave voltage, ε0 is the vacuum dielectric constant, and h is The power index, j is The power index, k is The power index of is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, Λ dp，12 ,Λ rp，12 ,Λ dp，22 ,Λ rdp，112 and Λ rp，22 is the nonlinear coefficient.

[0120] In the embodiment of the present invention, the determining of the nonlinear coefficients and boundary conditions of the initial radio frequency equation and the determining of the nonlinear coefficients and boundary conditions in the initial light wave equation are performed in the following manner:

[0121]

[0122]

[0123] Λ dp，12 =Λ rp，12

[0124] Λ dp，22 =Λ rp，22 =Λ drp，112

[0125] The boundary conditions are as follows:

[0126] A rf，R (L, t) = Γ L A rf，F (L, t)exp(2iβ rf L)

[0127]

[0128]

[0129]

[0130] Λ drp，hjk is the nonlinear coefficient, is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, h is The power index, j is The power index, k is The power index of h, j and k are positive integers and h+j+k-1≥2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δn eff，rf is the effective refractive index of radio frequency waves, k 0，rf is the vacuum wave vector of the radio frequency wave, A rf，R (L, t) is A at the coordinate (L, t) rf，R , A rf，F (L, t) is A at the coordinate (L, t) rf，F , β rf is the RF wave phase constant, L is the length of the modulator traveling wave region, ω rf is the angular frequency of the RF wave, V in (t) is the RF voltage signal input to the modulator, Z c is the impedance of the equivalent RF waveguide of the electrode, R g is the internal resistance of the equivalent RF waveguide of the electrode, R L is the load resistance of the equivalent RF waveguide of the electrode, Λ dp，12 ,Λ rp，12 ,Λdp，22 ,Λ rdp，112 ,Λ rp，22 ,Λ drp，hjk ,Λ dr，12 ,Λ dr，22 and Λ dr，32 is the nonlinear coefficient.

[0131] The boundary conditions in the above-mentioned embodiment of the present invention are only one method. There are other methods in other embodiments of the present invention, which will not be described in detail here.

[0132] In the embodiment of the present invention, the nonlinear coefficient can be calculated based on the change Δn of the effective refractive index of the radio frequency wave and the optical wave. eff，rf or Δn eff，opt With the static voltage V dc The specific calculation method is as described above. Compared with the definition of the nonlinear coefficient, this method is simpler and easy to combine with experimental measurement. In addition, some materials lack measured material parameters χ, such as PN junctions. The method of the embodiment of the present invention can circumvent this problem and is therefore more universal.

[0133] Figure 3 FIG. 1 is a schematic diagram of a radio frequency circuit model according to an embodiment of the present invention. Figure 3 As shown, based on this model, the boundary conditions can be derived. The electrodes are considered to have impedance Z c The RF waveguide is connected to a load resistor R L and an internal resistor R g The boundary conditions describe the reflection and transmission of RF waves at the interfaces at both ends of the RF waveguide.

[0134] In the embodiment of the present invention, in step 220, the RF signal is input into the RF time domain coupled mode equation to obtain the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude in the modulator in the following manner:

[0135]

[0136]

[0137] Δz rf =Δt / β 1，rf

[0138] Δz opt =Δt / β 1，opt

[0139]

[0140]

[0141] Δt is the time interval after the time t is discretized, α rf is the RF wave attenuation constant, V dc is the electrostatic field voltage, β 1，rf and β 1，opt is the first derivative of the propagation constant, v g，rf is the group velocity of RF, v g，opt is the group velocity of the light wave, Λ dr，12 ,Λ r，4 ,Λ dr，22 ,Λ dr，32 and Λ dr，14 is the nonlinear coefficient.

[0142] In the embodiment of the present invention, in step 240, the first lightwave amplitude and the second lightwave amplitude in the modulator are obtained according to the RF voltage distribution and the lightwave time-domain coupling mode equation in the following manner:

[0143] A opt，1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt，1 (z, t)

[0144] A opt，2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt，2 (z, t)

[0145] Δφ1(z, t)=-α opt Δz opt +2iΛ dp，12 (V rf，F (z,t)+V rf，R (z,t)+V dc )Δz opt +3iΛ dp，22 (V rf，F (z,t)+V rf，R (z,t)+V dc ) 2 Δz opt

[0146] Δφ2(z, t)=-α opt Δz opt +2iΛ dp，12 (-V rf，F (z, t)-V rf，R (z,t)+V dc )Δz opt

[0147] Among them, A opt，1 (z, t) is the amplitude of the first light wave, A opt，2(z, t) is the amplitude of the second light wave, Δz opt represents the discretization interval of the spatial coordinate z, α opt is the light wave attenuation constant, V rf，F (z, t) is the forward propagation RF wave voltage, V rf，R (z, t) is the reverse propagation RF wave voltage, V dc is the electrostatic field voltage, Λ dp，12 ,Λ dp，22 and Λ dp，12 is the nonlinear coefficient.

[0148] In the embodiment of the present invention, when solving the problem based on TCME, the difference form of the two-dimensional space-time grid is relied on. The TCME of radio frequency waves and light waves have their own two-dimensional space-time grids, where the grid size of t is the same for radio frequency waves and light waves. The grid size of z for radio frequency waves and light waves is different, which is Δz rf =Δt / β 1，rf and Δz opt =Δt / β 1，opt .

[0149] Figure 4 Demonstrated calculation A rf，F and A rf，R The numerical calculation process of TCME of radio frequency waves. The calculation is based on the spatial interval Δz rf and a two-dimensional space-time grid with a time interval of Δt, where each circle represents a storage A rf，F and A rf，R The initial value of all circles is zero. Black or gray arrows indicate the nodes where A is calculated. rf，F or A rf，R The boundary conditions are indicated by green arrows. The orange arrow marks the input voltage signal V in The algorithm starts with calculating A rf，F First, according to the boundary conditions on the left side of the RF circuit, V in Due to impedance mismatch, according to the boundary conditions on the right side of the RF circuit, A rf，R A rf，F In the simulation, Δt = 0.5ps was used to ensure numerical convergence. Finally, V on the grid was obtained according to the forward and reverse propagating RF wave voltages. rf，F and V rf，R distribution.

[0150] Figure 5 Demonstrated calculation A opt The numerical calculation process of optical wave TCME. The input optical wave A opt，in It is divided into two paths, namely the two arms of MZI (Mach-Zehnder interferometer, MZI for short), Aopt，1 and A opt，2 The black arrow indicates the calculation of A opt，1 and A opt，2 The dashed green arrows mark the nonlinear terms in the light-wave TCME. Figure 5 V rf V rf，F 、V rf，R and V dc . Due to Δz opt ≠Δz rf , we use interpolation to resample V along z = 0 to z = L at t rf Considering that optical waves oscillate faster than RF waves, we do not use a finer grid size Δz opt , but A opt，1 and A opt，2 The differential form is converted into the exponential differential form. The final modulator output light wave is obtained by coherent superposition A opt，1 and A opt，2 Get an A opt，out .

[0151] In the embodiment of the present invention, the method for obtaining the eye diagram may adopt the method in the prior art, which will not be described in detail here.

[0152] Figure 6 FIG. 4 is a comparison of the eye diagram obtained by the simulation method of the embodiment of the present invention and the eye diagram measured by the experiment. It can be seen that the eye diagram obtained by the simulation method of the embodiment of the present invention is as follows: Figure 6 (a) The eye diagram measured experimentally is shown in Figure 6 (b) Very close.

[0153] Figure 6 In the figure, the red and blue colors in the eye diagram are due to the device or system settings and do not affect the comparative structure of the eye diagram.

[0154] The method of the embodiment of the present invention can realize electrical simulation and optical simulation at the same time, so as to realize efficient and accurate design and simulation of the modulator, and promote the development of optical communication and optical computing systems. In addition, in the method of the present invention, the RF electrode is converted into an RF waveguide for processing instead of an equivalent circuit model, the silicon PN junction is regarded as an equivalent fourth-order nonlinear medium, and the nonlinear IV relationship of the PN junction is replaced by the self-phase modulation (SPM) and cross-phase modulation (XPM) of the RF wave. In the method of the embodiment of the present invention, the dispersion of the RF electrode is automatically solved by using the group velocity of the RF wave in the RF waveguide, avoiding the convolution calculation of the previous method.

[0155] like Figure 7 As shown, the present invention also provides a modulator optoelectronic joint simulation device, the device comprising:

[0156] A signal acquisition unit 710, configured to acquire a radio frequency signal input to the modulator;

[0157] The RF simulation unit 720 is used to input the RF signal into the RF time domain coupled mode equation to obtain the amplitude of the forward propagating RF wave and the amplitude of the reverse propagating RF wave in the modulator;

[0158] The RF simulation unit 720 is further used to obtain the RF voltage distribution in the modulator according to the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude;

[0159] A lightwave simulation unit 730, configured to obtain a first lightwave amplitude and a second lightwave amplitude in the modulator according to the RF voltage distribution and the lightwave time-domain coupling mode equation;

[0160] The eye diagram generating unit 740 is used to obtain the optical signal output by the modulator according to the first optical wave amplitude and the second optical wave amplitude, and obtain a simulated eye diagram according to the optical signal. In the embodiment of the present invention, the RF simulation unit 720 is also used to:

[0161] Get the initial RF equation,

[0162] Determine the nonlinear coefficients and boundary conditions of the initial RF equations,

[0163] The initial radio frequency equation for determining the nonlinear coefficient is used as the radio frequency time-domain coupled mode equation.

[0164] In the embodiment of the present invention, the initial radio frequency equation is:

[0165]

[0166]

[0167]

[0168]

[0169]

[0170]

[0171]

[0172]

[0173] V rf，F (z, t) = A rf，F (z,t)exp(iβ rf z-iω rft)+cc

[0174] V rf，R (z, t) = A rf，R (z,t)exp(-iβ rf z-iω rf t)+cc

[0175] Among them, α rf is the radio frequency wave attenuation constant, A rf，F (z, t) is the amplitude of the forward propagating RF wave, v g，rf is the RF wave group velocity, A rf，R (z, t) is the amplitude of the reverse propagating RF wave, V dc is the electrostatic field voltage, V rf，F (z, t) is the forward propagation RF wave voltage, V rf，R (z, t) is the reverse propagation RF wave voltage, β rf is the phase constant of the radio frequency wave, z is the spatial coordinate, t is the time, ω rf is the angular frequency of the radio frequency wave, cc represents the complex conjugate of the previous expression, is the second-order nonlinearity of the forward RF wave, is the third-order nonlinearity of the forward RF wave, is the fourth-order nonlinearity of the forward RF wave, is the second-order nonlinearity of the reflected RF wave, is the third-order nonlinearity of the reflected RF wave, is the fourth-order nonlinearity of the reflected RF wave, Λ dr，12 ,Λ r，4 ,Λ dr，22 ,Λ dr，22 ,Λ dr，14 and Λ dr，32 is the nonlinear coefficient.

[0176] In the embodiment of the present invention, the light wave simulation unit 730 is further used for:

[0177] Obtain the initial light wave equation;

[0178] Determining nonlinear coefficients and boundary conditions in the initial light wave equation;

[0179] The initial light wave equation for determining the nonlinear coefficient is used as the light wave time domain coupled mode equation.

[0180] In the embodiment of the present invention, the initial light wave equation is:

[0181]

[0182]

[0183]

[0184] Among them, α opt is the light wave attenuation constant, v g，opt is the light wave group velocity, A opt (z, t) is the light wave amplitude, is the second-order nonlinearity of the light wave, is the third-order nonlinearity of the light wave, V dc is the electrostatic field voltage, V rf，F is the forward propagating RF wave voltage, V rf，R is the reverse propagation RF wave voltage, ε0 is the vacuum dielectric constant, and h is The power index, j is The power index, k is The power index of is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, Λ dp，12 ,Λ rp，12 ,Λ dp，22 ,Λ rdp，112 and Λ rp，22 is the nonlinear coefficient.

[0185] In the embodiment of the present invention, the determining of the nonlinear coefficients and boundary conditions of the initial radio frequency equation and the determining of the nonlinear coefficients and boundary conditions in the initial light wave equation are performed in the following manner:

[0186]

[0187]

[0188] Λ dp，12 =Λ rp，12

[0189] Λ dp，22 =Λ rp，22 =Λ drp，112

[0190] The boundary conditions are as follows:

[0191] A rf，R (L, t) = Γ L A rf，F (L, t)exp(2iβ rf L)

[0192]

[0193]

[0194]

[0195] Λ drp，hjk is the nonlinear coefficient, is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, h is The power index, j is The power index, k is The power index of h, j and k are positive integers and h+j+k-1≥2, and d is used to refer to r is used to refer to p is used to refer to , χ is the nonlinear polarizability of the material, Δn eff，rf is the effective refractive index of radio frequency waves, k 0，rf is the vacuum wave vector of the radio frequency wave, A rf，R (L, t) is A at the coordinate (L, t) rf，R , A rf，F (L, t) is A at the coordinate (L, t) rf，F , β rf is the RF wave phase constant, L is the length of the modulator traveling wave region, ω rf is the angular frequency of the RF wave, V in (t) is the RF voltage signal input to the modulator, Z c is the impedance of the equivalent RF waveguide of the electrode, R g is the internal resistance of the equivalent RF waveguide of the electrode, R L is the load resistance of the equivalent RF waveguide of the electrode, Λ dp，12 ,Λ rp，12 ,Λ dp，22 ,Λ rdp，112 ,Λ rp，22 ,Λ drp，hjk ,Λ dr，12 ,Λ dr，22 and Λ dr，32 is the nonlinear coefficient.

[0196] In the embodiment of the present invention, the radio frequency simulation unit 720 is further used to obtain the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave in the following manner:

[0197]

[0198]

[0199] Δz rf =Δt / β 1，rf

[0200] Δz opt =Δt / β 1，opt

[0201]

[0202]

[0203] Δt is the time interval after the time t is discretized, α rf is the RF wave attenuation constant, V dc is the electrostatic field voltage, β 1，rf and β 1，opt is the first derivative of the propagation constant, v g，rf is the group velocity of RF, v g，opt is the group velocity of the light wave, Λ dr，12 ,Λ r，4 ,Λ dr，22 ,Λ dr，32 and Λ dr，14 is the nonlinear coefficient.

[0204] In the embodiment of the present invention, the light wave simulation unit 730 is further used to obtain the first light wave amplitude and the second light wave amplitude in the following manner:

[0205] A opt，1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt，1 (z, t)

[0206] A opt，2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt，2 (z, t)

[0207] Δφ1(z, t)=-α opt Δz opt +2iΛ dp，12 (V rf，F (z,t)+V rf，R (z,t)+V dc )Δz opt +3iΛ dp，22 (V rf , F(z, t)+V rf，R (z,t)+V dc ) 2 Δz opt

[0208] Δφ2(z, t)=-α opt Δz opt +2iΛ d，12 (-V rf，F (z, t)-V rf，R (z,t)+V dc )Δzopt

[0209] Among them, A opt，1 (z, t) is the amplitude of the first light wave, A opt，2 (z, t) is the amplitude of the second light wave, Δz opt represents the discretization interval of the spatial coordinate z, α opt is the light wave attenuation constant, V rf，F (z, t) is the forward propagation RF wave voltage, V rf，R (z, t) is the reverse propagation RF wave voltage, V dc is the electrostatic field voltage, Λ dp，12 ,Λ dp，22 and Λ dp，12 is the nonlinear coefficient.

[0210] The method of the embodiment of the present invention realizes electrical simulation and optical simulation simultaneously in the design process of the modulator, thereby realizing efficient and accurate design and simulation of the modulator, and promoting the development of optical communication and optical computing systems. In addition, the method of the embodiment of the present invention does not need to rely on specific materials, nor does it need to establish a specific circuit model, and has wider applicability and higher efficiency.

[0211] An embodiment of the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the following method when executing the computer program: obtaining an input RF signal of the modulator; inputting the RF signal into an RF time-domain coupling mode equation to obtain a forward propagation RF wave amplitude and a reverse propagation RF wave amplitude in the modulator; obtaining an RF voltage distribution in the modulator based on the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; obtaining a first light wave amplitude and a second light wave amplitude in the modulator based on the RF voltage distribution and the light wave time-domain coupling mode equation; obtaining the optical signal output by the modulator based on the first light wave amplitude and the second light wave amplitude, and obtaining a simulated eye diagram based on the optical signal.

[0212] An embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program implements the following method when executed by a processor: obtaining an input RF signal of the modulator; inputting the RF signal into an RF time-domain coupling mode equation to obtain a forward propagation RF wave amplitude and a reverse propagation RF wave amplitude in the modulator; obtaining an RF voltage distribution in the modulator based on the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; obtaining a first light wave amplitude and a second light wave amplitude in the modulator based on the RF voltage distribution and the light wave time-domain coupling mode equation; obtaining the optical signal output by the modulator based on the first light wave amplitude and the second light wave amplitude, and obtaining a simulated eye diagram based on the optical signal.

[0213] The above modulator optoelectronic joint simulation method achieves the beneficial effect of being able to solve the technical problems raised in the background technology.

[0214] Figure 2 FIG. 1 is a flow chart of a modulator optoelectronic co-simulation method in one embodiment. It should be understood that although Figure 2 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 2 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0215] Figure 8 The internal structure diagram of a computer device in one embodiment is shown. The computer device may specifically be Figure 1 The server 120 in FIG. Figure 8As shown, the computer device includes a processor, a memory, a network interface, an input device and a display screen connected through a system bus. Among them, the memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores an operating system and may also store a computer program. When the computer program is executed by the processor, the processor can implement the modulator optoelectronic co-simulation method. The internal memory may also store a computer program. When the computer program is executed by the processor, the processor can execute the modulator optoelectronic co-simulation method. The display screen of the computer device may be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device may be a touch layer covered on the display screen, or a button, trackball or touchpad provided on the housing of the computer device, or an external keyboard, touchpad or mouse, etc.

[0216] Those skilled in the art will understand that Figure 8 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present invention, and does not constitute a limitation on the computer device to which the solution of the present invention is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0217] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be completed by instructing the relevant hardware through a computer program, and the program can be stored in a non-volatile computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0218] It should be noted that, in this article, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0219] The foregoing is merely a specific embodiment of the present invention, which enables those skilled in the art to understand or implement the present invention. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but rather to the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A modulator optoelectronic joint simulation method, characterized in that: The modulator receives a radio frequency signal and outputs an optical signal, and the method comprises: Acquire a radio frequency signal input to the modulator; Inputting the radio frequency signal into the radio frequency time domain coupled mode equation to obtain the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave in the modulator; Acquiring a radio frequency voltage distribution in the modulator according to the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave; Obtaining a first light wave amplitude and a second light wave amplitude in the modulator according to the RF voltage distribution and the light wave time domain coupling mode equation; The optical signal output by the modulator is acquired according to the first optical wave amplitude and the second optical wave amplitude, and a simulated eye diagram is obtained according to the optical signal.

2. The method according to claim 1, characterized in that Before inputting the RF signal into the RF time-domain coupled mode equation to obtain the amplitude of the forward propagating RF wave and the amplitude of the reverse propagating RF wave in the modulator, the method further includes: Get the initial RF equation, Determine the nonlinear coefficients and boundary conditions of the initial RF equations, Using the initial RF equation for determining the nonlinear coefficient as the RF time-domain coupled mode equation; Before obtaining the first lightwave amplitude and the second lightwave amplitude in the modulator according to the RF voltage distribution and the lightwave time-domain coupling mode equation, the method further includes: Get the initial light wave equation, determining nonlinear coefficients and boundary conditions in the initial light wave equation, The initial light wave equation for determining the nonlinear coefficient is used as the light wave time domain coupled mode equation.

3. The method according to claim 2, characterized in that The initial RF equation is: V rf,F (z,t)=A rf,F (z,t)exp(iβ rf z-iω rf t)+c.c. V rf,R (z,t)=A rf,R (z,t)exp(-iβ rf z-iω rf t)+c.c. Among them, α rf is the radio frequency wave attenuation constant, A rf,F (z, t) is the amplitude of the forward propagating RF wave, v g,rf is the RF wave group velocity, A rf,R (z, t) is the amplitude of the reverse propagating RF wave, Λ dc is the electrostatic field voltage, Λ rf,F (z, t) is the voltage of the forward propagating RF wave, Λ rf,R (z, t) is the reverse propagation RF wave voltage, β rf is the phase constant of the radio frequency wave, Z is the spatial coordinate, t is the time, ω rf is the angular frequency of the radio frequency wave, cc represents the complex conjugate of the previous expression, is the second-order nonlinearity of the forward RF wave, is the third-order nonlinearity of the forward RF wave, is the fourth-order nonlinearity of the forward RF wave, is the second-order nonlinearity of the reflected RF wave, is the third-order nonlinearity of the reflected RF wave, is the fourth-order nonlinearity of the reflected RF wave, Λ dr,12 ,Λ r,4 ,Λ dr,22 ,Λ dr,22 ,Λ dr,14 and Λ dr,32 is the nonlinear coefficient.

4. The method according to claim 3, characterized in that The initial light wave equation is: Among them, α opt is the light wave attenuation constant, v g,opt is the light wave group velocity, A opt (z,t) is the light wave amplitude, is the second-order nonlinearity of the light wave, is the third-order nonlinearity of the light wave, V dc is the electrostatic field voltage, V rf,F is the forward propagating RF wave voltage, V rf,R is the reverse propagation RF wave voltage, ε0 is the vacuum dielectric constant, and h is The power index, j is The power index, k is The power index of is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, Λ dp,12 ,Λ rp,12 ,Λ dp,22 ,Λ rdp,112 and Λ rp,22 is the nonlinear coefficient.

5. The method according to claim 4, characterized in that The determining of the nonlinear coefficients and boundary conditions of the initial radio frequency equation and the determining of the nonlinear coefficients and boundary conditions in the initial light wave equation are carried out in the following manner: L dp,12 =L rp,12 L dp,22 =L rp,22 =L drp,112 The boundary conditions are as follows: From rf,R (L,t)=Γ L From rf,F (L,t)exp(2iβ rf L) Λ drp,hjk is the nonlinear coefficient, is the electrostatic field distribution, is the RF wave mode field distribution, is the light wave mode field distribution, h is The power index, j is The power index, k is The power index of h, j and k are positive integers and h+j+k-1≥2, and d is used to refer to r is used to refer to p is used to refer to χ is the nonlinear polarizability of the material, Δn eff,rf is the effective refractive index of radio frequency waves, k 0,rf is the vacuum wave vector of the radio frequency wave, A rf,R (L, t) is A at the coordinate (L, t) rf,R , A rf,F (L, t) is A at the coordinate (L, t) rf,F , β rf is the RF wave phase constant, L is the length of the modulator traveling wave region, ω rf is the angular frequency of the RF wave, V in (t) is the RF voltage signal input to the modulator, Z c is the impedance of the equivalent RF waveguide of the electrode, R g is the internal resistance of the equivalent RF waveguide of the electrode, R L is the load resistance of the equivalent RF waveguide of the electrode, Λ dp,12 ,Λ rp,12 ,Λ dp,22 ,Λ rdp,112 ,Λ rp,22 ,Λ drp,hjk ,Λ dr,12 ,Λ dr,22 and Λ dr,32 is the nonlinear coefficient.

6. The method according to claim 1, characterized in that The RF signal is input into the RF time domain coupled mode equation to obtain the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude in the modulator in the following manner: Δz rf =Δt / β 1,rf Δz opt =Δt / β 1,opt Δt is the time interval after the time t is discretized, α rf is the RF wave attenuation constant, V dc is the electrostatic field voltage, β 1,rf and β 1,opt is the first derivative of the propagation constant, v g,rf is the group velocity of RF, v g,opt is the group velocity of the light wave, Λ dr,12 ,Λ r,4 ,Λ dr,22 ,Λ dr,32 and Λ dr,14 is the nonlinear coefficient.

7. The method according to claim 1, characterized in that The first light wave amplitude and the second light wave amplitude in the modulator are obtained according to the RF voltage distribution and the light wave time domain coupling mode equation in the following manner: A opt,1 (z+Δz opt ,t+Δt)=exp(Δφ1(z,t))A opt,1 (z,t) A opt,2 (z+Δz opt ,t+Δt)=exp(Δφ2(z,t))A opt,2 (z,t) Δφ1(z,t)=-α opt Δz opt +2iΛ dp,12 (V rf,F (z,t)+V rf,R (z,t)+V dc )Δz opt +3iΛ dp,22 (V rf,F (z,t)+V rf,R (z,t)+V dc ) 2 Δz opt Δφ2(z,t)=-α opt Δz opt +2iΛ dp,12 (-V rf,F (z,t)-V rf,R (z,t)+V dc )Δz opt Among them, A opt,1 (z, t) is the amplitude of the first light wave, A opt,2 (z, t) is the amplitude of the second light wave, Δz opt represents the discretization interval of the spatial coordinate z, α opt is the light wave attenuation constant, V rf,F (z, t) is the forward propagation RF wave voltage, V rf,R (z, t) is the reverse propagation RF wave voltage, V dc is the electrostatic field voltage, Λ dp,12 ,Λ dp,22 and Λ dp,12 is the nonlinear coefficient.

8. A modulator optoelectronic joint simulation device, characterized in that: The modulator receives a radio frequency signal and outputs an optical signal. The device comprises: A signal acquisition unit, used to acquire the input radio frequency signal of the modulator; A radio frequency simulation unit, used for inputting the radio frequency signal into a radio frequency time domain coupled mode equation to obtain the amplitude of the forward propagating radio frequency wave and the amplitude of the reverse propagating radio frequency wave in the modulator; The RF simulation unit is further used to obtain the RF voltage distribution in the modulator according to the forward propagation RF wave amplitude and the reverse propagation RF wave amplitude; A light wave simulation unit, used for obtaining the first light wave amplitude and the second light wave amplitude in the modulator according to the RF voltage distribution and the light wave time domain coupling mode equation; The eye diagram generating unit is used to obtain the optical signal output by the modulator according to the first optical wave amplitude and the second optical wave amplitude, and obtain a simulated eye diagram according to the optical signal.

9. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images