CN114660537B - Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment - Google Patents
Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment Download PDFInfo
- Publication number
- CN114660537B CN114660537B CN202210396710.5A CN202210396710A CN114660537B CN 114660537 B CN114660537 B CN 114660537B CN 202210396710 A CN202210396710 A CN 202210396710A CN 114660537 B CN114660537 B CN 114660537B
- Authority
- CN
- China
- Prior art keywords
- layer
- echo
- frequency point
- pitch angle
- virtual height
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 69
- 239000005433 ionosphere Substances 0.000 title claims description 31
- 238000010586 diagram Methods 0.000 claims description 28
- 238000001514 detection method Methods 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 6
- 230000007547 defect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a method and equipment for estimating a sky wave pitch angle of an ionized layer F based on an Es layer. The method comprises the following steps: step 1 to step 4. According to the method, under the condition that the Es layer height is obtained by adopting the altimeter at the Es layer, the pitching angle of the sky wave of the F layer can be estimated only by erecting the transmitting station and the receiving station of a single antenna, the method is simple and high in precision, and the error of the pitching angle in a simulation experiment is within 1 degree.
Description
Technical Field
The embodiment of the invention relates to the technical field of ionosphere sky wave pitch angles, in particular to a method and equipment for estimating an ionosphere F layer sky wave pitch angle based on an Es layer.
Background
The ionosphere is a high-rise atmospheric region partially ionized from about 60km to about 1000km from the ground, and can be divided into 4 different regions, namely a D layer, an E layer, an F1 layer, and an F2 layer (collectively referred to as an F layer), according to electron density. F2 is layered above about 200km from the ground and is the main area for reflecting radio signals. There are sometimes very thin (several kilometers) of strongly ionized layer E in the E layer. Also called "burst E layer" (Es layer). Es is an inhomogeneity in the height of the ionosphere E layer, whose shape sometimes becomes a plate, completely masking the upper F layer, sometimes one by one, like a "flake". Es "flakes" can range up to tens or hundreds of kilometers and drift at speeds of about 500 m/s. Because the thickness of the Es layer is very thin, the false height h' Es of the Es layer echo obtained by detecting through a height measuring instrument is generally considered to be similar to the real height of the Es layer. The existing sky wave pitch angle estimation method mostly adopts array antennas to carry out beam forming for estimation, and has the main defect that a large array antenna is required to be used for receiving sky wave signals. Therefore, developing a method and a device for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer can effectively overcome the above-mentioned defects in the related art, and is a technical problem to be solved in the industry.
Disclosure of Invention
Aiming at the problems in the prior art, the embodiment of the invention provides a method and equipment for estimating a sky wave pitch angle of an ionized layer F based on an Es layer.
In a first aspect, an embodiment of the present invention provides a method for estimating a sky-wave pitch angle of an ionosphere F layer based on an Es layer, including: step 1: performing an oblique return detection experiment by using a transmitting station and a receiving station of a single antenna to obtain an echo diagram, detecting by using a vertical measuring instrument to obtain the height of an Es layer near the transmitting station or near an echo relay point, and calculating the minimum value and the maximum value of the pitching angle of the oblique return wave of a frequency point based on the echo of the Es layer in the echo diagram for the frequency point; and 2, step: obtaining the minimum value and the maximum value of an echo group path of a frequency point F layer by interpreting an echo diagram, obtaining the maximum value of the echo virtual height of the frequency point F layer through the maximum value of the pitch angle of the frequency point and the minimum value of the echo virtual height of the F layer, obtaining the minimum value of the echo virtual height of the frequency point F layer through the minimum value of the pitch angle of the frequency point and the maximum value of the echo virtual height of the F layer, and further obtaining a relational expression of the echo virtual height and the echo virtual height of the frequency point F layer; and step 3: for an echo of a frequency point F layer, obtaining the virtual height of the echo according to the relation between the group path and the virtual height of the F layer in the step 2, and obtaining the pitch angle of the echo by adopting a cosine formula; and 4, step 4: and constructing a relation between the group path and the virtual height of the F layer for the rest frequency points by using the same method, and realizing the estimation of the echo pitch angle of the F layer of each frequency point.
Based on the content of the above method embodiment, the method for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer provided in the embodiment of the present invention calculates the minimum value and the maximum value of the pitch angle of the one-frequency-point slant return wave based on the echo of the Es layer in the echo diagram in step 1, and includes:
wherein, theta i min Represents the minimum value of the pitch angle of the echo corresponding to the frequency point fi, theta i max Representing the maximum value of the echo pitch angle P of the frequency point fi iEs max Maximum value of Es layer echo group path, P, representing frequency point fi iEs min Represents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, h Es Indicating Es layer height.
On the basis of the content of the above method embodiment, the method for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer provided in the embodiment of the present invention obtains the maximum value of the false height of the echo of the F layer of a frequency point through the maximum value of the pitch angle of the frequency point and the minimum value of the group path of the F layer in step 2, and obtains the minimum value of the false height of the echo of the F layer of the frequency point through the minimum value of the pitch angle of the frequency point and the maximum value of the group path of the F layer, including:
hv i F max =sin(θ i max )·P i F min /2
hv i F min =sin(θ i min )·P i F max /2
wherein, P i F min Minimum value, P, of F-layer echo group path representing frequency point fi i F max Maximum value, hv, of the F-layer echo group path representing frequency point fi i F min Minimum value of F-layer echo virtual height, hv, representing frequency point fi i F max And the maximum value of the F-layer echo virtual height of the frequency point fi is represented.
On the basis of the content of the above method embodiment, in the method for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer provided in the embodiment of the present invention, the relation between the group path of the one-frequency F layer and the virtual height in step 2 includes:
wherein, P ij Group path, hv, of the j-th echo of the F layer representing frequency point fi ij And indicating the virtual height corresponding to the jth echo of the F layer of the frequency point fi.
Based on the content of the above method embodiment, the method for estimating the sky-wave pitch angle of the ionosphere F layer based on the Es layer according to the embodiment of the present invention obtains the virtual height of an echo according to the relationship between the group path of the F layer and the virtual height in step 2, and obtains the pitch angle of an echo by using a cosine formula, including:
wherein, theta ij And the pitch angle of the jth echo of the F layer at the frequency point fi is represented.
In a second aspect, an embodiment of the present invention provides an ionosphere F-layer sky-wave pitch angle estimation apparatus based on the Es layer, including: a first master module, configured to implement step 1: performing an oblique return detection experiment by using a transmitting station and a receiving station of a single antenna to obtain an echo diagram, detecting by using a vertical measuring instrument to obtain the height of an Es layer near the transmitting station or near an echo relay point, and calculating the minimum value and the maximum value of the pitching angle of the oblique return wave of a frequency point based on the echo of the Es layer in the echo diagram for the frequency point; a second master module, configured to implement step 2: obtaining the minimum value and the maximum value of an echo group path of a frequency point F layer by interpreting an echo diagram, obtaining the maximum value of the echo virtual height of the frequency point F layer through the maximum value of the pitch angle of the frequency point and the minimum value of the echo group path of the F layer, obtaining the minimum value of the echo virtual height of the frequency point F layer through the minimum value of the pitch angle of the frequency point and the maximum value of the echo virtual height of the F layer, and further obtaining a relational expression of the echo virtual height and the echo virtual height of the frequency point F layer; a third main module, configured to implement step 3: for an echo of a frequency point F layer, obtaining the virtual height of the echo according to the relation between the group path and the virtual height of the F layer in the step 2, and obtaining the pitch angle of the echo by adopting a cosine formula; a fourth master module, configured to implement step 4: and constructing a relation between the group path and the virtual height of the F layer for the rest frequency points by using the same method, and realizing the estimation of the echo pitch angle of the F layer of each frequency point.
In a third aspect, an embodiment of the present invention provides an electronic device, including:
at least one processor; and
at least one memory communicatively coupled to the processor, wherein:
the memory stores program instructions executable by the processor, and the processor calls the program instructions to execute the method for estimating the sky wave pitch angle of the ionosphere F based on the Es layer provided by any one of the various implementations of the first aspect.
In a fourth aspect, embodiments of the present invention provide a non-transitory computer-readable storage medium storing computer instructions that cause a computer to perform a method for estimating an ionospheric F-layer sky-wave pitch angle based on the Es-layer provided in any of the various implementations of the first aspect.
According to the method and the device for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer, the E layer is subjected to altimetry by adopting a altimeter, the sky wave pitch angle of the F layer can be estimated only by erecting a transmitting station and a receiving station of a single antenna, the method is simple and high in precision, and the pitch angle error of a simulation experiment is within 1 degree.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, a brief description will be given below to the drawings required for the description of the embodiments or the prior art, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a flowchart of a method for estimating a sky-wave pitch angle of an ionosphere F layer based on an Es layer according to an embodiment of the present invention;
fig. 2 is a schematic structural diagram of an ionosphere F-layer sky-wave pitch angle estimation device based on the Es layer according to an embodiment of the present invention;
fig. 3 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention;
fig. 4 is a comparison graph of the actual pitch angle and the estimated pitch angle at 15MHz provided by an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention. In addition, technical features of various embodiments or individual embodiments provided by the present invention may be arbitrarily combined with each other to form a feasible technical solution, and such combination is not limited by the sequence of steps and/or the structural composition mode, but must be realized by a person skilled in the art, and when the technical solution combination is contradictory or cannot be realized, such a technical solution combination should not be considered to exist and is not within the protection scope of the present invention.
The embodiment of the invention provides an ionosphere F layer sky wave pitch angle estimation method based on an Es layer, and referring to FIG. 1, the method comprises the following steps: step 1: performing an oblique return detection experiment by using a transmitting station and a receiving station of a single antenna to obtain an echo diagram, detecting by using a vertical measuring instrument to obtain the height of an Es layer near the transmitting station or near an echo relay point, and calculating the minimum value and the maximum value of the pitching angle of the oblique return wave of a frequency point based on the echo of the Es layer in the echo diagram for the frequency point; step 2: obtaining the minimum value and the maximum value of an echo group path of a frequency point F layer by interpreting an echo diagram, obtaining the maximum value of the echo virtual height of the frequency point F layer through the maximum value of the pitch angle of the frequency point and the minimum value of the echo group path of the F layer, obtaining the minimum value of the echo virtual height of the frequency point F layer through the minimum value of the pitch angle of the frequency point and the maximum value of the echo virtual height of the F layer, and further obtaining a relational expression of the echo virtual height and the echo virtual height of the frequency point F layer; and 3, step 3: for an echo of a frequency point F layer, obtaining the virtual height of the echo according to the relation between the group path and the virtual height of the F layer in the step 2, and obtaining the pitch angle of the echo by adopting a cosine formula; and 4, step 4: and constructing a relation between the group path and the virtual height of the F layer for the rest frequency points by using the same method, and realizing the estimation of the echo pitch angle of the F layer of each frequency point.
Based on the content of the above method embodiment, as an optional embodiment, the method for estimating the sky-wave pitch angle of the ionosphere F layer based on the Es layer provided in the embodiment of the present invention calculates the minimum value and the maximum value of the pitch angle of the one-frequency-point slant return wave based on the echo of the Es layer in the echo diagram in step 1, and includes:
wherein, theta i min Representing the minimum value of the pitch angle of the echo corresponding to the frequency point fi, theta i max Representing the maximum value of the echo pitch angle P of the frequency point fi iEs max Maximum value of Es layer echo group path, P, representing frequency point fi iEs min Represents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, h Es Indicating Es layer height.
In another embodiment, a numerical ray tracing is adopted to simulate an oblique return detection experiment, the longitude of a transmitting station of the experiment is 113 degrees E, the latitude of the transmitting station is 30 degrees N, firstly, an IRI model is utilized to generate a background electron density map of which the distance from 14 points in 12, 25 and 25 months in 2019 to the azimuth of the transmitting station is 180 degrees and the range is 2000km, an Es layer of which the range is 600km is added at the position of 110km of the height, and the adjacent frequency corresponding to the electron density of the Es layer is 9MHz. The detection experiment adopts sweep frequency detection, the frequency is from 15MHz to 16MHz, the frequency is stepped by 0.2MHz, the elevation angle range is from 14 degrees to 19 degrees, and the elevation angle step is 1 degree.
And judging the maximum value and the minimum value of the Es layer echo group path under 15MHz, and calculating the minimum value and the maximum value of the pitch angle of the frequency point according to a formula (1). Wherein, theta i min =13.83 degrees represents the minimum value of the echo pitch angle corresponding to the frequency point fi, theta i max =18.79 ° for frequency fi returnMaximum wave pitch angle, P iEs max =920.32km represents Es layer echo group path maximum value, P of frequency point fi iEs min =682.89km means Es layer echo group path minimum at frequency fi, R =6370km means earth radius, h Es =110km represents Es layer height.
Based on the content of the foregoing method embodiment, as an optional embodiment, the method for estimating an sky wave pitch angle of an ionosphere F layer based on an Es layer provided in the embodiment of the present invention, in step 2, the maximum value of an imaginary height of an echo in an F layer of a frequency point is obtained through the maximum value of the pitch angle of the frequency point and the minimum value of a group path in the F layer, and the minimum value of an imaginary height of an echo in an F layer of the frequency point is obtained through the minimum value of the pitch angle of the frequency point and the maximum value of the group path in the F layer, including:
wherein, P i F min Minimum value, P, of F-layer echo group path representing frequency point fi i F max Maximum value, hv, of the F-layer echo group path representing frequency point fi i F min Minimum value of F-layer echo virtual height, hv, representing frequency point fi i F max And the maximum value of the F-layer echo virtual height of the frequency point fi is represented.
In another embodiment, the minimum value and the maximum value of the group path of the F-layer echo at the frequency point of 15MHz are interpreted, and the maximum value and the minimum value of the virtual height of the F-layer echo are obtained by using the formula (2). Wherein P is i F min =1243.73km means minimum value of F-layer echo group path at frequency point 15MHz, P i F max =1360.41km denotes the maximum value of the F-layer echo group path at frequency point 15MHz, hv i F min =162.6km represents the minimum value of the F-layer echo virtual height at frequency point 15MHz, hv i F max And =200.3km represents the maximum value of the F-layer echo imaginary height at the frequency point of 15 MHz.
Based on the content of the foregoing method embodiment, as an optional embodiment, in the method for estimating an ionosphere F-layer sky-wave pitch angle based on the Es layer provided in the embodiment of the present invention, the relational expression between the one-frequency F-layer group path and the virtual height in step 2 includes:
wherein, P ij Group path, hv, of the j-th echo of the F layer representing frequency point fi ij And indicating the virtual height corresponding to the jth echo of the F layer of the frequency point fi.
In another embodiment, the relationship between the group path and the virtual height at the frequency point of 15MHz is constructed as shown in formula (3). Wherein P is ij Group path, hv, of the j-th echo of F layer representing frequency point 15MHz ij And the virtual height corresponding to the jth echo of the F layer at the frequency point of 15MHz is shown.
Based on the content of the foregoing method embodiment, as an optional embodiment, the method for estimating a sky-wave pitch angle of an ionosphere F layer based on an Es layer provided in the embodiment of the present invention obtains an imaginary height of an echo according to a relational expression between an F layer group path and the imaginary height in step 2, and obtains a pitch angle of an echo by using a cosine formula, including:
wherein, theta ij And the pitch angle of the jth echo of the F layer at the frequency point fi is represented.
Specifically, the 3 rd echo group path Pi3=1290.70km of the echo at 15MHz, for example, and the corresponding virtual height hvi =185.1km, the estimated value θ i3=16.03 ° is obtained according to the pitch angle calculation formula, and substantially coincides with the actual elevation angle 16.00 °, and the comparison between the actual pitch angle and the estimated pitch angle at 15MHz is shown in fig. 4.
According to the method for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer, provided by the embodiment of the invention, under the condition that the Es layer height is obtained by adopting a height finder on the Es layer, the sky wave pitch angle of the F layer can be estimated only by erecting a transmitting station and a receiving station of a single antenna, the method is simple and high in precision, and the pitch angle error of a simulation experiment is within 1 degree.
The implementation basis of the various embodiments of the present invention is realized by programmed processing performed by a device having a processor function. Therefore, in engineering practice, the technical solutions and functions thereof of the embodiments of the present invention can be packaged into various modules. Based on this reality, on the basis of the foregoing embodiments, embodiments of the present invention provide an ionospheric F-layer sky-wave pitch angle estimation apparatus based on the Es-layer, which is used to execute the ionospheric F-layer sky-wave pitch angle estimation method based on the Es-layer in the foregoing method embodiments. Referring to fig. 2, the apparatus includes: a first master module, configured to implement step 1: performing an oblique return detection experiment by using a transmitting station and a receiving station of a single antenna to obtain an echo diagram, detecting by using a vertical measuring instrument to obtain the height of an Es layer near the transmitting station or near an echo relay point, and calculating the minimum value and the maximum value of the pitching angle of the oblique return wave of a frequency point based on the echo of the Es layer in the echo diagram for the frequency point; a second master module, configured to implement step 2: obtaining the minimum value and the maximum value of an echo group path of a frequency point F layer by interpreting an echo diagram, obtaining the maximum value of the echo virtual height of the frequency point F layer through the maximum value of the pitch angle of the frequency point and the minimum value of the echo group path of the F layer, obtaining the minimum value of the echo virtual height of the frequency point F layer through the minimum value of the pitch angle of the frequency point and the maximum value of the echo virtual height of the F layer, and further obtaining a relational expression of the echo virtual height and the echo virtual height of the frequency point F layer; a third main module, configured to implement step 3: for an echo of a frequency point F layer, obtaining the virtual height of the echo according to the relation between the group path and the virtual height of the F layer in the step 2, and obtaining the pitch angle of the echo by adopting a cosine formula; a fourth master module, configured to implement step 4: and constructing a relation between the group path and the virtual height of the F layer for the rest frequency points by using the same method, and realizing the estimation of the echo pitch angle of the F layer of each frequency point.
According to the ionosphere F layer sky wave pitch angle estimation device based on the Es layer, provided by the embodiment of the invention, the modules in the figure 2 are adopted, and under the condition that the Es layer height is obtained by adopting a height indicator in the Es layer, the F layer sky wave pitch angle can be estimated only by erecting a transmitting station and a receiving station of a single antenna, so that the method is simple and high in precision, and the pitch angle error of a simulation experiment is within 1 degree.
It should be noted that, the apparatus in the apparatus embodiment provided by the present invention may be used for implementing methods in other method embodiments provided by the present invention, except that corresponding function modules are provided, and the principle of the apparatus embodiment provided by the present invention is basically the same as that of the apparatus embodiment provided by the present invention, so long as a person skilled in the art obtains corresponding technical means by combining technical features on the basis of the apparatus embodiment described above, and obtains a technical solution formed by these technical means, on the premise of ensuring that the technical solution has practicability, the apparatus in the apparatus embodiment described above may be modified, so as to obtain a corresponding apparatus class embodiment, which is used for implementing methods in other method class embodiments. For example:
based on the content of the above device embodiment, as an optional embodiment, the device for estimating a sky-wave pitch angle of an ionosphere F layer based on an Es layer provided in the embodiment of the present invention further includes: the first submodule is used for calculating the minimum value and the maximum value of a pitching angle of a frequency point of the slant return wave based on the echo of the Es layer in the echo diagram in the step 1, and comprises:
wherein, theta i min Representing the minimum value of the pitch angle of the echo corresponding to the frequency point fi, theta imax represents the maximum value of the echo pitch angle of the frequency point fi, P iEs max Maximum value of Es layer echo group path, P, representing frequency point fi iEs min Represents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, h Es Indicating Es layer height.
Based on the content of the above device embodiment, as an optional embodiment, the device for estimating a sky-wave pitch angle of an ionosphere F layer based on an Es layer provided in the embodiment of the present invention further includes: the second submodule is used for obtaining the maximum value of the virtual height of the echo of the F layer of the frequency point through the maximum value of the pitch angle of the frequency point and the minimum value of the group path of the F layer in the step 2, and obtaining the minimum value of the virtual height of the echo of the F layer of the frequency point through the minimum value of the pitch angle of the frequency point and the maximum value of the group path of the F layer, and comprises:
hv i F max =sin(θ i max )·P i F min /2
hv i F min =sin(θ i min )·P i F max /2
wherein, P i F min Minimum value, P, of F-layer echo group path representing frequency point fi i F max Maximum value, hv, of the F-layer echo group path representing frequency point fi i F min Minimum value of F-layer echo virtual height, hv, representing frequency point fi i F max And the maximum value of the F-layer echo virtual height of the frequency point fi is represented.
Based on the content of the foregoing device embodiment, as an optional embodiment, the device for estimating sky-wave pitch angle of ionosphere F layer based on Es layer provided in the embodiment of the present invention further includes: a third sub-module, configured to implement the relation between the group path of the one-frequency F layer and the virtual height in step 2, including:
wherein, P ij Group path, hv, of the j-th echo of the F layer representing frequency point fi ij And indicating the virtual height corresponding to the jth echo of the F layer of the frequency point fi.
Based on the content of the foregoing device embodiment, as an optional embodiment, the device for estimating sky-wave pitch angle of ionosphere F layer based on Es layer provided in the embodiment of the present invention further includes: the fourth sub-module is configured to obtain the virtual height of an echo according to the relationship between the group path of the F layer and the virtual height in step 2, and obtain the pitch angle of an echo by using a cosine formula, and includes:
wherein, theta ij And the pitch angle of the jth echo of the F layer at the frequency point fi is represented.
The method of the embodiment of the invention is realized by depending on the electronic equipment, so that the related electronic equipment is necessarily introduced. To this end, an embodiment of the present invention provides an electronic apparatus, as shown in fig. 3, including: the system comprises at least one processor (processor), a communication Interface (Communications Interface), at least one memory (memory) and a communication bus, wherein the at least one processor, the communication Interface and the at least one memory are communicated with each other through the communication bus. The at least one processor may invoke logic instructions in the at least one memory to perform all or a portion of the steps of the methods provided by the various method embodiments described above.
In addition, the logic instructions in the at least one memory may be implemented in software functional units and stored in a computer readable storage medium when sold or used as a stand-alone product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the method embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one position, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. Based on this recognition, each block in the flowchart or block diagrams may represent a module, a program segment, or a portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In this patent, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising … …" does not exclude the presence of another like element in a process, method, article, or apparatus that comprises the element.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, and not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (8)
1. An ionosphere F layer sky wave pitch angle estimation method based on an Es layer is characterized by comprising the following steps: step 1: performing an oblique return detection experiment by using a transmitting station and a receiving station of a single antenna to obtain an echo diagram, detecting by using a vertical measuring instrument to obtain the height of an Es layer near the transmitting station or near an echo relay point, and calculating the minimum value and the maximum value of the pitching angle of the oblique return wave of a frequency point based on the echo of the Es layer in the echo diagram for the frequency point; and 2, step: obtaining the minimum value and the maximum value of an echo group path of a frequency point F layer by interpreting an echo diagram, obtaining the maximum value of the echo virtual height of the frequency point F layer through the maximum value of the pitch angle of the frequency point and the minimum value of the echo group path of the F layer, obtaining the minimum value of the echo virtual height of the frequency point F layer through the minimum value of the pitch angle of the frequency point and the maximum value of the echo virtual height of the F layer, and further obtaining a relational expression of the echo virtual height and the echo virtual height of the frequency point F layer; and step 3: for an echo of a frequency point F layer, obtaining the virtual height of the echo according to the relation between the group path and the virtual height of the F layer in the step 2, and obtaining the pitch angle of the echo by adopting a cosine formula; and 4, step 4: and constructing a relation between the group path and the virtual height of the F layer for the rest frequency points by using the same method, and realizing the estimation of the echo pitch angle of the F layer of each frequency point.
2. The method for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer according to claim 1, wherein the step 1 of calculating the minimum value and the maximum value of the pitch angle of the one-frequency-point slant return wave based on the echo of the Es layer in the echo diagram comprises:
wherein, theta imin Represents the minimum value of the pitch angle of the echo corresponding to the frequency point fi, theta imax Representing the maximum value of the echo pitch angle P of the frequency point fi iEsmax Maximum value of Es layer echo group path, P, representing frequency point fi iEsmin Represents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, h Es Indicating Es layer height.
3. The method for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer according to claim 2, wherein the step 2 of obtaining the maximum virtual height of the echo of the F layer of a frequency point through the maximum value of the pitch angle of the frequency point and the minimum value of the group path of the F layer, and obtaining the minimum virtual height of the echo of the F layer of the frequency point through the minimum value of the pitch angle of the frequency point and the maximum value of the group path of the F layer comprises the following steps:
hv iFmax =sin(θ imax )·P iFmin /2
hv iFmin =sin(θ imin )·P iFmax /2
wherein, P iFmin Minimum value, P, of F-layer echo group path representing frequency point fi iFmax Maximum value, hv, of the F-layer echo group path representing frequency point fi iFmin Minimum value of F-layer echo virtual height, hv, representing frequency point fi iFmax And the maximum value of the F-layer echo virtual height of the frequency point fi is represented.
4. The method of claim 3, wherein the relationship between the group path of the first frequency point and the F-layer and the virtual height in step 2 comprises:
wherein, P ij Group path, hv, of the j-th echo of the F layer representing frequency point fi ij And indicating the virtual height corresponding to the jth echo of the F layer of the frequency point fi.
5. The method of claim 4, wherein the obtaining of the elevation of an echo according to the relation between the group path of the F layer and the elevation in step 2 and the obtaining of the elevation of an echo by the cosine formula comprises:
wherein, theta ij And the pitch angle of the jth echo of the F layer representing the frequency point fi.
6. An ionospheric F-layer sky-wave pitch angle estimation device based on an Es layer, comprising: a first main module, configured to implement step 1: performing an oblique return detection experiment by using a transmitting station and a receiving station of a single antenna to obtain an echo diagram, detecting by using a vertical measuring instrument to obtain the height of an Es layer near the transmitting station or near an echo relay point, and calculating the minimum value and the maximum value of the pitching angle of the oblique return wave of a frequency point based on the echo of the Es layer in the echo diagram for the frequency point; a second master module, configured to implement step 2: obtaining the minimum value and the maximum value of an echo group path of a frequency point F layer by interpreting an echo diagram, obtaining the maximum value of the echo virtual height of the frequency point F layer through the maximum value of the pitch angle of the frequency point and the minimum value of the echo group path of the F layer, obtaining the minimum value of the echo virtual height of the frequency point F layer through the minimum value of the pitch angle of the frequency point and the maximum value of the echo virtual height of the F layer, and further obtaining a relational expression of the echo virtual height and the echo virtual height of the frequency point F layer; a third main module, configured to implement step 3: for an echo of a frequency point F layer, obtaining the virtual height of the echo according to the relation between the group path and the virtual height of the F layer in the step 2, and obtaining the pitch angle of the echo by adopting a cosine formula; a fourth main module, configured to implement step 4: and constructing a relation between the group path and the virtual height of the F layer for the rest frequency points by using the same method, and realizing the estimation of the echo pitch angle of the F layer of each frequency point.
7. An electronic device, comprising:
at least one processor, at least one memory, and a communication interface; wherein the content of the first and second substances,
the processor, the memory and the communication interface are communicated with each other;
the memory stores program instructions executable by the processor, the processor invoking the program instructions to perform the method of any of claims 1 to 5.
8. A non-transitory computer-readable storage medium storing computer instructions for causing a computer to perform the method of any one of claims 1 to 5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210396710.5A CN114660537B (en) | 2022-04-15 | 2022-04-15 | Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210396710.5A CN114660537B (en) | 2022-04-15 | 2022-04-15 | Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114660537A CN114660537A (en) | 2022-06-24 |
CN114660537B true CN114660537B (en) | 2022-11-15 |
Family
ID=82036070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210396710.5A Expired - Fee Related CN114660537B (en) | 2022-04-15 | 2022-04-15 | Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114660537B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5943629A (en) * | 1996-04-01 | 1999-08-24 | Tci International, Inc. | Method and apparatus for real-time ionospheric mapping and dynamic forecasting |
CN101581782A (en) * | 2009-06-15 | 2009-11-18 | 武汉大学 | Method for inhibiting ionospheric clutter in portable high frequency groundwave radar |
CN108872948A (en) * | 2018-08-01 | 2018-11-23 | 哈尔滨工业大学 | A kind of high-frequency ground wave radar ionospheric clutter suppressing method |
CN111007489A (en) * | 2019-12-05 | 2020-04-14 | 武汉大学 | Buoy type ionized layer double-pass oblique return detection system and method |
CN111007490A (en) * | 2019-12-05 | 2020-04-14 | 武汉大学 | Sky wave over-the-horizon radar coordinate registration method based on buoy geographic information |
CN113376631A (en) * | 2021-08-11 | 2021-09-10 | 武汉大学 | Super-resolution vertical detection method for ionized layer Es (electron beams) |
-
2022
- 2022-04-15 CN CN202210396710.5A patent/CN114660537B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5943629A (en) * | 1996-04-01 | 1999-08-24 | Tci International, Inc. | Method and apparatus for real-time ionospheric mapping and dynamic forecasting |
CN101581782A (en) * | 2009-06-15 | 2009-11-18 | 武汉大学 | Method for inhibiting ionospheric clutter in portable high frequency groundwave radar |
CN108872948A (en) * | 2018-08-01 | 2018-11-23 | 哈尔滨工业大学 | A kind of high-frequency ground wave radar ionospheric clutter suppressing method |
CN111007489A (en) * | 2019-12-05 | 2020-04-14 | 武汉大学 | Buoy type ionized layer double-pass oblique return detection system and method |
CN111007490A (en) * | 2019-12-05 | 2020-04-14 | 武汉大学 | Sky wave over-the-horizon radar coordinate registration method based on buoy geographic information |
CN113376631A (en) * | 2021-08-11 | 2021-09-10 | 武汉大学 | Super-resolution vertical detection method for ionized layer Es (electron beams) |
Non-Patent Citations (4)
Title |
---|
Ionosonde observations of daytime spread F at low latitudes;Chunhua Jiang et al.;《Journal of Geophysical Research:Space Physics》;20161130;第121卷(第12期);第12093-12103页 * |
Large-Scale Ionospheric Irregularities Detected by Ionosonde and GNSS Receiver Network;Chunhua Jiang et al.;《IEEE Geoscience and Remote Sensing Letters》;20200508;第18卷(第6期);第940-943页 * |
电离层Es电波传播模型研究;张雅彬 等;《空间科学学报》;20190131;第39卷(第1期);第55-61页 * |
频高图Es层回波的自动识别;王胜 等;《空间科学学报》;20190731;第39卷(第4期);第432-441页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114660537A (en) | 2022-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10473791B1 (en) | GBAS integrity performance evaluation method based on pseudorange error distribution model | |
CN108508461B (en) | GNSS carrier phase based high-precision positioning integrity monitoring method | |
CN108061889B (en) | AIS and radar angle system deviation correlation method | |
US20080117106A1 (en) | Emitter Geolocation | |
Yoon et al. | Equatorial plasma bubble threat parameterization to support GBAS operations in the Brazilian region | |
CN112666588B (en) | Positioning method based on scene matching and machine learning in urban canyon environment | |
CN110889380B (en) | Ship identification method and device and computer storage medium | |
Naganawa et al. | ADS-B anti-spoofing performance of monopulse technique with sector antennas | |
CN110658494B (en) | Passive unmanned aerial vehicle positioning method and server | |
EP3889648A1 (en) | High-precision satellite positioning method, positioning terminal and positioning system | |
CN115982564A (en) | Method and device for calculating electron density of regional ionized layer and computer equipment | |
CN112083446A (en) | Method and device for positioning deception jamming source | |
CN115685260A (en) | Identification deception jamming identification method, device, platform and medium for GLS reference station | |
CN114235007B (en) | Positioning and integrity monitoring method and system for APNT service | |
CN114660537B (en) | Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment | |
CN111796313B (en) | Satellite positioning method and device, electronic equipment and storage medium | |
CN112198533B (en) | System and method for evaluating integrity of foundation enhancement system under multiple hypotheses | |
Jalloul et al. | DME/DME navigation using a single low-cost SDR and sequential operation | |
CN109490828B (en) | Positioning method based on homologous baseline array | |
CN111624589A (en) | Marine target data fusion system and method based on space-based radio monitoring | |
CN113640760B (en) | Radar discovery probability evaluation method and equipment based on air situation data | |
CN111276000A (en) | Positioning method and system based on vehicle-road cooperation, medium and vehicle-mounted equipment | |
CN110927765A (en) | Laser radar and satellite navigation fused target online positioning method | |
CN113534130B (en) | Multi-station radar multi-target data association method based on sight angle | |
CN105353385A (en) | ARAIM nominal bias evaluation method and device based on Beidou three frequency points |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221115 |