Atmospheric turbulence light intensity flicker model parameter estimation method, system, equipment and medium

Abstract

The invention discloses an atmospheric turbulence light intensity flicker model parameter estimation method, system and device and a medium, and relates to the technical field of free space optical communication and free space quantum communication, and the method comprises the steps: obtaining a plurality of channel sampling values of atmospheric turbulence, and obtaining initial channel parameters of the channel sampling values; generating a plurality of groups of light intensity flicker samples according to the estimated channel parameters, and obtaining a light intensity flicker estimation model of each group of light intensity flicker samples by using a # imgabs0 # nearest neighbor density estimation method; data which are not in the light intensity flicker sample range in the channel sampling values are eliminated, and the rest data are input into the light intensity flicker estimation model to obtain an average log-likelihood function value; and iterating the step of obtaining the average log-likelihood function value, and iterating the # imgabs1 # value and the channel parameter by using a gradient descent method to obtain the maximum average log-likelihood function value and the channel parameter value. According to the method, the balance between the estimation precision and the calculation complexity is realized, and the parameter estimation efficiency of different channel models is improved.

Description

Technical Field

[0001] The present invention relates to the fields of free-space optical communication technology and free-space quantum communication technology, and in particular to a method, system, device and medium for estimating parameters of an atmospheric turbulence light intensity scintillation model. Background Art

[0002] With the continuous development of society and the rapid advancement of science and technology, traditional ground communication technology has shortcomings in wide-area coverage, low latency and high-bandwidth transmission, which has led to the development of integrated space-ground communication. In recent years, communication systems based on optical carriers have developed rapidly on free-space channels. In satellite-to-ground laser or optical quantum communications, there is an atmosphere more than 20 kilometers thick. The random changes in temperature and pressure in space will produce atmospheric turbulence, which in turn causes random changes in the refractive index, and ultimately causes the power of light to propagate to the receiving end when it passes through the atmospheric channel. This fading effect is called scintillation. Due to the wave-particle duality of light, the interaction between quantum signals and atmospheric turbulence is consistent with classical light, so the model of the quantum light efficiency at the receiving end changing with time is consistent with classical light.

[0003] In order to effectively evaluate the performance of actual free-space channels and optical quantum communication systems, researchers have proposed many light intensity flicker models and compared and discussed the adaptability of light intensity flicker models in detail. According to the number of parameters used to determine the model, the existing light intensity flicker models can be divided into single-parameter models, dual-parameter models, and multi-parameter models. Among the single-parameter turbulence models, the lognormal model and the negative exponential model are commonly used, and their scope of application is limited to the weak turbulence and strong turbulence limit areas. Among the dual-parameter turbulence models, the light intensity flicker Gamma-Gamma model and the wireless channel Lognormal-Rician model are commonly used. They have the ability to adapt to the range of changes from weak turbulence to strong turbulence. In some channel scenarios, such as large aperture reception or when the laser transmission mode is spherical wave propagation, the Lognormal-Rician model is more accurate than the Gamma-Gamma model. The multi-parameter turbulence model includes the Malaga model (five-parameter model). Generally speaking, the more complex the light intensity flicker model is, the wider its adaptability and channel conditions are. Since the light intensity flicker model parameters need to be estimated when evaluating the actual system performance, the two-parameter turbulence model is the most studied in terms of trade-offs between complexity and model accuracy.

[0004] In the dual-parameter turbulence model, the parameter estimation of the Gamma-Gamma channel includes moment estimation, expected maximum algorithm, fractional order moment estimation, etc., and the estimation methods of the Lognormal-Rician include generalized order moment estimation method, expected maximum algorithm and saddle point approximation algorithm. However, the parameter estimation method used in the existing dual-parameter turbulence model involves complex functions or integral operations when obtaining a high-precision model, which will cause the estimation process to become complicated, which is not conducive to parameter estimation of different channel models through the model, and even cannot accurately and effectively evaluate the performance of free space channels and optical quantum communication systems. Summary of the invention

[0005] The purpose of the present invention is to address the deficiencies of the above-mentioned prior art and provide a method, system, device and medium for estimating the parameters of an atmospheric turbulence light intensity scintillation model, so as to solve the problem that the existing estimation methods in the prior art involve complex functions or integral operations, which will cause the estimation process to become complicated, making it difficult to estimate the parameters of different channel models through the model, and even making it impossible to accurately and effectively evaluate the performance of free space channels and optical quantum communication systems.

[0006] The present invention specifically provides the following technical solution: a method for estimating parameters of an atmospheric turbulence light intensity scintillation model, comprising the following steps: Acquire multiple channel sampling values ​​of atmospheric turbulence, and acquire initial channel parameters of the multiple channel sampling values; According to the initial channel parameters and the number of samples k Values ​​generate multiple sets of light intensity flicker samples, and use The nearest neighbor density estimation method is used to obtain the light intensity flicker estimation model for each group of light intensity flicker samples; Eliminate the data that are not in the range of the minimum value to the maximum value of each group of light intensity flicker samples from the channel sampling values, input the remaining channel sampling values ​​into the light intensity flicker estimation model, obtain the log-likelihood function value of each group of light intensity flicker samples, and average multiple log-likelihood function values; The step of obtaining the average log-likelihood function value is iterated to obtain the maximum average log-likelihood function value, and the channel parameter value of the light intensity flicker estimation model is represented by the maximum average log-likelihood function value.

[0007] Preferably, the obtaining of multiple channel sampling values ​​of atmospheric turbulence and obtaining initial channel parameters of the multiple channel sampling values ​​is specifically: Get channel sampling value parameters The initial estimate expression of , , specifically: ; in, The photodetector obtains atmospheric turbulence channel sampling values, satisfying , Indicates the total sample size, 1 and 2 indicate the first and second samples.

[0008] Preferably, the method according to the initial channel parameters and the number of samples k Generate multiple sets of light intensity flicker samples, including: The specific expressions of multiple groups of light intensity flicker samples are: ; in Subscript Indicates Group light intensity flicker samples, , Indicates Group light intensity flicker samples, Indicates samples, , each light intensity flicker sample obeys the parameter The wireless channel Lognormal-Rician model distribution, Indicates the number of each group of light intensity flicker samples.

[0009] Preferably, the use The nearest neighbor density estimation method obtains the light intensity flicker estimation model for each group of light intensity flicker samples, including: For any set of light intensity flicker sample data, Group as an example ,use Nearest neighbor density estimation method and current Worth getting sample points Density Estimation for: ; The variables , is the gamma function, Represents the dimension of each data sample. For the channel in the Lognormal-Rician model of wireless channels, , Representation sequence Zhongyu The nearest distance values; For data points Perform linear interpolation to obtain the light intensity flicker estimation model and obtain the log-likelihood function value of the light intensity flicker estimation model .

[0010] Preferably, the data that is not in the range of the minimum value to the maximum value of each group of light intensity flicker samples in the channel sampling values ​​are eliminated, the remaining channel sampling values ​​are input into the light intensity flicker estimation model, the logarithmic likelihood function value of each group of light intensity flicker samples is obtained, and multiple logarithmic likelihood function values ​​are averaged, specifically: For Group as an example , channel sampling sequence After elimination, we get , satisfy: ; in, M is the total sample size; for any set of light intensity flicker samples ,satisfy , Indicates the maximum number of light intensity flicker samples after culling. n is the sampling number; Will Substitute into And find the average log-likelihood function value of each group , the specific expression is: ; in, is the normalization factor, and the specific expression is: .

[0011] If the light intensity flicker sample is group, then the average log-likelihood function value for: .

[0012] Preferably, the step of obtaining the average log-likelihood function value is iterated to obtain the maximum average log-likelihood function value, and the channel parameter value of the light intensity flicker estimation model is represented by the maximum average log-likelihood function value, including: Iterate through the gradient descent method Values ​​and channel parameters are used to obtain the maximum average log-likelihood function value. The specific expression is: .

[0013] By taking the difference of the maximum average log-likelihood function value, we can approximate each parameter The first-order derivative of is: ; In the formula, Indicates The estimated results.

[0014] The present invention also includes an atmospheric turbulence light intensity scintillation model parameter estimation system, comprising: An acquisition module, used for acquiring multiple channel sampling values ​​of atmospheric turbulence and acquiring initial channel parameters of the multiple channel sampling values; The model acquisition module is used to obtain the initial channel parameters and sample number. k Values ​​generate multiple sets of light intensity flicker samples, and use The nearest neighbor density estimation method is used to obtain the light intensity flicker estimation model for each group of light intensity flicker samples; An average processing module is used to remove data from the channel sampling values ​​that are not in the range from the minimum value to the maximum value of each group of light intensity flicker samples, input the remaining channel sampling values ​​into the light intensity flicker estimation model, obtain the logarithmic likelihood function value of each group of light intensity flicker samples, and average multiple logarithmic likelihood function values; The estimation module is used to iterate the step of averaging the log-likelihood function value to obtain the maximum average log-likelihood function value, and to represent the channel parameter value of the light intensity flicker estimation model through the maximum average log-likelihood function value.

