Simulation method for propagation of inhomogeneous partially coherent light in turbulent air-sea links

Abstract

This invention discloses a simulation method for the transmission of non-uniform partially coherent light in turbulent air-sea links, comprising the following steps: S1: Modeling of non-uniform partially coherent light; S2: Modeling of air-sea turbulent links; S3: Modeling of the transmission of non-uniform partially coherent light in air-sea turbulent links; S4: Obtaining the light field reaching the receiving plane; S5: Calculating the light intensity on the receiving plane; S6: Calculating the aperture-average scintillation index; S7: Calculating the signal-to-noise ratio of the optical link. This invention primarily addresses the difficulty of constructing practical air-sea link experimental platforms in existing technologies. The simulation method of this invention can simulate the transmission of light beams in turbulent environments within air-sea links, and simultaneously enables autonomous control of the turbulent environment. This facilitates the use of controlled variables to study the influence of various parameters of the light beam and turbulence on the beam transmission characteristics, thereby exploring optical carriers with strong anti-turbulence capabilities and improving the performance of space optical communication systems.

Description

Technical Field

[0001] The invention belongs to the technical field of spatial optical transmission, and in particular relates to a simulation method for transmission of non-uniform partially coherent light in air-sea link turbulence. Background Art

[0002] At present, human exploration activities in space are increasing. The high-resolution space exploration platform and multi-sensor composite detection technology used have led to an explosive growth in the amount of communication data between the satellite and the ground. Satellite communication using radio frequency cannot meet the demand. This situation restricts the overall improvement of the level of informatization in modern society. The development of society needs to be supported by a communication system with fast data transmission rate, large channel capacity and wide coverage. Space optical communication is a communication technology that uses light waves as a carrier to transmit information in a vacuum or atmosphere. It has the advantages of large communication capacity and high data transmission rate. With the breakthrough of key technologies in satellite laser communication and the gradual manifestation of its advantages, it has become a general consensus in the research and industrial fields to use optical communication for satellite communication to meet the society's growing demand for high data rates and large communication capacity. The communication structure for the future world will be a communication system of integrated sea, land, air and space optical networks where satellite optical networks and ground optical fiber networks are intertwined.

[0003] When optical carriers are transmitted in the atmosphere and ocean channels, the turbulence present will cause optical carrier wavefront distortion, resulting in light spot dispersion, beam drift, and light intensity flicker on the receiving plane. These factors will lead to a decrease in the quality of the received optical signal and reduce the signal-to-noise ratio of the receiving end signal, thereby affecting the system's bit error rate, communication distance, and communication bandwidth, restricting the development of space optical communication. Therefore, space optical communication has its unique advantages, but it also faces some scientific problems that need to be solved in the actual application process. Among them, improving the anti-turbulence ability of optical carriers and thus improving the performance of space optical communication systems is an urgent problem to be solved.

[0004] Improving the anti-turbulence capability of optical carriers requires studying the transmission characteristics of light beams in turbulence and clarifying the interaction mechanism between light and turbulence. Among them, the platform for studying the transmission of light beams in turbulence is very important. Given that turbulence is an irregular and complex flow, the distance of air-sea links is thousands of meters, and the air experimental platform is not easy to build, it is of great significance to build a simulation platform that simulates the transmission of light beams in atmospheric turbulence in inclined links. Summary of the invention

[0005] In order to solve the difficulties in building actual air-sea links in the prior art, the present invention provides a simulation method for the transmission of non-uniform partially coherent light in turbulence of an air-sea link with scientific principles and simple operation.

[0006] In order to solve the above technical problems, the present invention adopts the following technical solution: a simulation method for the transmission of non-uniform partially coherent light in air-sea link turbulence comprises the following steps:

[0007] S1: Modeling of inhomogeneous partially coherent light;

[0008] S2: air-sea turbulence link modeling;

[0009] S3: Modeling of inhomogeneous partially coherent light transmission in air-sea turbulent links;

[0010] S4: Obtain the light field reaching the receiving plane;

[0011] S5: Calculate the light intensity on the receiving plane;

[0012] S6: Calculate the aperture average flicker index;

[0013] S7: Calculate the signal-to-noise ratio of the optical link.

[0014] Further, step S1 is specifically as follows: the light source can use any type of light beam, and the size of each parameter of the target light beam is determined according to the application environment; the function of modulating the coherence distribution can be independently selected according to the needs to achieve non-uniform distribution of spatial coherence; the spatial random transmission phase screen theory is used to complete the construction of the non-uniform partially coherent light beam, and it is emitted by the transmitting end;

