CN108573120A - A kind of simple algorithm of improved 350GHz or less electric waves zenith decaying - Google Patents
A kind of simple algorithm of improved 350GHz or less electric waves zenith decaying Download PDFInfo
- Publication number
- CN108573120A CN108573120A CN201810173309.9A CN201810173309A CN108573120A CN 108573120 A CN108573120 A CN 108573120A CN 201810173309 A CN201810173309 A CN 201810173309A CN 108573120 A CN108573120 A CN 108573120A
- Authority
- CN
- China
- Prior art keywords
- zenith
- decaying
- height equivlent
- dry air
- 350ghz
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of simple algorithms of improved 350GHz or less electric waves zenith decaying, introduce the concept of a height equivlent, height equivlent is a kind of equivalent distance, this regards that uniform and its path attenuation is equal to zenith decaying herein as apart from upper air, and zenith, which can conveniently be calculated, using height equivlent decays.The temperature and humidity that improved height equivlent model introduces earth surface keeps it more preferable to the applicability of different regions.The present invention keeps the estimation of the air-ground path attenuation of millimeter wave, submillimeter region more accurate and has real-time because in these frequency ranges, the decaying on air-ground path be communication link design in cannot ignore.
Description
Technical field
The present invention relates to a kind of easy evaluation methods of 350GHz or less the electric wave zeniths of simplicity decaying, introduce surface temperature
Height equivlent models of the ITU-R P.676-11 in recommendation is improved with humidity, the computational accuracy higher for making zenith decay and tool
There is real-time.
Background technology
In recent years, the development of millimeter wave, submillimeter wave wave band hardware makes it possible that its application range expands, for example, high
Fast satellite-based communication, radar detection and remote sensing.Therefore, in the applications such as communication link design, remote sensing signal analysis, earth-to-space link
Decaying has to consider.
Since last century the eighties, many propagation models propose on the basis of spectrum experiment:Liebe was in 1980
MPM (Millimeter-wave Propagation Model) model is proposed for latter stage;Urban was proposed in 2004
Moliere-5 (Microwave Observation LIne Estimation and REtrieval, version 5) model;
J Mendork propose SARTre ([Approximate] Spherical Atmospheric Radiative for 2006
Transfer model) model;Japanese communication research institute proposes AMATERASU (Model for Atmospheric within 2008
Terahertz Radiation Analysis and Simulation model) model;2005, P.Erikson and S.A
Buehler proposes ARTS (Atmospheric Radiative Transfer Simulator) this model up to now
It is improving always.But in above model, the link calculation of these frequency range upper pathways decaying needs complete path warm and humid
Pressure or molecule content profile information, these information can not easily obtain.So there is an urgent need for a kind of approximate models of simplicity to satisfy the need
Diameter decaying carries out accurate and convenient estimation in real time.
For air-ground path, Suen can settle steam to analyze the transmission characteristic of THz waves using air, and analyze the whole world
Air-ground path atmospheric attenuation.In ITU-R P.676-11 recommendation, height equivlent is used for estimating for air path decaying
It calculates.Wherein height equivlent is a kind of equivalent distance, this is equal to zenith decaying herein apart from upper path attenuation.The latter's
Input condition is easier to obtain, but only includes the atmospheric pressure of earth's surface in model, this is under different weather and seasonal conditions
It is not applicable.Because what it is to THz wave propagation effect biggers is humidity and temperature.So in improved model, we introduce
Temperature and humidity improves the precision of model, and it is made to have more universality.
Invention content
The technical problem to be solved by the present invention is to:A kind of height equivlent based in ITU-R P.676 recommendation is provided
Model utilizes surface temperature, humidity and the improved precision higher of pressure, the stronger model of universality.The present invention is directed more
It is carried out with the window frequency for communicating and detecting potentiality.
In order to solve the above technical problems, technical scheme of the present invention includes the following steps:
1) utilize ITU-R P.676 in recommendation atmospheric attenuation model and ITU-R P.835 in different latitude, Various Seasonal
Temperature and humidity pressure profile, obtain just pushing away height equivlent in the case of different profiles;
2) it will just push away height equivlent to be compared with original height equivlent model, obtain difference;
3) difference curve is fitted, records fitting coefficient;
4) dry air height equivlent coefficient is modified using temperature, obtains improved dry air height equivlent model;
5) steam height equivlent coefficient is modified using temperature and humidity, obtains improved steam height equivlent mould
Type.
The beneficial effects of the invention are as follows:The temperature and humidity that improved height equivlent model introduces earth surface keeps its right
The applicability of different regions is more preferable, keeps the estimation of the air-ground path attenuation of millimeter wave, submillimeter region more accurate and has real-time.
Description of the drawings
Fig. 1 is a kind of flow chart of the simple algorithm of improved 350GHz or less electric waves zenith decaying of the present invention;
Fig. 2 is absorption peak frequency and window frequency schematic diagram;
Fig. 3 is the relative error curve graph of dry air height equivlent difference atmospheric outline, wherein (a) is normal atmosphere, (b)
For low latitudes, it is (c) middle latitude summer, is (d) high latitude summer;
Fig. 4 is middle latitude winter dry air comparison of computational results figure, wherein (a) is height equivlent, it is (b) that zenith is decayed,
(c) it is relative error;
Fig. 5 is high latitude winter dry air comparison of computational results figure, wherein (a) is height equivlent, it is (b) that zenith is decayed,
(c) it is relative error;
Fig. 6 is the relative error curve of steam height equivlent difference atmospheric outline, wherein (a) is normal atmosphere, it is (b) low
Latitude is (c) middle latitude summer, is (d) high latitude summer;
Fig. 7 is middle latitude winter steam comparison of computational results figure, wherein (a) is height equivlent, it is (b) that zenith is decayed, (c)
For relative error;
Fig. 8 is high latitude winter steam comparison of computational results figure, wherein (a) is height equivlent, it is (b) that zenith is decayed, (c)
For relative error.
Specific implementation mode
With reference to attached drawing 1, using the specific implementation of the invention for calculating path attenuation, steps are as follows:
Step 1, determine that frequency be window frequency, detector less than 350GHz is zenith direction;
As shown in Fig. 2, in addition to the water vapor absorption peak indicated in figure:22GHz, 183GHz and 380GHz;Oxygen absorption peak:
All it is to absorb smaller window frequency except 50-70GHz frequency ranges and 118GHz these frequency ranges.It is improved etc. in these frequency ranges
The precision for imitating height model is higher.
Step 2, dry air height equivlent is calculated using approximate formula proposed by the present invention and Earth Surface Atmosphere temperature, pressure;
It is write based on surface temperature and the dry air height equivlent approximate formula of pressure:
Wherein:
t3=c0+c1f+c2f2+c3f3+c4f4+c5f5 (4)
Here:
c1=-0.001968, c2=5.353 × 10–5、c3=-3.707 × 10–7、c4=1.03 × 10–9、c5=-1.02 ×
10–1, c0Value be given by:
c0=-0.14+0.008t (5)
Wherein, t is temperature, and unit is DEG C.The limitation restrictive condition of this model:
HereWherein ptotRepresent the total gas pressure of earth surface.
Step 3, using ITU-R P.676-11 recommendation attachment 1 calculate dry air attenuation rate:
P.676-11, ITU-R can obtain the calculating of Atmospheric Absorption to be up to 1THz frequencies in (attachment 1 of recommendation)
On radio wave attenuation rate in an atmosphere, the attenuation by absorption rate of air is
κa=κo+κw=0.1820fN " (f) dB/km (7)
Wherein:γOUnit (dB/km) is attenuation by absorption rate under the conditions of dry air (only under Oxygen Condition, due to big
Nitrogen caused by atmospheric pressure and disresonance Debye decaying).γw(unit dB/km) is that feature under the conditions of certain water-vapo(u)r density declines
Subtract.F (unit GHz) is frequency.
N " (f)=∑ SiFi+N″D(f) (8)
Wherein, SiIt is the intensity of i-th line, FiBe the curve shape factor and summation to extend to institute wired (for being higher than
For the f frequencies of 118.750343GHz oxygen lines, the oxygen line that should only will be above 60GHz recombination rates is included in summary);N″D
(f) be nitrogen adsorption caused by atmospheric pressure and Debye frequency spectrums drying continuous band.It is calculated merely with oxygen absorption Spectral Overlapping
In addition N "D(f) and indicate dry air attenuation by absorption rate.
Step 4, dry air height equivlent is multiplied to obtain the decaying of dry air zenith with attenuation rate:
Ao=hoγo (9)
AoIt is the decaying of dry air zenith.
Step 5, steam etc. is calculated using approximate formula proposed by the present invention and Earth Surface Atmosphere temperature, pressure and water-vapo(u)r density
Effect height;
Water-vapo(u)r density height equivlent approximate formula writing based on surface temperature, water-vapo(u)r density and pressure:
hw=1.66 (B+C (h1+h2+h3)) (10)
Here
Wherein,
+ 0.030287 ρ (15) of B=42.01244-7.27769ln (t+273.15)
+ 0.043366 ρ (16) of C=10.35695-0.03411 (t+273.15)
Here, ρ is water-vapo(u)r density, unit g/m3。
Step 6, using ITU-R P.676-11 recommendation attachment 1 calculate Water vapor attenuation ratio:
The attenuation by absorption rate of steam is obtained using steam spectral line Spectral Overlapping in step 3.
Step 7, steam height equivlent is multiplied to obtain the decaying of steam zenith with attenuation rate;
Aw=hwγw (17)
AwIt is the decaying of dry air zenith.
Step 8, dry air is added to obtain the air complete attenuation in this frequency with the decaying of steam zenith.Total air zenith declines
It is kept to:
A=Ao+Aw (18)
The effect of the present invention can be further illustrated by following emulation:
1. simulated conditions and emulation content:
Using ITU-R P.835 in different latitude, Various Seasonal temperature and humidity pressure profile to improved height equivlent model
It is verified.
2. the simulation experiment result:
Based on ITU-R P.835-5 in profile verify two kinds of height equivlent models, as shown in Fig. 3-Fig. 8.From opposite
The curve of error can be seen that in the case of different atmospheric outlines, and the precision of new model has different degrees of improvement, especially
It it is the winter of middle latitude and high latitude, relative error is all even reduced to 10% or so from 20% or more close to 50%.So of the invention
It is more applicable for different atmospheric conditions, and inputs parameter and is easy to get, the more real-time for making zenith decay.
Claims (4)
1. a kind of simple algorithm of improved 350GHz or less electric waves zenith decaying, including it is following
Step:
1) determine that frequency be window frequency, detector less than 350GHz is zenith direction;
2) dry air height equivlent is calculated using based on surface temperature and the dry air height equivlent approximate formula of pressure;
3) dry air attenuation rate is calculated to the computational methods of Atmospheric Absorption using in the ITU-R P.676-11 attachment 1 of recommendation;
4) dry air height equivlent is multiplied to obtain the decaying of dry air zenith with attenuation rate;
5) water-vapo(u)r density height equivlent approximate formula and Earth Surface Atmosphere temperature based on surface temperature, water-vapo(u)r density and pressure are utilized
Degree, pressure and water-vapo(u)r density calculate steam height equivlent;
6) Water vapor attenuation ratio is calculated;
7) steam height equivlent is multiplied to obtain the decaying of steam zenith with attenuation rate;
8) dry air is added to obtain the air complete attenuation in this frequency with the decaying of steam zenith.
2. the simple algorithm of improved 350GHz or less electric waves zenith decaying as described in claim 1, which is characterized in that step
2) the dry air height equivlent approximate formula based on surface temperature and pressure described in is:
Wherein:
t3=c0+c1f+c2f2+c3f3+c4f4+c5f5 (4)
c1=-0.001968, c2=5.353 × 10–5、c3=-3.707 × 10–7、c4=1.03 × 10–9、c5=-1.02 × 10–1,
c0Value be given by:
c0=-0.14+0.008t (5)
Wherein, t is temperature, and unit is DEG C;
The restrictive condition of this model:
When f < 70GHz
Here f is the frequency of wave, unit GHz,Wherein ptotRepresent the total gas pressure of earth surface.
3. the simple algorithm of improved 350GHz or less electric waves zenith decaying as described in claim 1, which is characterized in that step
5) the water-vapo(u)r density height equivlent approximate formula based on surface temperature, water-vapo(u)r density and pressure described in is:
hw=1.66 (B+C (h1+h2+h3)) (7)
Here
Wherein,
+ 0.030287 ρ (12) of B=42.01244-7.27769ln (t+273.15)
+ 0.043366 ρ (13) of C=10.35695-0.03411 (t+273.15)
Here, ρ is water-vapo(u)r density, unit g/m3。
4. the simplicity of improved 350GHz or less electric waves zenith decaying as described in claim 1
Algorithm, which is characterized in that the height equivlent and zenith decaying between relationship be:
A=hoγo+hwγw (14)
Wherein, γoAnd γwRespectively oxygen and Water vapor attenuation ratio.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810173309.9A CN108573120B (en) | 2018-03-02 | 2018-03-02 | Improved simple algorithm for zenith attenuation of electric waves below 350GHz |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810173309.9A CN108573120B (en) | 2018-03-02 | 2018-03-02 | Improved simple algorithm for zenith attenuation of electric waves below 350GHz |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108573120A true CN108573120A (en) | 2018-09-25 |
CN108573120B CN108573120B (en) | 2020-04-28 |
Family
ID=63576737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810173309.9A Active CN108573120B (en) | 2018-03-02 | 2018-03-02 | Improved simple algorithm for zenith attenuation of electric waves below 350GHz |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108573120B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110568431A (en) * | 2019-08-07 | 2019-12-13 | 中国电子科技集团公司第二十九研究所 | Method for calculating detection distance of Blake-like radar |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4611929A (en) * | 1983-03-21 | 1986-09-16 | The United States Of America As Represented By The Secretary Of The Navy | Satellite method for measuring sea surface temperature |
US6489915B1 (en) * | 1999-11-11 | 2002-12-03 | Raytheon Company | Microwave icing avoidance system |
CN104348445A (en) * | 2013-07-26 | 2015-02-11 | 中兴通讯股份有限公司 | Attenuation device, system and attenuation method |
-
2018
- 2018-03-02 CN CN201810173309.9A patent/CN108573120B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4611929A (en) * | 1983-03-21 | 1986-09-16 | The United States Of America As Represented By The Secretary Of The Navy | Satellite method for measuring sea surface temperature |
US6489915B1 (en) * | 1999-11-11 | 2002-12-03 | Raytheon Company | Microwave icing avoidance system |
CN104348445A (en) * | 2013-07-26 | 2015-02-11 | 中兴通讯股份有限公司 | Attenuation device, system and attenuation method |
Non-Patent Citations (4)
Title |
---|
ITU: "Recommendation ITU-R P.676-11", 《ITU-R》 * |
KAMAL HARB ET AL: "Systems Adaptation for Satellite Signal under Dust, Sand and Gaseous Attenuations", 《JOURNAL OF WIRELESS NETWORKING AND COMMUNICATIONS》 * |
T.V.OMOTOSHO ET AL: "Total Atmospheric Absorption of Fixed Satellite Communication Signal Due to Oxygen and Water Vapor in Nigeria", 《IEEE ANTENNAS AND PROPAGATION MAGAZINE》 * |
卢昌胜等: "0.34GHz大气传播测试系统设计", 《微波学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110568431A (en) * | 2019-08-07 | 2019-12-13 | 中国电子科技集团公司第二十九研究所 | Method for calculating detection distance of Blake-like radar |
CN110568431B (en) * | 2019-08-07 | 2023-03-14 | 中国电子科技集团公司第二十九研究所 | Method for calculating detection distance of Blake-like radar |
Also Published As
Publication number | Publication date |
---|---|
CN108573120B (en) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109387487A (en) | Short-wave infrared high-spectral data atmospheric methane fast inversion method under the conditions of cirrus | |
CN109709558B (en) | Physical inversion method of space-borne microwave remote sensing land overhead PWV | |
Leinonen et al. | Improved retrieval of cloud liquid water from CloudSat and MODIS | |
CN103792538A (en) | Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer | |
Engelen et al. | Characterization of water‐vapour retrievals from TOVS/HIRS and SSM/T‐2 measurements | |
CN107656279A (en) | A kind of rain detection radar radiation transmission system based on two-parameter volume size distribution | |
Luini et al. | Attenuation in nonrainy conditions at millimeter wavelengths: Assessment of a procedure | |
CN115758667A (en) | Method for jointly detecting and inverting sea atmosphere temperature and humidity profile | |
CN108827963A (en) | Soil reference materials measuring method based on multiple dimensioned thermal inertia model inverting | |
CN103743679A (en) | Method for increasing inversion speed of concentration of greenhouse-gas vertical column | |
CN111665218A (en) | Method for improving inversion accuracy of carbon dioxide differential absorption laser radar | |
Dairo et al. | Radio refractivity gradients in the lowest 100m of the atmosphere over Lagos, Nigeria in the rainy-harmattan transition phase | |
CN113111529B (en) | Infrared bright temperature simulation method integrating numerical mode and satellite microwave cloud inversion data | |
CN108573120A (en) | A kind of simple algorithm of improved 350GHz or less electric waves zenith decaying | |
Wang et al. | Long-term over-the-horizon microwave channel measurements and statistical analysis in evaporation ducts over the Yellow Sea | |
Westwater et al. | Combined surface-and satellite-based microwave temperature profile retrieval | |
CN113534194A (en) | Troposphere temperature and humidity profile inversion method combining GNSS and wind lidar | |
Zhigang et al. | Retrieving atmospheric temperature profiles from AMSU-A data with neural networks | |
Smith et al. | Vertical distributions of atmospheric water vapor from satellite infrared spectrometer measurements | |
Lee et al. | All‐sky microwave humidity sounder assimilation in the Korean Integrated Model forecast system | |
CN115575914B (en) | Multi-band dual-polarization weather radar observed quantity error quantification method | |
CN115270068B (en) | Method for quickly estimating atmospheric moisture delay of offshore inclined path based on buoy platform | |
Martin et al. | Tropospheric water and temperature retrieval for ASMUWARA | |
Singh et al. | Prediction of radio wave attenuation due to clouds using ANN and its business aspects | |
Moreau et al. | Atmospheric liquid water retrieval using a gated experts neural network |
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 |