[0015] The present invention provides a storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps of the above-mentioned atmospheric turbulence light intensity scintillation model parameter estimation method are implemented.

[0016] The present invention provides a computer device, comprising a memory and a processor, wherein a program is stored in the memory, and when the program is executed by the processor, the processor executes the steps of the above-mentioned atmospheric turbulence light intensity scintillation model parameter estimation method.

[0017] Compared with the prior art, the present invention has the following significant advantages: The present invention estimates the parameters of the light intensity flicker model based on the idea of ​​light intensity flicker samples. During operation, the parameters are not directly estimated according to the specific expression of probability density, but multiple groups of light intensity flicker samples are generated according to the initial channel parameters and the sample number k value, and the data that is not within the light intensity flicker sample range in the channel sampling value is eliminated, so as to obtain high-quality channel sampling values, improve the accuracy of the estimated data, and input the high-quality channel sampling values ​​into the respectively selected light intensity flicker estimation models, output the log-likelihood function value of each group of light intensity flicker samples, and obtain the maximum average log-likelihood function value and channel parameter value through iteration, so as to achieve the effect of directly estimating the log-likelihood function value of each point based on each group of light intensity flicker sample data, avoid the complex functions or integral operations involved in the parameter estimation, achieve the trade-off between the estimation accuracy and the computational complexity, improve the efficiency of the atmospheric turbulence light intensity flicker model parameter estimation, and further improve the accuracy and effect of the free space channel and the optical quantum communication system performance evaluation. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] Figure 1 The present invention provides a Lognormal-Rician turbulent channel based on The maximum likelihood and nearest neighbor-based Simulation diagram of the minimum KL divergence estimation result of the nearest neighbor, the simulation parameters are , the number of trials is 30 to calculate the MSE performance; Figure 2 The present invention provides a Lognormal-Rician turbulent channel based on The maximum likelihood and nearest neighbor-based Simulation diagram of the minimum KL divergence estimation result of the nearest neighbor, the simulation parameters are , the number of trials is 30 to calculate the MSE performance; Figure 3 The present invention provides a Lognormal-Rician turbulent channel based on The maximum likelihood and nearest neighbor-based Simulation diagram of the minimum KL divergence estimation result of the nearest neighbor, the simulation parameters are , the number of trials is 30 to calculate the MSE performance; Figure 4 The present invention provides a Lognormal-Rician turbulent channel based on The maximum likelihood and nearest neighbor-based Simulation diagram of the minimum KL divergence estimation result of the nearest neighbor, the simulation parameters are , the number of trials is 30 to calculate the MSE performance; Figure 5 A flow chart of a method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to the present invention; Figure 6 This is a second flow chart of a method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to the present invention. DETAILED DESCRIPTION

[0019] The following is a clear and complete description of the technical solutions of the embodiments of the present invention in conjunction with the drawings in 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 should fall within the scope of protection of the present invention.

[0020] For the light intensity flicker model generated by the atmospheric turbulence effect, a two-parameter model is commonly used from the perspectives of model accuracy and model complexity, which are the light intensity flicker Gamma-Gamma model and the wireless channel Lognormal-Rician model. Since the Lognormal-Rician model is more adaptable to a wider range of channels than the Gamma-Gamma model, the implementation steps of the present invention are described below using the Lognormal-Rician model as an example.

[0021] The probability density function expression of the Lognormal-Rician light intensity flicker model is: .

[0022] There are currently three methods for estimating the parameters of the Lognormal-Rician light intensity scintillation model, which are described in detail as follows.

[0023] (1) Prior art solution 1, generalized moment estimation method: X. Song and J. Cheng, ``Jointestimation of the Lognormal-Rician atmospheric turbulence model by thegeneralized method of moments,'' Opt. Commun., vol. 285, no. 24, pp. 4727-4732, Nov. 2012. Step 1: Calculate the sample value Moments of ,Right now: .