[0015] (1)

[0016] In formula (1), represents the constructed partially coherent light; represents the selected fully coherent light; It represents the random transmission phase screen function for constructing partially coherent light, and its expression is:

[0017] (2)

[0018] (3)

[0019] In formula (3), The variance is The subscript “nn” represents a sampling grid point of size n×n; is the Gauss-Scherr mode, and its expression is:

[0020] (4)

[0021] in, is the width parameter; the coherence length The relationship is:

[0022] (5)

[0023] The mutual coherence function of partially coherent light is:

[0024] (6)

[0025] The coherence length is modulated to obtain non-uniform partially coherent light, and its mutual coherence function expression is:

[0026] (7)

[0027] Further, step S2 is specifically as follows: based on the power spectrum model of atmospheric and oceanic turbulence and the Hufnagel-Valley model describing the variation of the refractive index structure constant with height, a phase screen simulating atmospheric and oceanic turbulence is constructed using Monte-Carlo technology and subharmonic compensation theory;

[0028] (8)

[0029] In formula (8), Re represents the real part; F -1 represents inverse Fourier transform; C nn represents a complex random normal number with zero mean and unit variance, and the subscript "nn" represents the sampling points on a grid of size n×n; Represents a power spectrum model of the atmosphere or ocean;

[0030] The refractive index power spectrum of atmospheric turbulence adopts the Von Karman model:

[0031] (9)

[0032] In the above formula, , f is the spatial frequency on the phase screen plane; , ; The atmospheric turbulence intensity in formula (9) is is a function of height, and its expression is:

[0033] (10)

[0034] in, h Indicates the height from the ground in meters; wind speed The value is 21m / s; Represents the refractive index structure constant within 1 km from the ground;

[0035] The refractive index power spectrum model of ocean turbulence is:

[0036] (11)

[0037] (12)

[0038] (13)

[0039] (14)

[0040] (15)

[0041] (16)

[0042] In formulas (11)-(16), C 0 The value is 0.72; the turbulent kinetic energy dissipation rate per unit fluid mass The value range is 10 -1 m 2 / s 3 to 10 -10 m 2 / s 3 ; C 1 The value is 2.35; η is the Kolmogorov microscale; , ; and denote the dissipation rates of mean square temperature and salinity, respectively, express and relevance; ; Under constant mean temperature and salinity gradients, ω is the ratio of the temperature to salinity contributions to the spectrum, ranging from −5 (corresponding to temperature-dominated turbulence) to 0 (corresponding to salinity-dominated turbulence); and are the eddy thermal diffusivity and the diffusion coefficient of the salt, respectively;

[0043] A single turbulence phase screen is used to simulate the atmospheric turbulence effect at a certain moment in the communication channel, and two different turbulence phase screens are used to represent the atmospheric turbulence of the same intensity at different moments. C nn , so that different turbulence screens with the same turbulence intensity can simulate the changes in atmospheric turbulence; at the same time, by controlling the parameters of turbulence, a constant turbulence intensity in the transmission path can be achieved, which makes it easier to use the controlled variable method to study the influence of different parameters of the light beam on its transmission characteristics.

[0044] Furthermore, step S3 is specifically as follows: dividing the transmission distance into several segments, and completing the segmented transmission algorithm of the light beam in turbulence based on the extended Huygens-Fresnel theory; at the same time, using the turbulence phase screen generation method described in step S2, obtaining a phase screen that simulates the turbulence effect of each segment, and inserting the obtained phase screen of the segment into the corresponding segment; completing the modeling of the transmission of the non-uniform partially coherent light beam in the turbulent environment of the air-sea link;

[0045] (17)

[0046] The segmented transmission algorithm for light beam propagation in turbulence is shown in formula (17), where L is the transmission distance; the transmission distance is divided into n part, represents the position vector on the i-th receiving plane; F and F -1 represents the forward and inverse Fourier transform; and is the adjacent i and i+1 The light field on the receiving plane; Segmented phase screen for simulating turbulence; is the transmission function of the light beam in free space; by inputting the non-uniform partially coherent light on the source plane and using the segmented transmission algorithm, the light intensity on any segmented receiving plane can be obtained, thereby completing the modeling of the simulation platform for the transmission of the light beam in the turbulent environment of the air-sea link.

[0047] Furthermore, step S4 is specifically as follows: using a simulation platform for light beam transmission in a turbulent environment of an air-sea link, calculate the light field of the light beam emitted from the source plane and arriving at the receiving plane after passing through the air-sea turbulent link; input the non-uniform partially coherent light beam obtained in step S1 into the simulation platform for light beam transmission in a turbulent environment of an air-sea link to obtain the light field distribution on the receiving plane.

[0048] Further, step S5 is specifically as follows: based on the calculated light field, the light intensity on the receiving plane is obtained by using the light intensity equal to the product of the light field and its conjugate

[0049] (18) .

[0050] Further, step S6 is specifically as follows: based on the random number generated by the system each time in formula (8) of step 2 C nn Different, multiple realizations of random phase screens simulating turbulence of the same intensity are realized; on this basis, the traversal of the integral value of the light intensity and its square in the receiving aperture during multiple turbulence realizations is completed; the traversal average value of the integral value of the light intensity and its square in the receiving aperture corresponding to the multiple turbulence realizations is calculated; according to the calculation formula of the scintillation index, the corresponding values ​​obtained are substituted into the aperture average scintillation index:

[0051] (19)

[0052] In the above formula, R represents the radius of the receiving aperture; < > represents ergodic averaging.

[0053] Further, step S7 is specifically: based on the aperture average scintillation index obtained in step S6 and the signal-to-noise ratio model of the optical communication system in turbulence, the corresponding values ​​are substituted to obtain the signal-to-noise ratio of the optical link used in the model:

[0054] (20)

[0055] SNR 0 and P S0 are the output signal-to-noise ratio and signal power in free space respectively; P S > is the average signal power in atmospheric turbulence.

[0056] The above technical solution is very convenient for studying the transmission characteristics of light beams in air-sea link turbulence, which mainly solves the problem that it is difficult to build an actual air-sea link experimental platform in the prior art. The simulation method of the present invention can simulate the transmission of light beams in the turbulent environment of air-sea links, and can realize autonomous control of the turbulent environment. It is very convenient to use the method of controlling variables to study the influence of various parameters of light beams and turbulence on the transmission characteristics of light beams, and explore optical carriers with strong anti-turbulence capabilities, thereby improving the performance of space optical communication systems. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] The present invention will be further described below with reference to the accompanying drawings and embodiments, in which:

[0058] Figure 1 is a flow chart of the present invention;

[0059] Figure 2 It is a schematic diagram of the uplink and downlink transmission of light beams in the air-sea turbulent environment;

[0060] Figure 3 is the turbulent phase screen diagram of the first segment;

[0061] Figure 4 This is the turbulent phase screen diagram of the second segment;

[0062] Figure 5 This is the turbulent phase screen diagram of the third segment;

[0063] Figure 6 This is the turbulent phase screen diagram of the fourth segment;

[0064] Figure 7 This is the turbulent phase screen diagram of the fifth segment;

[0065] Figure 8 This is the turbulent phase screen diagram of the sixth segment;

[0066] Fig. 9 This is the turbulent phase screen image of the seventh segment;

[0067] Fig.10 This is the turbulent phase screen image of the eighth segment;

[0068] Fig.11 This is the turbulent phase screen image of the ninth segment;

[0069] Fig.12 This is the turbulent phase screen image of the tenth segment;

[0070] Fig.13 This is the turbulence phase screen of the eleventh segment;

[0071] Fig.14 This is the turbulent phase screen diagram of the twelfth segment;

[0072] Fig.15 This is the turbulent phase screen image of the thirteenth segment;

[0073] Fig.16 This is the turbulent phase screen diagram of the fourteenth segment;

[0074] Fig.17 This is the turbulent phase screen diagram of the fifteenth segment;

[0075] Fig.18 is the instantaneous light intensity diagram of the non-uniform partially coherent beam Bessel-sg1 uplink on the receiving plane;

[0076] Fig.19 It is the instantaneous intensity diagram of the non-uniform partially coherent beam Bessel-sg1 downlink on the receiving plane.

[0077] Table 1 shows the average aperture scintillation index and signal-to-noise ratio of the optical link of the non-uniform partially coherent beam Bessel-sg1 at different transmission distances in the uplink and downlink. DETAILED DESCRIPTION

[0078] The technical solution of the present invention will be further specifically described below through embodiments and in conjunction with the accompanying drawings. Figure 1 and Figure 2 As shown, the simulation method of the transmission of non-uniform partially coherent light in air-sea link turbulence of the present invention comprises the following steps:

[0079] S1: Modeling of inhomogeneous partially coherent light;

[0080] S2: air-sea turbulence link modeling;

[0081] S3: Modeling of inhomogeneous partially coherent light transmission in air-sea turbulent links;

[0082] S4: Obtain the light field reaching the receiving plane;

[0083] S5: Calculate the light intensity on the receiving plane;

[0084] S6: Calculate the aperture average flicker index;

[0085] S7: Calculate the signal-to-noise ratio of the optical link.

[0086] Step S1 specifically includes: using Bessel beam as the light source, using Gaussian Schell mode as the mutual coherence function, and using super Gaussian function as the function of modulating coherence distribution. Using a spatial random transmission phase screen to complete the construction of a non-uniform partially coherent Bessel beam, and emitting it from the transmitting end.

[0087] The mathematical expression of Bessel beam is:

[0088] (a)

[0089] in, represents the position vector on the source plane, , J n Indicates the first category n Bessel function of order, a is the width parameter;

[0090] The mutual coherence function of partially coherent light is:

[0091] (b)

[0092] l c is the coherence length, is the super Gaussian function of the position-based modulation coherence distribution, and its expression is as follows:

[0093] (c)

[0094] In formula (3), N is the spatial scale factor, is the width of the Gaussian function, d Indicates the rate of change of coherence for a specific modulation index.

[0095] Step S2 is specifically as follows: the inner and outer dimensions of the atmospheric turbulence are set to 5 mm and 5 m respectively, and based on the power spectrum model of atmospheric and oceanic turbulence and the Hufnagel-Valley refractive index structure constant model, a phase screen simulating atmospheric and oceanic turbulence is constructed using Monte-Carlo technology, turbulence spectrum inversion theory, and subharmonic compensation theory. The segmented turbulence phase screens in the uplink generated in this embodiment are shown in the attached figure. Figure 3-17 shown.

[0096] Step S3 is specifically as follows: dividing the transmission distance into several segments, inserting a phase screen simulating the turbulence effect in each segment, and the total phase screen effect is the same as the turbulence effect of the entire air-sea turbulent link; based on the extended Huygens-Fresnel theory, using a multi-step algorithm, the modeling of the transmission of non-uniform partially coherent Bessel light beams in the turbulent environment of the air-sea link is completed.

[0097] (d)

[0098] The segmented transmission algorithm for light beams in turbulence is shown in the above formula, where L is the transmission distance; the transmission distance is divided into n part, represents the position vector on the i-th receiving plane; F and F -1 represents the forward and inverse Fourier transform; and is the adjacent i and i+1 The light field on the receiving plane; Segmented phase screen for simulating turbulence; is the transmission function of the light beam in free space. By inputting the non-uniform partially coherent light on the source plane and using the segmented transmission algorithm, the light intensity on any segmented receiving plane can be obtained, thus completing the modeling of the simulation platform for the transmission of light beams in the turbulent environment of the air-sea link.

[0099] Step S4 is specifically: using a simulation platform for light beam transmission in a turbulent environment of an air-sea link, calculating the light field of a non-uniform partially coherent Bessel light beam emitted from a source plane and reaching a receiving plane after passing through the turbulence of the air-sea turbulent link.

[0100] Step S5 specifically includes: based on the calculated light field, the light intensity on the receiving plane is obtained by using the light intensity being equal to the product of the light field and its conjugate. In this embodiment, the instantaneous light intensity of the non-uniform partially coherent light beam Bessel-sg1 on the uplink and downlink receiving planes is shown in the attached figure. Fig.18 and 19 shown.

[0101] Step S6 is specifically: using the random number generated by the system each time C nnDifferent, multiple realizations of random phase screens simulating turbulence of the same intensity are realized. On this basis, the traversal of the integral value of the light intensity and its square in the receiving aperture during multiple turbulence realizations is completed. The integral values ​​of the light intensity and its square in the receiving aperture corresponding to the multiple turbulence realizations are processed to obtain their traversal average values. According to the calculation formula of the scintillation index, the corresponding values ​​are substituted to obtain the aperture average scintillation index. The aperture average scintillation index of the non-uniform partially coherent light beam Bessel-sg1 at different transmission distances in the uplink and downlink in this embodiment is shown in Table 1.

