CN114660537A - Method and equipment for estimating sky wave pitch angle of ionized layer F layer based on Es layer - Google Patents
Method and equipment for estimating sky wave pitch angle of ionized layer F layer based on Es layer Download PDFInfo
- Publication number
- CN114660537A CN114660537A CN202210396710.5A CN202210396710A CN114660537A CN 114660537 A CN114660537 A CN 114660537A CN 202210396710 A CN202210396710 A CN 202210396710A CN 114660537 A CN114660537 A CN 114660537A
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 239000005433 ionosphere Substances 0.000 claims description 29
- 238000010586 diagram Methods 0.000 claims description 28
- 238000002474 experimental method Methods 0.000 claims description 10
- 238000001514 detection 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 sky wave pitch angle estimation method and equipment for an ionized layer F layer 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 ionospheric sky-wave pitch angles, in particular to a method and equipment for estimating an ionospheric F-layer sky-wave pitch angle based on an Es layer.
Background
The ionosphere is a high-level atmospheric region partially ionized at about 60km to about 1000km from the ground, and can be divided into 4 different regions, namely a layer D, a layer E, a layer F1, and a layer F2 (collectively referred to as layer F) 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) strong ionosphere E layers 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 a speed 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, it is an urgent technical problem in the art to develop a method and a device for estimating the sky wave pitch angle of the ionosphere F layer based on the Es layer, which can effectively overcome the above-mentioned defects in the related art.
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; 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 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, thetai minRepresents the minimum value of the pitch angle of the echo corresponding to the frequency point fi, thetai maxRepresenting the maximum value of the echo pitch angle P of the frequency point fiiEs maxMaximum value of Es layer echo group path, P, representing frequency point fiiEs minRepresents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, hEsIndicating 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:
hvi F max=sin(θi max)·Pi F min/2
hvi F min=sin(θi min)·Pi F max/2
wherein, Pi F minMinimum value, P, of F-layer echo group path representing frequency point fii F maxMaximum value, hv, of the F-layer echo group path representing frequency point fii F minRepresents the minimum of the F-layer echo imaginary height at frequency point fi,hvi F maxand 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, 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 includes that the relational expression of the group path and the virtual height of the one-frequency point F layer in step 2 includes:
wherein, PijGroup path, hv, of the j-th echo of the F layer representing frequency point fiijAnd 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, thetaijAnd 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 of actual and estimated pitch angles 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 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, thetai minRepresents the minimum value of the pitch angle of the echo corresponding to the frequency point fi, thetai maxRepresenting the maximum value of the echo pitch angle P of the frequency point fiiEs maxMaximum value of Es layer echo group path, P, representing frequency point fiiEs minRepresents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, hEsIndicating 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 in the experiment is 113 degrees E, the latitude of the transmitting station in the experiment is 30 degrees N, firstly, an IRI model is utilized to generate a background electron density map of which the distance from 14 points 12, 25 and 25 in 2019 to the azimuth of the transmitting station is 180 degrees, an Es layer with the range of 600km is added at the position with the height of 110km, and the adjacent frequency corresponding to the electron density of the Es layer is 9 MHz. 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, thetai minThe minimum value of the echo pitch angle corresponding to the frequency point fi is represented as 13.83 degrees, thetai max18.79 degrees represents the maximum value of the echo pitch angle of the frequency point fi, PiEs max920.32km represents the Es layer echo group path maximum, P, of frequency point fiiEs min682.89km represents the minimum of the Es layer echo group path of the frequency point fi, R6370 km represents the radius of the earth, hEs110km indicates the Es layer height.
Based on the content of the above method embodiment, as an optional embodiment, 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 obtains the maximum value of the imaginary height of the F-layer echo at one frequency point through the maximum value of the pitch angle at one frequency point and the minimum value of the F-layer group path in step 2, and obtains the minimum value of the imaginary height of the F-layer echo at one frequency point through the minimum value of the pitch angle at one frequency point and the maximum value of the F-layer group path, including:
wherein, Pi F minMinimum value, P, of F-layer echo group path representing frequency point fii F maxMaximum value, hv, of the F-layer echo group path representing frequency point fii F minMinimum value of F-layer echo virtual height, hv, representing frequency point fii F maxAnd 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 isi F min1243.73km represents the minimum value of the F-layer echo group path at frequency point 15MHz, Pi F max1360.41km represents the maximum value of the F-layer echo group path at frequency 15MHz, hvi F min162.6km represents the minimum value of the F-layer echo imaginary height at the frequency point of 15MHz, hvi F maxThe maximum value of the F-layer echo imaginary height at the frequency point of 15MHz is represented by 200.3 km.
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, PijGroup path, hv, of the j-th echo of the F layer representing frequency point fiijAnd 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 isijRepresenting the group path, hv, of the j-th echo of the F layer at the frequency point of 15MHzijAnd the virtual height corresponding to the j-th echo of the F layer with 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, thetaijAnd the pitch angle of the jth echo of the F layer at the frequency point fi is represented.
Specifically, for example, the 3 rd echo group path Pi3 of the 15MHz echo is 1290.70km, the corresponding virtual height hvi3 is 185.1km, the estimated value θ i3 is 16.03 ° according to the pitch angle calculation formula, and substantially coincides with the actual pitch angle 16.00 °, and the comparison between the actual pitch angle at 15MHz and the estimated pitch angle 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 the above practical situation, on the basis of the foregoing embodiments, embodiments of the present invention provide an ionospheric F sky-wave pitch angle estimation apparatus based on the Es layer, where the apparatus is configured to execute the method for estimating the ionospheric F sky-wave pitch angle 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 detector 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 an 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 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; 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, thetai minRepresenting the minimum value of the pitch angle of the echo corresponding to the frequency point fi, theta imax represents the maximum value of the pitch angle of the echo of the frequency point fi, PiEs maxMaximum value of Es layer echo group path, P, representing frequency point fiiEs minRepresents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, hEsIndicating 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:
hvi F max=sin(θi max)·Pi F min/2
hvi F min=sin(θi min)·Pi F max/2
wherein, Pi F minMinimum value, P, of F-layer echo group path representing frequency point fii F maxMaximum value, hv, of the F-layer echo group path representing frequency point fii F minMinimum value of F-layer echo virtual height, hv, representing frequency point fii F maxAnd the maximum value of the F-layer echo virtual height of the frequency point fi is represented.
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: 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, PijGroup path, hv, of the j-th echo of the F layer representing frequency point fiijAnd 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 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 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, thetaijAnd 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. With this object in mind, an embodiment of the present invention provides an electronic device, 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 other identical elements 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, but 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; 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 F-layer group path and the virtual height 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, thetaiminRepresents the minimum value of the pitch angle of the echo corresponding to the frequency point fi, thetaimaxRepresenting the maximum value of the echo pitch angle P of the frequency point fiiEsmaxMaximum value of Es layer echo group path, P, representing frequency point fiiEsminRepresents the minimum value of the Es layer echo group path of the frequency point fi, R represents the earth radius, hEsIndicating 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:
hviFmax=sin(θimax)·PiFmin/2
hviFmin=sin(θimin)·PiFmax/2
wherein, PiFminMinimum value, P, of F-layer echo group path representing frequency point fiiFmaxMaximum value, hv, of the F-layer echo group path representing frequency point fiiFminMinimum value of F-layer echo virtual height, hv, representing frequency point fiiFmaxAnd 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, PijGroup path, hv, of the j-th echo of the F layer representing frequency point fiijAnd 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 an imaginary height of an echo is obtained according to a relation between the group path of the F layer and the imaginary height in step 2, and a pitch angle of an echo is obtained by using a cosine formula, the method comprising:
wherein, thetaijAnd the pitch angle of the jth echo of the F layer at the frequency point fi is represented.
6. An ionosphere F layer sky wave pitch angle estimation device based on an Es layer is characterized by comprising: 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.
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 true CN114660537A (en) | 2022-06-24 |
CN114660537B 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 |
---|
CHUNHUA JIANG ET AL.: "Ionosonde observations of daytime spread F at low latitudes", 《JOURNAL OF GEOPHYSICAL RESEARCH:SPACE PHYSICS》 * |
CHUNHUA JIANG ET AL.: "Large-Scale Ionospheric Irregularities Detected by Ionosonde and GNSS Receiver Network", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
张雅彬 等: "电离层Es电波传播模型研究", 《空间科学学报》 * |
王胜 等: "频高图Es层回波的自动识别", 《空间科学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114660537B (en) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10473791B1 (en) | GBAS integrity performance evaluation method based on pseudorange error distribution model | |
CN108061889B (en) | AIS and radar angle system deviation correlation method | |
CN101971047B (en) | Device and method for the real-time monitoring of the integrity of a satellite navigation system | |
Yoon et al. | Equatorial plasma bubble threat parameterization to support GBAS operations in the Brazilian region | |
CN110889380B (en) | Ship identification method and device and computer storage medium | |
EP1828799A1 (en) | Emitter geolocation | |
CN112666588B (en) | Positioning method based on scene matching and machine learning in urban canyon environment | |
Naganawa et al. | ADS-B anti-spoofing performance of monopulse technique with sector antennas | |
EP3889648A1 (en) | High-precision satellite positioning method, positioning terminal and positioning system | |
CN103529437A (en) | Method used for captive-balloon-borne phased array radar to distinguish open space targets under multi-target condition | |
CN114488230B (en) | Doppler positioning method and device, electronic equipment and storage medium | |
CN110658494A (en) | Passive unmanned aerial vehicle positioning method and server | |
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 | |
Margaria et al. | The need for GNSS position integrity and authentication in ITS: Conceptual and practical limitations in urban contexts | |
CN111796313B (en) | Satellite positioning method and device, electronic equipment and storage medium | |
CN104101869A (en) | Ground wave radar moving object simulation modeling method in polar coordinate | |
CN114660537B (en) | Es layer-based ionosphere F layer sky wave pitch angle estimation method and equipment | |
US11428802B2 (en) | Localization using particle filtering and image registration of radar against elevation datasets | |
CN111624589B (en) | Marine target data fusion system and method based on space-based radio monitoring | |
CN115728797B (en) | GNSS multi-system ionosphere disturbance analysis method for tsunami warning | |
CN111276000A (en) | Positioning method and system based on vehicle-road cooperation, medium and vehicle-mounted equipment | |
CN105353385A (en) | ARAIM nominal bias evaluation method and device based on Beidou three frequency points | |
CN113640760A (en) | Radar discovery probability evaluation method and equipment based on air situation data | |
Louart et al. | Detection of AIS messages falsifications and spoofing by checking messages compliance with TDMA protocol |
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 |