CN109342365A - Atmospheric turbulance refractive index structure parameter path profile approximate measure method - Google Patents
Atmospheric turbulance refractive index structure parameter path profile approximate measure method Download PDFInfo
- Publication number
- CN109342365A CN109342365A CN201811382166.9A CN201811382166A CN109342365A CN 109342365 A CN109342365 A CN 109342365A CN 201811382166 A CN201811382166 A CN 201811382166A CN 109342365 A CN109342365 A CN 109342365A
- Authority
- CN
- China
- Prior art keywords
- variance
- lens
- random phase
- light
- centroid position
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The present invention discloses a kind of atmospheric turbulance refractive index structure parameter path profile approximate measure method.Transmission path is divided into three sections by this method, i.e., close to one section of receiving and transmitting terminals A, it is one section and intermediate close to one section of receiving and transmitting terminals B, and with three discrete random phase screens come the continuous atmospheric turbulance of approximate description.This method obtains the effective refractive index structural constant of each random phase screen by measuring light wave light intensity fluctuation and the arrival angle fluctuation of bi-directional transmission channel, substantially variation characteristic of atmospheric turbulance refractive index structure parameter along transmission path is described with this, to provide reference for light wave propagation in atmosphere Design of Engineering Systems.
Description
Technical field
The invention belongs to atmospheric channel optical signal transmission technical fields, are related to a kind of atmospheric turbulance refractive index structure parameter road
Diameter profile approximate measure method.
Background technique
Phenomena such as light intensity flashing and arrival angle fluctuation, can occur when transmitting through atmospheric turbulance for laser, they would generally give laser
Propagation in atmosphere engineering system adversely affects.Before designing Laser Atmospheric Transmission engineering system, estimation atmosphere is generally required
Turbulent Flow Effects degree.Refractive index structure parameter is a key parameter for calculating atmospheric turbulance transmission light wave statistic.For ground
For horizontal transport near face, in default of the path profile data of refractive index structure parameter, very when carrying out theoretical calculation
It can only all assume that refractive index structure parameter is uniformly distributed and (does not change along transmission path) along transmission path when more.However
In real atmosphere environment, usual refractive index structure parameter can change along horizontal route.In order to study various reasons
By the reliability of calculated result, it is necessary to which the path profile of measurement of Atmospheric Turbulence refractive index structure parameter is (i.e. along transmission path
Change curve).Will accurately measurement of Atmospheric Turbulence refractive index structure parameter path profile, can be in the different location of transmission path
Series of temperature oscillometer is placed at place, and according to big at each position of micro-temperature sensor measurement result inverting at each position
Gas turbulent flow refractive index structure parameter.However, this method is on the one hand at high cost, on the other hand implement more troublesome.Very
In more situations, it is only necessary to which an approximate evaluation of atmospheric turbulance refractive index structure parameter path profile can answer many reality
With generating help.The present invention discloses a kind of very simple atmospheric turbulance refractive index structure parameter path profile measurement side in realization
Method can easily obtain the approximate evaluation of atmospheric turbulance refractive index structure parameter path profile.
As shown in Figure 1, receiving and transmitting terminals A is located at z=0, receiving and transmitting terminals B is located at z=L, receiving and transmitting terminals A and receiving and transmitting terminals
It is atmospheric turbulence channels between B.According to SPIE publishing house " the Numerical Simulation of published in 2010
Optical Wave Propagation with Examples in MATLAB " and Anhui science tech publishing house 2005 go out
" propagation of the light in turbulent atmosphere " of version is it is found that pass to the deviation of log-amplitude fluctuation of the spherical wave of z=L from z=0 are as follows:
Wherein L indicates transmission range, and k=2 π/λ indicates light wave wave number, and λ indicates optical wavelength,It indicates at the z of position
Atmospheric turbulance refractive index structure parameter；The deviation of log-amplitude fluctuation of the spherical wave of z=0 is passed to from z=L are as follows:
The arrival angle fluctuation variance of the spherical wave of z=L is passed to from z=0 are as follows:
The arrival angle fluctuation variance of the spherical wave of z=0 is passed to from z=L are as follows:
Wherein D indicates the diameter in circular reception aperture.
As shown in Figure 1, with three random phase screens (i.e. random phase screen 001, random phase screen 002 and random phase screen
003) carry out the influence that the continuous atmospheric turbulance of approximate description generates light wave transmissions.It is available as follows according to formula (1) to formula (4)
As a result:
WhereinWithIt is random phase screen 001, random phase screen 002 and random phase screen 003 respectively
Effective refractive index structural constant can be used to approximatively estimate atmospheric turbulance refractive index structure parameter along the change of transmission path
Change.IfWithApproximately equal, then atmospheric turbulance refractive index structure parameter is constant along transmission path approximation.IfSignificantly greater thanWithThen atmospheric turbulance refractive index structure parameter in a transmission path between position take the larger value.IfSignificantly greater thanWithThen atmospheric turbulance refractive index structure parameter takes the larger value in the position close to receiving and transmitting terminals A.Such as
FruitSignificantly greater thanWithThen atmospheric turbulance refractive index structure parameter takes the larger value in the position close to receiving and transmitting terminals B.
It can be in the form of wushu (5) be written as follow matrix and multiplication of vectors to formula (7):
Wherein
By formula (8) as long as it is found that measuringWithIt can easily find outWithIt is worth noting that, formula (1) and formula (2) are the deviation of log-amplitude fluctuation under a condition of acceptance, can only survey in practice
Measure the light-intensity oscillation variance under the condition of acceptance of finite aperture With the light-intensity oscillation variance under condition of acceptanceBetween relationship it is writeable are as follows:
Wherein A (D) indicates the aperture averaging factor of spherical wave light intensity fluctuation.Larry C.Andrews was delivered in 1992
It is provided in " the Journal of the Optical Society of America A " phase page 597 to 600 of volume 94 paper as follows
Calculation formula:
The relationship of light-intensity oscillation variance and the deviation of log-amplitude fluctuation under point condition of acceptance under point condition of acceptance is:
So when measuring the light-intensity oscillation variance for passing to the spherical wave of z=0 from z=LAfterwards, it can directly count
Calculate corresponding deviation of log-amplitude fluctuation are as follows:
In actual measurement, the wavelength meeting of the wavelength and the light wave that z=0 is passed to from z=L of the light wave of z=L is passed to from z=0
There is nuance；Due to the wavelength dependency very little of light wave atmospheric turbulance transmission statistic characteristic, small wavelength difference is to measurement
As a result influence can be ignored.When actual measurement calculates, the light wave wave number k in the formula of front can be set to 4 π/(λ1+
λ2), λ1For passed to from z=0 z=L light wave wavelength, λ2For passed to from z=L z=0 light wave wavelength.
Summary of the invention
The object of the present invention is to provide a kind of atmospheric turbulance refractive index structure parameter path profile approximate measure method,
With discrete random phase screen come the continuous atmospheric turbulance of approximate description, by measure bi-directional transmission channel light wave light intensity fluctuation and
Arrival angle fluctuation obtains the effective refractive index structural constant of each random phase screen, in this, as atmospheric turbulance refractive index structures
Constant along the entire change feature of transmission path approximate description, to provide ginseng for light wave propagation in atmosphere Design of Engineering Systems
It examines.
The technical solution of this method is achieved in that a kind of atmospheric turbulance refractive index structure parameter path profile approximation is surveyed
Amount method, it is characterised in that: the light-intensity oscillation variance and arrival angle fluctuation variance for using two-way lightwave atmospheric turbulance to transmit first
Measuring system measures the light-intensity oscillation variance in two transmission directions and arrival angle fluctuation variance；Then, it is transmitted according to two
It is random that light-intensity oscillation variance and arrival angle fluctuation variance on direction calculate three be located at z=L/6, z=L/2, z=5L/6
The effective refractive index structural constant of phase screen, wherein z-axis is overlapped with transmission path, and L is path-length.Two-way lightwave atmosphere
The light-intensity oscillation variance of turbulence transfer and the receiving and transmitting terminals A of arrival angle fluctuation variance measuring system are located at z=0, two-way lightwave
The light-intensity oscillation variance of atmospheric turbulance transmission and the receiving and transmitting terminals B of arrival angle fluctuation variance measuring system are located at z=L.Positioned at z
Random phase screen at=L/6 is referred to as random phase screen P1, and the random phase screen at z=L/2 is referred to as random phase
Shield P2, the random phase screen at z=5L/6 is referred to as random phase screen P3.Random phase screen P1, random phase screen P2, with
The effective refractive index structural constant approximate description of machine phase screen P3 atmospheric turbulance refractive index structure parameter is along the total of transmission path
Body variation characteristic.
As shown in Fig. 2, the light-intensity oscillation variance and arrival angle fluctuation variance measuring system of the transmission of two-way lightwave atmospheric turbulance
It is made of receiving and transmitting terminals A and receiving and transmitting terminals B；Receiving and transmitting terminals A includes first laser device (101), the first single mode optical fiber (102),
One lens (103), the first optical splitter (104), the second lens (105), the first ccd detector (106), the first computer (107)；
Receiving and transmitting terminals B include second laser (201), the second single mode optical fiber (202), the third lens (203), the second optical splitter (204),
4th lens (205), the second ccd detector (206), second computer (207)；The optical signal that first laser device (101) issues
A001 is introduced into the first single mode optical fiber (102) and is emitted to atmospheric turbulance using the first lens (103) and the first optical splitter (104)
In channel, optical signal A001 transmission reaches and is incident on the through the second optical splitter (204) and the 4th lens (205) after receiving and transmitting terminals B
On two ccd detectors (206)；The optical signal A002 that second laser (201) issues is introduced into the second single mode optical fiber (202) and passes through again
It crosses the third lens (203) and the second optical splitter (204) is emitted in atmospheric turbulence channels, optical signal A002 transmission reaches sending and receiving end
It is incident on the first ccd detector (106) after machine A through the first optical splitter (104) and the second lens (105)；First ccd detector
(106) image exported is acquired and is handled by the first computer (107), and the image of the second ccd detector (206) output is by second
Computer (207) acquisition and processing；First ccd detector (106) is located at the focal position of the second lens (105), and the 2nd CCD is visited
The focal position that device (206) are located at the 4th lens (205) is surveyed, the exit end of the first single mode optical fiber (102) is located at the first lens
(103) focal position, the exit end of the second single mode optical fiber (202) are located at the focal position of the third lens (203).Light wave passes through
Effective light passing diameter that first optical splitter (104) reaches when the second lens (105) penetrate the second lens (105) again is passed through equal to light wave
Effective light passing diameter when the 4th lens (205) of the second optical splitter (204) arrival penetrate the 4th lens (205) again is crossed, that is, is received and dispatched
Equivalent round receiving aperture diameter of the equivalent round receiving aperture diameter of terminal A equal to receiving and transmitting terminals B.
Specific step is as follows for the measurement method:
Step 101: the light-intensity oscillation variance and arrival angle fluctuation variance measuring system for transmitting two-way lightwave atmospheric turbulance
It works normally；
Step 102: the continuous N frame light spot image exported with the first computer (107) acquisition the first ccd detector (106),
Centroid position B001 is saved in disk text by the centroid position B001 and grey scale pixel value summation B002 for calculating every frame light spot image
In part CRLAOA, gray value summation B002 is saved in disk file CRLSCI；At the same time, with second computer (207)
The continuous N frame light spot image for acquiring the second ccd detector (206) output, calculates the centroid position B003 of every frame light spot image,
Centroid position B003 is saved in disk file CRRAOA；
Step 103: for all centroid position B001 data saved in disk file CRLAOA, each centroid position
B001 data are made of the centroid position data in the direction x and the centroid position data in the direction y, using statistical method according to all matter
The centroid position variance in the heart position B001 data calculating direction xWith the centroid position variance in the direction yZ is passed to from z=L
The arrival angle fluctuation variance of=0 light wavef2For the focal length of the second lens (105)；
Step 104: for all centroid position B003 data saved in disk file CRRAOA, each centroid position
B003 data are made of the centroid position data in the direction x and the centroid position data in the direction y, using statistical method according to all matter
The centroid position variance in the heart position B003 data calculating direction xWith the centroid position variance in the direction yZ is passed to from z=0
The arrival angle fluctuation variance of the light wave of=Lf4For the focal length of the 4th lens (205)；
Step 105: calculating the mean value of all gray value summation B002 data saved in disk file CRLSCIAnd varianceIt enablesIt enablesWherein
D indicates that the equivalent round receiving aperture diameter of receiving and transmitting terminals A, k indicate light wave wave number, k=4 π/(λ1+λ2), λ1It is first
The wavelength of laser (101), λ2For the wavelength of second laser (201)；
Step 106: according to the following formula calculate random phase screen P1, random phase screen P2, random phase screen P3 it is effective
Refractive index structure parameter
Wherein
Step 107: record random phase screen P1, random phase screen P2, the effective refractive index structure of random phase screen P3 are normal
Number
In order to enable the light wave of transmission to be approximately spherical wave, it is desirable that w0,A<[L/(10k)]1/2And w0,B<[L/(10k)]1/2,
Wherein w0,AAnd w0,BThe radius of respectively the first lens (103) and the outgoing beam of the third lens (203).
The positive effect of the present invention is can relatively easily to obtain atmospheric turbulance refractive index structure parameter along transmission path
Entire change feature.Specifically, transmission path is divided into three sections by the method for the present invention, i.e., close to receiving and transmitting terminals A one section (from
Z=0 to z=L/3), intermediate one section (from z=L/3 to z=2L/3) and close to one section of receiving and transmitting terminals B (from z=2L/3 to z
=L).The corresponding effective refractive index structural constant of three sections of transmission paths can be measured using the method for the present invention Pass throughThe size of value may determine that substantially variation of the refractive index structure parameter along transmission path
Feature, so that helping light wave propagation in atmosphere Design of Engineering Systems personnel to analyse whether can be Air Close To The Earth Surface turbulence transfer path
On refractive index structure parameter regard as and do not change along transmission path approximation.
Detailed description of the invention
Fig. 1 is the schematic diagram for placing three random phase screens on the transmit path.
Fig. 2 is that the light-intensity oscillation variance of two-way lightwave atmospheric turbulance transmission and arrival angle fluctuation variance measuring system are illustrated
Figure.
Fig. 3 is the relation schematic diagram in the direction x and the direction y and ccd image.
Specific embodiment
In order to which the feature and advantage of this method are more clearly understood, this method is made into one combined with specific embodiments below
The description of step.In the present embodiment, first laser device (101) selects wavelength for the semiconductor laser of 808nm, second laser
(201) select wavelength for the semiconductor laser of 793nm, the two wavelength belong to 800nm wave band；First ccd detector
(106) and the output frame of the second ccd detector (206) is set as 1000Hz, the first computer (107) and second computer
(207) 30000 frame light spot image of continuous acquisition realizes that light-intensity oscillation variance and arrival angle fluctuation variance statistic calculate respectively.
The technical solution of this method is achieved in that a kind of atmospheric turbulance refractive index structure parameter path profile approximation is surveyed
Amount method, it is characterised in that: the light-intensity oscillation variance and arrival angle fluctuation variance for using two-way lightwave atmospheric turbulance to transmit first
Measuring system measures the light-intensity oscillation variance in two transmission directions and arrival angle fluctuation variance；Then, it is transmitted according to two
It is random that light-intensity oscillation variance and arrival angle fluctuation variance on direction calculate three be located at z=L/6, z=L/2, z=5L/6
The effective refractive index structural constant of phase screen, wherein z-axis is overlapped with transmission path, and L is path-length.Two-way lightwave atmosphere
The light-intensity oscillation variance of turbulence transfer and the receiving and transmitting terminals A of arrival angle fluctuation variance measuring system are located at z=0, two-way lightwave
The light-intensity oscillation variance of atmospheric turbulance transmission and the receiving and transmitting terminals B of arrival angle fluctuation variance measuring system are located at z=L.Positioned at z
Random phase screen at=L/6 is referred to as random phase screen P1, and the random phase screen at z=L/2 is referred to as random phase
Shield P2, the random phase screen at z=5L/6 is referred to as random phase screen P3.Random phase screen P1, random phase screen P2, with
The effective refractive index structural constant approximate description of machine phase screen P3 atmospheric turbulance refractive index structure parameter is along the total of transmission path
Body variation characteristic.
As shown in Fig. 2, the light-intensity oscillation variance and arrival angle fluctuation variance measuring system of the transmission of two-way lightwave atmospheric turbulance
It is made of receiving and transmitting terminals A and receiving and transmitting terminals B；Receiving and transmitting terminals A includes first laser device (101), the first single mode optical fiber (102),
One lens (103), the first optical splitter (104), the second lens (105), the first ccd detector (106), the first computer (107)；
Receiving and transmitting terminals B include second laser (201), the second single mode optical fiber (202), the third lens (203), the second optical splitter (204),
4th lens (205), the second ccd detector (206), second computer (207)；The optical signal that first laser device (101) issues
A001 is introduced into the first single mode optical fiber (102) and is emitted to atmospheric turbulance using the first lens (103) and the first optical splitter (104)
In channel, optical signal A001 transmission reaches and is incident on the through the second optical splitter (204) and the 4th lens (205) after receiving and transmitting terminals B
On two ccd detectors (206)；The optical signal A002 that second laser (201) issues is introduced into the second single mode optical fiber (202) and passes through again
It crosses the third lens (203) and the second optical splitter (204) is emitted in atmospheric turbulence channels, optical signal A002 transmission reaches sending and receiving end
It is incident on the first ccd detector (106) after machine A through the first optical splitter (104) and the second lens (105)；First ccd detector
(106) image exported is acquired and is handled by the first computer (107), and the image of the second ccd detector (206) output is by second
Computer (207) acquisition and processing；First ccd detector (106) is located at the focal position of the second lens (105), and the 2nd CCD is visited
The focal position that device (206) are located at the 4th lens (205) is surveyed, the exit end of the first single mode optical fiber (102) is located at the first lens
(103) focal position, the exit end of the second single mode optical fiber (202) are located at the focal position of the third lens (203).Light wave passes through
Effective light passing diameter that first optical splitter (104) reaches when the second lens (105) penetrate the second lens (105) again is passed through equal to light wave
Effective light passing diameter when the 4th lens (205) of the second optical splitter (204) arrival penetrate the 4th lens (205) again is crossed, that is, is received and dispatched
Equivalent round receiving aperture diameter of the equivalent round receiving aperture diameter of terminal A equal to receiving and transmitting terminals B.
Specific step is as follows for the measurement method:
Step 101: the light-intensity oscillation variance and arrival angle fluctuation variance measuring system for transmitting two-way lightwave atmospheric turbulance
It works normally；
Step 102: the continuous N frame light spot image exported with the first computer (107) acquisition the first ccd detector (106),
Centroid position B001 is saved in disk text by the centroid position B001 and grey scale pixel value summation B002 for calculating every frame light spot image
In part CRLAOA, gray value summation B002 is saved in disk file CRLSCI；At the same time, with second computer (207)
The continuous N frame light spot image for acquiring the second ccd detector (206) output, calculates the centroid position B003 of every frame light spot image,
Centroid position B003 is saved in disk file CRRAOA；
Step 103: for all centroid position B001 data saved in disk file CRLAOA, each centroid position
B001 data are made of the centroid position data in the direction x and the centroid position data in the direction y, using statistical method according to all matter
The centroid position variance in the heart position B001 data calculating direction xWith the centroid position variance in the direction yZ is passed to from z=L
The arrival angle fluctuation variance of=0 light wavef2For the focal length of the second lens (105)；
Step 104: for all centroid position B003 data saved in disk file CRRAOA, each centroid position
B003 data are made of the centroid position data in the direction x and the centroid position data in the direction y, using statistical method according to all matter
The centroid position variance in the heart position B003 data calculating direction xWith the centroid position variance in the direction yZ is passed to from z=0
The arrival angle fluctuation variance of the light wave of=Lf4For the focal length of the 4th lens (205)；
Step 105: calculating the mean value of all gray value summation B002 data saved in disk file CRLSCIAnd varianceIt enablesIt enablesWherein
D indicates that the equivalent round receiving aperture diameter of receiving and transmitting terminals A, k indicate light wave wave number, k=4 π/(λ1+λ2), λ1It is first
The wavelength of laser (101), λ2For the wavelength of second laser (201)；
Step 106: according to the following formula calculate random phase screen P1, random phase screen P2, random phase screen P3 it is effective
Refractive index structure parameter
Wherein
Step 107: record random phase screen P1, random phase screen P2, the effective refractive index structure of random phase screen P3 are normal
Number
In order to enable the light wave of transmission to be approximately spherical wave, it is desirable that w0,A<[L/(10k)]1/2And w0,B<[L/(10k)]1/2,
Wherein w0,AAnd w0,BThe radius of respectively the first lens (103) and the outgoing beam of the third lens (203).
Grey scale pixel value summation B002 in step 102 refers to the total of the gray value of all pixels of every frame light spot image
With.
The definition in the direction x and the direction y in step 103 and step 104 is referring to Fig. 3.
Claims (1)
1. a kind of atmospheric turbulance refractive index structure parameter path profile approximate measure method, it is characterised in that: first using two-way
The light-intensity oscillation variance of light wave atmospheric turbulance transmission measures in two transmission directions with arrival angle fluctuation variance measuring system
Light-intensity oscillation variance and arrival angle fluctuation variance；Then, according to the light-intensity oscillation variance and angle of arrival in two transmission directions
The effective refractive index structural constant that variance calculates three random phase screens being located at z=L/6, z=L/2, z=5L/6 is lied prostrate,
Middle z-axis is overlapped with transmission path, and L is path-length；The light-intensity oscillation variance of two-way lightwave atmospheric turbulance transmission and arrival
The receiving and transmitting terminals A of angle fluctuating variance measuring system is located at z=0, two-way lightwave atmospheric turbulance transmission light-intensity oscillation variance with
The receiving and transmitting terminals B of arrival angle fluctuation variance measuring system is located at z=L；Random phase screen at z=L/6 be referred to as with
Machine phase screen P1, the random phase screen at z=L/2 are referred to as random phase screen P2, the random phase at z=5L/6
Screen is referred to as random phase screen P3；Random phase screen P1, random phase screen P2, the effective refractive index structure of random phase screen P3 are normal
Count approximate description entire change feature of the atmospheric turbulance refractive index structure parameter along transmission path；
The light-intensity oscillation variance and arrival angle fluctuation variance measuring system of two-way lightwave atmospheric turbulance transmission are by receiving and transmitting terminals A and receipts
Originator machine B composition；Receiving and transmitting terminals A include first laser device, the first single mode optical fiber, the first lens, the first optical splitter, second thoroughly
Mirror, the first ccd detector, the first computer；Receiving and transmitting terminals B includes second laser, the second single mode optical fiber, the third lens,
Two optical splitters, the 4th lens, the second ccd detector, second computer；The optical signal A001 that first laser device issues is introduced into the
One single mode optical fiber is emitted in atmospheric turbulence channels using the first lens and the first optical splitter, and optical signal A001 transmission, which reaches, receives
Through on the second optical splitter and the 4th lens entrance to the second ccd detector after originator machine B；The optical signal that second laser issues
A002 is introduced into the second single mode optical fiber and is emitted in atmospheric turbulence channels using the third lens and the second optical splitter, optical signal
Through on the first optical splitter and the second lens entrance to the first ccd detector after A002 transmission arrival receiving and transmitting terminals A；First CCD is visited
The image of device output is surveyed by the first computer acquisition and processing, the image of the second ccd detector output is acquired by second computer
And processing；First ccd detector is located at the focal position of the second lens, and the second ccd detector is located at the focus position of the 4th lens
It sets, the exit end of the first single mode optical fiber is located at the focal position of the first lens, and it is saturating that the exit end of the second single mode optical fiber is located at third
The focal position of mirror；Effective light passing diameter that light wave reaches when the second lens penetrate the second lens again by the first optical splitter is equal to
Light wave reaches effective light passing diameter when the 4th lens penetrate four lens again by the second optical splitter, i.e. receiving and transmitting terminals A etc.
Imitate equivalent round receiving aperture diameter of the circular reception aperture diameter equal to receiving and transmitting terminals B；
Specific step is as follows for the measurement method:
Step 101: the light-intensity oscillation variance and arrival angle fluctuation variance measuring system for transmitting two-way lightwave atmospheric turbulance are normal
Work；
Step 102: the continuous N frame light spot image exported with first the first ccd detector of computer acquisition calculates every frame hot spot figure
The centroid position B001 and grey scale pixel value summation B002 of picture, centroid position B001 are saved in disk file CRLAOA,
Gray value summation B002 is saved in disk file CRLSCI；At the same time, defeated with second computer the second ccd detector of acquisition
Continuous N frame light spot image out, calculates the centroid position B003 of every frame light spot image, and centroid position B003 is saved in disk text
In part CRRAOA；
Step 103: for all centroid position B001 data saved in disk file CRLAOA, each centroid position B001 number
It is formed according to the centroid position data of centroid position data and the direction y by the direction x, using statistical method according to all centroid positions
The centroid position variance in the B001 data calculating direction xWith the centroid position variance in the direction yThe light of z=0 is passed to from z=L
The arrival angle fluctuation variance of wavef2For the focal length of the second lens；
Step 104: for all centroid position B003 data saved in disk file CRRAOA, each centroid position B003 number
It is formed according to the centroid position data of centroid position data and the direction y by the direction x, using statistical method according to all centroid positions
The centroid position variance in the B003 data calculating direction xWith the centroid position variance in the direction yThe light of z=L is passed to from z=0
The arrival angle fluctuation variance of wavef4For the focal length of the 4th lens；
Step 105: calculating the mean value of all gray value summation B002 data saved in disk file CRLSCIAnd variance
It enablesIt enablesWherein
D indicates that the equivalent round receiving aperture diameter of receiving and transmitting terminals A, k indicate light wave wave number, k=4 π/(λ1+λ2), λ1It is first
The wavelength of laser, λ2For the wavelength of second laser；
Step 106: calculating effective refraction of random phase screen P1, random phase screen P2, random phase screen P3 according to the following formula
Rate structural constant
Wherein
Step 107: the effective refractive index structural constant of record random phase screen P1, random phase screen P2, random phase screen P3
In order to enable the light wave of transmission to be approximately spherical wave, it is desirable that w0,A<[L/(10k)]1/2And w0,B<[L/(10k)]1/2,
Middle w0,AAnd w0,BThe respectively radius of the first lens and the outgoing beam of the third lens.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811382166.9A CN109342365B (en) | 2018-11-20 | 2018-11-20 | Approximate measuring method for atmospheric turbulence refractive index structure constant path profile |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811382166.9A CN109342365B (en) | 2018-11-20 | 2018-11-20 | Approximate measuring method for atmospheric turbulence refractive index structure constant path profile |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109342365A true CN109342365A (en) | 2019-02-15 |
| CN109342365B CN109342365B (en) | 2020-12-29 |
Family
ID=65316521
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811382166.9A Active CN109342365B (en) | 2018-11-20 | 2018-11-20 | Approximate measuring method for atmospheric turbulence refractive index structure constant path profile |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109342365B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113686817A (en) * | 2021-08-24 | 2021-11-23 | 桂林电子科技大学 | Non-uniform path atmospheric phase screen distribution method in marine aerial environment |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1229576A1 (en) * | 1984-07-23 | 1986-05-07 | Войсковая часть 21109 | Method of measuring object displacements |
| CN1945355A (en) * | 2006-10-27 | 2007-04-11 | 中国科学院武汉物理与数学研究所 | Atmospheric turbulance detection laser rader using position-sensitive detector |
| CN101109702A (en) * | 2007-07-28 | 2008-01-23 | 中国科学院安徽光学精密机械研究所 | Large-caliber laser isotope scanner and method for measuring structure constant atmosphere index of refraction |
| CN101241190A (en) * | 2007-04-25 | 2008-08-13 | 中国科学院安徽光学精密机械研究所 | Zernike modal based laser atmospheric turbulence intensity profile measuring systems and method |
| CN101957317A (en) * | 2010-09-02 | 2011-01-26 | 长春理工大学 | Altitude distribution mode measurer of refractive index structural constants of atmospheric turbulence |
| CN102607798A (en) * | 2012-03-28 | 2012-07-25 | 中国人民解放军63655部队 | Device and method for measuring inner scale of atmospheric turbulence |
| CN106506068A (en) * | 2016-10-19 | 2017-03-15 | 长春理工大学 | The measurement of asymmetry method of two-way atmospheric turbulance optical transport channel |
| CN107942406A (en) * | 2017-11-09 | 2018-04-20 | 中国科学院光电技术研究所 | A kind of stratified atmosphere turbulence intensity measuring method based on the constraint of total atmospheric coherence length |
-
2018
- 2018-11-20 CN CN201811382166.9A patent/CN109342365B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1229576A1 (en) * | 1984-07-23 | 1986-05-07 | Войсковая часть 21109 | Method of measuring object displacements |
| CN1945355A (en) * | 2006-10-27 | 2007-04-11 | 中国科学院武汉物理与数学研究所 | Atmospheric turbulance detection laser rader using position-sensitive detector |
| CN101241190A (en) * | 2007-04-25 | 2008-08-13 | 中国科学院安徽光学精密机械研究所 | Zernike modal based laser atmospheric turbulence intensity profile measuring systems and method |
| CN101109702A (en) * | 2007-07-28 | 2008-01-23 | 中国科学院安徽光学精密机械研究所 | Large-caliber laser isotope scanner and method for measuring structure constant atmosphere index of refraction |
| CN101957317A (en) * | 2010-09-02 | 2011-01-26 | 长春理工大学 | Altitude distribution mode measurer of refractive index structural constants of atmospheric turbulence |
| CN102607798A (en) * | 2012-03-28 | 2012-07-25 | 中国人民解放军63655部队 | Device and method for measuring inner scale of atmospheric turbulence |
| CN106506068A (en) * | 2016-10-19 | 2017-03-15 | 长春理工大学 | The measurement of asymmetry method of two-way atmospheric turbulance optical transport channel |
| CN107942406A (en) * | 2017-11-09 | 2018-04-20 | 中国科学院光电技术研究所 | A kind of stratified atmosphere turbulence intensity measuring method based on the constraint of total atmospheric coherence length |
Non-Patent Citations (9)
| Title |
|---|
| JASON D. SCHMIDT: "《Numerical Simulation of Optical Wave Propagation with Examples in MATLAB》", 31 December 2010, SPIE * |
| LARRY C. ANDREWS: "Aperture-averaging factor for optical scintillations of plane and spherical waves in the atmosphere", 《JOURNAL OF OPTICAL SOCIETY OF AMERICA A》 * |
| MEI HAI-PING 等: "Evaluating refractive index structure constant and characterizing spectrum of atmospheric turbulence", 《强激光与粒子束》 * |
| 娄岩 等: "高斯涡旋光束在大气湍流传输中的特性研究", 《中国光学》 * |
| 张守川 等: "激光雷达测量大气湍流廓线", 《强激光与粒子束》 * |
| 陈纯毅 等: "激光大气折射率结构常数测量实验与分析", 《红外与激光工程》 * |
| 陈纯毅 等: "激光大气湍流传输数值实验建模与计算机模拟", 《系统仿真学报》 * |
| 饶瑞中: "《光在湍流大气中的传播》", 31 December 2005, 安徽科学技术出版 * |
| 黄剑锋 等: "大气湍流折射率结构常数的测试技术", 《光通信技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113686817A (en) * | 2021-08-24 | 2021-11-23 | 桂林电子科技大学 | Non-uniform path atmospheric phase screen distribution method in marine aerial environment |
| CN113686817B (en) * | 2021-08-24 | 2022-10-28 | 桂林电子科技大学 | Non-uniform path atmospheric phase screen distribution method in marine aerial environment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109342365B (en) | 2020-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101435700B (en) | Method and device for testing infrared laser lighting source diverging angle | |
| Mahalati et al. | Effect of fog on free-space optical links employing imaging receivers | |
| CN106506068B (en) | The measurement of asymmetry method of two-way atmospheric turbulance optical transport channel | |
| CN105334519B (en) | More atmospheric parameters based on triple channel F-P etalons while detected with high accuracy laser radar system | |
| CN101814952B (en) | Method for testing lightwave transmission characteristics in atmospheric channel | |
| CN100374875C (en) | Multi-purpose laser altimeter measuring device | |
| CN103472456A (en) | Active imaging system and method based on sparse aperture compressing calculation correlation | |
| CN108594257A (en) | Tachogenerator and its scaling method based on Doppler effect and measurement method | |
| Xian et al. | Parameter optimization of a visibility LiDAR for sea-fog early warnings | |
| CN109342365A (en) | Atmospheric turbulance refractive index structure parameter path profile approximate measure method | |
| CN109297685A (en) | A kind of spectral transmittance test device and method for heavy caliber parallel light tube | |
| Kokkalis | Using paraxial approximation to describe the optical setup of a typical EARLINET lidar system | |
| CN207423494U (en) | A kind of laser communication terminal vacuum test system | |
| CN106769738A (en) | A kind of reflection type optical fiber powder concentration measurement system | |
| David et al. | Cross-sectional measuring of optical beam | |
| US4475816A (en) | Method for determining in situ the absorption coefficient of particulate media using pulsed laser technique | |
| Song et al. | Influence of fog on the signal to interference plus noise ratio of the imaging laser radar using a 16-element APD array | |
| CN207051484U (en) | A kind of device for measuring the non-homogeneous horizontal air Aerosol Extinction in sea level | |
| Kim et al. | Effects of solar noise on the detection range performance of a laser spot tracker | |
| Raj et al. | Seasonal investigation on prediction accuracy of atmospheric turbulence strength with a new model at Punalkulam, Tamil Nadu | |
| US3514612A (en) | Astronomical seeing conditions monitor | |
| Hua et al. | Numerical investigation on the ranging performance of forward pulsed laser radar affected by shock waves | |
| Zeisse et al. | Measurement of low-altitude infrared propagation | |
| Mei et al. | Principle and feasibility of using a 3M reflective film based fold pass laser speckle imaging system for measuring atmospheric optical turbulence | |
| CN205898418U (en) | Optics efficiency testing arrangement of infrared optical fiber coherent fiber bundle |
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 |