CN106788815B - Short wave communication reliability assessment method based on multi-system detection data - Google Patents

Short wave communication reliability assessment method based on multi-system detection data Download PDF

Info

Publication number
CN106788815B
CN106788815B CN201611071301.9A CN201611071301A CN106788815B CN 106788815 B CN106788815 B CN 106788815B CN 201611071301 A CN201611071301 A CN 201611071301A CN 106788815 B CN106788815 B CN 106788815B
Authority
CN
China
Prior art keywords
short
wave communication
communication link
ionospheric
detection
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.)
Active
Application number
CN201611071301.9A
Other languages
Chinese (zh)
Other versions
CN106788815A (en
Inventor
张秀强
苏东林
谢树果
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201611071301.9A priority Critical patent/CN106788815B/en
Publication of CN106788815A publication Critical patent/CN106788815A/en
Application granted granted Critical
Publication of CN106788815B publication Critical patent/CN106788815B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 invention relates to a short wave communication reliability assessment method based on multi-system detection data, which comprises the steps of calculating the positions of reflection points of a short wave communication link and all detection links, and selecting detection link data close to the position of an ionosphere reflection point of the short wave communication link; forecasting to obtain ionospheric parameters of ionospheric reflection points of the detection link according to the determined ionospheric reflection point positions of the detection link and the detection link data; determining positions of ionospheric reflection points of the short-wave communication link and the detection link, predicting ionospheric parameters of the ionospheric reflection points of the detection link, and reconstructing to obtain ionospheric parameters of the ionospheric reflection points of the short-wave communication link; predicting the highest available frequency of the short-wave communication link by using a short-wave propagation prediction method according to the obtained ionized layer parameters; and calculating the receiving field intensity of the short-wave communication link when the short-wave communication link adopts the highest available frequency for communication according to the obtained highest available frequency, and analyzing the reliability of the short-wave communication link.

Description

Short wave communication reliability assessment method based on multi-system detection data
Technical Field
The invention relates to the field of short-wave communication, in particular to a short-wave communication reliability assessment method based on multi-system detection data.
Background
The ionosphere, as an important component of a near-earth space environment, has important influence (such as signal fading, error code and the like) which cannot be ignored on short-wave communication, and the accurate determination of the state characteristics of the ionosphere can effectively avoid or reduce the influence of ionosphere change on a short-wave communication link, thereby ensuring the reliable communication of the short-wave link. At present, the research of the ionosphere state characteristic prediction technology mainly focuses on predicting ionosphere parameters in a time domain and a space domain, and an engineering practical method for evaluating the communication quality of a short-wave communication link by combining measured data is lacked.
Disclosure of Invention
The invention solves the problems: in order to reasonably and effectively utilize short-wave detection data under different detection mechanisms, improve the utilization rate of the detection data and provide support for short-term frequency planning of a short-wave communication link, the invention provides a short-wave communication link characteristic analysis method based on vertical and oblique detection data by utilizing theories of nonlinear prediction, multidimensional space reconstruction and the like, so that the influence of ionosphere change on the short-wave communication link is effectively avoided or reduced, and stable and efficient communication of the short-wave communication link is ensured.
The technical scheme of the invention is as follows: a short wave communication reliability assessment method based on multi-system detection data comprises the following steps:
step A: calculating the positions of reflection points of the short-wave communication link and all the detection links, and selecting detection link data close to the positions of ionospheric reflection points of the short-wave communication link;
and B: b, forecasting and obtaining ionospheric parameters of the ionospheric reflection points of the detection link according to the ionospheric reflection point positions of the detection link and the detection link data determined in the step A;
and C: reconstructing to obtain ionospheric parameters of the ionospheric reflection points of the short-wave communication link according to the ionospheric reflection point positions of the short-wave communication link and the detection link determined in the step A and the ionospheric parameters of the ionospheric reflection points of the detection link obtained by prediction in the step B;
step D: c, predicting the highest available frequency of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link obtained in the step C;
step E: and D, calculating the receiving field intensity of the short-wave communication link when the short-wave communication link adopts the highest available frequency for communication according to the highest available frequency obtained in the step D, and analyzing the reliability of the short-wave communication link.
The step A comprises the following steps:
a1, calculating the longitude and latitude of the ionized layer reflection point of the short-wave communication link according to the short-wave communication physical model and the propagation ray theory;
step A2, acquiring the positions and detection types of all detection stations, calculating the longitude and latitude of the ionospheric reflection point of each detection link, and selecting a detection link with the position close to the position of the ionospheric reflection point of the short-wave communication link according to the longitude and latitude of the ionospheric reflection point of the short-wave communication link determined in the step A1;
step A3, according to the probe link selected in step A2, reading the data of the probe link of the last 5 days, and calculating the ionospheric parameters according to the data of the probe link.
The step C comprises the following steps:
step C1, calculating the ionospheric distance between the ionospheric reflection point of the detection link and the ionospheric reflection point of the short-wave communication link according to the position of the ionospheric reflection point of the detection link determined in the step A;
step C2, according to the ionosphere parameter value of the link ionosphere reflection point determined in the step B, analyzing and calculating the correlation quantity of the difference value between the ionosphere parameter value and the reference ionosphere model ionosphere parameter value;
step C3, according to the calculation results of the step C1 and the step C2, a weight coefficient W of the ionosphere reflection point of the detection link and the ionosphere reflection point of the short-wave communication link is obtained by utilizing a regional reconstruction equation set;
and C4, calculating the correlation quantity of the difference value between the ionospheric reflection point ionospheric parameter of the short-wave communication link and the reference ionospheric model value and the ionospheric reflection point ionospheric parameter value of the short-wave communication link according to the weight coefficient W determined in the step C3.
The step D comprises the following steps:
d1, calculating the highest available frequency E (D) MUF of the E layer of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link analyzed and obtained in the step C;
d2, calculating the F1 layer highest available frequency F1(D) MUF of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link obtained by analyzing in the step C;
d3, calculating the F2 layer highest available frequency F2(D) MUF of the short-wave communication link according to the ionized layer reflection point parameter value of the short-wave communication link ionized layer obtained by analyzing in the step C and by combining the propagation distance condition of the short-wave communication link;
and D4, calculating the highest available frequency MUF of the short-wave communication link according to the E (D) MUF, the F1(D) MUF and the F2(D) MUF analyzed and obtained in the steps D1, D2 and D3.
The step D3 is specifically as follows:
step D31, judging whether the propagation path distance D of the short wave communication link is larger than the maximum jump distance D of single jumpmax；
Step D32, if the propagation path distance D is less than or equal to the maximum hop distance D of single hopmaxCombining the ionosphere parameter values obtained by analysis in the step C, and calculating F2 layer basic MUF according to a single-hop mode; if the propagation path distance d is greater than the maximum hop distance d of a single hopmaxThe F2-layer basic MUF is calculated in accordance with the propagation of the radio wave as a multi-hop pattern.
The step E is specifically as follows:
e1, judging the propagation path distance of the short wave communication link, and if the propagation path distance is less than 7000km, executing the step E2; if the propagation path distance is greater than 9000km, executing a step E3; otherwise, performing steps E2, E3 and E4 simultaneously;
e2, if the distance of the propagation path is less than 7000km, analyzing and calculating various propagation modes and synthesized receiving field intensity when the shortwave communication link adopts the highest available frequency for communication;
step E3, if the distance of the propagation path is larger than 9000km, analyzing and calculating the receiving field intensity when the shortwave communication link adopts the highest available frequency for communication;
step E4, if the propagation path distance is between [7000,9000] km, field intensity interpolation is carried out on the basis of E2 and E3, and the receiving field intensity when the highest available frequency is adopted for communication on the propagation path distance is obtained;
and E5, analyzing the reliability of the short-wave communication link according to the field intensity calculated in the step E2, E3 or E4.
The step E5 is specifically as follows:
step E51, calculating the median value of the received power according to the field intensity calculated in the step E2, E3 or E4;
step E52, calculating the middle value of the signal-to-noise ratio, the upper ten value deviation of the signal-to-noise ratio and the lower ten value deviation of the signal-to-noise ratio according to the received power median value calculated by E51;
and E53, calculating the reliability of the short-wave communication link according to the signal-to-noise ratio intermediate value, the upper ten-point value deviation of the signal-to-noise ratio and the lower ten-point value deviation of the signal-to-noise ratio which are obtained by calculating in E52.
Compared with the prior art, the invention has the advantages that: the method has the greatest advantages that the established short-wave link reliability forecasting method can provide short-wave communication link quality assessment short-term forecasting service for a short-wave communication system, so that a user can know the receiving effect of short-wave communication at a future moment through an ionized layer in advance, a basis is provided for selecting other communication modes or frequency conversion when the short-wave reliability is not high, and economic loss caused by the fact that short-wave communication cannot be carried out can be effectively avoided. Meanwhile, under the condition that the short-wave communication link is initially connected or communication is suddenly interrupted, accurate working frequency parameters can be quickly obtained, communication can be quickly established or recovered, normal and efficient operation of the short-wave communication network is guaranteed, and the influence of ionosphere change on the optimal frequency of the short-wave channel is reduced.
(1) Two types of data of vertical and oblique detection are fused for short-wave communication link characteristic analysis;
(2) a nonlinear prediction theory is introduced for the short-term prediction of ionosphere parameters;
(3) multidimensional reconstruction theory is introduced for analysis of ionospheric spatial properties.
The invention solves the problem of poor fusion degree of multi-system detection data, and improves the accuracy of short-wave communication frequency prediction through ionosphere vertical measurement and oblique measurement data.
Drawings
FIG. 1 is a schematic flow diagram of the present invention;
FIG. 2 shows the absorption loss factor ATnoonA drawing;
FIG. 3 shows absorption penetration factor
Figure BDA0001163779430000041
A drawing;
FIG. 4 is a graph of daily absorption index p.
Detailed Description
To further explain the technical means and effects of the present invention adopted to achieve the intended purpose, the following detailed description of the technical solution proposed by the present invention is made with reference to the accompanying drawings and examples.
As shown in fig. 1, the short-wave communication reliability evaluation method based on multi-system probe data of the present invention includes the following specific steps:
step A: and calculating the positions of reflection points of the short-wave communication link and all the detection links, and selecting detection link data close to the positions of ionospheric reflection points of the short-wave communication link. The specific process comprises the following steps:
a1, calculating the longitude and latitude (theta) of the ionosphere reflection point of the short-wave communication link according to the short-wave communication physical model and the propagation ray theoryc,λc)。
Figure BDA0001163779430000051
Figure BDA0001163779430000052
Figure BDA0001163779430000053
Figure BDA0001163779430000054
Wherein d is the propagation path distance, a0Is the radius of the earth, thetat、λtLongitude and latitude of the emission point, lambdarIs the latitude of the receiving point. The equations (1) to (4) are described in detail in "radio wave propagation" (electronic industry Press) page 634.
A2, acquiring the positions and detection types of all the detection stations, calculating the longitude and latitude of the ionospheric reflection point of each detection link, and selecting a detection link with the position close to the position of the ionospheric reflection point of the short-wave communication link according to the longitude and latitude of the ionospheric reflection point of the short-wave communication link determined in the step A1.
A21, if the detection link is vertical detection, the longitude and latitude of the detection station are the longitude and latitude of the ionosphere reflection point of the detection link; if the detection link is oblique detection, the method for calculating the longitude and latitude of the ionosphere reflection point of the detection link is the same as the step A1;
and A22, selecting a detection link with the position close to the ionosphere reflection point of the short-wave communication link.
Presetting a distance threshold dlimThe numerical value is specified by a user, when the user does not specify the numerical value, the numerical value can be determined according to an empirical value, or a last historical set value is reserved, and accordingly N closest to the position of an ionospheric reflection point of the short wave communication link is obtainedlinkThe bars probe the links. N is a radical oflinkThe value is specified by the user, and when not specified by the user, can be determined based on empirical values, or the last historical setting is retained.
A3, according to the probe link selected in step A2, reading the data of the probe link of the last 5 days, and calculating the ionosphere parameters according to the data of the probe link.
If the probe link is vertical probe, the electric layer parameters such as E layer critical frequency (foE), F1 layer critical frequency (foF1), F2 layer critical frequency (foF2) and transmission factor (M (3000) F2) with distance of F2 layer of 3000km are directly obtained.
If the detection link is oblique detection, the ionosphere parameters such as foE, foF1, foF2 and M (3000) F2 are obtained according to the analysis of oblique detection data such as MUF, frequency-time delay and the like.
The parameters of the fo E and the fo F1 of the E layer and the F1 layer can be directly derived by utilizing the relation between basic MUF and critical frequency in the propagation ray theory, and the specific calculation method is as follows:
Figure BDA0001163779430000061
Figure BDA0001163779430000062
wherein E (d) MUF is E-layer MUF, and F1(d) MUF is F1-layer MUF.
In the formula (5), MEIs an E-layer conversion factor, expressed as:
ME＝3.94+2.80x-1.70x2-0.60x3+0.96x4(7)
x is a distance factor, expressed as:
x＝min(d/1150-1,0.74) (8)
in the formulas (5) and (6), d is the propagation path distance, MF1Is a conversion factor of F1 layer, expressed as:
MF1＝J0-0.01(J0-J100)R12(9)
in the formula, R12Mean flow value of Sun-Zi 12 months, J0And J100Distance factors for 0 and 100 sun black, respectively, are 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)
equations (5) through (11) are detailed on page 3 of ITU recommendation P.1240-1.
For the propagation condition of the F2 layer, M (3000) F2 is calculated by using MUF and corresponding time delay, and the calculation method is as follows:
Figure BDA0001163779430000071
in the formula, hrThe calculation method of the specular reflection height of the reflection point is as follows:
Figure BDA0001163779430000072
in the formula, a0The radius of the earth, d the propagation path distance, τ the time delay corresponding to MUF, and c the speed of light. The formulas (12) and (13) are described in detail in pages 638 of radio wave propagation (electronic industry press).
On the basis of the obtained M (3000) F2, a heuristic is used to calculate foF 2: the critical frequency foF2 is substituted into the formula (14) at [2,30] according to the step length of 0.01MHz to calculate MUF, and foF2 corresponding to the minimum error between the predicted value and the measured value is taken as the inversion result. The specific calculation 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)
Figure BDA0001163779430000073
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)
where d is the propagation path distance, dmaxMaximum hop distance, F, for layer F2 propagationHAt the midpoint of the propagation path, C3000Is C at 3000kmdThe value is obtained.
The formulas (14) to (19) are described in detail in "radio wave propagation" (electronic industry Press) p 646.
And B: forecasting ionospheric parameters of the link ionospheric reflection points according to the link ionospheric reflection point positions and link data determined in step a, wherein the method specifically comprises the following steps:
the invention adopts a weighting analysis method to realize 24-hour short-term prediction of the link ionosphere parameters:
Figure BDA0001163779430000081
in the formula InAnd (t) is a detection value corresponding to the t moment of the previous N days to be currently forecasted, t is 0,1,2, …, 23, N +1 is the maximum weighted analysis day number, and the maximum weighted analysis day number is 5 (namely N is 4).
And C: and B, reconstructing to obtain the ionospheric parameters of the ionospheric reflection points of the short-wave communication link according to the ionospheric reflection point positions of the short-wave communication link and the detection link determined in the step A and the ionospheric parameters of the ionospheric reflection points of the detection link obtained by prediction in the step B. The specific process comprises the following steps:
c1, according to the position of the ionospheric reflection point of the detection link determined in the step A, calculating the ionospheric distance between the ionospheric reflection point of the detection link and the ionospheric reflection point of the short-wave communication link. Define arbitrary two points (theta) in spacei,λi) And (theta)j,λj) Ionospheric distance d betweenijComprises the following steps:
Figure BDA0001163779430000082
where SF is a scale factor to account for differences in the ionospheric characteristic variable dependencies in the latitude and longitude directions, SF suggests a value of 2.
C2, according to the ionosphere parameter value I (theta, lambda) of the ionosphere reflection point of the detection link determined in the step B, analyzing and calculating the correlation quantity Z (theta, lambda) of the ionosphere parameter value difference with the reference ionosphere model:
Figure BDA0001163779430000083
in the formula (I), the compound is shown in the specification,
Figure BDA0001163779430000084
for reference ionosphere model ionosphere parameter values, (θ, λ) are longitude and latitude of the ionosphere reflection point of the probe link.
C3, according to the calculation results of the steps C1 and C2, calculating the weight coefficient W of the ionospheric reflection point of the detection link and the ionospheric reflection point of the short-wave communication link by using a regional reconstruction equation set:
Figure BDA0001163779430000091
where N is the number of probing links, dijIonospheric reflection points (theta) for the ith probe linki,λi) Ionospheric reflection point (theta) of j-th probe linkj,λj) In betweenDistance, di0For the ionospheric reflection point of the ith detection link and the ionospheric reflection point (theta) of the short-wave communication link0,λ0) And μ is the lagrange multiplier.
And C4, calculating the correlation quantity of the difference value between the ionospheric reflection point ionospheric parameter of the short-wave communication link and the reference ionospheric model value and the ionospheric reflection point ionospheric parameter value of the short-wave communication link according to the weight coefficient W determined in the step C3.
C41, calculating the related quantity Z (theta) of the difference value between the ionospheric reflection point ionospheric parameter of the short-wave communication link ionospheric reflection point and the reference ionospheric model value according to the weight coefficient determined in the step C30,λ0)。
Figure BDA0001163779430000092
C42, calculating ionospheric reflection point ionospheric parameter value I (theta) of short wave communication link according to the correlation quantity determined by C410,λ0)。
Figure BDA0001163779430000093
In the formula (I), the compound is shown in the specification,
Figure BDA0001163779430000094
and the ionospheric parameter values of the reference ionospheric model at the ionospheric reflection points of the short-wave communication link.
The formulas (21) to (25) are detailed in the journal of the institute of Chinese electronic sciences, 2013, 5 th stage 530.
Step D: and D, predicting the highest available frequency of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link obtained in the step C. The specific process comprises the following steps:
d1, calculating the E layer highest available frequency E (D) MUF of the short-wave communication link according to the ionized layer parameter value of the short-wave communication link ionized layer reflection point obtained by analyzing in the step C.
E(d)MUF＝ME·foE(MHz) (26)
In the formula:
ME＝3.94+2.80y-1.70y2-0.60y3+0.96y4；
Figure BDA0001163779430000101
wherein M isED is the propagation path distance, which is the E-layer conversion factor.
Equation (26) is detailed in ITU recommendation ITU-R P.1240-1, page 3.
D2, calculating the F1 layer highest available frequency F1(D) MUF of the short-wave communication link according to the ionized layer reflection point ionized layer parameter value of the short-wave communication link obtained by analyzing in the step C.
F1(d)MUF＝MF1·foF1(MHz) (27)
In the 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 M isF1Is a conversion factor of F1 layer, J0And J100Distance factors R for 0 and 100 sun black, respectively12The average flow value of Sun-Black 12 months, d is the propagation path distance.
The formula (27) is detailed on page 3 of ITU recommendation ITU-R P.1240-1.
D3, calculating the F2 layer highest available frequency F2(D) MUF of the short-wave communication link according to the ionized layer reflection point parameter value of the short-wave communication link obtained by analyzing in the step C and by combining the propagation path distance condition of the short-wave communication link.
D31, judging whether the short wave communication link propagation path distance D is larger than the single-hop maximum hop distance DmaxAccording to whether d is greater than dmaxThe highest available frequency F2(d) MUF at the F2 level is calculated from equations (14) and (28), respectively.
D32, when the propagation path distance D is less than or equal to DmaxElectric powerThe wave propagation is in a single-hop mode, the control point is the middle point of the propagation path, and the F2 layer basic MUF is calculated by the formula (14).
When propagation path distance d>dmaxThe radio wave propagation is a multi-hop mode, and the control point is a distance transmitting and receiving point d0At/2, the F2 layer basic MUF is calculated from the formula (28).
F2(d)MUF＝min[F2(dmax)MUF1,F2(dmax)MUF2](28)
In the formula (d)0Hop length for the least-hop mode, F2(d)max)MUF1And F2(d)max)MUF2F2(d) being the lowest hop patterns at the two control points, respectivelymax) MUF, calculated from the formula (14).
The formula (28) is described in detail in "radio wave propagation" (electronic industry Press) p 646.
And D4, calculating the highest available frequency MUF of the short-wave communication link according to the E (D) MUF, the F1(D) MUF and the F2(D) MUF analyzed and obtained in the steps D1, D2 and D3.
MUF＝max[E(d)MUF,F1(d)MUF,F2(d)MUF·Rop](MHz) (29)
In the formula, RopThe ratio of the working MUF to the basic MUF for the F2 layer is shown in table 1.
Table 1 is detailed on page 632 of radio wave propagation (electronic industry Press).
TABLE 1 ratio R of working MUF to basic MUF for layer F2op
Figure BDA0001163779430000111
Step E: and D, calculating the receiving field intensity of the short-wave communication link when the short-wave communication link adopts the highest available frequency for communication according to the highest available frequency obtained in the step D, and analyzing the reliability of the short-wave communication link. The specific process comprises the following steps:
e1, judging the propagation path distance d of the short wave communication link, and if the propagation path distance d is less than 7000km, executing the step E2; if the propagation path distance d is greater than 9000km, executing a step E3; otherwise, steps E2, E3 and E4 are performed simultaneously.
E2, if the distance d of the propagation path is less than 7000km, analyzing and calculating various propagation modes and composite receiving field intensity when the short wave communication link adopts the highest available frequency for communication.
Etw＝94.25+Pt+Gt-20lgp′-Lm-Lg-Lh-Li(30)
In the formula:
Etw-median received field strength, dB μ V/m;
Pt-radiation power, dBkW;
Gt-a transmit antenna gain factor, dB, in transmit azimuth and elevation relative to the isotropic antenna;
p' -skywave ray slant propagation distance, km, calculated using equation (31):
Figure BDA0001163779430000121
in the formula: Δ is given by equation (40).
dnTo use dnThe jump length of the n jump mode calculated as d/n is km;
the formulas (30) and (31) are detailed on page 10 of ITU recommendation ITU-R P.533-10.
LmFor transmission loss factors, dB, above the highest available frequency (MUF), the E-layer and F2-layer propagation modes are calculated using equations (32) and (33), respectively:
Figure BDA0001163779430000122
Figure BDA0001163779430000123
wherein n is the number of days, and t is 0,1,2 … 23.
LgFor the ground reflection loss factor, dB, for the n-hop mode, it is calculated using equation (34):
Lg＝2(n-1) (34)
equations (32) through (34) are detailed on page 11 of ITU recommendation ITU-R P.533-10.
LiFor the ionospheric absorption loss factor, dB, for the n-hop mode, calculated using equation (35):
Figure BDA0001163779430000124
wherein:
F(χ)＝max[cosp(0.881χ),0.02]；
fV＝f×cosi110；
in the formula: f. ofLThe lowest reference frequency, MHz (the average of the electron-revolution frequencies determined at the control points given in table 4, the longitudinal component of the earth's magnetic field at a height of about 100 km);
k is the number of control points (as listed in table 4);
i110is the inclination at 110 km;
R12the average flow value of sunset semen 12 months;
ATnoonto absorb the loss factor, dB, given by fig. 2;
χjthe vertex angle of the sun at the jth control point is 102 degrees, and the smaller is selected;
χjnoonis the x of the local noonjA value;
Figure BDA0001163779430000131
is the equivalent normal incident wave frequency f as the penetration factor of the absorption layerVA function of the ratio to FoE, given by FIG. 3;
p is the daily absorption index, given by figure 4.
LhIn dB, the polar region and other signal loss factors. Tables 2 and 3 show the geomagnetic latitude G according to the Earth's center dipole (78.5N, 68.2W)n(south-north equator) and local time thWhen corresponding value of Gn<At 42.5 deg., L h0. In northern hemisphere, 12 months to 2 months are winter; the spring and autumn seasons are 3 to 5 months and 9 to 11 months; the summer is 6 months to 8 months. In the southern hemisphere, winterThe months of the season and summer are interchanged.
TABLE 2 polar regions with propagation path distances of no more than 2500km and other signal loss factors
Figure BDA0001163779430000132
Figure BDA0001163779430000141
TABLE 3 polar regions with propagation path distances greater than 2500km and other signal loss factors
Figure BDA0001163779430000142
Tables 2 and 3 are detailed on pages 657 and 658 of radio wave transmission (electronic industry Press). Selecting the field intensity of two strong E-layer modes and three strong F2-layer modes from the calculated median value of the field intensity of each propagation mode
And the median is superposed with the power to obtain a composite receiving field intensity median, and the composite receiving field intensity median is calculated by using a formula (36):
Figure BDA0001163779430000151
in the formula:
Etsthe median received field strength, dB μ V/m,
Etw-the value of the field strength, dB μ V/m, for each mode.
The formulas (35) and (36) are detailed in ITU recommendation ITU-R P.533-10 pages 10 to 12.
E3, the distance of the propagation path is larger than 9000km, and the value of the receiving field intensity when the short-wave communication link adopts MUF for communication is calculated.
Dividing the distance of the propagation path into at least n sections according to the average jump distance of not more than 4000km, synthesizing the median value of the field intensity, and calculating by using a formula (37):
Figure BDA0001163779430000152
in the formula:
Etl-median synthetic field strength, dB μ V/m;
E0-free space field strength of 3MW equivalent isotropic radiated power, dB μ V/m, calculated using equation (38):
E0＝139.6-20lg p′ (38)
p' -skywave ray slant propagation distance, km, calculated using equations (39) and (40):
Figure BDA0001163779430000153
Figure BDA0001163779430000154
dnby dnThe jump length of the n jump mode calculated as d/n is km;
Gt1-maximum gain, dB, of the transmitting antenna in the range of 0-8 ° in elevation at the desired azimuth angle;
Gap-gain factor, dB, for field strength increase at distance due to focusing, calculated by equation (41) when d is π R0Taking 15dB when the multiple is obtained;
Figure BDA0001163779430000161
in the formula (41), R0The radius of the earth.
Ly-sky wave propagation effect factor, dB, recommended value-3.7 dB;
fM-highest reference frequency, MHz;
fL-lowest reference frequency, MHz;
fHthe average value of the electromagnetic rotation frequency, MHz, of the control point position of the reflection height of the ray path mirror in the table 4 is obtained.
Pt-transmit power, dBW.
TABLE 4 control point position of reflection height of ray road mirror
Figure BDA0001163779430000162
The formulas (37) and (38) are detailed in pages 658 of radio wave propagation (electronic industry Press); see table 4 on page 634 of radio wave propagation (electronic industry press).
E4, if the propagation path distance is between [7000,9000] km, carrying out field intensity interpolation on the basis of E2 and E3 to obtain the receiving field intensity when the highest available frequency is adopted for communication on the propagation path distance.
In this distance range, the mean value E of the intensity of the sky wave fieldtiBy using EtsAnd EtlAn interpolation calculation is performed, using equation (42):
Eti＝100lg10Xi(42)
in the formula:
Eti-median sky wave field strength, dB μ V/m;
Figure BDA0001163779430000171
Xs＝100.01Ets
Xl＝100.01Etl
Ets-receiving the median field strength, dB μ V/m;
Etl-median synthetic field strength, dB μ V/m.
The formula (42) is described in detail in "radio wave propagation" (electronic industry press) page 660.
E5, analyzing the reliability of the short-wave communication link according to the field intensity calculated in the step E2, E3 or E4.
E51, calculating the median value of the received power according to the field intensity calculated in the step E2, E3 or E4.
S＝E+G-20logf-107.2 (43)
Where E is the reception field strength, G represents the reception antenna gain in the direction of incidence, and f is the communication signal frequency.
E52, calculating the S/N of the middle value of the signal-to-noise ratio and the upper tenth value deviation D of the signal-to-noise ratio according to the received power median value calculated by E51uLower tenth deviation D of SN and SNRlSN。
Figure BDA0001163779430000172
Wherein, FaA、FaM、FaG is the middle value of the atmospheric noise, the artificial noise and the noise coefficient of the Galaxy noise respectively, and b is the bandwidth.
Figure BDA0001163779430000173
Wherein D isuSd、DuShThe deviation of the signal in terms of the number of tenths (daily), the deviation of the signal in terms of the number of tenths (hourly), DlA、DlM、DlG is the ten-point value deviation under atmospheric noise, artificial noise and galaxy noise respectively.
Figure BDA0001163779430000181
Wherein D islSd、DlShThe deviation of the signal in the tenth value (every day), the deviation of the signal in the tenth value (every hour), DuA、DuM、DuG is the upper ten-degree deviation of atmospheric noise, artificial noise and galaxy noise respectively.
And E53, calculating the reliability of the short-wave communication link according to the signal-to-noise ratio middle value, the signal-to-noise ratio upper ten-tenth value deviation and the signal-to-noise ratio lower ten-tenth value deviation which are obtained by calculating in E52.
Figure BDA0001163779430000182
Wherein, S/NrThe signal-to-noise ratio required by the user.
The formula (44) is described in detail in page 663 of radio wave propagation (electronic industry press); the formulas (45) to (47) are described in detail in radio wave propagation (electronic industry Press) 666 and 669.
The above examples are provided only for the purpose of describing the present invention, and are not intended to limit the scope of the present invention. The scope of the invention is defined by the appended claims. Various equivalent substitutions and modifications can be made without departing from the spirit and principles of the invention, and are intended to be within the scope of the invention.

Claims (1)

1. A short wave communication reliability assessment method based on multi-system detection data is characterized by comprising the following steps:
step A: calculating the positions of reflection points of the short-wave communication link and all the detection links, and selecting detection link data close to the positions of ionospheric reflection points of the short-wave communication link, wherein the specific process comprises the following steps:
a1, calculating the longitude and latitude of the ionosphere reflection point of the short-wave communication link according to the short-wave communication physical model and the propagation ray theory;
a2, acquiring the positions and detection types of all detection stations, calculating the longitude and latitude of the ionospheric reflection point of each detection link, and selecting a detection link with the position close to the position of the ionospheric reflection point of the short-wave communication link according to the longitude and latitude of the ionospheric reflection point of the short-wave communication link determined in the step A1, wherein the specific process comprises the following steps:
a21, if the detection link is vertical detection, the longitude and latitude of the detection station are the longitude and latitude of the ionosphere reflection point of the detection link; if the detection link is oblique detection, the method for calculating the longitude and latitude of the ionosphere reflection point of the detection link is the same as the step A1;
a22, passing through the preset distance threshold dlimSelecting N detection links with similar positions of ionosphere reflection points from the short-wave communication link;
a3, reading the data of the link according to the link selected in the step A2, and calculating the ionosphere parameters according to the link data;
if the detection link is vertical detection, directly acquiring transmission factor M (3000) F2 ionosphere parameters of E layer critical frequency foE, F1 layer critical frequency foF1, F2 layer critical frequency foF2 and F2 layer distance of 3000km obtained by detection;
if the detection link is oblique detection, analyzing and obtaining foE, foF1 and M (3000) F2 ionosphere parameters according to the oblique detection data of the highest available frequency MUF and the frequency-time delay, and calculating to obtain foF2 by using a heuristic method on the basis of the M (3000) F2 parameters;
and B: b, forecasting ionospheric parameters of the ionospheric reflection points of the detection link by adopting a weighting analysis method according to the ionospheric reflection point positions of the detection link and the detection link data determined in the step A;
and C: reconstructing to obtain ionospheric parameters of ionospheric reflection points of the short-wave communication link according to the ionospheric reflection point positions of the short-wave communication link and the detection link determined in the step A and the ionospheric parameters of the ionospheric reflection points of the detection link obtained by prediction in the step B, wherein the specific process comprises the following steps:
c1, calculating the ionospheric distance between the ionospheric reflection point of the detection link and the ionospheric reflection point of the short-wave communication link according to the position of the ionospheric reflection point of the detection link determined in the step A;
c2, analyzing and calculating the correlation quantity of the difference value between the ionospheric parameter value of the probe link ionospheric reflection point and the ionospheric parameter value of the reference ionospheric model according to the ionospheric parameter value of the probe link ionospheric reflection point determined in the step B;
c3, according to the calculation results of the steps C1 and C2, the weight coefficients of the ionosphere reflection points of the detection link and the ionosphere reflection points of the short-wave communication link are obtained by utilizing a regional reconstruction equation set;
c4, calculating the correlation quantity of the difference value between the ionospheric reflection point ionospheric parameter of the short-wave communication link and the reference ionospheric model value and the ionospheric reflection point ionospheric parameter value of the short-wave communication link according to the weight coefficient determined in the step C3, wherein the specific process comprises the following steps:
c41, calculating the correlation quantity of the difference value between the ionospheric reflection point ionospheric parameter of the short-wave communication link and the reference ionospheric model value according to the weight coefficient determined in the step C3;
c42, calculating ionospheric reflection point ionospheric parameter values of the short-wave communication link according to the correlation quantity determined by C41;
step D: and C, predicting the highest available frequency of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link obtained in the step C, wherein the specific process comprises the following steps:
d1, calculating the highest available frequency E (D) MUF of the E layer of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link analyzed and obtained in the step C;
d2, calculating the F1 layer highest available frequency F1(D) MUF of the short-wave communication link according to the ionized layer parameter value of the ionized layer reflection point of the short-wave communication link analyzed and obtained in the step C;
d3, calculating the F2 layer highest available frequency F2(D) MUF of the short-wave communication link according to the ionized layer reflection point parameter value of the short-wave communication link obtained by analyzing in the step C and by combining the propagation path distance condition of the short-wave communication link, wherein the specific process comprises the following steps:
d31, judging whether the short wave communication link propagation path distance D is larger than the single-hop maximum hop distance DmaxAccording to whether d is greater than dmaxRespectively calculating the highest available frequency F2(d) MUF of the F2 layer;
d32, if the propagation path distance D is less than or equal to the maximum hop distance D of one hopmaxCombining the ionosphere parameter values obtained by analysis in the step C, and calculating F2 layer basic MUF according to a single-hop mode; if the propagation path distance d is greater than the maximum hop distance d of a single hopmaxCalculating F2 layer basic MUF according to the multi-hop mode of radio wave propagation;
step D4, calculating the highest available frequency MUF of the short-wave communication link according to the E (D) MUF, F1(D) MUF and F2(D) MUF analyzed and obtained in the steps D1, D2 and D3;
step E: and D, calculating the receiving field intensity of the short-wave communication link when the short-wave communication link adopts the highest available frequency for communication according to the highest available frequency obtained in the step D, and analyzing the reliability of the short-wave communication link, wherein the specific process comprises the following steps:
e1, judging the propagation path distance d of the short wave communication link, and if the propagation path distance d is less than 7000km, executing the step E2; if the propagation path distance d is greater than 9000km, executing a step E3; otherwise, executing step E4;
e2, if the distance d of the propagation path is less than 7000km, analyzing and calculating various propagation modes and synthesized receiving field intensity when the shortwave communication link adopts the highest available frequency for communication;
e3, if the distance d of the propagation path is greater than 9000km, calculating a receiving field intensity value when the shortwave communication link adopts the highest available frequency for communication;
e4, if the propagation path distance d is between [7000,9000] km, carrying out field intensity interpolation on the basis of executing E2 and E3 to obtain the receiving field intensity when the highest available frequency is adopted for communication on the propagation path distance;
e5, analyzing the reliability of the short-wave communication link according to the field intensity calculated in the step E2, E3 or E4, wherein the specific process comprises the following steps:
e51, calculating the median value of the received power according to the field intensity calculated in the step E2, E3 or E4;
e52, calculating the signal-to-noise ratio middle value, the upper ten value deviation of the signal-to-noise ratio and the lower ten value deviation of the signal-to-noise ratio according to the received power median value calculated by E51;
and E53, calculating the reliability of the short-wave communication link according to the signal-to-noise ratio middle value, the signal-to-noise ratio upper ten-tenth value deviation and the signal-to-noise ratio lower ten-tenth value deviation which are obtained by calculating in E52.
CN201611071301.9A 2016-11-28 2016-11-28 Short wave communication reliability assessment method based on multi-system detection data Active CN106788815B (en)

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 CN106788815A (en) 2017-05-31
CN106788815B true 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)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
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
CN110288117B (en) * 2019-05-10 2021-06-29 中国人民解放军31007部队 Regional reconstruction method for critical frequency of ionosphere parameters
CN110245316A (en) * 2019-05-10 2019-09-17 中国人民解放军31007部队 A kind of inversion method of Ionospheric Parameters
CN111148145B (en) * 2019-12-20 2022-03-22 中国电波传播研究所（中国电子科技集团公司第二十二研究所） Short wave sky wave channel maximum accessible bandwidth evaluation method
CN111881584B (en) * 2020-07-30 2023-02-28 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 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
CN112272067B (en) * 2020-10-15 2022-04-08 天津大学 Short wave broadcast frequency scheduling method based on multi-source data processing
CN113109632B (en) * 2021-04-08 2023-04-14 三门峡职业技术学院 Method for inverting F2 layer parameters by using oblique ionogram

Citations (1)

* Cited by examiner, † Cited by third party
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

Patent Citations (1)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
ITU-R methods of basic MUF, operational MUF and ray-path prediction;国际电信联盟;《Rec.ITU-R P.1240-1》;20070201;第3页 *
Method for the prediction of the;国家电信联盟;《Recommendation ITU-R P.533-10》;20091001;第10页 *
中国地区电离层f oF2重构方法及其在短期预报中的应用;刘瑞源，张北辰等;《地球物理学报》;20080331;第2-3页 *
付炜，王健等.高频通信可用频率现报与预报方法.《中国电子科学研究院学报》.2013, *
高频通信可用频率现报与预报方法;付炜，王健等;《中国电子科学研究院学报》;20131031;第2页 *

Also Published As

Publication number Publication date
CN106788815A (en) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106788815B (en) Short wave communication reliability assessment method based on multi-system detection data
Matolak et al. Air–ground channel characterization for unmanned aircraft systems—Part I: Methods, measurements, and models for over-water settings
CN109964422B (en) Method and apparatus for analyzing communication channel and designing wireless network considering information about actual environment
Willink et al. Measurement and characterization of low-altitude air-to-ground MIMO channels
CN102843750B (en) Short-wave frequency searching device and control method thereof
US11032791B2 (en) Method and system for dynamically adjusting communication characteristics of satellite communication links based on weather nowcasting
CN106231621A (en) A kind of many scene adaptives optimization method of propagation model in FDD LTE system
Francis et al. Performance prediction for design of a network of skywave over-the-horizon radars
Bhardwaj et al. Geometrical-empirical channel propagation model for human presence at 60 GHz
CN109217947A (en) The method that test data chain radiation of equipment signal is trapped probability performance
Chou et al. Fast position-aided MIMO beam training via noisy tensor completion
Oestges et al. Inclusive radio communications for 5G and beyond
Yussuff et al. Modified ITU-R rain attenuation prediction model for a tropical station
CN106850036A (en) A kind of removable spot beam dispatching method of middle rail satellite system based on priority
CN117082536A (en) Air-ground network collaborative coverage method based on reinforcement learning
CN107395245B (en) A kind of radio astronomy and the dual-purpose Science Popularization System of amateur radio communication
WO2007064953A2 (en) Hf radio network frequency management
Abdulrahman et al. An improved slant path attenuation prediction method in tropical climates
Biscarini et al. Statistical modeling of atmospheric propagation channel at W-band through sun-tracking microwave radiometric measurements for non-geostationary satellite links
Laspougeas et al. Radio propagation in urban small cells environment at 2 GHz: experimental spatio-temporal characterization and spatial wideband channel model
Luini et al. On the relationship between the spatial correlation of point rain rate and of rain attenuation on earth-space radio links
Paulson Estimating 1 min rain rate distributions from numerical weather prediction
Matolak et al. Ag channel measurement and modeling results for over-water and hilly terrain conditions
US20240323719A1 (en) Radio frequency radiance field models for communication system control
Prakash et al. Digital receiver based Ka Band beacon receiver for improved beacon power estimation

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