Sky wave propagation time delay correction method for improving ionized layer equivalent height

Abstract

The invention discloses a sky wave propagation time delay correction method for improving the equivalent height of an ionized layer. The method comprises the following steps: converting obtained geographical location information of a long-wave transmitting station and a signal receiver into latitude and longitude information; the latitude and the longitude corresponding to the position information of the sky wave signal ionosphere reflection point are obtained by combining a long wave transmitting station; obtaining an annual day and an hour moment corresponding to the emission time of the long-wave emission station for emitting the time service signal; inputting the solar activity F10.7 index and the geomagnetic activity planet index corresponding to the annual day, the hour moment and the emission time into an ionized layer equivalent height model to obtain an ionized layer equivalent height value; obtaining the sight distance extreme value of the radio signal through the ionosphere equivalent height value, and obtaining the actual signal propagation path length of the sky wave signal through the sight distance extreme value; and further acquiring the propagation time delay of the sky wave signal, and acquiring the complete propagation time of the sky wave signal from the long-wave transmitting station to the signal receiver through the propagation time delay. According to the invention, the propagation time delay precision of sky wave signals and the time service precision of a long-wave time service system are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of time and frequency, and relates to, but is not limited to, a method for correcting the skywave propagation delay of the ionospheric equivalent height. Background Art

[0002] High-precision time service technology is one of the most important space-time infrastructure, which is directly related to national security and social development. At present, time service means are mainly divided into land-based time service and satellite-based time service. Among them, satellite-based time service systems have problems such as being vulnerable to interference, being vulnerable to attacks, and having weak signals, and there are huge security risks. As a kind of land-based time service means, the long-wave time service system has the advantages of large transmission power, high reliability, strong anti-interference ability, etc., and can effectively make up for the defects existing in the satellite-based time service system.

[0003] The skywave propagation path of the long-wave time service system signal refers to the propagation mode in which the time service signal emitted by the antenna is reflected by the ionosphere and then received by the receiver. The skywave signal can effectively expand the application range of the long-wave time service system. The propagation delay of the skywave signal refers to the time experienced by the skywave signal from the transmitting station antenna to the receiving antenna. The propagation delay of the skywave signal is a key parameter that needs to be considered in the monitoring and evaluation of the long-wave time service system. Since the skywave propagation mode relies on the reflection of the signal by the ionosphere, changes in the ionospheric morphology and structure will affect the amplitude and phase of the time service signal, thereby further affecting the skywave propagation delay. Considering the refined changes in the ionosphere helps to improve the determination of the skywave propagation delay, and thus improve the time service accuracy of the long-wave time service system.

[0004] In the related technology, when determining the propagation delay of the skywave signal, the empirical value obtained by looking up the table is used to calculate the ionospheric equivalent height value that affects the total length of the skywave path. On the one hand, this method ignores the time and space variation characteristics of the ionospheric equivalent height itself. On the other hand, solar activity and geomagnetic activity will also affect the ionospheric equivalent height. Therefore, the method of determining the ionospheric equivalent height reduces the accuracy of the skywave signal propagation delay.

[0005] Therefore, how to accurately improve the propagation delay accuracy of the skywave signal and the time service accuracy of the long-wave time service system has become an urgent problem to be solved. Summary of the Invention

[0006] In view of this, an embodiment of the present invention provides a method for correcting the skywave propagation delay by improving the ionospheric equivalent height, which at least solves the problem that the related technology cannot accurately improve the accuracy of the skywave propagation delay and the time service accuracy of the long-wave time service system.

[0007] According to the first aspect of the embodiment of the present invention, a method for correcting the skywave propagation delay by improving the ionospheric equivalent height is provided, including: Obtain the geographical location information of the long-wave transmitting station and the signal receiver respectively, and convert the geographical location information into longitude and latitude information respectively; Obtain the position information of the ionospheric reflection point of the skywave signal through the long-wave transmitting station and the longitude and latitude information, and obtain the latitude of the reflection point and the longitude of the reflection point based on the position information; Obtain the transmission time when the long-wave transmitting station transmits the time signal, and obtain the corresponding day of the year and hour based on the transmission time; Input the obtained parameters into the ionospheric equivalent height model to obtain the ionospheric equivalent height value; the parameters include the latitude, the longitude, the day of the year, the hour, the solar activity F10.7 index corresponding to the obtained transmission time, and the geomagnetic activity planetary index; Obtain the line-of-sight extreme value of the radio signal based on the ionospheric equivalent height value, and obtain the actual signal propagation path length of the skywave signal based on the line-of-sight extreme value; Obtain the propagation delay of the skywave signal based on the signal propagation path length, and obtain the complete propagation time of the skywave signal from the long-wave transmitting station to the signal receiver based on the propagation delay.