[0102] Step S7 is specifically as follows: based on the obtained aperture average scintillation index and the signal-to-noise ratio model of the optical communication system in turbulence, the signal-to-noise ratio of the optical link used in the model is obtained. The average signal-to-noise ratio of the non-uniform partially coherent light beam Bessel-sg1 at different transmission distances in the uplink and downlink in this embodiment is shown in Table 2.

[0103] In this embodiment, the aperture radius of the transmitting end is 25 mm, and the aperture radius of the receiving end is 50 mm. The wavelength used is 532 nm, and the inner and outer dimensions of the turbulence are 5 mm and 5 m respectively. The number of realizations of the transmission phase screen for partially coherent light is 100, and the number of realizations of the turbulent phase screen is 500. The coherence modulation index is set to 0.5, which means the coherence change rate under a specific coherence modulation index. d Take 1. The transmission distance in the ocean is 50 m. The transmission distance in the atmosphere is divided into 10 segments, and the transmission distance in the ocean is divided into 5 segments. After the above steps 1 to 7, the light intensity distribution, aperture average scintillation index and signal-to-noise ratio of the Bessel-sg1 beam at different turbulence intensities and transmission distances in the uplink and downlink can be obtained. Specific data Fig.18 , 19 , as shown in Table 1 and Table 2.

[0104] Table 1 shows the average aperture scintillation index and signal-to-noise ratio of the optical link of the non-uniform partially coherent beam Bessel-sg1 at different transmission distances in the uplink and downlink:

[0105] Table 1 Aperture average flicker index

[0106]

[0107] Table 2 Signal-to-noise ratio

[0108]

[0109] The above description is only a preferred embodiment of the present invention. Although the present invention has been described in detail with reference to the aforementioned embodiments, it is still possible for those skilled in the art to modify the technical solutions described in the aforementioned embodiments, or to make equivalent substitutions for some of the technical features therein. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A simulation method for the transmission of non-uniform partially coherent light in turbulence of an air-sea link, characterized in that: The following steps are involved: S1: Modeling of inhomogeneous partially coherent light; S2: air-sea turbulence link modeling; S3: Modeling of inhomogeneous partially coherent light transmission in air-sea turbulent links; S4: Obtain the light field reaching the receiving plane; S5: Calculate the light intensity on the receiving plane; S6: Calculate the aperture average flicker index; S7: Calculate the signal-to-noise ratio of the optical link; Step S1 is specifically as follows: the light source can use any type of light beam, and the size of each parameter of the target light beam is determined according to the application environment; the function of modulating the coherence distribution can be independently selected according to the needs to achieve non-uniform distribution of spatial coherence; the spatial random transmission phase screen theory is used to complete the construction of the non-uniform partially coherent light beam, and it is emitted by the transmitting end; U(x,y)=U in (x, y)P t (x, y) (1) In formula (1), U(x, y) represents the constructed partially coherent light; U in (x, y) represents the selected fully coherent light; P t (x, y) represents the random transmission phase screen function for constructing partially coherent light, and its expression is: In formula (3), R nn (x, y) is the variance The subscript "nn" represents the sampling grid points of size n×n; f(x, y) is a Gauss-Scherr pattern, and its expression is: in, is the width parameter; the coherence length l c The relationship is: The mutual coherence function of partially coherent light is: The coherence length is modulated to obtain non-uniform partially coherent light, and its mutual coherence function expression is: in, Step S2 specifically includes: constructing a phase screen simulating atmospheric and oceanic turbulence based on a power spectrum model of atmospheric and oceanic turbulence and a Hufnagel-Valley model describing the variation of the refractive index structure constant with altitude by using Monte-Carlo technology and subharmonic compensation theory; In formula (8), Re represents the real part; F -1 represents inverse Fourier transform; C nn represents a complex random normal number with zero mean and unit variance, and the subscript "nn" represents the sampling points on an n×n grid; Φ φ (κ) represents the power spectrum model of the atmosphere or ocean; The refractive index power spectrum of atmospheric turbulence adopts the Von Karman model: In the above formula, κ = 2πf, f is the spatial frequency on the phase screen plane; κ m =5.92 / l0, κ0=2π / L0; the formula (9) represents the atmospheric turbulence intensity is a function of height, and its expression is: Wherein, h represents the height from the ground, in meters; the wind speed v is 21m / s; Represents the refractive index structure constant within 1 km from the ground; The refractive index power spectrum model of ocean turbulence is: x n =α2χ T +β2x S -2abx TS (12) δ=8.824(κη) 4 / 3 +12.978(ki) 2 (16) In formulas (11)-(16), C0 is taken as 0.72; the turbulent kinetic energy dissipation rate ε per unit fluid mass ranges from 10 - 1 m 2 / s 3 to 10 -10 m 2 / s 3 ; C1 is 2.35; η is the Kolmogorov microscale; α = 2.6 × 10 -4 L / deg,β=1.75×10 -4 L / g; T and χ S denote the dissipation rate of mean square temperature and salinity, respectively, TS Represents χ T and χ S The correlation of T =1.683×10 -2 , A S =1.9×10 -4 , A TS =9.41×10 -3 ; When the mean temperature and salinity gradients are constant, ω is the ratio of the temperature to salinity contribution to the spectrum, ranging from -5 to 0; K T and K S are the eddy thermal diffusivity and the diffusion coefficient of the salt, respectively; A single turbulence phase screen is used to simulate the atmospheric turbulence effect at a certain moment in the communication channel, and two different turbulence phase screens are used to represent the atmospheric turbulence of the same intensity at different moments. nn , so that different turbulence phase screens with the same turbulence intensity can simulate the changes in atmospheric turbulence; at the same time, by controlling the turbulence parameters, a constant turbulence intensity in the transmission path can be achieved, making it easier to use the control variable method to study the influence of different light beam parameters on its transmission characteristics; Step S3 is specifically as follows: the transmission distance is divided into several segments, and based on the extended Huygens-Fresnel theory, a segmented transmission algorithm for the transmission of a light beam in turbulence is completed; at the same time, the turbulence phase screen generation method described in step S2 is used to obtain a phase screen that simulates the turbulence effect of each segment, and the obtained phase screen of the segment is inserted into the corresponding segment; the modeling of the transmission of a non-uniform partially coherent light beam in a turbulent environment of an air-sea link is completed; The segmented transmission algorithm for light beam transmission in turbulence is shown in formula (17), where L is the transmission distance; the transmission distance is divided into n segments, represents the position vector on the i-th receiving plane; F and F -1 represents the forward and inverse Fourier transform; U i and U i+1 is the light field on the adjacent i-th and i+1-th receiving planes in the propagation path; Segmented phase screen for simulating turbulence; is the transmission function of the light beam in free space; by inputting the non-uniform partially coherent light on the source plane and using the segmented transmission algorithm, the light intensity on any segmented receiving plane can be obtained, thereby completing the modeling of the simulation platform for the transmission of the light beam in the turbulent environment of the air-sea link; Step S4 is specifically as follows: using a simulation platform for light beam transmission in a turbulent air-sea link environment, calculate the light field of the light beam emitted from the source plane and arriving at the receiving plane after passing through the air-sea turbulent link; input the non-uniform partially coherent light beam obtained in step S1 into the simulation platform for light beam transmission in a turbulent air-sea link environment to obtain the light field distribution on the receiving plane.

