CN109800543A - Optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen - Google Patents

Optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen Download PDF

Info

Publication number
CN109800543A
CN109800543A CN201910234014.2A CN201910234014A CN109800543A CN 109800543 A CN109800543 A CN 109800543A CN 201910234014 A CN201910234014 A CN 201910234014A CN 109800543 A CN109800543 A CN 109800543A
Authority
CN
China
Prior art keywords
optimal
overtone order
phase screen
comprehensive power
primary condition
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
Application number
CN201910234014.2A
Other languages
Chinese (zh)
Other versions
CN109800543B (en
Inventor
刘涛
张景芝
雷艳旭
朱聪
孙春阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Publication of CN109800543A publication Critical patent/CN109800543A/en
Application granted granted Critical
Publication of CN109800543B publication Critical patent/CN109800543B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Abstract

A kind of optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen, the accuracy chosen for improving optimal overtone order.Its technical solution is, primary condition is arranged first and chooses comprehensive power ratio for the method, then preliminary optimal overtone order is calculated according to primary condition and comprehensive power radiometer, approximate rounding finally is carried out to preliminary optimal overtone order according to the actual situation, obtains final optimal overtone order.The present invention determines optimal overtone order using primary condition and comprehensive power ratio, it is short not only to calculate the time, but also can rapidly and accurately choose required optimal overtone order according to the actual situation.The method of the present invention is not only suitable for von K á rm á n phase power spectral-density model, is also applied for modified von K á rm á n phase power spectral-density model, can provide convenience for the research of atmospheric turbulence effect.

Description

Optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen
Technical field
The present invention relates to a kind of atmospheric turbulance technical fields, especially a kind of to utilize subharmonic method simulated atmosphere turbulent flow phase Method using when screen, for choosing optimal overtone order, belongs to field of measuring technique.
Background technique
Free space optical communication due to big bandwidth, without frequency license, that transmission rate is fast, is easily installed erection etc. is excellent Point, so that it is widely used in the fields such as wireless sensing, navigation space flight, municipal LAN.But due to atmosphere rapids Stream can cause the fluctuating of the light intensity and phase of laser, and free space optical communication performance degradation is caused even to cause communication link It interrupts, it is therefore necessary to which atmospheric turbulence effect is studied.
Currently, the method for research and simulated atmosphere turbulent flow mainly has test analysis mathed, theoretical analysis and Method for Numerical, Wherein Method for Numerical is easily realized due to simple, is most widely used.The core of Method for Numerical research atmospheric turbulence effect The heart is such as to be based on Fast Fourier Transform (FFT) (Fast Fourier using random phase screen come simulated atmosphere turbulence effect Transform, FFT) the spectrum method of inversion, the Zernike polynomial method based on orthogonal polynomial, based on Random Mid-point Displacement point Shape method, based on random data member extension method of formation and based on the Method of Multiple Scales of wavelet analysis.Since the spectrum method of inversion has meter It calculates that speed is fast, the advantage applicable to various forms of atmospheric turbulance spectrum models, becomes most common phase screen simulation side One of method, subharmonic method especially therein.The accuracy for generating phase screen using subharmonic method is related with overtone order, with The increase of overtone order, the phase screen of generation are also more and more accurate.But after overtone order is more than certain numerical value, phase screen Accuracy no longer significantly improve, continuing growing overtone order only will increase simulation time, or even numerical value change can occur, therefore There are optimal values for overtone order.Currently, also seldom on how to accurately determine the research of optimal overtone order, early simulation is big When gas turbulent flow phase screen, add that how many subharmonic have no basis, such as when simulating Kolmogorov phase screen, Lane et al. addition 5 subharmonic, and Sedmak is added to 10 subharmonic.Later, Sedmak had studied overtone order and phase screen dimensions and big Relationship between gas turbulent flow external measurement gives the practical evaluation criteria of some non-quantitations.Then, Marcel Carbillet couple It is improved, and two kinds of qualitative assessment standards are defined --- comprehensive power ratio and structure function ratio, but comprehensive power is than only Suitable for a kind of this spectrum model of von K á rm á n, and using structure function than obtaining the result is that a big probable value, not accurate enough.
From the studies above as can be seen that the feasibility and validity of subharmonic simulated atmosphere turbulent flow phase screen are unquestionable, But when generating phase screen, how accurately to choose optimal overtone order and still have problem, accuracy is slightly inadequate.Therefore have The evaluation criteria that the new optimal overtone order of necessity design is chosen, improves the accuracy that optimal overtone order is chosen, so as to Generate accurate atmospheric turbulance phase screen.
Summary of the invention
It is an object of the invention to aiming at the disadvantages of the prior art, provide it is a kind of for simulated atmosphere turbulent flow phase screen most Excellent overtone order choosing method, to improve the accuracy that optimal overtone order is chosen.
Problem of the present invention is solved with following technical proposals:
A kind of optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen, the method are arranged initially first Condition simultaneously chooses comprehensive power ratio, then calculates preliminary optimal harmonic wave according to primary condition and comprehensive power radiometer Number finally carries out approximate rounding to preliminary optimal overtone order according to the actual situation, obtains final optimal overtone order.
The above-mentioned optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen, the specific steps of the method Are as follows:
A., primary condition is set, comprehensive power ratio cc is selectedΦ:
1. primary condition includes: the dimension D of phase screen, atmospheric turbulence outer scale L0
2. selecting comprehensive power ratio, αΦValue range be [0,1];
B. according to the comprehensive power ratio cc of the primary condition of setting and selectionΦ, calculate preliminary optimal overtone order pin_optimal:
C. according to the actual situation to pin_optimalApproximate rounding is carried out, final optimal overtone order p is obtainedoptimal:
1. when the requirement to phase screen accuracy is not less than the requirement to the speed of service, final optimal overtone order Are as follows:
2. when the requirement to phase screen accuracy is lower than the requirement to the speed of service, final optimal overtone order:
The above-mentioned optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen, the comprehensive power ratio setting For αΦ=0.99.
The present invention determines optimal overtone order using primary condition and comprehensive power ratio, and not only the calculating time is short, but also Required optimal overtone order can be rapidly and accurately chosen according to the actual situation.The method of the present invention is not only suitable for von K á rm á n Phase power spectral-density model is also applied for modified von K á rm á n phase power spectral-density model, can imitate for atmospheric turbulance The research answered provides convenience.
Detailed description of the invention
The invention will be further described with reference to the accompanying drawings and examples.
Fig. 1 is the flow chart that optimal overtone order is chosen using method of the invention;
Fig. 2 is two kinds of phase power spectral-density models of von K á rm á n and modified von K á rm á n;
Fig. 3 (a) is the phase indirect detection of modified von K á rm á n model with the variation of overtone order;
Fig. 3 (b) is the relative error function of modified von K á rm á n model with the variation of overtone order.
Text neutralizes symbol used in figure are as follows: ΦφIt (f) is phase power spectrum density, f is spatial frequency, DφIt (r) is phase knot Structure function, r are the distance between two sampled points, and p is overtone order, and α Φ is comprehensive power ratio, and D is the size of phase screen, L0It is atmospheric turbulence outer scale, pin_optimalIt is preliminary optimal overtone order, poptimalIt is final optimal overtone order.
Specific embodiment
It is used when using subharmonic method simulated atmosphere turbulent flow phase screen the present invention provides a kind of, it is optimal for choosing The method of overtone order, the method for the present invention calculating time is short, is suitable for two kinds of realities of von K á rm á n and modified von K á rm á n Phase power spectral-density model, and required optimal overtone order can be accurately chosen according to the actual situation.
The basic ideas of technical solution of the present invention are:
(1) comprehensive power ratio cc is selectedΦ.(2) according to comprehensive power ratio ccΦObtain preliminary optimal overtone order pin_optimal.(3) according to the actual situation to pin_optimalApproximate rounding is carried out, final optimal overtone order p is obtainedoptimal
Specific steps of the present invention are described in detail as follows:
Primary condition is arranged in step 1, selects comprehensive power ratio ccΦ:
1. primary condition includes: the dimension D of phase screen, atmospheric turbulence outer scale L0;Wherein, atmospheric turbulence outer scale L0: it is The main parameters for referring to description optical turbulence characteristic and its being applied in light propagation effect analysis.It counts and manages in Kolmogorov turbulent flow In, external measurement refers to the out to out in inertia area.
2. selecting comprehensive power ratio, it is traditionally arranged to be αΦ=0.99;
Step 2, by primary condition and comprehensive power ratio ccΦSubstitute into equation
, obtain preliminary optimal overtone order pin_optimal
Step 3, according to the actual situation to pin_optimalApproximate rounding is carried out, final optimal overtone order p is obtainedoptimal, As shown in Fig. 1 flow chart, following below scheme is carried out:
1. when phase screen needs high accuracy, to pin_optimalApproximate rounding upwards is carried out, i.e., final is optimal Overtone order
2. when the requirement to phase screen accuracy is lower than the requirement to the speed of service, to pin_optimalIt carries out to lower aprons It is rounded, i.e., final optimal overtone order
Technical solution for a better understanding of the present invention makees further the present invention below with reference to calculated examples It is bright.
Calculated examples:
1. Fig. 2 is the practical phase power spectral-density model figure of von K á rm á n and two kinds of n of modified von K á rm á, I For the modified von K á rm á n model shown in the dotted line, utilize method of the invention to choose optimal overtone order.Setting Dimension D=2m, the atmospheric turbulence outer scale L of phase screen0=10m, comprehensive power ratio are αΦ=0.99.
2. according to comprehensive power ratio ccΦ, the step as shown in Fig. 1 flow chart is carried out, preliminary optimal overtone order is obtained For pin_optimal=2.8419.
3. according to the actual situation to poptimalApproximate rounding is carried out, is obtained:
1. final optimal overtone order is when phase screen needs high accuracy (not less than to the speed of service) poptimal=3;
2. when the requirement to phase screen accuracy is lower than the requirement to the speed of service, final optimal overtone order is poptimal=2.
4. in order to verify poptimalAccuracy, when will use modified von K á rm á n model, different overtone orders are raw It is compared at the phase indirect detection and phase indirect detection theoretical value of phase screen, as a result as shown in Fig. 3 (a);Different harmonic waves Number generates relative error function result such as Fig. 3 between the phase indirect detection and phase indirect detection theoretical value of phase screen (b) shown in., it is apparent that the phase indirect detection for generating phase screen is humorous twice with increasing from Fig. 3 (a), Fig. 3 (b) Wave fast approaching phase indirect detection theoretical value, the relative error between them are also reduced rapidly;Increase more overtone orders Seem the further improvement that will not cause phase indirect detection and relative error function, but can from the enlarged drawing of Fig. 3 (b) Out, increase the corresponding relative error of triple-frequency harmonics relative error more corresponding than two subharmonic of increase and reduce about 0.2%, i.e., still have Certain improvement;But increasing more overtone orders but can not almost see improvement.So when phase screen needs When wanting high accuracy, should choose optimal overtone order is poptimal=3;When the requirement to phase screen accuracy is lower than to fortune When the requirement of scanning frequency degree, should choose optimal overtone order is poptimal=2.This is consistent with the result that the invention is calculated, because This demonstrates the accuracy of the method for the present invention.

Claims (3)

1. a kind of optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen, characterized in that the method is first Setting primary condition simultaneously chooses comprehensive power ratio, then calculates tentatively optimal according to primary condition and comprehensive power radiometer Overtone order finally carries out approximate rounding to preliminary optimal overtone order according to the actual situation, obtains final optimal harmonic wave Number.
2. a kind of optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen according to claim 1, It is characterized in, the specific steps of the choosing method are as follows:
A., primary condition is set, comprehensive power ratio cc is selectedΦ:
1. primary condition includes: the dimension D of phase screen, atmospheric turbulence outer scale L0
2. selecting comprehensive power ratio, αΦValue range be [0,1];
B. according to the comprehensive power ratio cc of the primary condition of setting and selectionΦ, calculate preliminary optimal overtone order pin_optimal:
C. according to the actual situation to pin_optimalApproximate rounding is carried out, final optimal overtone order p is obtainedoptimal:
1. when the requirement to phase screen accuracy is not less than the requirement to the speed of service, final optimal overtone order are as follows:
2. when the requirement to phase screen accuracy is lower than the requirement to the speed of service, final optimal overtone order:
3. a kind of optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen according to claim 1 or 2, It is characterized in that the comprehensive power ratio is set as αΦ=0.99.
CN201910234014.2A 2018-11-01 2019-03-26 Optimal harmonic frequency selection method for simulating atmospheric turbulence phase screen Active CN109800543B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201811297218 2018-11-01
CN2018112972182 2018-11-01

Publications (2)

Publication Number Publication Date
CN109800543A true CN109800543A (en) 2019-05-24
CN109800543B CN109800543B (en) 2023-06-16

Family

ID=66563881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910234014.2A Active CN109800543B (en) 2018-11-01 2019-03-26 Optimal harmonic frequency selection method for simulating atmospheric turbulence phase screen

Country Status (1)

Country Link
CN (1) CN109800543B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6023977A (en) * 1997-08-01 2000-02-15 Acuson Corporation Ultrasonic imaging aberration correction system and method
CN102164004A (en) * 2011-03-02 2011-08-24 中国科学院上海光学精密机械研究所 Method for fast simulating phase screen distorted by time-varying turbulence and used for light-spot positioning
CN103605860A (en) * 2013-11-28 2014-02-26 中国科学院软件研究所 Non-ideal light source simulation method based on equivalent phase screen method
CN106644104A (en) * 2016-10-13 2017-05-10 哈尔滨工业大学 Phase screen modeling method for discrete raindrop media based on spectral inversion method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6023977A (en) * 1997-08-01 2000-02-15 Acuson Corporation Ultrasonic imaging aberration correction system and method
CN102164004A (en) * 2011-03-02 2011-08-24 中国科学院上海光学精密机械研究所 Method for fast simulating phase screen distorted by time-varying turbulence and used for light-spot positioning
CN103605860A (en) * 2013-11-28 2014-02-26 中国科学院软件研究所 Non-ideal light source simulation method based on equivalent phase screen method
CN106644104A (en) * 2016-10-13 2017-05-10 哈尔滨工业大学 Phase screen modeling method for discrete raindrop media based on spectral inversion method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MARCEL CARBILLET ET AL.: "Numerical modeling of atmospherically perturbed phase screens: new solutions for classical fast Fourier transform and Zernike methods", 《APPLIED OPTICS》 *

Also Published As

Publication number Publication date
CN109800543B (en) 2023-06-16

Similar Documents

Publication Publication Date Title
Sobash et al. The impact of covariance localization for radar data on EnKF analyses of a developing MCS: Observing system simulation experiments
CN106557828A (en) A kind of long time scale photovoltaic is exerted oneself time series modeling method and apparatus
CN103874090B (en) A kind of gsm communication base station electromagnetic radiation Forecasting Methodology
CN102914623B (en) Fusing method of temperature compensation of humidity sensor
CN105572617B (en) A kind of calibration system and method for pulse power standard
Salomons Reduction of the performance of a noise screen due to screen-induced wind-speed gradients. Numerical computations and wind-tunnel experiments
CN104318593A (en) Simulation method and system of radar sea clusters
CN104486000B (en) Optical pulse signal simulation generation method influenced by atmospheric turbulence
CN104850752A (en) Parameter search range determining method based on adaptive random resonance
CN106772295A (en) A kind of distance and Doppler's matching process for the measurement of moving-target radar cross section
CN105373708B (en) A kind of Time-Frequency Analysis Method of the improvement generalized S-transform based on parameter optimization
CN106559146A (en) A kind of signal generator and signal generating method
CN105510706A (en) High-precision undersampling frequency measurement method
Tamura et al. Current-induced modulation of the ocean wave spectrum and the role of nonlinear energy transfer
CN107490464A (en) The back wave separation method of Nonlinear Wave based on addition of waveforms principle
CN109800543A (en) Optimal overtone order choosing method for simulated atmosphere turbulent flow phase screen
Zhang et al. Derivation and verification of three-dimensional wake model of multiple wind turbines based on super-Gaussian function
CN109655775A (en) A kind of multiple dimensioned calibration method of arbitrary waveform generator amplitude frequency sweep and device
CN106644104B (en) A kind of phase screen modeling method of the discrete raindrop medium based on the spectrum method of inversion
CN108763610A (en) A kind of atmospheric turbulance phase screen emulation mode based on the spectrum method of inversion
CN105334506B (en) A kind of method and apparatus that ocean surface wind speed is estimated based on radar return center line spectral intensity
CN107292031A (en) A kind of wave slope inversion method based on improved quasi- mirror image surface scattering model
CN110186533A (en) A kind of short-term tide prediction method in high-precision river mouth
CN106155980A (en) Overhead transmission line field wire coupling modeling method based on inhomogeneous excitation electric field
CN115542327A (en) Millimeter wave radar rain measurement method and system

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
GR01 Patent grant
GR01 Patent grant