GNSS amplitude flicker index construction method and system

Abstract

The invention belongs to the field of satellite positioning technology application, and particularly discloses a GNSS amplitude flicker index construction method and system. The construction method comprises the following steps: acquiring observation data of an ionospheric scintillation monitoring receiver, extracting an in-phase component I and an orthogonal phase component Q in the data according to the observation data, and calculating GNSS signal intensity according to the in-phase component I and the orthogonal phase component Q; carrying out detrending processing on the GNSS signal intensity by adopting the average value; calculating an amplitude flicker index # imgabs0 # and a signal attenuation speed # imgabs1 # according to the strength of the GNSS signal after the de-trending processing; and calculating the amplitude flicker index # imgabs3 # of the ionized layer by using a cosine function of # imgabs2 #. According to the method, the amplitude flicker index # imgabs4 # and the signal attenuation speed # imgabs5 # are combined to construct the amplitude flicker index # imgabs6 # considering the satellite signal attenuation speed and amplitude, the GNSS amplitude flicker condition can be accurately reflected, and the method has good practicability.

Description

Technical Field

[0001] This application belongs to the field of satellite positioning technology applications, and more specifically, relates to a method and system for constructing a GNSS amplitude scintillation index. Background Art

[0002] Ionospheric scintillation refers to the physical phenomenon in which the amplitude and phase of radio signals passing through ionospheric irregularities change rapidly, randomly, and fluctuating. Ionospheric scintillation can reduce the quality of Global Navigation Satellite System (GNSS) observations, and can even cause satellite signal interruption, which has become one of the difficult problems that need to be solved urgently in current GNSS precise positioning applications. Currently, the sun is in the high year of the 25th solar activity cycle, and ionospheric scintillation events closely related to the solar activity level have increased significantly. How to accurately detect ionospheric scintillation is the key to reliably determining the intensity of ionospheric scintillation.

[0003] The index is currently the most widely used indicator to characterize the severity of the impact of GNSS signals on ionospheric scintillation. Existing research has shown that when > 0.6, the error of GPS L1 carrier phase observations can exceed 0.3 m, while the errors of GNSS precise point positioning (PPP) and real-time differential positioning can reach several meters. In

[0004] Although the index has been widely used to characterize the intensity of ionospheric scintillation, it still has the problem of insufficient availability. Existing scholars have conducted experiments using station data in many regions such as China, Brazil, India, and Kenya, and shown that a strong ionospheric scintillation environment characterized by > 0.6 does not lead to an increase in the dynamic PPP positioning error of GNSS. In fact, the existing index can only reflect the amplitude of signal attenuation when GNSS signals pass through ionospheric irregularities, and cannot characterize the signal attenuation speed. The GNSS signal attenuation speed is related to the depletion degree and drift speed of the electron density in ionospheric irregularities, and can be measured by the index . Compared with the single scintillation index

[0005] Therefore, how to construct an ionospheric scintillation intensity index that takes into account the GNSS signal attenuation speed using GNSS receiver data is a difficult problem that needs to be solved currently. Summary of the invention

[0006] In view of the defects of the prior art, the purpose of this application is to provide a GNSS amplitude scintillation index construction method and system, which can accurately monitor ionospheric scintillation based on different types of GNSS receivers around the world.

[0007] To achieve the above objectives, in a first aspect, the present application provides a method for constructing a GNSS amplitude scintillation index, comprising the following steps: S10, obtaining and extracting the in-phase component in the data according to the observation data of the ionospheric scintillation monitoring receiver I With orthogonal phase component Q , based on which the GNSS signal strength is calculated; S20, detrending the GNSS signal strength using the average value; S30, calculate the amplitude scintillation index based on the GNSS signal strength after detrending and signal attenuation rate ; S40, using The ionospheric amplitude scintillation index is calculated by using the cosine function of .

[0008] The beneficial effects of the present application are as follows: the GNSS amplitude scintillation index construction method provided by the present application utilizes The ionospheric amplitude scintillation index is calculated by using the cosine function of , through the ionospheric amplitude scintillation index To reflect the ionospheric amplitude scintillation, that is When the GNSS signal attenuation rate is slow, reduce The improved weighting of the index makes Numerical and Equivalent, ensuring that the signal decays slowly The accuracy of detecting ionospheric scintillation; when the GNSS signal decays quickly, Increase the pair The index is adjusted to make the satellite signals more seriously affected by ionospheric irregularities Larger, improved The disadvantage of the value being underestimated is to ensure that the signal decays quickly. Detection of ionospheric scintillation is more accurate, effectively improving the accuracy of ionospheric scintillation monitoring.

[0009] As a further preferred embodiment, step S10 is specifically as follows: Convert the observation data and select the required satellite system and corresponding time; Retrieve the in-phase component in the transformed observation data I and the quadrature-phase component Q in the corresponding columns, so as to extract the in-phase component I and the quadrature-phase component Q data; Use the in-phase component I and the quadrature-phase component Q data to calculate the GNSS signal strength , and the calculation formula is: .

[0010] As a further preference, step S20 is specifically: Determine the GNSS signal strength data of the current observation epoch; Determine the GNSS signal strength data for a period of time before and after the current observation epoch; Divide the GNSS signal strength of the current observation epoch by the corresponding signal strength mean value to obtain the detrended GNSS signal strength .

[0011] As a further preference, in step S30, the steps of calculating the amplitude scintillation index are specifically: Calculate the mean value of the square of the detrended GNSS signal strength ; Based on the detrended GNSS signal strength calculate the square of its mean value; According to the mean value of the square of the detrended GNSS signal strength and the square of the mean value of the detrended GNSS signal strength , calculate the amplitude scintillation index , and the calculation formula is:

[0012] where the angular brackets represent taking the mean value for calculation.

[0013] As a further preference, in step S30, the steps of calculating the signal attenuation rate are specifically: According to the detrended GNSS signal strength time series within the time period, calculate its autocorrelation function , and the calculation formula is:

[0014] In the formula, represents Variance; Denotes the expectation operation; Calculate During the time period Autocorrelation function of the time series Drop to Time delay when, that is, the signal attenuation rate , The calculation formula is:

[0015] Wherein, Is The maximum value during the time period.

[0016] Before step S40, it also includes: Statistical analysis of the amplitude scintillation index And the signal attenuation rate Relationship.

[0017] As a further preference, in step S30, the statistical analysis of the amplitude scintillation index And the signal attenuation rate The specific steps are: Use the ionospheric monitoring equipment to determine the real ionospheric scintillation event; Based on a large number of ionospheric scintillation events, statistical analysis , And the distribution of the detrended signal strength change.

[0018] As a further preference, step S40 is specifically: Use the cosine function to adjust the value of the signal attenuation rate So that the adjusted Is between 0 and 1; Based on the statistical analysis method, determine the joint And Construct Fine-tuning factor of ; Calculate the GNSS amplitude scintillation index , The calculation formula is:

[0019] Wherein, Is the data sampling rate; For calculating The time period used.

[0020] In a second aspect, the present application provides a GNSS amplitude scintillation index construction system, including: Signal strength calculation module, used to obtain and extract the in-phase component from the original data of the ionospheric scintillation monitoring receiver according to the ionospheric scintillation monitoring receiver original dataI With the quadrature component Q , and calculate the GNSS signal strength accordingly; A signal strength processing module for detrending the GNSS signal strength using the average value; and A solution module for calculating the amplitude scintillation index and the signal attenuation rate ; An index construction module for calculating the ionospheric amplitude scintillation index using the cosine function of .

[0021] It can be understood that the beneficial effects of the second aspect above can be referred to the relevant descriptions in the first aspect above, and will not be repeated here. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] Figure 1 is a flowchart of the GNSS amplitude scintillation index construction method provided by an embodiment of the present application; Figure 2 is a time series comparison diagram of the signal attenuation rate index calculated by different solutions provided by a specific embodiment of the present application ; where (a) is a schematic diagram of the GNSS detrended signal strength time series, (b) is a time series diagram obtained by the scheme of subtracting the mean value, (c) is a time series diagram obtained by the scheme of not subtracting the mean value, and (d) is a GNSS amplitude scintillation index time series diagram; Figure 3 is a time series diagram of the GNSS signal strength, signal attenuation rate index , and provided by a specific embodiment of the present application; Figure 4 is a comparative analysis diagram of the in-phase component , quadrature component and the GNSS signal strength provided by a specific embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0023] In order to make the objectives, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0024] It should be understood that in the description of the present application, the term "a number of" means at least one, such as one, two, etc., unless otherwise specifically defined; the term "a plurality of" means two or more, unless otherwise specifically defined; the terms "first" and "second" etc. are used to distinguish different objects, rather than to describe a specific order of the objects; the term "and / or" includes any and all combinations of one or more of the related listed items.

[0025] In addition, throughout the description of this specification, the reference to "an embodiment"; the language such as "an embodiment", "an example" or the like means that the specific features, structures or characteristics described in connection with that embodiment are included in at least one embodiment of the present application. Therefore, the appearance of the phrase "in one embodiment;" throughout the specification, "in one embodiment" and similar language may or may not all refer to the same embodiment.

[0026] The main technical problem to be solved by the present application is how to construct a GNSS amplitude scintillation index that takes into account both the attenuation amplitude and speed of satellite signals, so as to accurately monitor ionospheric scintillation based on different types of GNSS receivers globally.

[0027] As Figure 1 shown, the present application provides a method for constructing a GNSS amplitude scintillation index considering the attenuation speed of satellite signals, which is applicable to real-time monitoring of GNSS ionospheric scintillation, including steps S10 to S40, described in detail as follows: Step S10, obtain and extract the in-phase component I and the quadrature component Q from the observation data of the ionospheric scintillation monitoring receiver with high-frequency sampling, and calculate the GNSS signal strength accordingly.

[0028] Specifically, step S10 may include the following sub-steps: Step S11, transform the observation data and select the required satellite system and corresponding time; Step S12, retrieve the column numbers of the in-phase component I and the quadrature component Q in the transformed observation data, so as to extract the data of the in-phase component I and the quadrature component Q ; Step S13, calculate the GNSS signal strength I using the data of the in-phase component Q and the quadrature component , and the calculation formula is: .

[0029] Step S20, use the average value method for the GNSS signal strength Perform detrending to eliminate the long-term change trend and extract short-term fluctuations.

[0030] It should be noted that traditional signal detrending methods usually use a low-pass Butterworth filter for processing. Among them, the filter order is usually set to 6, and the cut-off frequency is set to 0.1. However, when using the Butterworth filter to detrend the signal strength, its effect is affected by the cut-off frequency. In contrast, the average value method adopted in this application is not affected by the cut-off frequency, and the program implementation is very simple.

[0031] Specifically, step S20 may include the following sub-steps: Step S21, determine the GNSS signal strength data of the current observation epoch; Step S22, determine the GNSS signal strength data for a period of time (usually 30 seconds) before and after the current observation epoch; Step S23, divide the GNSS signal strength of the current observation epoch by the corresponding signal strength mean value, so as to eliminate the long-term change trend of the signal strength and extract the short-term fluctuations affected by ionospheric irregularities. The specific calculation formula is:

[0032] Among them, ; Indicates the mean value operation.

[0033] Step S30, calculate the amplitude scintillation index and the signal attenuation rate and .

[0034] It should be noted that is an index of the severity of signal amplitude fading, is an index of the signal attenuation rate.

[0035] In step S30, the steps to calculate the amplitude scintillation index can be specifically as follows: Step S30a, calculate the mean value of squared; Step S30b, calculate the square of its mean value based on ; Step S30c, calculate the GNSS amplitude scintillation index according to the mean value of squared and the square of the mean value of , the specific formula is:

[0036] Among them, the angular brackets Indicates the calculation of taking the mean.

[0037] In step S30, calculate the signal attenuation rate , that is, by calculating the autocorrelation function, find the time delay when the autocorrelation value drops to . The specific steps can be: Step S30A, according to within the time period (usually taking 1 minute) of time series, calculate its autocorrelation function , and the specific formula is:

[0038] where, represents the mean and variance of; represents the expectation operation.

[0039] Step S30B, calculate within the time period of time series of the autocorrelation function when it drops to . The specific formula is:

[0040] where, is the maximum value within the time period of.

[0041] Compared with the traditional GNSS signal strength calculated by subtracting the mean , which will cause the situation of having outliers, the present application uses the detrended signal strength without subtracting the mean for autocorrelation calculation to obtain , which can ensure the characteristics of the ionospheric scintillation time period , and also makes the values of the non-scintillation time period normal. In addition, this processing also lays a foundation for constructing a scintillation intensity index considering the signal attenuation rate based on the data of geodetic receivers.

