CN108280520A - Atmospheric outline computational methods and device - Google Patents
Atmospheric outline computational methods and device Download PDFInfo
- Publication number
- CN108280520A CN108280520A CN201810157606.4A CN201810157606A CN108280520A CN 108280520 A CN108280520 A CN 108280520A CN 201810157606 A CN201810157606 A CN 201810157606A CN 108280520 A CN108280520 A CN 108280520A
- Authority
- CN
- China
- Prior art keywords
- chromosome
- value
- population
- fitness
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
- G01K13/024—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow of moving gases
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- General Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Physiology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides atmospheric outline computational methods and devices, are related to meteorological detection field.Atmospheric outline computational methods provided by the invention, use genetic algorithm, and it first generates initially for population, then, calculate the fitness value corresponding to each chromosome in current population, the then iterative mode for being selected, being intersected and being made a variation generates preferably target chromosome to integrate, and finally generates atmosphere temperature profile using target chromosome so that the atmosphere temperature profile of generation is more accurate.
Description
Technical field
The present invention relates to meteorological detection fields, in particular to atmospheric outline computational methods and device.
Background technology
Meteorological rocket is for detecting the meteorologic parameters such as the temperature of middle atmosphere parameter, pressure, density and wind speed, wind direction
Rocket in all kinds of sounding rockets, meteorological rocket is a kind of rocket that launch activities are most frequent, transmitting quantity is most, uses rocket
It is to collect a kind of important means of atmosphere data to carry out meteorological detection.
Later, with the rapid development of spacecraft and missile armament, there is an urgent need to the meteorologies under Middle and upper atmosphere environment to protect
Barrier, this has greatly pushed the progress of meteorological rocket sounding technology, wherein two achievements more outstanding are that thermosensitive resistance measurement is big
Temperature degree and falling ball method measure atmospheric density and wind field, this is also that technology is the most ripe so far, practical application is widest
Two methods.Thermosensitive resistance measurement atmospheric temperature technology is influenced by its principle, more than measurement error after 60km height drastically on
It rises, is only applicable to 60km detections below.And falling sphere Detection Techniques can measure the meteorologies such as 30~100km density, wind field, temperature and want
Element is a kind of detection method of great application value, and rocket sounding means most commonly used in the world at present.
Invention content
The purpose of the present invention is to provide atmospheric outline computational methods and devices.
In a first aspect, an embodiment of the present invention provides a kind of atmospheric outline computational methods, including：
It generates initially for population, includes initially the chromosome of predetermined number for population, chromosome is temperature value；
Calculate the fitness value corresponding to each chromosome in current population；
According to fitness value, selection operation is carried out to chromosome；
Crossover operation is carried out to chromosome, to generate new chromosome；
Mutation operation is carried out to chromosome, to generate new chromosome；
Judge whether present chromosome population meets stop condition；If so, according to fitness in present chromosome population
It is worth highest target chromosome and generates atmosphere temperature profile；It is each dyed in current population if it is not, then re-executing step and calculating
Fitness value corresponding to body.
With reference to first aspect, an embodiment of the present invention provides the first possible embodiments of first aspect, wherein step
Suddenly judge whether present chromosome population meets stop condition and include：
Judge whether maximum genetic algebra Ngen is more than the predetermined threshold value of maximum genetic algebra, and judges current chromosome
Whether the absolute value of the fitness of the highest target chromosome of fitness value is less than 1e-5 in population；
If maximum genetic algebra Ngen is more than the predetermined threshold value of maximum genetic algebra, and fitness in present chromosome population
The absolute value for being worth the fitness of highest target chromosome is less than 1e-5, then meets stop condition；If or, maximum genetic algebra
Ngen be more than maximum genetic algebra predetermined threshold value or the highest target chromosome of fitness value in present chromosome population it is suitable
The absolute value of response is less than 1e-5, then meets stop condition.
With reference to first aspect, an embodiment of the present invention provides second of possible embodiments of first aspect, wherein presses
The fitness value corresponding to each chromosome is calculated according to following formula：
Wherein T100Represent the calculating temperature at 100km height, T99.5The calculating temperature of height at 99.5km is represented, j is represented
Chromosome serial number, ξ represent the temperature gradient between 100km height and 99.5km height；
ξ values are calculate by the following formula acquisition：
Wherein, HmRepresent middle layer heights of roofs.
With reference to first aspect, an embodiment of the present invention provides the third possible embodiments of first aspect, wherein Hm
The height of temperature lowest point between 85km -100km.
With reference to first aspect, an embodiment of the present invention provides the 4th kind of possible embodiments of first aspect, wherein presses
It is by select probability according to what following formula calculated each individual：
Wherein, j is chromosome numbers (j=1,2 ..., 20), and Fitness is chromosome fitness value.
With reference to first aspect, an embodiment of the present invention provides the 5th kind of possible embodiments of first aspect, wherein hands over
The value of crossover probability in fork operation is calculated according to following formula：
Wherein, Pcross represents crossover probability, and iteration indicates that iterations, Ngen indicate the greatest iteration of setting
Number.
With reference to first aspect, an embodiment of the present invention provides the 6th kind of possible embodiments of first aspect, wherein hands over
In fork operation, child chromosome is generated according to following formula：
Wherein chrom1 and chrom2 are represented carries out crossover operation parent chromosome by selection, newchrom1 with
Newchrom2 represents two new child chromosomes, and rand represents the random number between 0 to 1.
With reference to first aspect, an embodiment of the present invention provides the 7th kind of possible embodiments of first aspect, wherein becomes
The value of mutation probability in ETTHER-OR operation is calculated according to following formula：
Wherein, Pmutation indicates that mutation probability, iteration indicate that iterations, Ngen indicate that the maximum of setting changes
Generation number.
With reference to first aspect, an embodiment of the present invention provides the 8th kind of possible embodiments of first aspect, wherein becomes
The numerical value of offspring individual in ETTHER-OR operation is calculated according to following formula：
Wherein, chrom and newchrom respectively represents parent individuality and offspring individual, and rand1 and rand2 represents two 0
Random number between to 1, Ub and Lb respectively represent the up-and-down boundary of temperature value at elemental height 30km.
Second aspect, the embodiment of the present invention additionally provide a kind of atmospheric outline computing device, including：
Generation module includes initially the chromosome of predetermined number for population for generating initially for population, and chromosome is
Temperature value；
Computing module, for calculating the fitness value in current population corresponding to each chromosome；
Selecting module, for according to fitness value, selection operation to be carried out to chromosome；
Cross module, for carrying out crossover operation to chromosome, to generate new chromosome；
Make a variation module, for carrying out mutation operation to chromosome, to generate new chromosome；
Judgment module, for judging whether present chromosome population meets stop condition；If so, according to current chromosome
The highest target chromosome of fitness value generates atmosphere temperature profile in population；It reworks if it is not, then triggering computing module.
Atmospheric outline computational methods provided in an embodiment of the present invention, use genetic algorithm, and first generate initial generation kind
Group, then, calculates the fitness value corresponding to each chromosome in current population, then iterative select, intersects
Preferably target chromosome is generated to integrate with the mode of variation, and finally generates atmospheric temperature exterior feature using target chromosome
Line so that the atmosphere temperature profile of generation is more accurate.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 is shown in the related technology, generates the acquisition flow of required parameter when atmospheric outline；
Fig. 2 shows the basic flow charts for the atmospheric outline computational methods that the embodiment of the present invention is provided；
Fig. 3 shows the basic module figure for the atmospheric outline computing device that the embodiment of the present invention is provided；
Fig. 4 shows the details module map of the judgment module for the atmospheric outline computing device that the embodiment of the present invention is provided.
Specific implementation mode
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause
This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below
Range, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing
The every other embodiment obtained under the premise of going out creative work, shall fall within the protection scope of the present invention.
Such as explanation hereinbefore, occur meteorological rocket falling sphere Detection Techniques in the related technology.
Near space meteorological rocket falling sphere Detection Techniques principle：Falling sphere is during declining mainly by gravity, resistance and floating
Force effect, according to Newton's second law, the equation of motion is represented by as shown in formula 1：
Wherein：
M is falling sphere quality,For acceleration of gravity, ρ is atmospheric density, CdFor damped coefficient, AsFor the cross section of falling sphere
Product,Speed for falling sphere relative to air, VsFor falling sphere volume；
It solves the component equation for solving the equation of motion (formula 1) in three-dimensional system of coordinate and arranges the table that density and wind can be obtained
Up to formula, the expression formula such as density is as shown in formula 2, and the expression formula of wind is as shown in formula 3-5：
Wz=0 (formula 5)；
Wherein, x indicates that east-west direction coordinate (eastwards for just), y indicate that North and South direction coordinate (northwards for just) and z indicate edge
Position coordinates perpendicular to terrestrial reference direction, Wx、WyAnd WzIndicate that component of the wind speed in earth axes, h are height, r represents ground
The radius of a ball.The speed extracted from track data is represented, andIndicate that acceleration, ρ are atmospheric density, m is to fall
Ball quality, CdFor damped coefficient, AsFor the cross-sectional area of falling sphere, g0Indicate sea level acceleration of gravity.
The damped coefficient C of falling spheredIt is most important for the calculating of atmospheric parameter.Damped coefficient is mainly by Mach number and Reynolds
Number determines.Mach number can be calculated by following formula, as shown in formula 6：
V is speed of the falling sphere relative to airɑ represents the velocity of sound, R
For dry air gas constant (R=287.05m2/s2K), T indicates that temperature, Ma are that Mach number (is defined as gas flow rate and locality
The ratio between velocity of sound, velocity of sound are related to temperature), γ indicates that specific heat ratio, specific heat ratio are represented by formula 7：
γ=A0+A1T+A2T2+A3T3+A4T4+A5T5(formula 7)；
Wherein, A0=1.371, A1=2.460 × 10-4, A2=-6.436 × 10-7, A3=5.200 × 10-10, A4
=-1.796 × 10-13, A5=2.182 × 10-17, T indicate temperature.
Reynolds number is the dimensionless number for characterizing fluid flow state, is defined as formula 8：
Wherein L is characterized length, indicates falling sphere diameter herein, Re is Reynolds number, and ρ is atmospheric density, and V is that falling sphere is opposite
In the speed of air, μ is Dynamic Viscosity, is determined by temperature：
μ=(1.458 × 10-6×T3/2)/(T+110.4) (formula 9)；
It can be seen that the calculating of Mach number and Reynolds number depends on temperature, T to indicate temperature.
On the basis of the speed of known falling sphere and acceleration, density, temperature and wind profile can pass through iteration
Method, which calculates, to be obtained, and detailed process is shown in Fig. 1, includes the following steps：
Step 1, obtaining mode density p0With model temperature T0；
Step 2, according to pattern density ρ0With model temperature T0Calculate Mach number Ma and reynolds number Re；
Step 3, according to pattern density ρ0, model temperature T0, Mach number Ma and reynolds number Re solve Nonlinear System of Equations, with
Obtain density p, temperature T, Horizontal Winds WxAnd Wy；
Step 4, judge whether to meet | ρ-ρ0|<10-5ρ0；If being judged as NO, 5 are thened follow the steps；If being judged as YES, hold
Row step 6；
Step 5, assignment as follows, and trigger step 2 and step 3 executes：ρ0=ρ0+(ρ-ρ0)×λ；0<λ<1；
T0=T；
Step 6, output density ρ, temperature T, Horizontal Winds WxAnd Wy。
In meteorological rocket falling sphere detection data processing system, temperature T0It is the call parameter for calculating Mach number and Reynolds number,
Meanwhile atmospheric temperature T is also the important environmental parameters for needing to know, it usually utilizes static(al) under conditions of known atmospheric density
Equilibrium equation and equation for ideal gases inverting obtain.The calculation formula of atmospheric temperature T such as formula 10 and 11：
Wherein P represents atmospheric pressure, and ρ represents atmospheric density, and R represents gas constant (R=287.0389Jkg-1K),
G represents vertical gravitational acceleration, and z represents height above sea level, and subscript i represents height above sea level sequence number (being divided into N layer heights).Existing skill
In art, the calculating of atmospheric temperature is general to be all from top to bottom iterated calculating to formula 10, i.e. elemental height is highest point
100km, it is therefore desirable to first give the temperature conjecture value on elemental height, in addition the integral and calculating in formula 10 is completed by following formula：
Therefore, under conditions of known atmospheric density profile, the calculation process of atmosphere temperature profile is as follows：
Step 1, it generates and initially guesses constant temperature angle value T0；
Step 2, the calculating temperature value of each height is calculated according to following formula
Wherein, TiWith ρiIndicate the calculating temperature value (atmospheric temperature value) and density of the i-th layer height, Ti-1With ρi-1Indicate the
The atmospheric temperature and density of i-1 layer heights, g indicate acceleration of gravity, and vertical separation is 500m in this programme.
Step 3, according to temperature value is calculated, atmosphere temperature profile is generated.
It can derive that above-mentioned technology calculates the error of atmosphere temperature profile according to error transfer function, first to formula 11
Total differential is unfolded, and obtains formula：
On the basis of the error of atmospheric density profile is zero, the variance for calculating temperature is represented by：
Atmospheric density is with exponential damping：Wherein H is about 7km, and substituting into formula 14 can
：
From the conjecture temperature value T known to formula 15 on elemental height0Error will with height decline and continuously decrease, pass through
Inventor verifies, and after elemental height declines 10-15km, the error for calculating temperature will drop to 1-2%.However due to 100km high
Initial temperature at degree is difficult to determine so that the calculating temperature under elemental height at least in 10km altitude ranges is unreliable, also
It is in the related technology, to calculate numerical value (T of the temperature within the scope of 90km-100kmi) not accurate enough.
For the problem that the numerical value of above-mentioned calculating temperature is not accurate enough.According to the error propagation of formula 15 it is found that if by
Under supreme iterative calculation atmosphere temperature profile, temperature error will rise with height and gradually increase, and generate in highest point very big
Error.Therefore most important, the provided herein technical solution of the selection of the initial guess temperature value at elemental height 30km
It is middle that optimal initial temperature value is determined using genetic algorithm.
Genetic algorithm is a kind of self-organizing, adaptive simulated nature biotechnology evolutionary process and mechanism and solve extreme-value problem
Artificial intelligence technology, basic thought be simulate nature genetic mechanism and theory of biological evolution and a kind of process searches for being formed most
The algorithm of excellent solution；It is suitable for, without any class letter for expressing or having expression, having achievable parallel computation behavior；It can be solved
Any type practical problem, is with a wide range of applications.Therefore, genetic algorithm be widely used in automatically controlling, computational science,
The fields such as pattern-recognition are suitable for solving complicated non-linear and hyperspace optimization problem.Although genetic algorithm is in many necks
Have successful application in domain, but its own there is also deficiencies, as local search ability is poor, there are immature oils and random trip
Phenomena such as walking causes convergence poor, the problems such as taking a long time just find optimal solution.To overcome basic genetic
The deficiency of algorithm, this programme uses dynamic self-adapting genetic algorithm, during algorithm carries out, is adjusted according to real-time results dynamic
Crossover operator and mutation operator.
The actual conditions that temperature computation is detected in conjunction with meteorological rocket falling sphere, using real coding mode, dynamic self-adapting is lost
Propagation algorithm specific implementation step is as follows, as shown in Figure 2：
Step S101：It generates initially for population, includes initially the chromosome of predetermined number for population, chromosome is temperature
Value.Genetic algorithm it needs to be determined that unique parameters be elemental height 30km place temperature value, value range be 150K (Kai Er
Text) to 350K (preferably 200K-300K), generate the initial temperature value of predetermined number (such as 20) at random in this value range
(chromosome), as genetic algorithm initially for population.
Step S102：Calculate the fitness value corresponding to each chromosome in current population.The initial temperature at known 30km
Under the condition (initial temperature value generated in step S101) of angle value, each initial temperature value institute can be iterated to calculate using formula 11
Corresponding atmosphere temperature profile.Then, (30km-100km's is big for the atmosphere temperature profile corresponding to each initial temperature value
Gas temperature profile) and the Atmospheric Characteristics of middle layer top and low-heat layer establish fitness and calculate function, and then be calculated each
Fitness value corresponding to initial temperature value.It calculates used in the fitness value corresponding to each chromosome (initial temperature value)
It is as follows that fitness calculates function：
Wherein T100Represent the calculating temperature at 100km height, T99.5The calculating temperature of height at 99.5km is represented, j is represented
Chromosome serial number (totally 20 group chromosome), ξ represent the temperature gradient between 100km height and 99.5km height.Middle layer top one
As be defined as temperature lowest point between 85km and 100km height, at this time temperature gradient with height increase and switched to by negative gradient
Positive gradient, therefore ξ should be positive number and change with the variation of middle layer heights of roofs, ξ values can be calculate by the following formula acquisition：
Wherein, HmRepresent middle layer heights of roofs (height of temperature lowest point between 85km -100km).
Step S103：According to fitness value, selection operation is carried out to chromosome.Selection operation is selected high-quality by fitness
Individual and abandon worst individual, embody the principle of " survival of the fittest ", wherein pass through atmospheric outline calculate obtain fitness value
It is high-quality individual when smaller, and when fitness value is larger is worst individual.Using sequencing selection method, ideal adaptation angle value is taken just
Value or the numerical value difference degree between negative value and individual adaptation degree are without particular/special requirement, to all individuals in group by its adaptation
Degree size is ranked up, and is distributed each individual selected probability according to sequence.Each individual is by select probability：
Wherein, j is chromosome numbers (j=1,2 ..., 20), and Fitness is chromosome fitness value；
It can be seen that when the inverse of fitness value is smaller, the corresponding selected probability of chromosome is smaller, and fitness
Value inverse it is larger when, the corresponding selected probability of chromosome is larger, that is, according to the selected probability of chromosome with
Being proportionate property of fitness value corresponding to chromosome.
Step S104：Crossover operation is carried out to chromosome, to generate new chromosome.Intersection refers to intersecting to two
Chromosome be exchanged with each other its portion gene in some way, to formed two new chromosomes, it be generate new individual
Main method, it is resolved that the ability of searching optimum of genetic algorithm plays a crucial role in genetic algorithm.It is introduced in crossover operation
Crossover probability, for judging whether two chromosomes are intersected.The selection of Different Individual crossover probability directly influences algorithm
Convergence rate and convergence precision, so the selection of crossover probability is also critically important.Adaptive crossover mutation is with evolutionary process
Progress adaptively adjust, in the incipient stage of evolution, crossover probability to be selected larger, such rough search process has
It conducive to the diversity of holding population, and in the later stage, then needs to carry out careful search, prevents from brokenly changing optimal solution, accelerate convergence rate.
In the present solution, the value of crossover probability is determined according to formula 19：
Wherein, Pcross represents crossover probability, and iteration indicates that iterations, Ngen indicate the greatest iteration of setting
Number.Since only there are one parameters for chromosome in this problem, when selected two chromosome need to carry out crossover operation, meter
Calculation method is represented by：
Wherein chrom1 and chrom2 are represented carries out crossover operation parent chromosome by selection, newchrom1 with
Newchrom2 represents two new filial generations, and rand represents the random number between 0 to 1.
Step S105：Mutation operation is carried out to chromosome, to generate new chromosome.Variation refers to compiling individual chromosome
Other allele of the genic value locus on certain locus in sequence are replaced, to forming new
Body, it is the householder method for generating new individual, determines the local search ability of genetic algorithm.It needs to introduce in mutation operation and become
Different probability, for deciding whether into row variation.General mutation probability is all smaller, and mutation probability is excessive, to the destructive of solution
It is bigger, it is easy the optimal solution made loss, mutation probability is too small, is susceptible to precocious phenomenon.So adaptive variation
Probability can take variation pattern from big to small.In this way, the various of population can be kept by carrying out search extensively in the incipient stage
Property, will at the end of carry out careful search and optimal solution can be prevented to be destroyed.In turn, preferred embodiment provided herein
In, the value of mutation probability can be calculated according to following formula：
Wherein, Pmutation indicates that mutation probability, iteration indicate that iterations, Ngen indicate that the maximum of setting changes
Generation number.When that need to carry out mutation operation by selection parent individuality for one, its offspring individual can be calculated as：
Wherein, chrom and newchrom respectively represents parent individuality and offspring individual, and rand1 and rand2 represents two 0
Random number between to 1, Ub and Lb respectively represent the up-and-down boundary of temperature value at elemental height 30km.
Step S106, judges whether present chromosome population meets stop condition, if so, by present chromosome population
The highest target chromosome output of fitness value；If it is not, then re-executing step S102.
Specifically, judging whether stop condition includes present chromosome population：Judge whether maximum genetic algebra Ngen surpasses
The predetermined threshold value (Ngen=200) of maximum genetic algebra is crossed, and judges the highest mesh of fitness value in present chromosome population
Whether the absolute value for marking the fitness of chromosome is less than 1e-5.If above-mentioned two Rule of judgment meets at least one or two judgements
Condition is satisfied by, then present chromosome population meets stop condition.
After genetic algorithm end of run, the chromosome for having minimum fitness value is exactly the solution of genetic algorithm,
Using the temperature at its corresponding 30km height, accurate complete atmosphere temperature profile (30km-100km) can be obtained.
It is, after by the highest target chromosome output of fitness value in present chromosome population, can also hold
Row following steps：
According to target chromosome, atmosphere temperature profile is generated using following formula：
Wherein, TiWith ρiIndicate calculating temperature value (atmospheric temperature value, the i.e. target exported in step S106 of the i-th layer height
Chromosome) and density, Ti-1With ρi-1Indicate that the atmospheric temperature and density of the (i-1)-th layer height, g indicate acceleration of gravity, this programme
Middle vertical separation is 500m.
Corresponding with preceding method, present invention also provides a kind of atmospheric outline computing devices, as shown in Figures 3 and 4, packet
It includes：
Generation module 301 includes initially the chromosome of predetermined number for population for generating initially for population, dyeing
Body is temperature value；
Computing module 302, for calculating the fitness value in current population corresponding to each chromosome；
Selecting module 303, for according to fitness value, selection operation to be carried out to chromosome；
Cross module 304, for carrying out crossover operation to chromosome, to generate new chromosome；
Make a variation module 305, for carrying out mutation operation to chromosome, to generate new chromosome；
Judgment module 306, for judging whether present chromosome population meets stop condition；If so, according to current dye
The highest target chromosome of fitness value generates atmosphere temperature profile in colour solid population；If it is not, then triggering computing module 302 again
Work.
Preferably, judgment module includes：
First judging unit 402, for judge maximum genetic algebra Ngen whether be more than maximum genetic algebra default threshold
Value；
Second judgment unit 403, for judging the suitable of the highest target chromosome of fitness value in present chromosome population
Whether the absolute value of response is less than 1e-5；
Trigger element 401, for when the first judging unit 402 and second judgment unit 403 are judged as YES, then judging
Meet stop condition；And for when one in the first judging unit 402 or second judgment unit 403 is judged as YES, sentencing
It is disconnected to meet stop condition.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, it can also
It is that each unit physically exists alone, it can also be during two or more units be integrated in one unit.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer read/write memory medium.Based on this understanding, technical scheme of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be expressed in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes：USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic disc or CD.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (10)
1. a kind of atmospheric outline computational methods, which is characterized in that including：
It generates initially for population, includes initially the chromosome of predetermined number for population, chromosome is temperature value；
Calculate the fitness value corresponding to each chromosome in current population；
According to fitness value, selection operation is carried out to chromosome；
Crossover operation is carried out to chromosome, to generate new chromosome；
Mutation operation is carried out to chromosome, to generate new chromosome；
Judge whether present chromosome population meets stop condition；If so, most according to fitness value in present chromosome population
High target chromosome generates atmosphere temperature profile；If it is not, then re-executing step calculates each chromosome institute in current population
Corresponding fitness value.
2. according to the method described in claim 1, it is characterized in that, step judges whether present chromosome population meets stopping item
Part includes：
Judge whether maximum genetic algebra Ngen is more than the predetermined threshold value of maximum genetic algebra, and judges present chromosome population
Whether the absolute value of the fitness of the middle highest target chromosome of fitness value is less than 1e-5；
If maximum genetic algebra Ngen is more than the predetermined threshold value of maximum genetic algebra, and fitness value is most in present chromosome population
The absolute value of the fitness of high target chromosome is less than 1e-5, then meets stop condition；Or, if maximum genetic algebra Ngen is super
Cross the predetermined threshold value of maximum genetic algebra or the fitness of the highest target chromosome of fitness value in present chromosome population
Absolute value is less than 1e-5, then meets stop condition.
3. according to the method described in claim 1, it is characterized in that, being calculated according to following formula suitable corresponding to each chromosome
Answer angle value：
Wherein T100Represent the calculating temperature at 100km height, T99.5The calculating temperature of height at 99.5km is represented, j represents dyeing
Body serial number, ξ represent the temperature gradient between 100km height and 99.5km height；
ξ values are calculate by the following formula acquisition：
Wherein, HmRepresent middle layer heights of roofs.
4. according to the method described in claim 3, it is characterized in that, HmThe height of temperature lowest point between 85km -100km.
5. according to the method described in claim 1, it is characterized in that, according to each individual of following formula calculating by select probability
For：
Wherein, j is chromosome numbers (j=1,2 ..., 20), and Fitness is chromosome fitness value.
6. according to the method described in claim 1, it is characterized in that, the value of the crossover probability in crossover operation is according to following public affairs
Formula calculates：
Wherein, Pcross represents crossover probability, and iteration indicates that iterations, Ngen indicate the maximum iteration of setting.
7. according to the method described in claim 6, it is characterized in that, in crossover operation, child chromosome is given birth to according to following formula
At：
Wherein chrom1 and chrom2 is represented carries out crossover operation parent chromosome, newchrom1 and newchrom2 generations by selection
Two new child chromosomes of table, rand represent the random number between 0 to 1.
8. according to the method described in claim 1, it is characterized in that, the value of the mutation probability in mutation operation is according to following public affairs
Formula calculates：
Wherein, Pmutation indicates that mutation probability, iteration indicate that iterations, Ngen indicate the greatest iteration time of setting
Number.
9. according to the method described in claim 8, it is characterized in that, the numerical value of the offspring individual in mutation operation is according to following public affairs
Formula calculates：
Wherein, chrom and newchrom respectively represent parent individuality and offspring individual, rand1 and rand2 represent two 0 to 1 it
Between random number, Ub and Lb respectively represent the up-and-down boundary of temperature value at elemental height 30km.
10. a kind of atmospheric outline computing device, which is characterized in that including：
Generation module includes initially the chromosome of predetermined number for population, chromosome is temperature for generating initially for population
Value；
Computing module, for calculating the fitness value in current population corresponding to each chromosome；
Selecting module, for according to fitness value, selection operation to be carried out to chromosome；
Cross module, for carrying out crossover operation to chromosome, to generate new chromosome；
Make a variation module, for carrying out mutation operation to chromosome, to generate new chromosome；
Judgment module, for judging whether present chromosome population meets stop condition；If so, according to present chromosome population
The middle highest target chromosome of fitness value generates atmosphere temperature profile；It reworks if it is not, then triggering computing module.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810157606.4A CN108280520B (en) | 2018-02-24 | 2018-02-24 | Atmosphere profile calculation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810157606.4A CN108280520B (en) | 2018-02-24 | 2018-02-24 | Atmosphere profile calculation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108280520A true CN108280520A (en) | 2018-07-13 |
| CN108280520B CN108280520B (en) | 2020-07-17 |
Family
ID=62808695
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810157606.4A Active CN108280520B (en) | 2018-02-24 | 2018-02-24 | Atmosphere profile calculation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108280520B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114154333A (en) * | 2021-12-06 | 2022-03-08 | 西安电子科技大学 | Atmospheric temperature profile prediction method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101587545A (en) * | 2009-06-19 | 2009-11-25 | 中国农业大学 | Method and system for selecting feature of cotton heterosexual fiber target image |
| US20100082264A1 (en) * | 2008-09-29 | 2010-04-01 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | MLD-Modified Synthetic Ocean Profiles |
| CN101706883A (en) * | 2009-11-09 | 2010-05-12 | 北京航空航天大学 | Data mining method and device |
| CN102663499A (en) * | 2012-03-11 | 2012-09-12 | 西安电子科技大学 | Network community division method based on simulated annealing genetic algorithm |
| CN103792538A (en) * | 2014-02-24 | 2014-05-14 | 北京航空航天大学 | Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer |
-
2018
- 2018-02-24 CN CN201810157606.4A patent/CN108280520B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100082264A1 (en) * | 2008-09-29 | 2010-04-01 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | MLD-Modified Synthetic Ocean Profiles |
| CN101587545A (en) * | 2009-06-19 | 2009-11-25 | 中国农业大学 | Method and system for selecting feature of cotton heterosexual fiber target image |
| CN101706883A (en) * | 2009-11-09 | 2010-05-12 | 北京航空航天大学 | Data mining method and device |
| CN102663499A (en) * | 2012-03-11 | 2012-09-12 | 西安电子科技大学 | Network community division method based on simulated annealing genetic algorithm |
| CN103792538A (en) * | 2014-02-24 | 2014-05-14 | 北京航空航天大学 | Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer |
Non-Patent Citations (1)
| Title |
|---|
| 刘亚亚 等: "地基微波辐射计遥感大气廓线的BP神经网络反演方法研究", 《高原气象》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114154333A (en) * | 2021-12-06 | 2022-03-08 | 西安电子科技大学 | Atmospheric temperature profile prediction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108280520B (en) | 2020-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111914481B (en) | Slope stability prediction method based on improved PSO-RBF algorithm | |
| Sharma et al. | Identification of photovoltaic module parameters by implementing a novel teaching learning based optimization with unique exemplar generation scheme (TLBO-UEGS) | |
| CA2392063A1 (en) | Method and program for simulating a physical system using object-oriented programming | |
| CN107194040B (en) | Water quality monitoring network multi-objective optimization deployment method based on bidirectional water flow | |
| CN109933067A (en) | A kind of unmanned boat collision prevention method based on genetic algorithm and particle swarm algorithm | |
| Siddiqui et al. | Lab-scale, closed-loop experimental characterization, model refinement, and validation of a hydrokinetic energy-harvesting ocean kite | |
| CN114168796B (en) | Method for establishing high-altitude aerodynamic database of aircraft | |
| CN108062585A (en) | A kind of method that Function Extreme value is calculated based on a flying moth darts into the fire algorithm | |
| CN106503365B (en) | A kind of sector search method for SPH algorithm | |
| Ma et al. | Location of natural gas leakage sources on offshore platform by a multi-robot system using particle swarm optimization algorithm | |
| CN106886657B (en) | A kind of finite element model method for building up based on kriging function | |
| CN111880412A (en) | Reconfigurable robot zero and neural optimal control method based on single evaluation network | |
| Wong et al. | Can transfer neuroevolution tractably solve your differential equations? | |
| CN109871934A (en) | Feature selection approach based on the distributed parallel binary of Spark a flying moth darts into the fire algorithm | |
| CN108280520A (en) | Atmospheric outline computational methods and device | |
| CN108805287A (en) | The Gyro Random error compensating method and system of genetic wavelet neural network | |
| CN116151102A (en) | Intelligent determination method for space target ultra-short arc initial orbit | |
| Wang et al. | Feed-Forward neural network soft-sensor modeling of flotation process based on particle swarm optimization and gravitational search algorithm | |
| Rastgoo et al. | A novel study on forecasting the airfoil self-noise, using a hybrid model based on the combination of CatBoost and Arithmetic Optimization Algorithm | |
| Ito et al. | Design space exploration using Self-Organizing Map based adaptive sampling | |
| CN113312832B (en) | Quick forecasting method for dynamic response of hull plate frame structure under underwater explosion load based on machine learning | |
| Hoyos et al. | Airfoil shape optimization: Comparative study of meta-heuristic algorithms, airfoil parameterization methods and reynolds number impact | |
| CN114912331A (en) | Cabin reinforcing rib optimization method, device, equipment and medium | |
| CN114996947A (en) | Three-dimensional oil reservoir numerical simulation method | |
| Hassan et al. | A Data-Driven Approach for Predicting the Onset of Entrainment in Two-Fluid Stratified Systems During Selective Withdrawal Process Using Machine Learning Techniques |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190916 Address after: 100000 No. 26 North Qing Road, Beijing, Haidian District Applicant after: Chen Shu Chi Applicant after: Beijing Aierda Electronic Equipment Co.,Ltd. Applicant after: SHANGHAI CHANGWANG METEOTECH CO., LTD. Address before: 100000 room 1013, East Tower, Cui Wei building, 33 Fuxing Road, Haidian District, Beijing Applicant before: Beijing Aierda Electronic Equipment Co.,Ltd. Applicant before: SHANGHAI CHANGWANG METEOTECH CO., LTD. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |