阅读说明：本技术 一种基于遗传算法的非相干散射雷达数据反演方法 (Genetic algorithm-based incoherent scattering radar data inversion method ) 是由 王俊逸 乐新安 赵必强 丁锋 万卫星 于 2020-05-07 设计创作，主要内容包括：本发明涉及一种基于遗传算法的非相干散射雷达数据反演方法,根据非相干散射雷达探测得到的电离层散射信号自相关函数,结合信号的调制与滤波包络,对系统进行积分离散化,完成电离层多参量反演系统建模,确定适应度评价函数；将一个电离层参量作为一个元素,依次排列组成一个电离层参量样本,根据遗传算法设定的样本个数,产生电离层参量样本的初始集合,并根据适应度评价完成父代筛选；通过对父代电离层参量样本集合的交叉和变异计算,循环得到适应度最好的电离层参量样本即为最终电离层各个目标参量的反演结果。本发明具有结果准确且可靠性高,输入信息少、鲁棒性好等优点,可用于非相干散射雷达系统的数据处理、空间物理和空间天气研究等。(The invention relates to a genetic algorithm-based incoherent scattering radar data inversion method, which is characterized in that integral discretization is carried out on a system according to an ionospheric scattering signal autocorrelation function obtained by incoherent scattering radar detection in combination with signal modulation and filtering envelope, modeling of an ionospheric multi-parameter inversion system is completed, and a fitness evaluation function is determined; taking an ionospheric parameter as an element, sequentially arranging the ionospheric parameter to form an ionospheric parameter sample, generating an initial set of ionospheric parameter samples according to the number of samples set by a genetic algorithm, and completing parent screening according to fitness evaluation; and (3) circularly obtaining ionized layer parameter samples with the best fitness through the intersection and variation calculation of the parent ionized layer parameter sample set, namely the inversion result of each target parameter of the final ionized layer. The method has the advantages of accurate result, high reliability, less input information, good robustness and the like, and can be used for data processing, space physics and space weather research and the like of an incoherent scattering radar system.)

1. A genetic algorithm-based incoherent scattering radar data inversion method is characterized by comprising the following steps:

firstly, obtaining an autocorrelation function of an ionospheric scattering signal and modulation and filtering of the signal according to incoherent scattering radar detection, carrying out integral discretization on an ionospheric multi-parameter inversion system, completing system modeling and determining a fitness evaluation function;

secondly, target ionospheric parameters in the ionospheric multi-parameter inversion system are used as elements to form parameter sample vectors of the ionospheric, an initial parameter sample set representing ionospheric characteristics is generated according to the number of samples set by a genetic algorithm, and parent screening is completed according to fitness evaluation; and finally, circularly calculating the ionospheric parameter sample with the highest fitness through crossing and variation of the parent ionospheric parameter sample set to obtain a final inversion result of a plurality of ionospheric parameters.

2. The genetic algorithm-based incoherent scattering radar data inversion method of claim 1, characterized in that: the first step is specifically realized as follows:

(1) according to the detection principle of the incoherent scattering radar, when a modulation signal is received after being scattered by an ionized layer, the radar signal obtained through filtering treatment has the expression as follows:

where p (t) is the filter envelope of the signal receiver,for the received signal, env (t) is the envelope of the modulated signal, e (t, d)³r) is the randomly scattered signal of electrons per unit volume at time t, r is the distance of the target of detection from the radar, and s (r) is the propagation time of the signal from transmission to reception; the autocorrelation function of the ionospheric scattering signal obtained by the incoherent scattering radar detection is:

wherein the content of the first and second substances,are parameters related to radar hardware and geometry (R is receiver impedance,Is a single electron scattering cross section, and χ is a polarization angle, P_tTo transmit power, G₀(r),G₁(R) transmitting and receiving antenna gain at R, R₀,R₁For the distance of the probe target to the transmitting and receiving antennas, λ is the radar wavelength), σ (ω, r, a) is the ionospheric plasma scattering power spectrum;the time delay distance fuzzy function is obtained by performing Fourier transform calculation on the time delay distance fuzzy function, and the calculation formula of the time delay distance fuzzy function is as follows:

wherein, W_t ^A(μ, s (r)) p (t- μ) env (μ -s (r)) is a distance amplitude ambiguity function whose fourier transform product at different sampling times is equal to the spectral distance ambiguity function;