[0008] According to the second aspect of the embodiments of the present invention, an electronic device is provided, including: a processor, a memory, a communication interface, and a communication bus. The processor, the memory, and the communication interface complete communication with each other through the communication bus; the memory is used to store at least one executable instruction, and the executable instruction causes the processor to execute the operations corresponding to the method described in the first aspect.

[0009] According to the third aspect of the embodiments of the present invention, a computer storage medium is provided, on which a computer program is stored, and when the program is executed by a processor, the method described in the first aspect is implemented.

[0010] According to the solution provided by the embodiments of the present invention, the geographical location information of the long-wave transmitting station and the signal receiver is obtained respectively, and the geographical location information is respectively converted into longitude and latitude information; the position information of the ionospheric reflection point of the skywave signal is obtained through the long-wave transmitting station and the longitude and latitude information, and the latitude of the reflection point and the longitude of the reflection point are obtained based on the position information; the emission time of the timing signal emitted by the long-wave transmitting station is obtained, and the corresponding day of the year and hour moment are obtained based on the emission time; the obtained parameters are input into the ionospheric equivalent height model to obtain the ionospheric equivalent height value; the parameters include the latitude, the longitude, the day of the year, the hour moment, the solar activity F10.7 index corresponding to the obtained emission time, and the geomagnetic activity planetary index; the line-of-sight extreme value of the radio signal is obtained based on the ionospheric equivalent height value, and the actual signal propagation path length of the skywave signal is obtained based on the line-of-sight extreme value; the propagation delay of the skywave signal is obtained based on the signal propagation path length, and the complete propagation time of the skywave signal from the long-wave transmitting station to the signal receiver is obtained based on the propagation delay. In this process, various parameters such as the day of the year, the hour moment, the solar activity F10.7 index corresponding to the obtained emission time, and the geomagnetic activity planetary index are calculated and obtained, fully considering the variation characteristics of the ionospheric equivalent height in time and space, and combining the physical indices of solar activity and geomagnetic activity to perform spatio-temporal refinement modeling on the ionospheric equivalent height that affects the propagation delay of the skywave signal, and using the obtained model to calculate the equivalent height at the ionospheric reflection point of the skywave signal, so as to improve the timing accuracy of the long-wave timing system. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings, where: Figure 1 It is a schematic flow chart of a skywave propagation delay correction method for improving the ionospheric equivalent height provided by the embodiments of the present invention; Figure 2 It is a schematic diagram of the effect of an ionospheric equivalent height model provided by the embodiments of the present invention; Figure 3 It is a schematic diagram of the effect of the propagation path of the timing signal forming a skywave signal through the ionosphere provided by the embodiments of the present invention; Figure 4 It is a schematic diagram of the effect of the propagation time from the timing signal to the signal receiver provided by the embodiments of the present invention; Figure 5 It is a schematic diagram of the structure of an electronic device provided by the embodiments of the present invention. Detailed implementation manners

[0012] To make the objectives, 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 accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, rather than all, of the embodiments of the present invention. The following embodiments are used to illustrate the present invention, but are not intended to limit the scope of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0013] In the following description, reference is made to "some embodiments", which describe subsets of all possible embodiments. However, it can be understood that "some embodiments" can be the same or different subsets of all possible embodiments and can be combined with each other without conflict.

[0014] It should be noted that the terms "first", "second", and "third" involved in the embodiments of the present invention are only used to distinguish similar objects and do not represent a specific order for the objects. It can be understood that "first", "second", and "third" can be interchanged with a specific order or sequence when allowed, so that the embodiments of the present invention described herein can be implemented in an order other than that illustrated or described herein.