[0024] Its theory The error between the order central moments is expressed as: ,in The expression of the order central moment is: .

[0025] Step 2: Let the weight matrix ( dimensional identity matrix), solve the following optimization function to obtain the channel parameter estimation The initial value of: .

[0026] Step 3: Using Initial Estimates Calculate the residual vector, which is expressed as: .

[0027] And based on the residual vector, calculate the autocovariance matrix: .

[0028] Thus we get The covariance matrix of can be estimated as: ; Let the weight coefficient , update the weight matrix .

[0029] Step 4: Solve the following optimization function get Update , repeat steps 3 and 4, and finally get Estimates.

[0030] (2) Existing technology solution 2: L. Yang, J. Cheng, and JF Holzman, ``Maximumlikelihood estimation of the Lognormal-Rician FSO channel model,'' IEEE Photon. Technol. Lett., vol. 27, no. 15, pp. 1656-1659, Aug. 1, 2015. Step 1: Initialize parameter estimates ,in The expression is: .

[0031] Step 2: Calculate the mean of the following function: ; In the formula, .

[0032] Step 3: Update parameters ,in The expressions are: .

[0033] Repeat steps 2 and 3, and finally get Estimates.

[0034] (3) Prior Art Solution 3: MK Miao, and XF Li, `` Parameter Estimation of the Lognormal-Rician Channel Model Using Saddlepoint Approximation,'' IEEE Access, vol. 8, pp. 152924-152931, 2020. Step 1: Initialize parameter estimates ,in The expression is: .

[0035] Step 2: Sampling value Sort from small to large and record it as a sequence , and for each sample Find the saddle point of the following function ,Right now .

[0036] The variables The expression is: .

[0037] Step 3: , Saddle Point Substitute it into the following function to get the density approximation : .

[0038] Then use linear interpolation to get the normalization factor according to the following formula: : .

[0039] Step 4: Calculate the approximate log-likelihood function , .

[0040] Step 5: Solving Forward and backward differences , and use the gradient descent method to iterate the parameters.

[0041] like Figure 5 As shown, the present invention provides a method for estimating parameters of an atmospheric turbulence light intensity scintillation model, specifically: Step S1: Acquire multiple channel sampling values ​​of atmospheric turbulence, and acquire initial channel parameters of the multiple channel sampling values.

[0042] Step S1 is specifically as follows: after obtaining the atmospheric turbulence channel sampling values, they are sorted from small to large, and the initial channel parameters are obtained by using the moment estimation method and the light intensity flicker theoretical model, wherein the parameters in the channel sampling values ​​are obtained. The initial estimate expression of , , specifically: ; in, The photodetector obtains atmospheric turbulence channel sampling values, satisfying , Indicates the total sample size, 1 and 2 indicate the first and second samples.

[0043] Step S2: Based on the initial channel parameters and sample number k Values ​​generate multiple sets of light intensity flicker samples, and use The nearest neighbor density estimation method obtains the light intensity flicker estimation model for each group of light intensity flicker samples.

[0044] According to the initial channel parameters and the number of samples k The values ​​generate multiple sets of light intensity flicker samples, including: According to the estimated channel parameters, the receiving end processor generates multiple groups of light intensity flicker samples. Each group of generated light intensity flicker samples is sorted in ascending order from small to large. The specific expressions of the multiple groups of light intensity flicker samples are: ; in Subscript Indicates Group light intensity flicker samples, , Indicates Group light intensity flicker samples, Indicates samples, , each light intensity flicker sample obeys the parameter The wireless channel Lognormal-Rician model distribution, Indicates the number of each group of light intensity flicker samples.

[0045] use The nearest neighbor density estimation method obtains the light intensity flicker estimation model for each group of light intensity flicker samples, including: For any set of light intensity flicker sample data, Group as an example ,use Nearest neighbor density estimation method and current Worth getting sample points Density Estimation for: ; The variables , is the gamma function, Represents the dimension of each data sample. For the channel in the Lognormal-Rician model of wireless channels, , Representation sequence Zhongyu The nearest distance value.

[0046] For data points Perform linear interpolation to obtain the light intensity flicker estimation model and obtain the log-likelihood function value of the light intensity flicker estimation model , and integrate the interpolated curve to obtain the current normalization factor using the following formula: , specifically: .