(2) spectrum recording fuzzy function W_t,t'(ω)＝∫_rd³rW_t,t'(ω, s (r)) is taken into formula (2), the frequency integral is divided into n sections according to an integral interval, the integral value is approximately equal to the sum of n divided rectangular areas with the function value in each section as the length and the section as the width, and the integral discretization expression of formula (2) is obtained as follows:

let Y_i＝z(t)z^*(t ') denotes the autocorrelation function of the received signal, t' ═ t + i △ t, i ═ 1, 2.. times, m, △ t are unit delays, m is the total number of delays of the autocorrelation function of the signal, H_ij＝P_0i×W_t,t+i△t(ω_j) Representing a system matrix calculated from radar coefficients and a spectral ambiguity function, S (a)_j＝σ(ω_jN represents a plasma scattering power spectrum of a certain distance gate, ω is a frequency coordinate of the spectrum, n is the number of frequency points of the spectrum, and a ═ Ne [ Ne, Ti, Te, V ], Te, and Te]The ion density is a vector composed of parameters of an ionized layer, Ne is electron density, Ti is ion temperature, Te is electron temperature, and V is plasma drift velocity, and the four parameters are inversion target parameters; the plasma scattering power spectrum is determined by the state of an ionized layer, a group of ionized layer parameter values are calculated to obtain a corresponding plasma scattering theoretical spectrum, and the specific calculation formula is as follows:

wherein the constant physical quantity λ_D0.0037798 is the electron debye length, σ₀＝4πr₀ ²≈10^-28m²(ii) a k is the radar wave vector, which is determined by the radar detection performance; p_iIndicates the ratio of the ionic components, coll_inPlasma collision frequencies are given by a theoretical model according to the conditions of different heights of an ionized layer; in addition, the constant coefficient for auxiliary calculation is set to be T by self₀,N₀,m₀,κ＝1.380658×10^-23J/K is Boltzmann constant; therefore, the ionosphere multi-parameter inversion system model of the incoherent scattering radar is obtained as follows:

Y_m＝H_m×nS(a)_n(6)

the system inverts parameters of an ionized layer through a signal autocorrelation function obtained by detection, wherein target parameters comprise electron density, ion temperature, electron temperature and plasma drift velocity;

(3) the inversion process is a nonlinear optimization problem, according to Bayesian theorem and Gaussian probability distribution hypothesis, the minimum residual error of the autocorrelation function is used as the target of the optimization problem, and the target is used as a fitness evaluation function of the ionosphere parameter sample, wherein the smaller the fitness evaluation function value is, the higher the fitness of the ionosphere parameter sample is, and the calculation formula of the fitness evaluation function is as follows:

F(a)＝∑(Y-HS(a))²(7)。

3. the genetic algorithm-based incoherent scattering radar data inversion method of claim 1, characterized in that: in the second step, the ionospheric multi-parameter inversion system model is solved by adopting a Genetic Algorithm (GA), the method does not need to input the initial values of all the inverted ionospheric target parameters, but randomly searches the optimal solution of the ionospheric parameters by parallel computing, and the method is specifically realized as follows:

(1) firstly, an ionosphere parameter is taken as an element to be sequentially arranged to form a parameter sample vector of the ionosphere, and the parameter sample vector is obtained according to the set number n of samples_pAnd the value limit of the ionospheric parameters: [ min, max]_4×1Expressing four target ionospheric parameters in the ionospheric multi-parameter inversion system, namely upper and lower limits of values of electron density, ion temperature, electron temperature and drift velocity, and calculating a sample set initial value representing ionospheric parameter informationCalculating the fitness value of each group of ionized layer parameter samples according to the fitness evaluation function, namely the sum of residual errors between the observed signal autocorrelation function and the signal autocorrelation function obtained theoretically;

(2) reserving n number of excellent ionospheric parameter samples according to setting_grAnd cross sample ratio p_crAnd respectively calculating the number of cross and variation samples:

wherein n is_crNumber of new samples of ionospheric parameters, n, generated for interleaving_chGenerating the number of new samples of ionospheric parameters for the variation; selecting parent samples for generating crossover and variant filial generations from initial samples carrying ionospheric parameter information, wherein the number of the parent samples is as follows:

n_f＝2×n_cr+n_ch(9)

finally, generating a parent from ionospheric parameter samples sorted according to fitness function values by adopting a gambling method, and specifically calculating as follows:

wherein p is₀Is the fitness probability, p, of the ionospheric parameter samples_FIs an accumulated value, p, of ionospheric parameter sample fitness probabilities_cTo select a parentThe fitness probability of the ionospheric parameter sample, rand, is in the range of [0, 1%]A generated random number;

(3) new ionospheric parameter information is obtained by cross-over and variance calculations on parent ionospheric parameter samples, wherein the cross calculation is to randomly select a certain ionospheric parameter for the screened parent ionospheric parameter sample, then interchanging the ionospheric parameter values of the two samples, and selecting an ionospheric parameter to be varied randomly in the original parent for variation calculation, the ionospheric information in the original sample is updated by changing the value of the ionospheric parameter, finally, the collection of the ionospheric parameter of the descendant generated by the inheritance is obtained, and determining whether the new ionospheric parameter sample continues to be inherited or not through the set inheritance ending constraint, if the new ionospheric parameter sample does not meet the inheritance ending condition, entering next filial generation inheritance circulation until the inheritance ending search obtains a sample with the best fitness in the ionospheric parameter sample set, wherein the value of each parameter is the inversion result of each parameter of the ionospheric layer.

Technical Field

The invention relates to a genetic algorithm-based incoherent scattering radar data inversion method which can be used for data processing of an incoherent scattering radar system, such as ionosphere multi-parameter monitoring, space weather research and the like.

Background

The ionized layer is a partial ionized plasma region 60km to 1000km above the ground surface, is connected with the magnetic layer and is connected with the atmosphere below, and is an important composition of the space of the day and the earth. The ionosphere is closely coupled with the upper and lower layers of the ring through various energy processes, and various energy outbreak processes of the magnetic layer and even the surface of the sun can respond in the ionosphere, so that the ionosphere is called a display screen of space weather. Meanwhile, because the ionosphere has influence on electric wave propagation, the ionosphere is also the key layer which is most closely related to human activities, and the time-space continuous detection of the ionosphere state has important significance on space weather research, satellite navigation and positioning, radio communication and the like.

The ground incoherent scattering radar is used as the most advanced ionosphere detection equipment at present, can realize ionosphere multi-parameter simultaneous detection, has the advantages of large height range coverage, high resolution and the like, and has a dominant position in ionosphere detection. In order to realize ionospheric multi-parameter detection based on incoherent scattering radar measurement, a data inversion method is one of the technical cores for realizing the function. The traditional least square method needs parameter initial information input from other sources, and simultaneously needs first-order Markov approximation on a nonlinear system, so that errors are easily introduced to influence inversion accuracy. Genetic Algorithm (GA) is an evolutionary algorithm simulating the evolution rule of 'survival of suitable persons' in the biological world, solves the limitations of the traditional method by hiding parallelism and the capability of global search for an optimal solution, and is widely applied to the solution of an optimization problem.

Disclosure of Invention

The technical problem to be solved by the invention is as follows: aiming at the data inversion problem of the incoherent scattering radar, the defects of the prior art are overcome, the incoherent scattering radar data inversion method based on the genetic algorithm is provided, the problem of ionospheric multi-parameter inversion by using incoherent scattering radar data is solved, and the ionospheric parameter result with high robustness and reliability is obtained.

The technical solution of the invention is as follows: a genetic algorithm-based incoherent scattering radar data inversion method is characterized by comprising the following steps:

firstly, obtaining an autocorrelation function of an ionospheric scattering signal according to incoherent scattering radar detection, and carrying out integral discretization on an ionospheric multi-parameter inversion system by using a known signal modulation and filtering method to complete system modeling and determine a fitness evaluation function; secondly, arranging a plurality of parameters of the ionized layer as elements to form an ionized layer parameter sample vector, generating an initial set of ionized layer parameter samples according to the number of samples set by a genetic algorithm, and completing parent screening by fitness evaluation; and finally, circularly screening the ionized layer parameter samples with high fitness through the intersection and variation calculation of the ionized layer parameter sample set to obtain the inversion results of a plurality of parameters of the ionized layer. The method comprises the following specific steps:

firstly, constructing a system model and a fitness evaluation function:

according to the detection principle of the incoherent scattering radar, when a modulation signal is received after being scattered by an ionized layer, the radar signal obtained through filtering treatment has the expression as follows:

where p (t) is the filter envelope of the signal receiver,for the received signal, env (t) is the envelope of the modulated signal, e (t, d)³r) is the randomly scattered signal of electrons per unit volume at time t, r is the distance of the target of detection from the radar, and s (r) is the propagation time of the signal from transmission to reception. Autocorrelation of ionospheric scatter signals detected by incoherent scatter radarThe function is:

wherein the content of the first and second substances,are parameters related to radar hardware and geometry (R is receiver impedance,Is a single electron scattering cross section, and χ is a polarization angle, P_tTo transmit power, G₀(r),G₁(R) transmitting and receiving antenna gain at R, R₀,R₁For the distance of the probe target to the transmitting and receiving antennas, λ is the radar wavelength), σ (ω, r, a) is the ionospheric plasma scattering power spectrum;the time delay distance fuzzy function is obtained by performing Fourier transform calculation on the time delay distance fuzzy function, and the calculation formula of the time delay distance fuzzy function is as follows:

wherein the content of the first and second substances,is a range-amplitude ambiguity function whose fourier transform product at different sampling times is equal to the spectral range ambiguity function. According to the definition of integral operation, the integral value in a certain interval is equal to the area value enclosed by the function curve and the variable interval in the rectangular coordinate. Spectrum recording fuzzy functionTaking the formula (2) into the system, dividing the frequency integral into n sections according to the integral interval, wherein the integral value is approximately equal to the sum of n divided rectangular areas with the function value in each section as the length and the section as the width to obtain the integral of the formula (2)The expression after discretization is:

let Y_i＝z(t)z^*(t ') denotes the autocorrelation function of the received signal, t' ═ t + i △ t, i ═ 1, 2.. times, m, △ t are unit delays, m is the total number of delays of the autocorrelation function of the signal, H_ij＝P_0i×W_t,t+i△t(ω_j) Representing a system matrix calculated from radar coefficients and a spectral ambiguity function, S (a)_j＝σ(ω_jN represents a plasma scattering power spectrum of a certain distance gate, ω is a frequency coordinate of the spectrum, n is the number of frequency points of the spectrum, and a ═ Ne [ Ne, Ti, Te, V ], Te, and Te]The ion density is a vector composed of parameters of an ionized layer, Ne is electron density, Ti is ion temperature, Te is electron temperature, and V is plasma drift velocity, and the four parameters are inversion target parameters; the plasma scattering power spectrum is determined by the state of an ionized layer, a group of ionized layer parameter values are calculated to obtain a corresponding plasma scattering theoretical spectrum, and the specific calculation formula is as follows:

wherein the constant physical quantity λ_D0.0037798 is the electron debye length,k is the radar wave vector, which is determined by the radar detection performance; p_iIndicates the ratio of the ionic components, coll_inPlasma collision frequencies are given by a theoretical model according to the conditions of different heights of an ionized layer; in addition, the constant value coefficient for auxiliary calculation is set to beκ＝1.380658×10^-23J/K is the Boltz constant. Therefore, the ionosphere multi-parameter inversion system model of the incoherent scattering radar is obtained as follows:

Y_m＝H_m×nS(a)_n(6)

the data inversion of the incoherent scattering radar is a nonlinear optimization problem for solving a plurality of parameters of an ionosphere through a signal autocorrelation function. According to Bayes theorem and Gaussian probability distribution hypothesis, the residual error of the autocorrelation function is used as the design target of the optimization problem, and is used as the fitness evaluation function of the ionospheric parameter sample, the smaller the fitness evaluation function value is, the higher the fitness of the ionospheric parameter sample is, and the calculation formula of the fitness evaluation function is as follows:

F(a)＝∑(Y-HS(a))²(7)

and step two, generating a primary ionospheric parameter sample and completing parent screening:

taking an ionized layer parameter as an element, connecting all parameters to form an ionized layer parameter sample vector, setting the number of initial samples, and calculating an initial sample set representing ionized layer parameter information according to the upper limit and the lower limit of the ionized layer parameterAnd (4) obtaining the fitness of each ionospheric parameter sample by the fitness evaluation function selected by the formula (6).

Then, the number n of excellent ionospheric parameter samples to be reserved is set_grAnd the sample ratio p to be cross-operated_crTo obtain cross offspring n respectively_crAnd the number n of variant progeny_ch：

And determining the number of parent ionospheric parameter samples used to generate crossover and variant offspring:

n_f＝2×n_cr+n_ch(9)

then according to the genetic principle of high-out-of-the-victory, in order to enable the ionospheric parameter information with high adaptability to be reserved through heredity, the ionospheric parameter samples are sequenced by utilizing the adaptability function values, and the ionospheric parameter samples with high adaptability exist more in the parent by adopting a gambling method, the specific calculation is as follows:

wherein p is₀Is the fitness probability, p, of the ionospheric parameter samples_FIs an accumulated value, p, of ionospheric parameter sample fitness probabilities_cTo select the fitness probability of the parent ionospheric parameter samples, rand is in the range of [0, 1%]The generated random number.

And step three, completing the calculation of a new generation of ionospheric parameter sample set:

in the genetic process, new ionospheric parameter information is obtained through the cross calculation and the variation calculation of parent ionospheric parameter samples, wherein the cross calculation is based on the random truncation of the screened parent ionospheric parameter samples, then the ionospheric parameter information of the two samples at the positions is exchanged, and the variation calculation is to replace the original ionospheric parameter value by introducing a new ionospheric parameter value, and finally obtain a filial ionospheric parameter sample set generated by the inheritance. And finally, determining whether a new ionospheric parameter sample set continues to be inherited or not through set inheritance ending constraints, if the set ionospheric parameter sample set does not meet inheritance ending conditions, entering next filial generation inheritance circulation until inheritance ending searching to obtain an ionospheric parameter sample with the best fitness in the ionospheric parameter sample set, wherein the sample consists of ionospheric parameters which are arranged, and the value of each parameter is the inversion result of a plurality of ionospheric parameters.

Compared with the prior art, the invention has the advantages that:

(1) compared with the traditional least square method, the genetic algorithm adopted by the invention does not need to carry out first-order Markov approximation on a nonlinear system any more, and does not need to calculate the initial value of the target ionospheric parameter through a theoretical model, so that the error introduced by model approximation is reduced, the input information of calculation is simplified, and the inversion precision and the calculation simplicity of incoherent scattering radar data are improved.

(2) The ionosphere multi-parameter synchronous inversion adopts a genetic algorithm, and the global optimal solution is searched based on the parallel computing principle instead of the local optimal solution based on the initial value information. The genetic algorithm has inherent implicit parallelism and better global optimization capability, adopts a probabilistic optimization method, can automatically acquire and guide an optimized search space, adaptively adjusts the search direction, does not need a determined rule, reduces the influence of the accuracy of external input information on the accuracy of an inversion result to a certain extent, and improves the reliability of the result.

Drawings

FIG. 1 is a processing flow chart of the incoherent scattering radar data inversion method based on the genetic algorithm.

Detailed Description

The present invention will be described in detail below with reference to the accompanying drawings and examples.

As shown in fig. 1, the method of the invention comprises the following steps:

according to the detection principle of the incoherent scattering radar, when a modulation signal is received after being scattered by an ionized layer, the radar signal obtained through filtering treatment has the expression as follows:

where p (t) is the filter envelope of the signal receiver,for the received signal, env (t) is the envelope of the modulated signal, e (t, d)³r) is the randomly scattered signal of electrons per unit volume at time t, r is the distance of the target of detection from the radar, and s (r) is the propagation time of the signal from transmission to reception. The autocorrelation function of the ionospheric scattering signal obtained by the incoherent scattering radar detection is:

wherein the content of the first and second substances,are parameters related to radar hardware and geometry (R is receiver impedance,Is a single electron scattering cross section, and χ is a polarization angle, P_tTo transmit power, G₀(r),G₁(R) transmitting and receiving antenna gain at R, R₀,R₁For the distance of the probe target to the transmitting and receiving antennas, λ is the radar wavelength), σ (ω, r, a) is the ionospheric plasma scattering power spectrum;the time delay distance fuzzy function is obtained by performing Fourier transform calculation on the time delay distance fuzzy function, and the calculation formula of the time delay distance fuzzy function is as follows:

wherein the content of the first and second substances,is a distance amplitude ambiguity function whose Fourier transform product at different sampling times is equal to the spectral distance ambiguity function; according to the definition of integral operation, the integral value in a certain interval is equal to the area value enclosed by the function curve and the variable interval in the rectangular coordinate. Spectrum recording fuzzy function W_t,t'(ω)＝∫_rd³rW_t,t'(ω, s (r)) is taken into formula (2), the frequency integral is divided into n sections according to an integral interval, the integral value is approximately equal to the sum of n divided rectangular areas with the function value in each section as the length and the section as the width, and the integral discretization expression of formula (2) is obtained as follows:

let Y_i＝z(t)z^*(t ') denotes the autocorrelation function of the received signal, t' ═ t + i △ t, i ═ 1, 2.. times, m, △ t are unit time delays, and m is the autocorrelation function of the signalTotal number of delays, H_ij＝P_0i×W_t,t+i△t(ω_j) Representing a system matrix calculated from radar coefficients and a spectral ambiguity function, S (a)_j＝σ(ω_jN represents a plasma scattering power spectrum of a certain distance gate, ω is a frequency coordinate of the spectrum, n is the number of frequency points of the spectrum, and a ═ Ne [ Ne, Ti, Te, V ], Te, and Te]The ion density is a vector composed of parameters of an ionized layer, Ne is electron density, Ti is ion temperature, Te is electron temperature, and V is plasma drift velocity, and the four parameters are inversion target parameters; the plasma scattering power spectrum is determined by the state of an ionized layer, a group of ionized layer parameter values are calculated to obtain a corresponding plasma scattering theoretical spectrum, and the specific calculation formula is as follows:

wherein the constant physical quantity λ_D0.0037798 is the electron debye length,taking a backscatter radar with a frequency of 500MHz as an example, the radar wave vector k is 20.958, and the ionospheric ion component ratio P_iAnd the plasma collision frequency coll_inGiven by the international ionospheric model (IRI). In addition, the coefficient for auxiliary calculation is set to be T by self₀＝300K、N₀＝10¹¹/m³、m₀＝30.5、ω₀8476.2 Hz. In conclusion, the ionospheric multi-parametric inversion system model of the incoherent scattering radar is as follows:

Y_m＝H_m×nS(a)_n(6)

the incoherent scattering radar data inversion is a nonlinear optimization problem for solving a plurality of parameters of an ionosphere through a signal autocorrelation function. According to Bayes theorem and Gaussian probability distribution hypothesis, the residual error of the autocorrelation function is used as the design target of the optimization problem, and is used as the fitness evaluation function of the ionospheric parameter sample, the smaller the fitness evaluation function value is, the higher the fitness of the ionospheric parameter sample is, namely the fitness evaluation function is:

F(a)＝∑(Y-HS(a))²(7)

then, the four ionized layer parameters are used as four elements to form a vector, called an ionized layer parameter sample, in the order of the electron density, the ion temperature, the electron temperature and the plasma drift velocity from front to back, the total number of the samples is set to be 150, the initial value of the sample set representing the ionized layer parameter information is calculated according to the value limit of the ionized layer parameters, and the fitness evaluation value of each ionized layer parameter sample is calculated according to the fitness evaluation function of the formula (7).

The number of the excellent ionospheric parameter samples to be reserved is set to be 8, the sample proportion of the cross operation is 80 percent, and the number n of the cross samples can be respectively obtained_crAnd the number n of variant samples_ch：

And determining the number of parent ionospheric parameter samples used to generate crossover and variation samples:

n_f＝2×n_cr+n_ch(9)

then according to the genetic principle of 'out-of-good' and according to the fitness evaluation function value ascending, the ionospheric parameter sample sets are arranged, and a betting method is adopted to enable ionospheric parameter samples with low fitness evaluation function values to exist more in the parent, and the specific calculation is as follows:

wherein p is₀Is the fitness probability, p, of the ionospheric parameter samples_FIs an accumulated value, p, of ionospheric parameter sample fitness probabilities_cTo select the fitness probability of the parent ionospheric parameter samples, rand is in the range of [0, 1%]The generated random number.

And finally, in the genetic process, new ionospheric parameter information is obtained through the crossing and variation calculation of parent ionospheric parameter samples, wherein the crossing calculation is to randomly select ionospheric parameters to be crossed based on the screened parent ionospheric parameter samples, then the parameter values of the two samples are exchanged, and the variation calculation is to introduce new ionospheric parameter values to replace the ionospheric parameter values of the original positions, so that a filial ionospheric parameter sample set generated by the heredity is finally obtained. The genetic termination conditions set were: the ionospheric parameter sample set is inherited to the maximum algebra, the ionospheric sample stagnation maximum algebra with optimal fitness and the minimum threshold value of the change of the evaluation function value of the optimal fitness, the inheritance ending condition is judged according to the ionospheric parameter sample set generated by the new generation of inheritance, whether to continue the new inheritance is determined, if the ionospheric parameter sample set does not accord with the inheritance ending condition, the next filial generation inheritance loop is entered, until the inheritance ending search is carried out, the ionospheric parameter sample with the optimal fitness in the ionospheric parameter sample set is obtained, and the value of each ionospheric parameter in the sample is the inversion result of a plurality of ionospheric parameters.

Those skilled in the art will appreciate that the invention may be practiced without these specific details.