[0015] Those skilled in the art of the present technology can understand that, unless otherwise defined, all terms (including technical terms and scientific terms) used herein have the same meaning as that generally understood by those of ordinary skill in the field to which the embodiments of the present invention belong. It should also be understood that terms such as those defined in a general dictionary should be understood to have a meaning consistent with the meaning in the context of the prior art, and will not be interpreted in an idealized or overly formal sense unless specifically defined as herein.

[0016] Figure 1 As shown in the flowchart of a method for correcting the skywave propagation delay by improving the ionospheric equivalent height provided for the embodiments of the present invention, the method for correcting the skywave propagation delay by improving the ionospheric equivalent height provided for the embodiments of the present invention can be executed by an electronic device, and the electronic device can be, for example, a computer, a server, or the like.

[0017] As Figure 1 shown, the method for correcting the skywave propagation delay by improving the ionospheric equivalent height includes: S101. Obtain the geographical location information of the longwave transmitting station and the signal receiver respectively, and convert the geographical location information into longitude and latitude information respectively.

[0018] In an embodiment of the present invention, the geographical location information includes first geographical location information and second geographical location information, and the longitude and latitude information includes first longitude and latitude information and second longitude and latitude information. The first geographical location information of the long-wave transmitting station and the second geographical location information of the signal receiver can be obtained respectively through technical detection means or map services and the like. Then, the first geographical location information is converted into first longitude and latitude information, and the second geographical location information is converted into second longitude and latitude information.

[0019] S102. Obtain the position information of the ionospheric reflection point of the sky wave signal through the long-wave transmitting station and the longitude and latitude information, and obtain the latitude of the reflection point and the longitude of the reflection point based on the position information.

[0020] In an embodiment of the present invention, after the long-wave transmitting station transmits the time service signal, it propagates in the form of an electromagnetic wave. After encountering the ionosphere, it is reflected to form a sky wave signal. The ionospheric reflection point of the sky wave signal is the specific position where the time service signal encounters the ionosphere to form a sky wave signal during the propagation process. The position information of the ionospheric reflection point of the sky wave signal can be calculated through the long-wave transmitting station, the first longitude and latitude information, and the second longitude and latitude information. Based on the position information, the latitude information of the reflection point and the longitude information of the reflection point are further obtained.

[0021] S103. Obtain the transmission time when the long-wave transmitting station transmits the time service signal, and obtain the corresponding day of the year and hour of the day based on the transmission time.

[0022] In an embodiment of the present invention, the transmission time when the long-wave transmitting station transmits the time service signal can be obtained through the official website of the long-wave time service station, or a specially designed long-wave receiver is used to receive the long-wave time service signal, and the transmission time and characteristics of the time service signal are observed. The real-time transmission time information can be directly obtained, and the day of the year and hour of the day corresponding to the transmission time are recorded. Among them, the day of the year and hour of the day corresponding to the transmission time can be calculated according to the specific transmission date and time, and the hour of the day is directly recorded as the specific time of transmission.

[0023] S104. Input the obtained parameters into the ionospheric equivalent height model to obtain the equivalent ionospheric equivalent height value; the parameters include latitude, longitude, day of the year, hour of the day, the F10.7 index of solar activity corresponding to the obtained occurrence time, and the planetary index of geomagnetic activity.

[0024] In an embodiment of the present invention, the ionospheric equivalent height model includes an input layer, three hidden layers, and an output layer. The nodes of each hidden layer are set according to (32, 64, 128). The latitude, longitude, day of the year, hour of the day, the F10.7 index of solar activity corresponding to the obtained occurrence time, and the planetary index of geomagnetic activity are input into the input layer and the three hidden layers of the ionospheric equivalent height model for processing, and the output layer outputs the equivalent ionospheric equivalent height value.

[0025] Furthermore, the ionospheric equivalent height model can be expressed by the following formula: ; In the formula, is the ionospheric equivalent height value, is the latitude of the reflection point of the skywave signal, is the longitude of the reflection point of the skywave signal, is the day of the year when the skywave signal is located, is the approximate hour when the skywave signal is located, F10.7 is the solar activity index, and Kp is the geomagnetic activity index value.

[0026] Among them, the model function can be illustrated by the following calculation process: ; In the formula, is the output of the hidden layer, is the model weight parameter, is the bias parameter, is the input parameter, is the non-linear activation function, is the number of hidden layers, is the model output result, that is, the ionospheric equivalent height value .

[0027] As Figure 2 shown, Figure 2 is a schematic diagram of the effect of an ionospheric equivalent height model provided by an embodiment of the present invention. In Figure 2 , the ionospheric equivalent height model includes an input layer, a hidden layer, and an output layer. The input layer inputs parameters such as the longitude of the reflection point, the latitude of the reflection point, the day of the year DOY, the hour HR, the solar activity F10.7 index corresponding to the obtained occurrence time, and the geomagnetic activity planetary index Kp. After being processed by the hidden layer and the output layer, the ionospheric equivalent height value is obtained.

[0028] S105. Obtain the line-of-sight extreme value of the radio signal based on the ionospheric equivalent height value, and obtain the actual signal propagation path length of the skywave signal based on the line-of-sight extreme value.

[0029] In an embodiment of the present invention, the timing signal transmitted by the long-wave transmitting station reaches the ground receiver after one reflection from the ionosphere. The actual propagation path of the skywave signal is the superposition of the paths before and after the signal emission point, that is, the timing signal includes not only the ascending section from the emission point to the ionosphere, but also the descending section from the ionosphere back to the ground after the skywave signal is formed. Therefore, the actual signal propagation path length of the skywave signal can be achieved through this step. First, substitute the equivalent height value of the ionosphere into the line-of-sight extreme value formula to obtain the line-of-sight extreme value, where the line-of-sight extreme value formula is as follows: ; In the formula, is the equivalent earth radius, calculated using , , is the equivalent height of the ionosphere, is the angle formed by the connection lines between the long-wave transmitting station and the signal receiver and the center of the earth, is the line-of-sight extreme value of the radio signal. The line-of-sight extreme value of the radio signal is the great circle distance between the transmitting station and the signal receiver when the angle formed by the skywave signal and the earth's tangent plane is zero.

[0030] Furthermore, substitute the line-of-sight extreme value into the signal propagation path length formula to obtain the actual signal propagation path length of the skywave signal. Among them, the signal propagation path length formula is as follows: ; In the above formula, d s is the signal propagation path length.

[0031] In an embodiment of the present invention, the extreme value of the skywave signal propagation path corresponding to the line-of-sight extreme value of the radio signal can be calculated by the following formula: ; It can be seen from the extreme value of the skywave signal propagation path that when the emission angle of the skywave signal is not zero, the actual skywave signal propagation path length is less than the extreme value of the skywave signal propagation path. When the calculated actual skywave signal propagation path length is not less than the extreme value of the skywave signal propagation path, recalculate the actual skywave signal propagation path length.

[0032] Figure 3 is the effect schematic diagram of the propagation path of the timing signal passing through the ionosphere to form the skywave signal provided by the embodiment of the present invention. In Figure 3In the process that the time service signal transmitted by a long-wave transmitting station (such as the BPL long-wave time service station) reaches the signal receiver on the ground after being reflected by the ionosphere, after the time service signal is emitted from the long-wave transmitting station, it first propagates towards the sky in an almost straight line. The time service signal passes through the atmosphere and enters the ionosphere. In the ionosphere, due to the existence of electrons and ions, the time service signal will undergo refraction and partial reflection to form a skywave signal. The skywave signal reflected by the ionosphere passes through the atmosphere again and returns to the Earth's surface, and is finally received by the ground receiver (signal receiver). In the whole process, there are multiple parameters, which are respectively the signal propagation path length d s , the latitude of the reflection point , longitude , the equivalent height value h of the ionosphere, the great circle distance (the extreme value of the radio signal line-of-sight), the equivalent Earth radius and the included angle formed by the connection lines between the long-wave transmitting station and the signal receiver and the center of the Earth .

[0033] S106. Obtain the propagation delay of the skywave signal based on the signal propagation path length, and obtain the complete propagation time of the skywave signal from the long-wave transmitting station to the signal receiver based on the propagation delay of the skywave signal.

[0034] In the embodiment of the present invention, the correct propagation delay of the skywave signal is calculated based on the signal propagation path length, and the complete time is further corrected through the calculated propagation delay of the skywave signal, that is, the complete propagation time of the skywave signal from the long-wave transmitting station to the signal receiver is calculated according to the propagation delay of the skywave signal.

[0035] It can be understood that in the embodiments of the present invention, the geographical location information of the long-wave transmitting station and the signal receiver is respectively obtained, and the geographical location information is respectively converted into longitude and latitude information; the position information of the ionospheric reflection point of the skywave signal is obtained through the long-wave transmitting station and the longitude and latitude information, and the latitude of the reflection point and the longitude of the reflection point are obtained based on the position information; the transmission time of the timing signal transmitted by the long-wave transmitting station is obtained, and the corresponding day of the year and hour are obtained based on the transmission time; the obtained parameters are input into the ionospheric equivalent height model to obtain the ionospheric equivalent height value; the parameters include latitude, longitude, day of the year, hour, the solar activity F10.7 index corresponding to the obtained transmission time, and the geomagnetic activity planetary index; the line-of-sight extreme value of the radio signal is obtained based on the ionospheric equivalent height value, and the actual signal propagation path length of the skywave signal is obtained based on the line-of-sight extreme value; the propagation delay of the skywave signal is obtained based on the signal propagation path length, and the complete propagation time of the skywave signal from the long-wave transmitting station to the signal receiver is obtained based on the propagation delay. In this process, various parameters such as the day of the year, hour, the solar activity F10.7 index corresponding to the obtained transmission time, and the geomagnetic activity planetary index are calculated and obtained, fully considering the variation characteristics of the ionospheric equivalent height in time and space, and combining the physical indices of solar activity and geomagnetic activity to perform spatio-temporal refinement modeling on the ionospheric equivalent height that affects the propagation delay of the skywave signal, and using the obtained model to calculate the equivalent height at the ionospheric reflection point of the skywave signal, so as to improve the timing accuracy of the long-wave timing system.

[0036] In the embodiments of the present invention, obtaining the propagation delay of the skywave signal based on the signal propagation path length in S106 can be implemented through S106A, and the following steps are used for illustration.

[0037] S106A. Obtain the propagation delay according to the signal propagation path length and the propagation delay formula.

[0038] In the embodiments of the present invention, the propagation delay formula is as follows: ; In the above formula, is the propagation delay, are respectively the phase correction coefficients of the skywave signal, which can be calculated by an empirical model, is the speed of light.

[0039] Among them, the signal propagation path length can be substituted into the above propagation delay formula to obtain the propagation delay.

[0040] In the embodiments of the present invention, obtaining the complete propagation time of the skywave signal from the long-wave transmitting station to the signal receiver based on the propagation delay in S106 can be implemented through S1061 to S1062, and the following steps are used for illustration.

[0041] S1061. Obtain the first delay time when the skywave signal arrives at the receiving antenna, the second delay time when the skywave signal is transmitted through the cable, the third delay time when the skywave signal is transmitted inside the signal receiver, and the fourth delay time of an additional tracking point, respectively.

[0042] In some embodiments of the present invention, the first delay time when the skywave signal arrives at the receiving antenna, that is, the delay generated when the skywave signal first touches the receiving antenna, can be achieved through measurement. The second delay time when the skywave signal is transmitted through the cable, that is, the delay generated by the skywave signal passing through the cable during the process of transmitting from the receiving antenna to the signal receiver, can be calculated according to the cable length and characteristics. The third delay time when the skywave signal is transmitted inside the signal receiver, that is, the time consumed by the skywave signal after entering the receiver through various processing processes (such as filtering, amplification, down-conversion, etc.), can be obtained by referring to relevant documents. And the fourth delay time of an additional tracking point, that is, the 30-microsecond delay of the tracking point.

[0043] S1062. Determine the complete propagation time based on the first delay time, the second delay time, the third delay time, the fourth delay time, and the propagation delay of the skywave signal.