[0042] Preferably, since and have different magnitudes under different ionospheric scintillation intensities, therefore, before step S40, it also includes the step of statistically analyzing the relationship between the amplitude scintillation index and the signal attenuation rate . The specific steps of this step are: based on the ionospheric scintillation monitoring device, determine different ionospheric scintillation intensity events; compare and analyze under different ionospheric scintillation intensities, , and Using a large number of ionospheric scintillation events, statistical analysis and The scattered distribution characteristics of .

[0043] Step S40, using The ionospheric amplitude scintillation index is calculated by using the cosine function of .

[0044] Specifically, step S40 may include the following sub-steps: Step S41, using the cosine function The value of Between 0-1; Step S42, and Based on the scattered distribution of and Build Fine-tuning factor ; Step S43, calculating the improved GNSS amplitude flicker index , the specific formula is:

[0045] in, is the data sampling rate; For calculation the time period used; and The maximum value of is equal to that of , thus ensuring that the cosine function value is between 0 and 1; is the fine-tuning factor, and its value is 1.2.

[0046] This embodiment utilizes The improved ionospheric amplitude scintillation index is calculated by using the cosine function of Compared with the widely used index, It can detect ionospheric scintillation more accurately. The specific working principle is: When the GNSS signal attenuation rate is slow, reduce The improved weighting of the index makes Numerical and It can ensure that the signal decays slowly. The accuracy of detecting ionospheric scintillation; when the GNSS signal decays quickly, Increase the pair The index is adjusted to make the satellite signals more seriously affected by ionospheric irregularities Larger in magnitude, can be improved the disadvantage of underestimated values, thus ensuring that in the case of fast signal attenuation ionospheric scintillation can be detected more accurately.

[0047] The beneficial effects of this application are as follows: The GNSS amplitude scintillation index construction method provided by this application uses the cosine function to calculate the ionospheric amplitude scintillation index , and reflects the ionospheric amplitude scintillation situation through the ionospheric amplitude scintillation index , that is when the GNSS signal attenuation speed is slow, by reducing the improved weight of the index, making the value equivalent to , ensuring the accuracy of detecting ionospheric scintillation in the case of slow signal attenuation when the GNSS signal attenuation speed is fast, by increasing the adjustment weight of the index, making the magnitude of the satellite signal more severely affected by ionospheric irregularities larger, improving the disadvantage of underestimated values, thus ensuring more accurate detection of ionospheric scintillation in the case of fast signal attenuation and effectively improving the accuracy of ionospheric scintillation situation monitoring.

[0048] Based on the same inventive concept, this application also provides a GNSS amplitude scintillation index construction system considering the satellite signal attenuation speed, which is applicable to real-time monitoring of GNSS ionospheric scintillation, including a signal strength calculation module, a signal strength processing module, and a solution module and an index construction module.

[0049] Among them, the signal strength calculation module is used to obtain and extract the in-phase component I and the quadrature component Q from the observation data of the ionospheric scintillation monitoring receiver according to the data.

[0050] The signal strength processing module is used to perform detrending processing on the GNSS signal strength by using the average value; and the solution module is used to calculate the amplitude scintillation index and the signal attenuation speed according to the detrended GNSS signal strength; The index construction module is used to combine and Construct an index that takes into account both the amplitude attenuation and speed of the signal , that is, use the cosine function to calculate the ionospheric amplitude scintillation index .

[0051] The following combines specific embodiments to illustrate the GNSS amplitude scintillation index construction method provided by this application.

[0052] Embodiment 1, as shown in Figure 1 , is a schematic flow diagram of the GNSS amplitude scintillation index construction method that takes into account the satellite signal attenuation speed of this application, including the following processing steps: Step 1: Extract the in-phase component I and the quadrature component Q from the raw data of the ionospheric scintillation monitoring receiver.

[0053] Step 2: Calculate the signal strength according to the in-phase component I and the quadrature component Q . The specific formula is:

[0054] where , are the narrowband power and the broadband power respectively, is an empirical constant, generally taking a value of 20; and respectively represent the in-phase component and the quadrature component of the th epoch.

[0055] It should be noted that when the empirical constant takes a value of 20, the sampling frequencies of the in-phase component and the quadrature component are 1000 Hz. For the commonly used ionospheric scintillation monitoring receiver with a sampling frequency of 50 Hz, this application uses the following formula to calculate the signal strength:

[0056] Step 3: Perform detrending processing on the signal strength. The existing signal detrending processing methods usually use a low-pass Butterworth filter for processing, where the filter order is set to 6 and the cut-off frequency is set to 0.1. However, when using the Butterworth filter to perform detrending processing on , its effect is affected by the cut-off frequency; in addition, the Butterworth filtering algorithm is relatively complex and difficult to implement. In this regard, this application uses the method of taking the mean to perform detrending processing on , that is

[0057] in, Indicates the mean operation.

[0058] Step 4: Calculate the amplitude flicker index based on the detrended signal strength , the specific formula is:

[0059] The angle brackets Indicates the mean calculation.

[0060] Step 5: Based on the detrending signal strength Calculate signal attenuation rate GNSS signal attenuation speed indicator Defined as During the period Autocorrelation function of time series Drop to a certain threshold Time delay:

[0061] Among them, for data with a sampling rate of 50 Hz, Generally, 1 minute is taken, that is, 3000 epochs; the threshold Generally take . Autocorrelation function The calculation formula is:

[0062] in, and Respectively The mean and variance of Indicates the expectation operation.

[0063] In the case of ionospheric scintillation, the above formula is used to subtract the mean of Calculated The value range is basically between 0-2s, however, during the non-blinking period There are outliers. Figure 2 As shown in the figure, the output of the commonly used Septentrio brand flicker monitoring receiver (model: PolaRx5S; sampling rate: 50 Hz) and The GPS satellite G19 data was calculated on February 18, 2023. (See Figure 2 Part (b) of the study). Preliminary results show that The value range is basically between 0 and 2 s, and the strong flash period is smaller in general than that in the weak scintillation period . However, there are outliers during the non-scintillation period (for example, at local time 20:49, it is 7.420), and at this time there is no obvious fluctuation (please see Figure 2 part (a) in ), and this situation also exists in the time series of other satellites . It should be noted that the existing literature usually only analyzes the characteristics of the scintillation period . In this regard, the present application uses the without subtracting the mean value to perform autocorrelation calculation to obtain the Figure 2 time series (please see part (c) in ), which can ensure the characteristics of the scintillation period

[0064] Step 6. After calculating the GNSS signal attenuation rate index and the scintillation index , the present application uses the cosine function to adjust the value of so that the adjusted is between 0 and 1; and when the value is large (the signal attenuation rate is slow), the improvement weight for the index is reduced; while when the value is small (the signal attenuation rate is fast), the adjustment weight for the index is increased, so as to ensure that the scintillation intensity of satellite signals more severely affected by low-latitude ionospheric irregularities is greater, and finally construct a low-latitude scintillation intensity index considering the signal attenuation rate . According to the above principle, the present application constructs the following formula to calculate :

[0065] where is the data sampling rate; is the time period used to calculate ; is equivalent to the maximum value of and is the fine-tuning factor, which can generally be taken as 1.2.

[0066] Example 2, as Figure 3 shown, using the in-phase component of GPS satellite G06 at Ledong Station, Hainan on February 18, 2023, with a frequency of 50 Hz and the quadrature-phase component data, demonstrating the detrended signal strength and signal attenuation rate indicators of satellite G06 , and preliminary results. It can be seen that, compared with the output by ISMR, the index has a larger magnitude when the signal attenuation rate is faster, while when the signal attenuation rate is slower, the magnitude is smaller. For example, at local times 21:51 and 21:55, the indices are very close, being 0.904 and 0.890 respectively, but the signal attenuation rate indicators at these two moments are significantly different, being 23.480 and 6.700 respectively.

[0067] Example 3, as Figure 4 shown, respectively presents the time series of the in-phase component , the quadrature-phase component and the detrended signal strength at 21:51 and 21:55. It can be found that although the values at these two moments are approximately equal, the change in the detrended signal strength at local time 21:55 is more obvious. In particular, data interruptions occurred during the two time periods from 21:55:27 to 21:55:31 and from 21:55:44 to 21:55:47, indicating that the index underestimates the ionospheric scintillation intensity at local time 21:55, while the index considering the signal attenuation rate can better characterize the scintillation intensity at this moment.

