Solar radio storm intensity judgment method and system based on utility function comprehensive evaluation
阅读说明：本技术 基于效用函数综合评价的太阳射电暴强度判定方法及系统 (Solar radio storm intensity judgment method and system based on utility function comprehensive evaluation ) 是由 陈新 祝雪芬 于 2020-05-07 设计创作，主要内容包括：本发明提供了一种基于效用函数综合评价的太阳射电暴强度判定方法及系统,包括：步骤1：获取观测地卫星的载噪比、定位误差、几何因子和卫星失锁数目,进行数据预处理,计算出观测地载噪比下降值和观测地总定位误差；步骤2：进行相对无量纲化处理,把数据预处理后的实测值转化成单项评价指数；步骤3：对太阳射电暴强度进行判定划分并进行无量纲化处理,确定指标的权重；步骤4：根据指标的权重,采用加权合成准则将单项评价指数合成对太阳射电暴强度进行判定的综合评价指数。该判定方法能全天候计算各个地区太阳射电暴强度,效率高,过程简单,并且不依赖于射电望远镜,成本低。(The invention provides a solar radio storm intensity judgment method and system based on utility function comprehensive evaluation, which comprises the following steps: step 1: acquiring a carrier-to-noise ratio, a positioning error, a geometric factor and a satellite lock losing number of an observation ground satellite, performing data preprocessing, and calculating a carrier-to-noise ratio reduction value of the observation ground and a total positioning error of the observation ground; step 2: carrying out relative dimensionless treatment, and converting the measured value after data pretreatment into a single evaluation index; and step 3: judging and dividing the intensity of the solar radio storm, carrying out dimensionless processing, and determining the weight of an index; and 4, step 4: and synthesizing the single evaluation index into a comprehensive evaluation index for judging the intensity of the solar radio storm by adopting a weighted synthesis rule according to the weight of the index. The method can calculate the solar radio storm intensity of each area in all weather, has high efficiency and simple process, does not depend on a radio telescope, and has low cost.)
1. A solar radio storm intensity judgment method based on utility function comprehensive evaluation is characterized by comprising the following steps:
step 1: acquiring a carrier-to-noise ratio, a positioning error, a geometric factor and a satellite lock losing number of an observation ground satellite, performing data preprocessing, and calculating a carrier-to-noise ratio reduction value of the observation ground and a total positioning error of the observation ground;
step 2: carrying out relative dimensionless treatment, and converting the measured value after data pretreatment into a single evaluation index;
and step 3: judging and dividing the intensity of the solar radio storm, carrying out dimensionless processing, and determining the weight of an index;
and 4, step 4: and synthesizing the single evaluation index into a comprehensive evaluation index for judging the intensity of the solar radio storm by adopting a weighted synthesis rule according to the weight of the index.
2. The method for determining the intensity of a solar radio storm based on the utility function comprehensive evaluation as claimed in claim 1, wherein the carrier-to-noise ratio is the carrier-to-noise ratio of each satellite, when no lock is lost in an observation site, the carrier-to-noise ratio is measured by at least four satellites, and the average value of the reduction of the carrier-to-noise ratios of the four satellites is taken as the reduction of the carrier-to-noise ratio of the station and is recorded as x1Then, there are:
wherein, snr isiThe unit is dBHz, which is the reduction of the carrier-to-noise ratio of each satellite; i represents: the ith satellite, i ═ 1,2,3, 4.
3. The method for determining solar radiostorm intensity based on utility function comprehensive evaluation according to claim 1, wherein the observation place positioning errors include positioning errors in X, Y, Z three directions in a northeast coordinate system, and an average value of absolute values of the three positioning errors is taken as a positioning error of the station and is recorded as x2Then, there are:
wherein r isi(i is 1,2, and 3) represents a positioning error in the direction X, Y, Z, respectively, and the unit is m.
4. The method for determining solar radio storm intensity based on utility function comprehensive evaluation as claimed in claim 1, wherein the geometric factor and the number of satellite lost locks are expressed by specific numbers respectively marked as x3、x4The unit is particle.
5. The method for determining the intensity of a solar radiostorm based on the utility function comprehensive evaluation as claimed in claim 1, wherein the step 2 comprises:
the more the carrier-to-noise ratio is reduced, the larger the positioning error, the larger the geometric factor, the larger the number of lost locks of the satellite and the larger the intensity of the solar radio storm;
adopting an increasing S-shaped function as a relative dimensionless function;
when the intensity of the solar radio storm is lower than a preset threshold value, the influence of the intensity of the solar radio storm on each measured value is small;
when the solar radio storm is gradually increased to 0.5, each measured value is gradually increased under the influence of the intensity of the solar radio storm;
when the intensity of the solar radio storm exceeds 0.5 and is gradually increased, the influence of the solar radio storm on each measured value is gradually reduced;
the increasing sigmoid function is represented as:
wherein a is the lower bound of the fluctuation of the index but not caused by the solar radio storm, and the dimensionless value is 0 when the lower bound is less than the index; b is an actual measurement value corresponding to the dimensionless index value of 0.5, namely, from the point (b,0.5), the influence of the solar radiation storm on each actual measurement value is gradually reduced; c is an actual measurement value corresponding to a dimensionless value of 1; and x is a carrier-to-noise ratio reduction value of the observation ground, a total positioning error of the observation ground, a geometric factor and the number of lost locks of the satellites, which are obtained after data preprocessing.
6. The method for determining the intensity of a solar radiostorm based on the utility function comprehensive evaluation as claimed in claim 1, wherein the step 3 comprises:
when the intensity of the solar radio storm is higher than a preset threshold value, the phenomenon of earth lock losing occurs, the number of locked satellites is less than 4, and the intensity of the solar radio storm is judged as follows: no lock loss occurred observably and lock loss occurred observably.
7. The method for determining solar radio storm intensity based on utility function comprehensive evaluation as claimed in claim 6, wherein when no lock loss occurs in observation ground, the weighting q is decreased for carrier-to-noise ratio12; for positioning error and geometric factor, the attached weight is q2＝q31 is ═ 1; setting the sum of the weights of the three indexes asThen there are:
wherein, w1,w2,w3The weights of three indexes, namely, carrier-to-noise ratio reduction, positioning error and geometric factor, are respectively calculated as follows: w is a1＝0.4,w2＝w30.2; j represents: the carrier-to-noise ratio is reduced, the positioning error and the geometric factor correspond to j being 1,2 and 3; w is ajRepresenting the weight corresponding to each index; q. q.sjRepresenting the attached weight of each index;
when the observation place is unlocked, the positioning cannot be carried out, the geometric factor cannot be detected, the carrier-to-noise ratio can not accurately reflect the change of the intensity of the solar radio storm, the relative dimensionless values corresponding to the three indexes are 1, and w4Weight of the index, w, for the number of out-of-lock satellites4＝0.2。
8. The method for determining solar radiostorm intensity based on utility function comprehensive evaluation according to claim 1, wherein the step 4 comprises:
synthesizing the intensity of the solar radiostorm by adopting a weighted average method synthetic model, and recording as S:
wherein, y1,y2,y3,y4Respectively representing relative dimensionless values corresponding to four indexes of observation ground carrier-to-noise ratio reduction, positioning error, geometric factor and satellite lock losing number; w is aiRepresenting the weight corresponding to four indexes of the reduction of the carrier-to-noise ratio, the positioning error, the GDOP factor and the number of the out-of-lock satellites, i is 1,2,3，4；yirepresents: the four indexes correspond to relatively dimensionless values.
9. The method for determining solar radiostorm intensity based on utility function comprehensive evaluation according to claim 8, wherein when no unlocking occurs in an observation place:
S＝0.4y1+0.2y2+0.2y3；
when the observation place is unlocked:
S＝0.8+0.2y4。
10. a solar radio storm intensity judgment system based on utility function comprehensive evaluation is characterized by comprising the following components:
module M1: inputting a carrier-to-noise ratio, a positioning error, a geometric factor and a satellite unlocking number of an observation ground satellite, performing data preprocessing, and calculating a carrier-to-noise ratio reduction value of the observation ground and a total positioning error of the observation ground;
module M2: carrying out relative dimensionless treatment, and converting the measured value after data pretreatment into a single evaluation index;
module M3: judging and dividing the intensity of the solar radio storm, carrying out dimensionless processing, and determining the weight of an index;
module M4: and synthesizing the single evaluation index into a comprehensive evaluation index for judging the intensity of the solar radio storm by adopting a weighted synthesis rule according to the weight of the index.
Technical Field
The invention relates to the technical field of wireless communication, in particular to a method and a system for judging the intensity of a solar radio storm based on comprehensive evaluation of a utility function.
Background
With the increasingly wide application of satellite technologies such as GNSS in modern society, the influence of solar radio storm intensity on GNSS signals becomes a significant part. The input of the sun to the earth's material and energy directly determines the maintenance and transition of the earth and its surrounding spatial environment. The observation of the solar radio storm plays a vital role in researching solar physics and monitoring space weather environment, and has important research value of space physics science. The solar radiation explosion is a phenomenon of ray enhancement and sudden increase of radio noise generated when strong disturbance suddenly appears on the sun, and commonly occurs with the phenomena of flare, X-ray explosion, even proton explosion or cosmic ray explosion and the like in the solar activity area. During the past multiple solar radio outbreaks, the satellite carrier-to-noise ratio is reduced due to the solar radio storm, the positioning error is increased, the geometric accuracy factor is increased, the satellite is unlocked, the number of visible stars is greatly reduced, and the like, so that the fact that a noise signal generated by the solar radio outbreak is one of the influence factors of a navigation signal is proved.
The influence of solar radio explosion on the GNSS signal is manifold, so that the navigation system can normally work when a solar radio explosion occurs in order to better deal with the interference of the solar radio explosion on the GNSS signal, the judgment on the intensity of the solar radio explosion is very important, and the navigation system has great significance for maintaining the normal operation of the satellite navigation system.
Patent document CN106771653B (application number: 201611063812.6) discloses a fast radiostorm real-time detection device, system and method, wherein the fast radiostorm real-time detection device includes a data acquisition and preprocessing unit, a real-time search and storage unit, and a search result fast publishing unit.
The traditional method for detecting the solar radio activity mainly depends on a radio telescope, but is expensive in manufacturing cost, sparse in distribution and incapable of carrying out all-weather real-time monitoring on the solar radio.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a solar radio storm intensity judgment method and system based on utility function comprehensive evaluation.
The solar radio storm intensity judgment method based on utility function comprehensive evaluation provided by the invention comprises the following steps:
step 1: acquiring a carrier-to-noise ratio, a positioning error, a geometric factor and a satellite lock losing number of an observation ground satellite, performing data preprocessing, and calculating a carrier-to-noise ratio reduction value of the observation ground and a total positioning error of the observation ground;
step 2: carrying out relative dimensionless treatment, and converting the measured value after data pretreatment into a single evaluation index;
and step 3: judging and dividing the intensity of the solar radio storm, carrying out dimensionless processing, and determining the weight of an index;
and 4, step 4: and synthesizing the single evaluation index into a comprehensive evaluation index for judging the intensity of the solar radio storm by adopting a weighted synthesis rule according to the weight of the index.
Preferably, the carrier-to-noise ratio is the carrier-to-noise ratio of each satellite, when the observation ground is not unlocked, at least four satellites measure the carrier-to-noise ratio, and the average value of the reduction of the carrier-to-noise ratios of the four satellites is taken as the reduction of the carrier-to-noise ratio of the station and is marked as x1Then, there are:
wherein, snr isiThe unit is dBHz, which is the reduction of the carrier-to-noise ratio of each satellite; i represents: the ith satellite, i ═ 1,2,3, 4.
Preferably, the observation location error includes X, Y, Z location errors in three directions in the northeast coordinate system, and the average value of the absolute values of the three location errors is taken as the location error of the station and recorded as x2Then, there are:
wherein r isi(i is 1,2, and 3) represents a positioning error in the direction X, Y, Z, respectively, and the unit is m.
Preferably, the geometric factor and the number of lost locks of the satellite are expressed by specific numbers respectively marked as x3、x4The unit is particle.
Preferably, the step 2 comprises:
the more the carrier-to-noise ratio is reduced, the larger the positioning error, the larger the geometric factor, the larger the number of lost locks of the satellite and the larger the intensity of the solar radio storm;
adopting an increasing S-shaped function as a relative dimensionless function;
when the intensity of the solar radio storm is lower than a preset threshold value, the influence of the intensity of the solar radio storm on each measured value is small;
when the solar radio storm is gradually increased to 0.5, each measured value is gradually increased under the influence of the intensity of the solar radio storm;
when the intensity of the solar radio storm exceeds 0.5 and is gradually increased, the influence of the solar radio storm on each measured value is gradually reduced;
the increasing sigmoid function is represented as:
wherein a is the lower bound of the fluctuation of the index but not caused by the solar radio storm, and the dimensionless value is 0 when the lower bound is less than the index; b is an actual measurement value corresponding to the dimensionless index value of 0.5, namely, from the point (b,0.5), the influence of the solar radiation storm on each actual measurement value is gradually reduced; c is an actual measurement value corresponding to a dimensionless value of 1; and x is a carrier-to-noise ratio reduction value of the observation ground, a total positioning error of the observation ground, a geometric factor and the number of lost locks of the satellites, which are obtained after data preprocessing.
Preferably, the step 3 comprises:
when the intensity of the solar radio storm is higher than a preset threshold value, the phenomenon of earth lock losing occurs, the number of locked satellites is less than 4, and the intensity of the solar radio storm is judged as follows: no lock loss occurred observably and lock loss occurred observably.
Preferably, when no lock loss occurs in an observation place, the carrier-to-noise ratio is reduced by an additional weight q12; for positioning error and geometric factor, the attached weight is q2＝q31 is ═ 1; setting the sum of the weights of the three indexes asThen there are:
wherein, w1,w2,w3Respectively representing carrier-to-noise ratio drop, positioning error and geometric factorCalculating the weights of the three indexes: w is a1＝0.4,w2＝w30.2; j represents: the carrier-to-noise ratio is reduced, the positioning error and the geometric factor correspond to j being 1,2 and 3; w is ajRepresenting the weight corresponding to each index; q. q.sjRepresenting the attached weight of each index;
when the observation place is unlocked, the positioning cannot be carried out, the geometric factor cannot be detected, the carrier-to-noise ratio can not accurately reflect the change of the intensity of the solar radio storm, the relative dimensionless values corresponding to the three indexes are 1, and w4Weight of the index, w, for the number of out-of-lock satellites4＝0.2。
Preferably, the step 4 comprises:
synthesizing the intensity of the solar radiostorm by adopting a weighted average method synthetic model, and recording as S:
wherein, y1,y2,y3,y4Respectively representing relative dimensionless values corresponding to four indexes of observation ground carrier-to-noise ratio reduction, positioning error, geometric factor and satellite lock losing number; w is aiRepresenting weights corresponding to four indexes of carrier-to-noise ratio reduction, positioning error, GDOP factor and number of unlocked satellites, wherein i is 1,2,3 and 4; y isiRepresents: the four indexes correspond to relatively dimensionless values.
Preferably, when no loss of lock has occurred at the observation site:
S＝0.4y1+0.2y2+0.2y3；
when the observation place is unlocked:
S＝0.8+0.2y4。
the solar radio storm intensity judgment system based on utility function comprehensive evaluation provided by the invention comprises:
module M1: inputting a carrier-to-noise ratio, a positioning error, a geometric factor and a satellite unlocking number of an observation ground satellite, performing data preprocessing, and calculating a carrier-to-noise ratio reduction value of the observation ground and a total positioning error of the observation ground;
module M2: carrying out relative dimensionless treatment, and converting the measured value after data pretreatment into a single evaluation index;
module M3: judging and dividing the intensity of the solar radio storm, carrying out dimensionless processing, and determining the weight of an index;
module M4: and synthesizing the single evaluation index into a comprehensive evaluation index for judging the solar radio storm intensity by adopting a weighted synthesis rule according to the weight of the index.
Compared with the prior art, the invention has the following beneficial effects:
1. the judging method can calculate the solar radio storm intensity of each area in all weather, has high efficiency and simple process, does not depend on a radio telescope, and has low cost;
2. the method has feasibility for evaluating the intensity of the solar radio storm, and quantitative evaluation results not only can show the change of the intensity of the solar radio storm of a single station, but also can simultaneously give the influence of the solar radio storm on a plurality of stations;
3. compared with the traditional method, the method has the advantages of low cost, simple algorithm, combination of multiple factors and relatively improved identification accuracy and efficiency.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic flow diagram of one embodiment of the present invention;
FIG. 2 is a diagram of a relatively non-dimensionalized "S" type function.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.