[0047] Step S3: Eliminate the data in the channel sampling values ​​that are not in the range of the minimum value to the maximum value of each group of light intensity flicker samples, input the remaining channel sampling values ​​into the light intensity flicker estimation model, obtain the logarithmic likelihood function value of each group of light intensity flicker samples, and average multiple logarithmic likelihood function values. Specifically: For Group as an example , channel sampling sequence After elimination, we get , satisfy: ; in, M is the total sample size; for any set of light intensity flicker samples ,satisfy , Indicates the maximum number of light intensity flicker samples after culling. n is the number of samples.

[0048] Will Substitute into And find the average log-likelihood function value of each group , the specific expression is: .

[0049] Therefore, if the light intensity flicker sample is group, the average log-likelihood function is: ; in, is the log-likelihood function value.

[0050] Step S4: iterate the step of averaging the log-likelihood function values ​​to obtain the maximum average log-likelihood function value, and use the maximum average log-likelihood function value to represent the channel parameter value of the light intensity flicker estimation model.

[0051] S4 includes: Iterate through gradient descent Values ​​and channel parameters are used to obtain the maximum average log-likelihood function value. The specific expression is: .

[0052] By taking the difference of the maximum average log-likelihood function value, we can approximate each parameter The first-order derivative of is: ; In the formula, Indicates The estimated results.

[0053] Figure 1~Figure 4 By depicting the parameter estimation mean square error (MSE) performance in different scenarios, we can draw the following two important conclusions: 1) The two data-driven parameter estimation methods proposed in this invention are feasible. With the increase of the number of sampled samples and the number of light intensity flicker samples, the parameter estimation MSE is improved. For the case of low sampled samples, for example , only when the signal amplitude or variance When is small, the deviation between the estimated result and the true value will not be large.

[0054] 2) For data volume , it can be found based on The minimum KL divergence estimation method for the nearest neighbor is The estimated performance is slightly better than that based on The maximum likelihood estimation method for the nearest neighbor. , two methods Performance in Slight up and down fluctuation.

[0055] in, Figure 1~Figure 4 The diamonds in the figure represent the variance condition ( The minimum KL divergence estimate of the nearest neighbor), the rectangle represents the signal amplitude condition ( The minimum KL divergence estimate of the nearest neighbor), the circle represents the variance condition ( The maximum likelihood estimate of the nearest neighbor), asterisks indicate the signal amplitude condition ( The maximum likelihood estimate of the nearest neighbor).

[0056] Table 1 compares in detail several estimation methods of the Lognormal-Rician light intensity flicker model.

[0057] Table 1 Comparison results between different estimation methods of Lognormal-Rician light intensity flicker model

[0058] The present invention also includes an atmospheric turbulence light intensity flicker model parameter estimation system, comprising: a collection module, a model acquisition module, an average processing module and an estimation module.

[0059] The acquisition module is used to obtain multiple channel sampling values ​​of atmospheric turbulence and obtain the initial channel parameters of the multiple channel sampling values; the sorting module is used to sort the initial channel parameters and the number of samples. k Values ​​generate multiple sets of light intensity flicker samples, and use The nearest neighbor density estimation method is used to obtain the light intensity flicker estimation model for each group of light intensity flicker samples; the average processing module is used to eliminate the data in the channel sampling values ​​that are not in the range of the minimum value to the maximum value of each group of light intensity flicker samples, and the remaining channel sampling values ​​are input into the light intensity flicker estimation model to obtain the log-likelihood function value of each group of light intensity flicker samples, and the multiple log-likelihood function values ​​are averaged; the estimation module is used to iterate the steps of obtaining the average log-likelihood function value to obtain the maximum average log-likelihood function value, and the channel parameter value of the light intensity flicker estimation model is represented by the maximum average log-likelihood function value.

[0060] The present invention provides a storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps of a method for estimating parameters of an atmospheric turbulence light intensity scintillation model are realized.

[0061] According to the disclosed embodiments, the storage medium may be a non-volatile computer-readable storage medium, such as but not limited to: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, the storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, apparatus, or device.

[0062] The present invention provides a computer device, comprising a memory and a processor. The memory stores a program, and when the program is executed by the processor, the processor executes the steps of a method for estimating parameters of an atmospheric turbulence light intensity flicker model.

