CN106788815A - A kind of short wave communication reliability estimation method based on many system detection datas - Google Patents
A kind of short wave communication reliability estimation method based on many system detection datas Download PDFInfo
- Publication number
- CN106788815A CN106788815A CN201611071301.9A CN201611071301A CN106788815A CN 106788815 A CN106788815 A CN 106788815A CN 201611071301 A CN201611071301 A CN 201611071301A CN 106788815 A CN106788815 A CN 106788815A
- Authority
- CN
- China
- Prior art keywords
- wave communication
- short wave
- communication link
- link
- ionosphere
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/373—Predicting channel quality or other radio frequency [RF] parameters
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/22—Scatter propagation systems, e.g. ionospheric, tropospheric or meteor scatter
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Quality & Reliability (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
The present invention relates to a kind of short wave communication reliability estimation method based on many system detection datas, the reflection point position of short wave communication link and all detecting links is calculated, choose the detecting link data close with short wave communication link ionosphere reflection point position；According to the detecting link ionospheric reflection point position and detecting link data that determine, forecast obtains the Ionospheric Parameters of detecting link ionospheric reflection point；The short wave communication link of determination and detecting link ionospheric reflection point position, and the Ionospheric Parameters of detecting link ionospheric reflection point that forecast is obtained, reconstruct obtain the Ionospheric Parameters of short wave communication link ionosphere pip；According to the Ionospheric Parameters for obtaining, the maximum usable frequency of short wave communication link is predicted using shortwave propagation Forecasting Methodology；According to the maximum usable frequency for obtaining, received field strength when short wave communication link is communicated using maximum usable frequency is calculated, analyze the reliability of short wave communication link.
Description
Technical field
The present invention relates to short wave communication field, more particularly to a kind of short wave communication reliability based on many system detection datas
Appraisal procedure.
Background technology
Ionosphere has the material impact that can not be ignored as the important component of near earth space environment to short wave communication
Electricity can effectively be evaded or be reduced to (such as causing signal fadeout, error code), the accurate judgement to the state feature in ionosphere
Absciss layer changes the influence to short wave communication link, it is ensured that the reliable communication of short-wave link.The prediction of ionization layer state feature is pre- at present
Report technical research is concentrated mainly in time-domain and spatial domain the forecast for carrying out Ionospheric Parameters, still lacks to enter with reference to measured data
The practical engineering method of row short wave communication link communication quality assessment.
The content of the invention
The technology of the present invention solve problem：In order to rationally and effectively utilize the shortwave detection data under different detection mechanism, carry
The utilization rate of detection data high, plans for the short-term frequency of short wave communication link and provides support, present invention utilization nonlinear prediction,
Hyperspace reconstructs scheduling theory, there is provided a kind of short wave communication link property analysis method based on vertical and oblique incidence souding data,
Effectively evade or reduce influence of the Ionospheric variability to short wave communication link, it is ensured that short wave communication link stability and high efficiency it is logical
Letter.
The technology of the present invention solution：A kind of short wave communication reliability estimation method based on many system detection datas, bag
Include following steps：
Step A：The reflection point position of short wave communication link and all detecting links is calculated, is chosen and short wave communication link electrical
The close detecting link data of absciss layer reflection point position；
Step B：The detecting link ionospheric reflection point position determined according to step A and detecting link data, forecast are obtained
The Ionospheric Parameters of detecting link ionospheric reflection point；
Step C：The short wave communication link determined according to step A and detecting link ionospheric reflection point position, and step B
The Ionospheric Parameters of the detecting link ionospheric reflection point that forecast is obtained, reconstruct obtains short wave communication link ionosphere pip
Ionospheric Parameters；
Step D：According to the short wave communication link ionosphere pip Ionospheric Parameters value that step C is obtained, short-pass is predicted
Believe the maximum usable frequency of link；
Step E：According to the maximum usable frequency that step D is obtained, calculate short wave communication link and entered using maximum usable frequency
Received field strength during row communication, analyzes the reliability of short wave communication link.
The step A is comprised the following steps：
Step A1, it is theoretical according to short wave communication physical model and propagation rays, calculate short wave communication link ionospheric reflection
Point longitude and latitude；
Step A2, obtains all detection station locations and detection type, calculates each detecting link ionospheric reflection point longitude and latitude,
And the short wave communication link ionosphere pip longitude and latitude determined according to step A1, choose and short wave communication link ionospheric reflection
The close detecting link in point position；
Step A3, the detecting link according to selected by step A2 reads its detecting link data of nearest 5 days, and according to spy
Surveyor's chain circuit-switched data calculates Ionospheric Parameters.
The step C is comprised the following steps：
Step C1, according to the detecting link ionospheric reflection point position that step A determines, calculates detecting link ionospheric reflection
Put the ionosphere distance with short wave communication link ionosphere pip；
Step C2, according to step B determine detecting link ionospheric reflection point Ionospheric Parameters value, analytical calculation its with
The correlative of Reference ionosphere model Ionospheric Parameters value difference value；
Step C3, according to step C1 and C2 result of calculation, asks for detecting link ionosphere anti-using regional restructuring equation group
The weight coefficient W of exit point and short wave communication link ionosphere pip；
Step C4, according to the weight coefficient W that step C3 determines, calculates short wave communication link ionosphere pip ionosphere ginseng
The correlative of number and Reference ionosphere model value difference, and short wave communication link ionosphere pip Ionospheric Parameters value.
The step D is comprised the following steps：
Step D1, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, calculates short
E layers of maximum usable frequency E (d) MUF of wave communication link；
Step D2, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, calculates short
F1 layers of maximum usable frequency F1 (d) MUF of wave communication link；
Step D3, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, with reference to short
The propagation distance situation of wave communication link, calculates F2 layers of maximum usable frequency F2 (d) MUF of short wave communication link；
Step D4, according to E (d) MUF, F1 (d) MUF and F2 (d) MUF that step D1, D2, D3 analyses are obtained, calculates shortwave
Communication link maximum usable frequency MUF.
The step D3 is specific as follows：
Step D31, judges whether the propagation path of short wave communication link is jumped away from d apart from d more than single-hop maximummax；
Step D32, if propagation path is jumped away from d apart from d less than or equal to single-hop maximummax, with reference to the electricity that step C analyses are obtained
Absciss layer parameter value, F2 layers of basic MUF is calculated according to single-hop mode；If propagation path is apart from d more than the maximum jump of single-hop away from dmax, press
It is F2 layers of basic MUF of multi-hop mode computation according to radio wave propagation.
The step E is specific as follows：
Step E1, judges the propagation path distance of short wave communication link, if propagation path distance is less than 7000km, performs step
Rapid E2；If propagation path distance is more than 9000km, step E3 is performed；Step E2, E3 and E4 are otherwise performed simultaneously；
Step E2, if propagation path distance is less than 7000km, analytical calculation short wave communication link uses maximum usable frequency
Various communication modes and synthesis received field strength when being communicated；
Step E3, if propagation path distance is more than 9000km, analytical calculation short wave communication link uses maximum usable frequency
Received field strength when being communicated；
Step E4, if propagation path distance is between [7000,9000] km, field is carried out on the basis of E2 and E3 is performed
Intrusion value, obtain the propagation path apart from upper use maximum usable frequency communicated when received field strength；
Step E5, according to the field intensity that step E2 or E3 or E4 are calculated, analyzes the reliability of short wave communication link.
The step E5 is specific as follows：
Step E51, according to the field intensity that step E2 or E3 or E4 are calculated, calculates receiving power intermediate value；
Step E52, according to the receiving power intermediate value that E51 is calculated, calculate signal to noise ratio median, signal to noise ratio it is upper very
The lower ten score values deviation of value deviation and signal to noise ratio；
Step E53, the signal to noise ratio median being calculated according to E52, the upper ten score values deviation of signal to noise ratio and signal to noise ratio
Lower ten score value deviations, calculate the reliability of short wave communication link.
Present invention advantage compared with prior art is：Invent the short-wave link reliability that maximum advantage is that foundation
Property forecasting procedure, can for Shortwave Communication System provide short wave communication link-quality assessment short-period forecast service, allow users to
The reception that future time instance shortwave is communicated via ionosphere is known in advance, selects other to lead to when being shortwave reliability not high
Letter mode changes frequency offer foundation, and can effectively evade cannot carry out the economic loss that short wave communication causes.Simultaneously in short-pass
In the case that letter link sets up connection or communication interruption suddenly for the first time, accurate working frequency parameter can be faster obtained,
It is quick to set up or recover communication, normal, the Efficient Operation of short wave communication network are ensured, Ionospheric variability is reduced to short wave channel
The influence of optimized frequency.
(1) the vertical and class data of oblique incidence souding two have been merged to be analyzed for short wave communication link property；
(2) refer to the theoretical short-period forecast for Ionospheric Parameters of nonlinear prediction；
(3) multidimensional re-construction theory is refer to for the analysis of ionosphere spatial character.
The present invention solves the problems, such as that many system detection data degrees of fusion are bad, is hung down by ionosphere and surveys, tiltedly surveys data, carries
The degree of accuracy of short wave communication frequency predication high.
Brief description of the drawings
Fig. 1 is schematic flow sheet of the invention；
Fig. 2 is absorption loss factors A TnoonFigure；
Fig. 3 is absorption penetration factorFigure；
Fig. 4 schemes for day absorption index p.
Specific embodiment
Further to illustrate the present invention to reach technological means and effect that predetermined purpose is taken, below in conjunction with accompanying drawing
And example, to technical scheme proposed by the present invention, describe in detail as after.
As shown in figure 1, short wave communication reliability estimation method of the present invention based on many system detection datas, including following tool
Body step：
Step A：The reflection point position of short wave communication link and all detecting links is calculated, is chosen and short wave communication link electrical
The close detecting link data of absciss layer reflection point position.Its detailed process includes：
A1, it is theoretical according to short wave communication physical model and propagation rays, calculate short wave communication link ionosphere pip warp
Latitude (θc,λc)。
In formula, d is propagation path distance, a0It is earth radius, θt、λtWarp, latitude for launch point, λrIt is receiving point
Latitude.Formula (1) to (4) is referred to《Radio propagation》(Electronic Industry Press) page 634.
A2, obtains all detection station locations and detection type, calculates each detecting link ionospheric reflection point longitude and latitude, and root
According to the short wave communication link ionosphere pip longitude and latitude that step A1 determines, choose and short wave communication link ionosphere pip position
Put close detecting link.
A21, if detecting link is vertical sounding, acquisition station longitude and latitude is detecting link ionospheric reflection point longitude and latitude
Degree；If detecting link is oblique incidence souding, detecting link ionospheric reflection point calculation of longitude ＆ latitude method is with step A1；
A22, chooses the detecting link close with short wave communication link ionosphere reflection point position.
Preset distance thresholding dlim, numerical value is specified by user, when user does not specify, can based on experience value come true
It is fixed, or retain last time history setting value, accordingly, obtain nearest apart from short wave communication link ionosphere reflection point position
NlinkBar detecting link.NlinkValue is specified by user, when user does not specify, can be determined based on experience value, or retain most
A history setting value afterwards.
A3, the detecting link according to selected by step A2 reads its detecting link data of nearest 5 days, and according to detection chain
Circuit-switched data calculates Ionospheric Parameters.
If detecting link is vertical sounding, direct access detection gained E layers of critical frequency (foE), F1 layers of critical frequency
(foF1), F2 layers of critical frequency (foF2) and F2 layers of distance are the Ionospheric Parameters such as the transmission factor (M (3000) F2) of 3000km.
If detecting link is oblique incidence souding, according to the oblique incidence souding data analyses such as MUF, frequency-time delay obtain foE,
The Ionospheric Parameters such as foF1, foF2 and M (3000) F2.
E layers and F1 layers of foE, foF1 parameter can directly using the relation of basic MUF and critical frequency in propagation rays theory
Derive, its circular is as follows：
In formula, E (d) MUF are E layers of MUF, and F1 (d) MUF are F1 layers of MUF.
In formula (5), MEIt is E layers of conversion factor, is expressed as：
ME=3.94+2.80x-1.70x2-0.60x3+0.96x4 (7)
X is, apart from the factor, to be expressed as：
X=min (d/1150-1,0.74) (8)
In formula (5), (6), d is propagation path distance, MF1It is F1 layers of conversion factor, is expressed as：
MF1=J0-0.01(J0-J100)R12 (9)
In formula, R12It is the monthly average slamp value of sunspot 12, J0And J100Distance when respectively sunspot is 0 and 100
The factor, is expressed as：
J0=0.16+2.64 × 10-3d-0.40×10-6d2 (10)
J100=-0.52+2.69 × 10-3d-0.39×10-6d2 (11)
Formula (5) to (11) refers to P.1240-1 page 3 of ITU recommendations.
To F2 Es-region propagations situations, first with MUF and its correspondence time-delay calculation M (3000) F2, its computational methods is as follows：
In formula, hrFor pip mirror-reflection highly, its computational methods is as follows：
In formula, a0It is earth radius, d is propagation path distance, and τ is the corresponding time delays of MUF, and c is the light velocity.Formula (12),
(13) refer to《Radio propagation》(Electronic Industry Press) page 638.
On the basis of M (3000) F2 is obtained, foF2 is calculated using heuristic：Make critical frequency foF2 in [2,30] by step-length
0.01MHz substitutes into (14) formula and calculates MUF, and corresponding foF2 is inversion result when taking predicted value with measured value minimal error.Its tool
Body computing formula is as follows：
F2 (d) MUF=[1+Cd/C3000(B-1)]foF2+0.5·fH(1-d/dmax) (14)
Cd=0.74-0.591Z-0.424Z2-0.090Z3+0.088Z4+0.181Z5+0.096Z6 (15)
Z=1-2d/dmax (16)
B=M (3000) F2-0.124+ { [M (3000) F2]2-4}·[0.0215+0.005·sin(7.854/x-
1.9635)]
(18)
X=max (foF2/foE, 2) (19)
Wherein, d is propagation path distance, dmaxJumped away from f for F2 Es-region propagations are maximumHFor the magnetic rotation of propagation path midpoint frequently
Rate, C3000It is the C at 3000kmdValue.
Formula (14) to (19) is referred to《Radio propagation》(Electronic Industry Press) page 646.
Step B：The detecting link ionospheric reflection point position determined according to step A and detecting link data, forecast detection
The Ionospheric Parameters of link ionosphere pip, its detailed process includes：
The present invention realizes short-period forecasts in 24 hours of detecting link Ionospheric Parameters using weighted analysis method：
In formula, InT () is the first n days probe values of correspondence t currently to be forecast, t is 0,1,2 ..., and 23, N+1 for most
Big weighted analysis number of days, it is 5 (i.e. N=4) that the present invention takes maximum weighted analysis number of days.
Step C：The short wave communication link determined according to step A and detecting link ionospheric reflection point position, and step B
The Ionospheric Parameters of the detecting link ionospheric reflection point that forecast is obtained, reconstruct obtains short wave communication link ionosphere pip
Ionospheric Parameters.Its detailed process includes：
C1, according to step A determine detecting link ionospheric reflection point position, calculate detecting link ionospheric reflection point with
The ionosphere distance of short wave communication link ionosphere pip.Definition space any two points (θi,λi) and (θj,λj) between electricity
Absciss layer is apart from dijFor：
In formula, SF is scale factor, is used to consider characteristic variable's correlation in ionosphere in latitude and the difference of longitudinal
Not, SF recommended values 2.
C2, according to step B determine detecting link ionospheric reflection point Ionospheric Parameters value I (θ, λ), analytical calculation its
With the correlative Z (θ, λ) of Reference ionosphere model Ionospheric Parameters value difference value：
In formula,It is Reference ionosphere model Ionospheric Parameters value, (θ, λ) is warp, the latitude of detecting link ionospheric reflection point
Degree.
C3, according to step C1 and C2 result of calculation, detecting link ionospheric reflection point is asked for using regional restructuring equation group
With the weight coefficient W of short wave communication link ionosphere pip：
In formula, N is the bar number of detecting link, dijIt is i-th detecting link ionospheric reflection point (θi,λi) visited with j-th strip
Surveyor's chain road ionospheric reflection point (θj,λjThe distance between), di0It is i-th detecting link ionospheric reflection point and short wave communication chain
Road ionospheric reflection point (θ0,λ0The distance between), μ is Lagrange's multiplier.
C4, according to step C3 determine weight coefficient W, calculate short wave communication link ionosphere pip Ionospheric Parameters and
The correlative of Reference ionosphere model value difference, and short wave communication link ionosphere pip Ionospheric Parameters value.
C41, according to step C3 determine weight coefficient, calculate short wave communication link ionosphere pip Ionospheric Parameters and
Correlative Z (the θ of Reference ionosphere model value difference0,λ0)。
C42, according to the correlative that C41 determines, calculates short wave communication link ionosphere pip Ionospheric Parameters value I (θ0,
λ0)。
In formula,It is the Reference ionosphere model Ionospheric Parameters at the pip of short wave communication link ionosphere
Value.
Formula (21) to (25) is referred to《Research institute of China Electronics journal》5th phase page 530 in 2013.
Step D：According to the short wave communication link ionosphere pip Ionospheric Parameters value that step C is obtained, short-pass is predicted
Believe the maximum usable frequency of link.Its detailed process includes：
D1, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, calculates short-pass
Believe E layers of maximum usable frequency E (d) MUF of link.
E (d) MUF=ME·foE(MHz) (26)
In formula：
ME=3.94+2.80y-1.70y2-0.60y3+0.96y4；
Wherein, MEIt is E layers of conversion factor, d is propagation path distance.
Formula (26) refers to P.1240-1 page 3 of ITU recommendations ITU-R.
D2, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, calculates short-pass
Believe F1 layers of maximum usable frequency F1 (d) MUF of link.
F1 (d) MUF=MF1·foF1(MHz) (27)
In formula：
MF1=J0-0.01(J0-J100)R12；
J0=0.16+2.64 × 10-3d-0.40×10-6d2；
J100=-0.52+2.69 × 10-3d-0.39×10-6d2。
Wherein, MF1It is F1 layers of conversion factor, J0And J100When respectively sunspot is 0 and 100 apart from the factor, R12For
The monthly average slamp value of sunspot 12, d is propagation path distance.
Formula (27) refers to P.1240-1 page 3 of ITU recommendations ITU-R.
D3, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, with reference to short-pass
Believe that the propagation path of link, apart from situation, calculates F2 layers of maximum usable frequency F2 (d) MUF of short wave communication link.
D31, judges whether short wave communication link propagation path distance d is maximum more than single-hop and jumps away from dmax, according to whether being more than
dmax, calculate F2 layers of maximum usable frequency F2 (d) MUF by formula (14) and (28) respectively.
D32, when propagation path is apart from d≤dmax, radio wave propagation is single-hop mode, and control point is propagation path midpoint, F2 layers
Basic MUF is calculated by (14) formula.
When propagation path is apart from d>dmax, radio wave propagation is many hop patterns, and control point is away from hair sink d0At/2, F2 layers of base
This MUF is calculated by (28) formula.
F2 (d) MUF=min [F2 (dmax)MUF1,F2(dmax)MUF2] (28)
In formula, d0At least to hop the jump of pattern away from F2 (dmax)MUF1With F2 (dmax)MUF2At respectively two control points
F2 (the d of minimum hop patternmax) MUF, it is calculated by (14) formula.
Formula (28) is referred to《Radio propagation》(Electronic Industry Press) page 646.
Step D4, according to E (d) MUF, F1 (d) MUF and F2 (d) MUF that step D1, D2, D3 analyses are obtained, calculates shortwave
Communication link maximum usable frequency MUF.
MUF=max [E (d) MUF, F1 (d) MUF, F2 (d) MUFRop](MHz) (29)
In formula, RopIt is the ratio of F2 layers of work MUF and basic MUF, as listed in table 1.
Table 1 is referred to《Radio propagation》(Electronic Industry Press) page 632.
The ratio R of the F2 layers of work MUF of table 1 and basic MUFop
Step E：According to the maximum usable frequency that step D is obtained, calculate short wave communication link and entered using maximum usable frequency
Received field strength during row communication, analyzes the reliability of short wave communication link.Its detailed process includes：
E1, judges that the propagation path of short wave communication link, apart from d, if propagation path is less than 7000km apart from d, performs step
E2；If propagation path is more than 9000km apart from d, step E3 is performed；Step E2, E3 and E4 are otherwise performed simultaneously.
E2, if propagation path is less than 7000km apart from d, analytical calculation short wave communication link is carried out using maximum usable frequency
Various communication modes and synthesis received field strength during communication.
Etw=94.25+Pt+Gt-20lgp′-Lm-Lg-Lh-Li (30)
In formula：
Etw- received field strength intermediate value, dB μ V/m；
Pt- radiant power, dBkW；
GtRelative to the transmitter antenna gain (dBi) factor of isotropic antenna, dB on-place launch azimuth and the elevation angle；
P'-day wave ray oblique propagation distance, km, are calculated with formula (31)：
In formula：Δ is given by formula (40).
dnIt is to use dnThe jump length of the n mode hoppings that=d/n is calculated, km；
Formula (30), (31) refer to P.533-10 page 10 of ITU recommendations ITU-R.
LmIt is that, higher than the loss factor of maximum usable frequency (MUF), dB, for E layers with F2 Es-region propagations patterns, distinguishes
Calculated with formula (32) and (33)：
In formula, n is number of days, and t is 0,1,2 ... 23.
LgIt is the ground reflection loss factor, dB, for n mode hoppings, is calculated with formula (34)：
Lg=2 (n-1) (34)
Formula (32) to (34) refers to P.533-10 page 11 of ITU recommendations ITU-R.
LiIt is ionospheric absorption fissipation factor, dB, for n mode hoppings, is calculated with formula (35)：
Wherein：
F (χ)=max [cosp(0.881χ),0.02]；
fV=f × cosi110；
In formula：fLIt is minimum reference frequency, (the electronics gyrofrequency of determination is average at the control point given by table 4 for MHz
Value, the about longitudinal component in the magnetic field of the earth of 100km height)；
K is the quantity (as listed by table 4) at control point；
i110It is the inclination angle at 110km；
R12It is the monthly average slamp value of sunspot 12；
ATnoonIt is the absorption loss factor, dB is given by Fig. 2；
χjIt is the sun drift angle at j-th control point or 102 °, takes smaller；
χjnoonIt is the χ at local noonjValue；
It is absorbed layer penetration factor, is equivalent vertical incidence wave frequency rate, fVWith the function of the ratio between FoE, by Fig. 3
Be given；
P is day absorption index, is given by Fig. 4.
LhIt is auroral region and other loss of signal factors, dB.Table 2 and table 3 are given according to earth center dipole
(78.5 ° of N, 68.2 ° of W) geomagnetic latitude Gn(equator north and south) and local time t can try to achieve LhAnalog value, work as Gn<At 42.5 °, Lh
=0.On the Northern Hemisphere, December was winter to 2 months；March to May and September to November are season in spring and autumn point；June to August is summer.
The month of the Southern Hemisphere, winter and summer exchanges.
The propagation path of table 2 the distance no more than auroral region of 2500km and other loss of signal factors
Auroral region of the propagation path of table 3 distance more than 2500km and other loss of signal factors
Table 2, table 3 are referred to《Radio propagation》(Electronic Industry Press) page 657,658.From each communication mode of calculating
Persistent erection value selects the field intensity of two strong E layer models and three strong F2 layer models
Intermediate value carries out power and is superimposed that received field strength intermediate value must be synthesized, and is calculated with formula (36)：
In formula：
Ets- received field strength intermediate value, dB μ V/m,
EtwThe corresponding field intensity value of-each pattern, dB μ V/m.
Formula (35), (36) refer to P.533-10 page 10 to 12 of ITU recommendations ITU-R.
E3, propagation path distance is more than 9000km, calculates received field strength when short wave communication link is communicated using MUF
Value.
Propagation path distance is divided into minimum n sections, its formate field intensity intermediate value according to the jump anomaly no more than 4000km, is used
Formula (37) is calculated：
In formula：
Etl- formate field intensity intermediate value, dB μ V/m；
E0The free-space field strength of -3MW equivalent isotropically radiated power, dB μ V/m are calculated with formula (38)：
E0=139.6-20lg p ' (38)
P'-day wave ray oblique propagation distance, km, are calculated with formula (39) and (40)：
dn- use dnThe jump length of the n mode hoppings that=d/n is calculated, km；
Gt1Maximum gain of-the transmitting antenna in required azimuth vertical angle in the range of 0 ° -8 °, dB；
Gap- remote because focusing makes the increased gain factor of field intensity, dB be calculated with formula (41), when d is π R0Times
During number, 15dB is taken；
In formula (41), R0It is earth radius.
Ly- sky wave propagation effector, dB, it is proposed that be worth for -3.7dB；
fM- highest reference frequency, MHz；
fL- minimum reference frequency, MHz；
fH- average value of its electronics gyro-frequency is tried to achieve with the control point position of ray path mirror reflection height in table 4,
MHz。
Pt- transmission power, dBW.
The control point position of the ray Lu Jingjing reflection heights of table 4
Formula (37), (38) are referred to《Radio propagation》(Electronic Industry Press) page 658；Table 4 is shown in《Radio wave is passed
Broadcast》(Electronic Industry Press) page 634.
E4, if propagation path distance is between [7000,9000] km, field intensity is carried out on the basis of E2 and E3 is performed and is inserted
Value, obtain the propagation path apart from upper use maximum usable frequency communicated when received field strength.
In this distance range, sky wave median of field strength EtiUsing EtsAnd EtlInterpolation calculating is carried out, is calculated with formula (42)：
Eti=100lg10Xi (42)
In formula：
Eti- sky wave median of field strength, dB μ V/m；
Xs=100.01Ets
Xl=100.01Etl
Ets--- received field strength intermediate value, dB μ V/m；
Etl--- formate field intensity intermediate value, dB μ V/m.
Formula (42) is referred to《Radio propagation》(Electronic Industry Press) page 660.
E5, according to the field intensity that step E2 or E3 or E4 are calculated, analyzes the reliability of short wave communication link.
E51, according to the field intensity that step E2 or E3 or E4 are calculated, calculates receiving power intermediate value.
S=E+G-20logf-107.2 (43)
Wherein, E is received field strength, and G represents receiving antenna gain on the direction of incidence, and f is communication signal frequency.
E52, according to the receiving power intermediate value that E51 is calculated, calculate signal to noise ratio median S/N, signal to noise ratio it is upper very
Value deviation DuThe lower ten score values deviation D of SN and signal to noise ratiolSN。
Wherein, FaA、FaM、FaG is respectively atmospheric noise, man-made noise, galactic noise noise coefficient median, and b is band
It is wide.
Wherein, DuSd、DuShTen score value deviations on ten score value deviations (daily), signal respectively on signal (per when), DlA、
DlM、DlG is respectively ten score value deviations under atmospheric noise, man-made noise, galactic noise.
Wherein, DlSd、DlShTen score value deviations under ten score value deviations (daily), signal respectively under signal (per when), DuA、
DuM、DuG is respectively atmospheric noise, man-made noise, ten score value deviations on galactic noise.
E53, lower the ten of the signal to noise ratio median being calculated according to E52, the upper ten score values deviation of signal to noise ratio and signal to noise ratio
Score value deviation, calculates the reliability of short wave communication link.
Wherein, S/NrFor the signal to noise ratio that user requires.
Formula (44) is referred to《Radio propagation》(Electronic Industry Press) page 663；Formula (45) is detailed to formula (47)
See《Radio propagation》666-669 pages of (Electronic Industry Press).
Above example is provided just for the sake of the description purpose of the present invention, and is not intended to limit the scope of the present invention.This
The scope of invention is defined by the following claims.The various equivalents that do not depart from spirit and principles of the present invention and make and repair
Change, all should cover within the scope of the present invention.
Claims (7)
1. a kind of short wave communication reliability estimation method based on many system detection datas, it is characterised in that comprise the following steps：
Step A：The reflection point position of short wave communication link and all detecting links is calculated, is chosen and short wave communication link ionosphere
The close detecting link data of reflection point position；
Step B：The detecting link ionospheric reflection point position determined according to step A and detecting link data, forecast are detected
The Ionospheric Parameters of link ionosphere pip；
Step C：The short wave communication link determined according to step A and detecting link ionospheric reflection point position, and step B forecast
The Ionospheric Parameters of the detecting link ionospheric reflection point for obtaining, reconstruct obtains the ionization of short wave communication link ionosphere pip
Layer parameter；
Step D：According to the short wave communication link ionosphere pip Ionospheric Parameters value that step C is obtained, short wave communication chain is predicted
The maximum usable frequency on road；
Step E：According to the maximum usable frequency that step D is obtained, calculate short wave communication link and led to using maximum usable frequency
Received field strength during letter, analyzes the reliability of short wave communication link.
2. the short wave communication reliability estimation method based on many system detection datas according to claim 1, its feature exists
In：The step A is comprised the following steps：
Step A1, it is theoretical according to short wave communication physical model and propagation rays, calculate short wave communication link ionosphere pip warp
Latitude；
Step A2, obtains all detection station locations and detection type, calculates each detecting link ionospheric reflection point longitude and latitude, and root
According to the short wave communication link ionosphere pip longitude and latitude that step A1 determines, choose and short wave communication link ionosphere pip position
Put close detecting link；
Step A3, the detecting link according to selected by step A2 reads its detecting link data of nearest 5 days, and according to detection chain
Circuit-switched data calculates Ionospheric Parameters.
3. the short wave communication reliability estimation method based on many system detection datas according to claim 1, its feature exists
In：The step C is comprised the following steps：
Step C1, according to step A determine detecting link ionospheric reflection point position, calculate detecting link ionospheric reflection point with
The ionosphere distance of short wave communication link ionosphere pip；
Step C2, the Ionospheric Parameters value of the detecting link ionospheric reflection point determined according to step B, analytical calculation its with reference
The correlative of ionospheric model Ionospheric Parameters value difference value；
Step C3, according to step C1 and C2 result of calculation, detecting link ionospheric reflection point is asked for using regional restructuring equation group
With the weight coefficient W of short wave communication link ionosphere pip；
Step C4, according to step C3 determine weight coefficient W, calculate short wave communication link ionosphere pip Ionospheric Parameters and
The correlative of Reference ionosphere model value difference, and short wave communication link ionosphere pip Ionospheric Parameters value.
4. the short wave communication reliability estimation method based on many system detection datas according to claim 1, its feature exists
In：The step D is comprised the following steps：
Step D1, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, calculates short-pass
Believe E layers of maximum usable frequency E (d) MUF of link；
Step D2, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, calculates short-pass
Believe F1 layers of maximum usable frequency F1 (d) MUF of link；
Step D3, according to the short wave communication link ionosphere pip Ionospheric Parameters value that step C analyses are obtained, with reference to short-pass
Believe that the propagation path of link, apart from situation, calculates F2 layers of maximum usable frequency F2 (d) MUF of short wave communication link；
Step D4, according to E (d) MUF, F1 (d) MUF and F2 (d) MUF that step D1, D2, D3 analyses are obtained, calculates short wave communication
Link maximum usable frequency MUF.
5. the short wave communication reliability estimation method based on many system detection datas according to claim 4, its feature exists
In：The step D3 is specific as follows：
Step D31, judges whether the propagation path of short wave communication link is jumped away from d apart from d more than single-hop maximummax；
Step D32, if propagation path is jumped away from d apart from d less than or equal to single-hop maximummax, with reference to the ionosphere that step C analyses are obtained
Parameter value, F2 layers of basic MUF is calculated according to single-hop mode；If propagation path is apart from d more than the maximum jump of single-hop away from dmax, according to electricity
It is F2 layers of basic MUF of multi-hop mode computation that ripple is propagated.
6. the short wave communication reliability estimation method based on many system detection datas according to claim 1, its feature exists
In：The step E is specific as follows：
Step E1, judges the propagation path distance of short wave communication link, if propagation path distance is less than 7000km, performs step
E2；If propagation path distance is more than 9000km, step E3 is performed；Step E2, E3 and E4 are otherwise performed simultaneously；
Step E2, if propagation path distance is less than 7000km, analytical calculation short wave communication link is carried out using maximum usable frequency
Various communication modes and synthesis received field strength during communication；
Step E3, if propagation path distance is more than 9000km, analytical calculation short wave communication link is carried out using maximum usable frequency
Received field strength during communication；
Step E4, if propagation path distance is between [7000,9000] km, field intensity is carried out on the basis of E2 and E3 is performed and is inserted
Value, obtain the propagation path apart from upper use maximum usable frequency communicated when received field strength；
Step E5, according to the field intensity that step E2 or E3 or E4 are calculated, analyzes the reliability of short wave communication link.
7. the short wave communication reliability estimation method based on many system detection datas according to claim 6, its feature exists
In：The step E5 is specific as follows：
Step E51, according to the field intensity that step E2 or E3 or E4 are calculated, calculates receiving power intermediate value；
Step E52, according to the receiving power intermediate value that E51 is calculated, calculate signal to noise ratio median, signal to noise ratio upper ten score value it is inclined
The lower ten score values deviation of difference and signal to noise ratio；
Step E53, lower the ten of the signal to noise ratio median being calculated according to E52, the upper ten score values deviation of signal to noise ratio and signal to noise ratio
Score value deviation, calculates the reliability of short wave communication link.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611071301.9A CN106788815B (en) | 2016-11-28 | 2016-11-28 | Short wave communication reliability assessment method based on multi-system detection data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611071301.9A CN106788815B (en) | 2016-11-28 | 2016-11-28 | Short wave communication reliability assessment method based on multi-system detection data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106788815A true CN106788815A (en) | 2017-05-31 |
| CN106788815B CN106788815B (en) | 2020-09-22 |
Family
ID=58904254
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611071301.9A Active CN106788815B (en) | 2016-11-28 | 2016-11-28 | Short wave communication reliability assessment method based on multi-system detection data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106788815B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107911185A (en) * | 2017-11-03 | 2018-04-13 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | A kind of maximum usable frequency computational methods suitable for short-wave link during ionospheric storm |
| CN108462545A (en) * | 2018-01-29 | 2018-08-28 | 武汉小石科技有限公司 | A kind of ionosphere foF based on single receiving station2Parameter reconstructing method |
| CN109490641A (en) * | 2019-01-05 | 2019-03-19 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | A kind of calculation method of middle latitude area Sporadic E layer short wave field strength |
| WO2019070857A2 (en) | 2017-10-04 | 2019-04-11 | Skywave Networks Llc | Adjusting transmissions based on direct sensing of the ionosphere |
| CN110245316A (en) * | 2019-05-10 | 2019-09-17 | 中国人民解放军31007部队 | A kind of inversion method of Ionospheric Parameters |
| CN110288117A (en) * | 2019-05-10 | 2019-09-27 | 中国人民解放军31007部队 | A kind of regional restructuring method of Ionospheric Parameters critical frequency |
| CN111148145A (en) * | 2019-12-20 | 2020-05-12 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Short wave sky wave channel maximum accessible bandwidth evaluation method |
| CN111881584A (en) * | 2020-07-30 | 2020-11-03 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Random generation method of short-wave link library |
| CN112036011A (en) * | 2020-08-05 | 2020-12-04 | 中国人民解放军海军工程大学 | Very low frequency wave communication transmission analysis method and system for underwater vehicle |
| CN112272067A (en) * | 2020-10-15 | 2021-01-26 | 天津大学 | Short wave broadcast frequency scheduling method based on multi-source data processing |
| CN113109632A (en) * | 2021-04-08 | 2021-07-13 | 三门峡职业技术学院 | Method for inverting F2 layer parameters by using oblique ionogram |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101719802A (en) * | 2009-10-30 | 2010-06-02 | 北京航空航天大学 | Device and calculation method for predicting maximum usable frequency (MUF) of short-wave communication |
-
2016
- 2016-11-28 CN CN201611071301.9A patent/CN106788815B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101719802A (en) * | 2009-10-30 | 2010-06-02 | 北京航空航天大学 | Device and calculation method for predicting maximum usable frequency (MUF) of short-wave communication |
Non-Patent Citations (5)
| Title |
|---|
| 付炜，王健等: "高频通信可用频率现报与预报方法", 《中国电子科学研究院学报》 * |
| 刘瑞源，张北辰等: "中国地区电离层f oF2重构方法及其在短期预报中的应用", 《地球物理学报》 * |
| 国家电信联盟: "Method for the prediction of the", 《RECOMMENDATION ITU-R P.533-10》 * |
| 国际电信联盟: "ITU-R methods of basic MUF, operational MUF and ray-path prediction", 《REC.ITU-R P.1240-1》 * |
| 熊皓等: "《无线电波传播》", 29 February 2000 * |
Cited By (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020537404A (en) * | 2017-10-04 | 2020-12-17 | スカイウェイブ・ネットワークス・エルエルシー | Adjustment of transmission based on direct detection of the ionosphere |
| GB2581699B (en) * | 2017-10-04 | 2022-08-31 | Skywave Networks Llc | Adjusting transmissions based on direct sensing of the ionosphere |
| WO2019070857A2 (en) | 2017-10-04 | 2019-04-11 | Skywave Networks Llc | Adjusting transmissions based on direct sensing of the ionosphere |
| JP7142687B2 (en) | 2017-10-04 | 2022-09-27 | スカイウェイブ・ネットワークス・エルエルシー | Coordination of transmissions based on direct sensing of the ionosphere |
| US11496210B2 (en) | 2017-10-04 | 2022-11-08 | Skywave Networks Llc | Adjusting transmissions based on direct sensing of the ionosphere |
| EP3692652A4 (en) * | 2017-10-04 | 2021-06-30 | Skywave Networks LLC | Adjusting transmissions based on direct sensing of the ionosphere |
| CN107911185B (en) * | 2017-11-03 | 2020-12-04 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Method for calculating highest available frequency of short wave link in ionosphere storm period |
| CN107911185A (en) * | 2017-11-03 | 2018-04-13 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | A kind of maximum usable frequency computational methods suitable for short-wave link during ionospheric storm |
| CN108462545A (en) * | 2018-01-29 | 2018-08-28 | 武汉小石科技有限公司 | A kind of ionosphere foF based on single receiving station2Parameter reconstructing method |
| CN108462545B (en) * | 2018-01-29 | 2021-12-28 | 武汉小石科技有限公司 | Ionosphere foF based on single receiving station2Parameter reconstruction method |
| CN109490641B (en) * | 2019-01-05 | 2020-12-08 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Method for calculating field intensity of short wave of sporadic E layer in mid-latitude area |
| CN109490641A (en) * | 2019-01-05 | 2019-03-19 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | A kind of calculation method of middle latitude area Sporadic E layer short wave field strength |
| CN110288117B (en) * | 2019-05-10 | 2021-06-29 | 中国人民解放军31007部队 | Regional reconstruction method for critical frequency of ionosphere parameters |
| CN110288117A (en) * | 2019-05-10 | 2019-09-27 | 中国人民解放军31007部队 | A kind of regional restructuring method of Ionospheric Parameters critical frequency |
| CN110245316A (en) * | 2019-05-10 | 2019-09-17 | 中国人民解放军31007部队 | A kind of inversion method of Ionospheric Parameters |
| CN111148145A (en) * | 2019-12-20 | 2020-05-12 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Short wave sky wave channel maximum accessible bandwidth evaluation method |
| CN111148145B (en) * | 2019-12-20 | 2022-03-22 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Short wave sky wave channel maximum accessible bandwidth evaluation method |
| CN111881584A (en) * | 2020-07-30 | 2020-11-03 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Random generation method of short-wave link library |
| CN112036011B (en) * | 2020-08-05 | 2022-09-20 | 中国人民解放军海军工程大学 | Very low frequency wave communication transmission analysis method and system for underwater vehicle |
| CN112036011A (en) * | 2020-08-05 | 2020-12-04 | 中国人民解放军海军工程大学 | Very low frequency wave communication transmission analysis method and system for underwater vehicle |
| CN112272067A (en) * | 2020-10-15 | 2021-01-26 | 天津大学 | Short wave broadcast frequency scheduling method based on multi-source data processing |
| CN112272067B (en) * | 2020-10-15 | 2022-04-08 | 天津大学 | Short wave broadcast frequency scheduling method based on multi-source data processing |
| CN113109632A (en) * | 2021-04-08 | 2021-07-13 | 三门峡职业技术学院 | Method for inverting F2 layer parameters by using oblique ionogram |
| CN113109632B (en) * | 2021-04-08 | 2023-04-14 | 三门峡职业技术学院 | Method for inverting F2 layer parameters by using oblique ionogram |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106788815B (en) | 2020-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106788815A (en) | A kind of short wave communication reliability estimation method based on many system detection datas | |
| EP1312234B1 (en) | Systems and methods for determining signal coverage | |
| Said et al. | Long‐range lightning geolocation using a VLF radio atmospheric waveform bank | |
| JP5123897B2 (en) | Indoor / outdoor determination apparatus and indoor / outdoor determination method | |
| Chen et al. | A statistical method based on remote sensing for the estimation of air temperature in China. | |
| US20230206545A1 (en) | Iterative ray-tracing for autoscaling of oblique ionograms | |
| CN107733515B (en) | Satellite communication link analysis method under in-orbit complex environment | |
| US20060199546A1 (en) | Electro-magnetic propagation modeling | |
| CN106793082A (en) | A kind of positioning of mobile equipment method in WLAN/ bluetooth heterogeneous network environments | |
| CN103698750A (en) | HY-2 satellite scatterometer sea surface wind field retrieval method and device | |
| US8909242B2 (en) | Cooperative spectrum sensing method and system for locationing primary transmitters in a cognitive radio system | |
| CN109946235A (en) | The multi layer cloud inversion method of wind and cloud 4A meteorological satellite Multichannel Scan Imagery Radiometer | |
| CN105182322A (en) | Passive positioning method based on reflected signal phase difference | |
| Wang et al. | Regional refined long-term predictions method of usable frequency for HF communication based on machine learning over Asia | |
| CN108462545A (en) | A kind of ionosphere foF based on single receiving station2Parameter reconstructing method | |
| Béniguel et al. | Ionospheric scintillation monitoring and modelling | |
| Popleteev | Improving ambient FM indoor localization using multipath-induced amplitude modulation effect: A year-long experiment | |
| Türkyılmaz | Environment aware location estimation in cellular networks | |
| Lwas et al. | Effects of wind velocity on slant path rain-attenuation for satellite application in Malaysia | |
| WO2007064953A2 (en) | Hf radio network frequency management | |
| Abdulrahman et al. | An improved slant path attenuation prediction method in tropical climates | |
| Valtr et al. | A terrestrial multiple-receiver radio link experiment at 10.7 GHz–comparisons of results with parabolic equation calculations | |
| Basere et al. | Cell coverage area estimation from receive signal level (RSL) measurements | |
| Matolak et al. | Ag channel measurement and modeling results for over-water and hilly terrain conditions | |
| WenHui et al. | The Test of Wireless Data Transmission System in Sea-Area Seismic |
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 |