[0068] The effects after the implementation of this application are mainly as follows: (1) Initially solved the problem that the commonly used ionospheric scintillation index cannot reliably reflect the severity of the impact of GNSS signals by ionospheric irregularities; (2) Unified the method for constructing the signal attenuation rate indicator based on data from different types of GNSS receivers, and solved the problem of outliers in the constructed by traditional methods; (3) Can detect ionospheric scintillation more accurately and can be applied to geodetic GNSS receivers widely distributed globally, reducing the cost of global ionospheric scintillation monitoring and early warning.

[0069] Those skilled in the art can easily understand that the above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for constructing a GNSS amplitude scintillation index, characterized in that It includes the following steps: S10, Obtain and extract the in-phase component and the quadrature-phase component from the observation data of the ionospheric scintillation monitoring receiver according to the data, and calculate the GNSS signal strength based on this; I and the quadrature-phase component Q , and calculate the GNSS signal strength accordingly; S20, perform detrending processing on the GNSS signal strength using the average value; S30. Calculate the amplitude scintillation index and the signal attenuation rate based on the detrended GNSS signal strength. and the signal attenuation rate ; S40, calculate the ionospheric amplitude scintillation index using 's cosine function .

2. The GNSS amplitude scintillation index construction method according to claim 1, wherein Specifically, step S10 is as follows: Convert the said observation data, and select the required satellite system and corresponding time; Retrieve the in-phase component in the transformed observed data I and the quadrature component Q in the corresponding column numbers, so as to extract the in-phase component I and the quadrature component Q data; Using the in-phase component I and the quadrature-phase component Q to calculate the GNSS signal strength , the calculation formula is: 。 3. The GNSS amplitude scintillation index construction method according to claim 1, wherein Specifically, step S20 is as follows: Determine the GNSS signal strength data of the current observation epoch; Determine the GNSS signal strength data for a period of time before and after the current observation epoch; Divide the GNSS signal strength at the current observation epoch by the corresponding mean signal strength to obtain the detrended GNSS signal strength .

4. The GNSS amplitude scintillation index construction method according to claim 1, wherein In step S30, the step of calculating the amplitude scintillation index is specifically as follows: Calculate the detrended GNSS signal strength Mean of the squares; Based on the detrended GNSS signal strength Calculate the square of its mean value; According to the detrended GNSS signal strength The mean of the square and the detrended GNSS signal strength The square of the mean, calculate the amplitude scintillation index , and the calculation formula is: Among them, the angled brackets represent the calculation of taking the mean value.

5. The GNSS amplitude scintillation index construction method according to claim 1, wherein In step S30, calculate the signal attenuation rate The specific steps are as follows: According to the detrended GNSS signal strength within a time period time series, calculate its autocorrelation function , and the calculation formula is: In the formula, represents the variance of; represents the expectation operation; Calculation Within a time period Autocorrelation function of a time series Drops to The time delay at which, i.e., the signal attenuation rate , and the calculation formula is: Among them, is the maximum value within the time period.

6. The GNSS amplitude scintillation index construction method according to claim 1, characterized in that Before step S40, it also includes: Statistical analysis of the amplitude scintillation index and the signal attenuation rate relationship 7. The GNSS amplitude scintillation index construction method according to claim 1, characterized in that Specifically, step S40 is as follows: Use the cosine function to adjust the value of the signal attenuation rate so that the adjusted is between 0 and 1; Determine the combination based on statistical analysis methods and Construct fine-tuning factor ; Calculate the GNSS amplitude scintillation index , and the calculation formula is as follows: Among them, is the data sampling rate; is for calculation the time period used.

8. A GNSS amplitude scintillation index construction system, characterized in that It includes: A signal strength calculation module, configured to obtain and, based on the observation data of an ionospheric scintillation monitoring receiver, extract the in-phase component in the data I and the quadrature-phase component Q , and calculate the GNSS signal strength accordingly; A signal strength processing module for performing detrending processing on the GNSS signal strength using the average value; and a solution module, configured to calculate an amplitude scintillation index based on the detrended GNSS signal strength and a signal attenuation rate ; An index construction module for calculating the ionospheric amplitude scintillation index using the cosine function of .

Images