[0063] According to the disclosed embodiments, a computing device may communicate with one or more external devices (e.g., keyboards, pointing devices, Bluetooth communications, etc.), or with any device (e.g., routers, modems, etc.) that enables a computing device to communicate with one or more other computing devices.

[0064] like Figure 6 As shown in the following, a second estimation method for the atmospheric turbulence intensity scintillation model parameters is introduced, based on The nearest neighbor minimum KL divergence estimation method includes (1) checking the theoretical model of light intensity flicker to see whether it satisfies four conditional expressions; (2) data operations between atmospheric turbulence channel sampling values ​​(data); (3) calculating the initial (channel) parameter estimation and k value according to the theoretical model; (4) data operations between light intensity flicker samples and atmospheric turbulence channel sampling values; (5) calculating the average KL divergence between atmospheric turbulence channel sampling values ​​and generated light intensity flicker samples. (6) parameter iteration. The following is a detailed description of these six steps. Specifically: Step S1: Check if it exists The constructed light intensity flicker theoretical model satisfies four conditional expressions.

[0065] Among them, check whether there is The constructed light intensity flicker theoretical model satisfies four conditions, specifically: ; In the formula They represent the light intensity flicker sampling model and the light intensity flicker generation model respectively. In the subsequent simulation results, they are both Lognormal-Rician models.

[0066] Step S2: Calculate the distance between each sampling value and other sampling values ​​according to the atmospheric turbulence channel sampling values, sort them from small to large, and keep the first indivual.

[0067] S2 is specifically: Assumptions The photodetector acquires atmospheric turbulence channel sampling value. Data ,calculate ,right Sort from small to large and take the The distance size is obtained as follows: .

[0068] Step S3: using the moment estimation method and the light intensity scintillation theoretical model to obtain the initial channel parameter estimation of the atmospheric turbulence channel sampling values.

[0069] S3 is specifically: .

[0070] Step S4: Based on the estimated channel parameters, the receiving end processor generates multiple groups of light intensity flicker samples, calculates the distance between each atmospheric turbulence channel sampling value and each group of generated light intensity flicker samples (or other sampling values), sorts them from small to large, and retains the first sample each time. A numerical value.

[0071] S4 is specifically: ; in Subscript Indicates Group light intensity flicker samples, subject to the parameters Lognormal-Rician distribution. Group light intensity flicker samples , for each sample value calculate ,right Sort from small to large and take the The distance size is denoted as .

[0072] Step S5: Calculate the average KL divergence between the atmospheric turbulence channel sampling values ​​and the generated light intensity scintillation samples using a KL divergence estimator.

[0073] S5 is as follows: No. Group light intensity flicker samples and The KL divergence between is estimated as: ; in represents the dimension of each data sample. The average KL divergence expression of a group of light intensity flicker sample sequences is: .

[0074] Step S6: Repeat steps S4-S5, iterating using the gradient descent method value and channel parameters so that the average KL divergence is minimized.

[0075] S6 is specifically: ; The above optimization problem can be solved using the gradient descent method. The first-order derivatives of can be approximated by their respective differences, which are: ; In the formula Indicates The estimated results.

[0076] Table 2 compares in detail several estimation methods of the Lognormal-Rician light intensity flicker model.

[0077] Table 2 Comparison results between different estimation methods of Lognormal-Rician light intensity flicker model

[0078] The above content is a further detailed description of the present invention in combination with a specific preferred embodiment. For technicians in the technical field to which the present invention belongs, several simple deductions or substitutions can be made without departing from the concept of the present invention, which should be regarded as belonging to the protection scope of the present invention.

Claims

1. A method for estimating parameters of an atmospheric turbulence light intensity scintillation model, characterized in that: The steps include: Acquire multiple channel sampling values ​​of atmospheric turbulence, and acquire initial channel parameters of the multiple channel sampling values; According to the initial channel parameters and the number of samples k Values ​​generate multiple sets of light intensity flicker samples, and use The nearest neighbor density estimation method is used to obtain the light intensity flicker estimation model for each group of light intensity flicker samples; Eliminate the data that are not in the range of the minimum value to the maximum value of each group of light intensity flicker samples from the channel sampling values, input the remaining channel sampling values ​​into the light intensity flicker estimation model, obtain the log-likelihood function value of each group of light intensity flicker samples, and average multiple log-likelihood function values; The step of obtaining the average log-likelihood function value is iterated to obtain the maximum average log-likelihood function value, and the channel parameter value of the light intensity flicker estimation model is represented by the maximum average log-likelihood function value.

2. The method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to claim 1, characterized in that: The method of obtaining a plurality of channel sampling values ​​of atmospheric turbulence and obtaining initial channel parameters of the plurality of channel sampling values ​​is specifically as follows: Get channel sampling value parameters The initial estimate expression of , , specifically: ； in, The photodetector obtains atmospheric turbulence channel sampling values, satisfying , Indicates the total sample size, 1 and 2 indicate the first and second samples.

3. The method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to claim 2, characterized in that: According to the initial channel parameters and the number of samples k Generate multiple sets of light intensity flicker samples, including: The specific expressions of multiple groups of light intensity flicker samples are: ； in Subscript Indicates Group light intensity flicker samples, , Indicates Group light intensity flicker samples, Indicates samples, , each light intensity flicker sample obeys the parameter The wireless channel Lognormal-Rician model distribution, Indicates the number of each group of light intensity flicker samples.

4. The method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to claim 3, characterized in that: The utilization The nearest neighbor density estimation method obtains the light intensity flicker estimation model for each group of light intensity flicker samples, including: For any set of light intensity flicker sample data, Group as an example ,use Nearest neighbor density estimation method and current Worth getting sample points Density Estimation for: ； The variables , is the gamma function, Represents the dimension of each data sample. For the channel in the Lognormal-Rician model of wireless channels, , Representation sequence Zhongyu The nearest distance values; For data points Perform linear interpolation to obtain the light intensity flicker estimation model and obtain the log-likelihood function value of the light intensity flicker estimation model .

5. The method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to claim 4, characterized in that: The data that is not in the range of the minimum value to the maximum value of each group of light intensity flicker samples in the channel sampling values ​​are removed, and the remaining channel sampling values ​​are input into the light intensity flicker estimation model to obtain the logarithmic likelihood function value of each group of light intensity flicker samples, and the multiple logarithmic likelihood function values ​​are averaged, specifically: For Group as an example , channel sampling sequence After elimination, we get , satisfy: ； in, M is the total sample size; for any set of light intensity flicker samples ,satisfy , Indicates the maximum number of light intensity flicker samples after culling. n is the sampling number; Will Substitute into And find the average log-likelihood function value of each group , the specific expression is: ； in, is the normalization factor, and the specific expression is: ； If the light intensity flicker sample is group, then the average log-likelihood function value for: 。 6. The method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to claim 5, characterized in that: The step of iterating the average log-likelihood function value to obtain the maximum average log-likelihood function value, and characterizing the channel parameter value of the light intensity flicker estimation model by the maximum average log-likelihood function value, comprises: Iterate through the gradient descent method Values ​​and channel parameters are used to obtain the maximum average log-likelihood function value. The specific expression is: ； By taking the difference of the maximum average log-likelihood function value, we can approximate each parameter The first-order derivative of is: ； In the formula, Indicates The estimated results.

7. An atmospheric turbulence light intensity scintillation model parameter estimation system, characterized in that: include: An acquisition module, used for acquiring multiple channel sampling values ​​of atmospheric turbulence and acquiring initial channel parameters of the multiple channel sampling values; The model acquisition module is used to obtain the initial channel parameters and sample number. k Values ​​generate multiple sets of light intensity flicker samples, and use The nearest neighbor density estimation method is used to obtain the light intensity flicker estimation model for each group of light intensity flicker samples; An average processing module is used to remove data from the channel sampling values ​​that are not in the range from the minimum value to the maximum value of each group of light intensity flicker samples, input the remaining channel sampling values ​​into the light intensity flicker estimation model, obtain the logarithmic likelihood function value of each group of light intensity flicker samples, and average multiple logarithmic likelihood function values; The estimation module is used to iterate the step of averaging the log-likelihood function value to obtain the maximum average log-likelihood function value, and to represent the channel parameter value of the light intensity flicker estimation model through the maximum average log-likelihood function value.

8. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of a method for estimating parameters of an atmospheric turbulence light intensity scintillation model according to any one of claims 1 to 6 are implemented.

9. A computer device, characterized in that: It comprises a memory and a processor, wherein a program is stored in the memory, and when the program is executed by the processor, the processor executes the steps of a method for estimating parameters of an atmospheric turbulence light intensity scintillation model as claimed in any one of claims 1 to 6.

Images