[0044] In some embodiments of the present invention, substitute the first delay time, the second delay time, the third delay time, the fourth delay time, and the propagation delay of the skywave signal into the summation formula for summation to obtain the complete propagation time. The summation formula is as follows: ; In the above formula, TOA is the complete propagation time, T 1 is the propagation delay of the skywave signal, and T 2 is the sum of the delays of the first delay time, the second delay time, the third delay time, and the fourth delay time.

[0045] As Figure 4 shown, Figure 4 is a schematic diagram of the effect of the propagation time of the timing signal to the receiver reception provided by the embodiment of the present invention. In Figure 4 , the propagation delay T 1 is calculated during the process from the moment when the timing signal is transmitted by the long-wave transmitting station to the moment when the skywave signal formed by the reflection of the timing signal by the ionosphere reaches the signal receiver to receive the skywave signal. The internal delay T 2 (sum of delays) is calculated during the process from the moment when the signal receiver receives the skywave signal to the moment when the signal receiver locks the skywave signal. Finally, the complete propagation time TOA is calculated through the propagation delay T 1 and the internal delay T 2 .

[0046] It can be understood that in some embodiments of the present invention, the first delay time for the skywave signal to reach the receiving antenna, the second delay time when the skywave signal is transmitted through the cable, the third delay time for the skywave signal to be transmitted inside the signal receiver, and the fourth delay time for an additional tracking point are respectively obtained, and the first delay time, the second delay time, the third delay time, the fourth delay time, and the skywave signal propagation delay are used to determine the complete propagation time. In this process, by calculating and correcting the total delay of the skywave signal through multiple delay times, the time synchronization accuracy can be significantly improved, the communication reliability can be enhanced, the system performance can be optimized, complex application scenarios can be supported, and maintenance and debugging are facilitated.

[0047] Referring to Figure 5 , a schematic structural diagram of an electronic device according to an embodiment of the present invention is shown. The specific implementation of the electronic device is not limited in the specific embodiments of the present invention.

[0048] As Figure 5 shown, the electronic device may include: a processor 502, a communication interface 504, a memory 506, and a communication bus 508.

[0049] Wherein: The processor 502, the communication interface 504, and the memory 506 communicate with each other through the communication bus 508.

[0050] The communication interface 504 is used to communicate with other electronic devices or servers.

[0051] The processor 502 is used to execute the program 510, and specifically can execute the relevant steps in the above method embodiments.

[0052] Specifically, the program 510 may include program code, and the program code includes computer operation instructions.

[0053] The processor 502 may be a central processing unit CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present invention. One or more processors included in the intelligent device may be of the same type of processor, such as one or more CPUs; or may be of different types of processors, such as one or more CPUs and one or more ASICs.

[0054] The memory 506 is used to store the program 510. The memory 506 may include a high-speed RAM memory, and may also include a non-volatile memory, such as at least one disk memory.

[0055] The program 510 can specifically be used to cause the processor 502 to perform the operations corresponding to the methods described in the above method embodiments.

[0056] For the specific implementation of each step in the program 510, reference may be made to the corresponding descriptions in the corresponding steps and units in the above method embodiments, which will not be elaborated here. Those skilled in the art can clearly understand that for the convenience and conciseness of description, the specific working processes of the devices and modules described above may refer to the corresponding process descriptions in the foregoing method embodiments, which will not be repeated here.

[0057] It should be noted that according to the needs of implementation, each component / step described in the embodiments of the present invention can be split into more components / steps, or two or more components / steps or partial operations of components / steps can be combined into new components / steps to achieve the purpose of the embodiments of the present invention.

[0058] The method according to the embodiments of the present invention described above can be implemented in hardware, firmware, or be implemented as software or computer code that can be stored in a recording medium (such as a CD ROM, RAM, floppy disk, hard disk, or magneto-optical disk), or be implemented as computer code originally stored in a remote recording medium or a non-transitory machine-readable medium and downloaded through a network and to be stored in a local recording medium, so that the method described herein can be stored on such a software process on a recording medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware (such as an ASIC or FPGA). It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component (such as a RAM, a ROM, a flash memory, etc.) that can store or receive software or computer code, and when the software or computer code is accessed and executed by the computer, the processor, or the hardware, the method described herein is implemented. In addition, when a general-purpose computer accesses the code for implementing the method shown herein, the execution of the code converts the general-purpose computer into a dedicated computer for implementing the method shown herein.