2. The simulation method for transmission of non-uniform partially coherent light in turbulence of air-sea links according to claim 1, characterized in that: Step S5 is specifically as follows: based on the calculated light field, the light intensity on the receiving plane is obtained by using the light intensity being equal to the product of the light field and its conjugate 3. The simulation method for transmission of non-uniform partially coherent light in turbulence of air-sea links according to claim 2, characterized in that: Step S6 is specifically as follows: based on the random number C generated by the system each time in formula (8) in step 2 nn Different, multiple realizations of random phase screens simulating turbulence of the same intensity are realized; on this basis, the traversal of the integral value of the light intensity and its square in the receiving aperture during multiple turbulence realizations is completed; the traversal average value of the integral value of the light intensity and its square in the receiving aperture corresponding to the multiple turbulence realizations is calculated; according to the calculation formula of the scintillation index, the corresponding values ​​obtained are substituted into the aperture average scintillation index: In the above formula, R represents the radius of the receiving aperture; <> represents the ergodic average; represents the position vector on the receiving plane.

4. The simulation method for transmission of non-uniform partially coherent light in air-sea link turbulence according to claim 2, characterized in that: Step S7 specifically includes: based on the aperture average scintillation index obtained in step S6 and the signal-to-noise ratio model of the optical communication system in turbulence, substituting corresponding values ​​into the signal-to-noise ratio of the optical link used in the model: SNR0 and P S0 are the output signal-to-noise ratio and signal power in free space, respectively; <P S > is the average signal power in atmospheric turbulence.

Images