[0059] Those of ordinary skill in the art can realize that the units and method steps of each example described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Skilled professionals can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the embodiments of the present invention.

[0060] The above embodiments are only used to illustrate the embodiments of the present invention, rather than to limit the embodiments of the present invention. Those of ordinary skill in the relevant technical field can also make various changes and modifications without departing from the spirit and scope of the embodiments of the present invention. Therefore, all equivalent technical solutions also belong to the scope of the embodiments of the present invention. The patent protection scope of the embodiments of the present invention shall be defined by the claims.

Claims

1. A method for improving the sky wave propagation delay correction of the ionospheric equivalent height, characterized in that: include: Respectively obtaining geographical location information of the long-wave transmitter and the signal receiver, and respectively converting the geographical location information into longitude and latitude information; Acquire the position information of the ionospheric reflection point of the sky wave signal through the long-wave transmitting station and the longitude and latitude information, and acquire the latitude and longitude of the reflection point based on the position information; Obtaining the transmission time of the timing signal transmitted by the long-wave transmitting station, and obtaining the corresponding annual day and hour based on the transmission time; Inputting the acquired parameters into an ionosphere equivalent height model to obtain an ionosphere equivalent height value; the parameters include the latitude, the longitude, the annual cumulative day, the hour, the solar activity F10.7 index corresponding to the acquired launch time, and the geomagnetic activity planet index; Acquire a line-of-sight extreme value of a radio signal based on the ionosphere equivalent height value, and acquire an actual signal propagation path length of a sky wave signal based on the line-of-sight extreme value; The propagation delay of the sky wave signal is obtained based on the length of the signal propagation path, and the complete propagation time of the sky wave signal from the long wave transmitting station to the signal receiver is obtained based on the propagation delay.

2. The method according to claim 1, characterized in that The ionospheric equivalent height model includes an input layer, three hidden layers and an output layer, and the nodes of each hidden layer are set according to (32, 64, 128).

3. The method according to claim 1, characterized in that The obtaining of the line-of-sight extreme value of the radio signal based on the ionosphere equivalent height value comprises: The line-of-sight extreme value is obtained based on the ionosphere equivalent height value and the line-of-sight extreme value formula, wherein the line-of-sight extreme value formula is as follows: ; In the above formula, is the equivalent earth radius, calculate, , is the equivalent ionosphere height value, is the angle formed by the line connecting the long-wave transmitting station, the signal receiver and the center of the earth, is the extreme value of the line of sight, which is the great circle distance between the long wave transmitting station and the signal receiver when the angle formed by the sky wave signal and the tangent plane of the earth is zero.

4. The method according to claim 3, characterized in that The obtaining the actual signal propagation path length of the sky wave signal based on the line-of-sight extreme value comprises: The signal propagation path length is obtained based on the line-of-sight extreme value and the signal propagation path length formula, wherein the signal propagation path length formula is as follows: ; In the above formula, d s is the signal propagation path length.

5. The method according to claim 4, characterized in that The acquiring the propagation delay of the sky wave signal based on the signal propagation path length includes: The propagation delay is obtained according to the signal propagation path length and the propagation delay formula, wherein the propagation delay formula is as follows: ; In the above formula, is the propagation delay, are the phase correction coefficients of the sky wave signal, The speed of light.

6. The method according to any one of claims 1 to 5, characterized in that: The acquiring the complete propagation time of the sky wave signal from the long wave transmitting station to the signal receiver based on the propagation delay comprises: Respectively acquiring a first delay time of the sky wave signal reaching a receiving antenna, a second delay time of the sky wave signal when being transmitted via a cable, a third delay time of the sky wave signal being transmitted inside the signal receiver, and a fourth delay time of an additional tracking point; The first delay time, the second delay time, the third delay time, the fourth delay time and the sky wave signal propagation delay are summed to obtain the complete propagation time.

Images