Multi-GNSS Jamming Source Deployment Optimization Method and System

Abstract

The present invention relates to the field of satellite navigation technology, and particularly relates to a method and system for optimizing the deployment of multi-GNSS jamming sources. Constraint conditions for generating jamming sources are set according to the types of area protection, where the types of area protection include key suppression area types, key protection area types, and general control area types. Based on the constraint conditions for generating jamming sources, when the number of jamming sources is fixed, the GNSS signals in the area are jointly suppressed and interfered with by using various jamming source types, or when the types of jamming sources are fixed, the least number of jamming sources are used to suppress and interfere with the GNSS signals in the area. By refining the different purposes of area jamming operations, the present invention distinguishes three types of areas, namely key suppression, key protection, and general control areas, and more pertinently implements the control of GNSS area signal suppression and interference, ensuring that the deployment plan has good practical application value, capable of automatically generating a multi-GNSS jamming source deployment plan, reducing manual intervention, and improving the efficiency of generating the deployment plan.

Description

Technical Field

[0001] The present invention relates to the technical field of satellite navigation, and particularly to a method and system for optimizing the deployment of multi-GNSS jamming sources. Background Art

[0002] Currently, the Global Navigation Satellite System (GNSS) has been deeply applied to many important fields such as economy, national defense, and social security. The wide application of GNSS makes it easy for people to obtain spatio-temporal information, which also provides opportunities for some illegal or even terrorist activities. Therefore, it is of great significance to effectively suppress and control GNSS signals in a specific area for the security of spatio-temporal information applications. However, the current technology mainly uses a single jamming source to directly implement jamming. However, this method is easily affected by terrain occlusion, making it difficult to achieve full coverage jamming of the area, and it is also easy to expose its own position. Summary of the Invention

[0003] Therefore, the present invention provides a method and system for optimizing the deployment of multi-GNSS jamming sources, which can achieve more effective regional jamming control through the cooperation of multiple jamming sources considering the influence of terrain conditions, and improve the security of spatio-temporal information applications.

[0004] According to the design scheme provided by the present invention, a method for optimizing the deployment of multi-GNSS jamming sources is provided, including:

[0005] Setting the generation constraint conditions of the jamming source according to the regional protection type, where the regional protection type includes the key jamming area type, the key protection area type, and the general control area type;

[0006] Based on the generation constraint conditions of the jamming source, and using each jamming source type to jointly jam the regional GNSS signals under the condition of a fixed number of jamming sources, or using the minimum number of jamming sources to jam the regional GNSS signals under the condition of a fixed jamming source type.

[0007] As the method for optimizing the deployment of multi-GNSS jamming sources of the present invention, further, the generation constraint conditions of the jamming source set according to the regional protection type include: for the key jamming area, setting the area greater than X-fold jamming coverage according to the application requirements; for the key protection area, setting the area less than Y-fold jamming coverage according to the anti-jamming performance of the protected target objects in the area; for the general control area, setting single-fold jamming for the area outside the key jamming area and the key protection area, where X and Y are user-defined parameters.

[0008] As the multi-GNSS jamming source deployment optimization method of the present invention, further, when the number of jamming sources is fixed, the GNSS signals in the region are collaboratively jammed using various jamming source types, including the following:

[0009] First, for the mission target area, the fixed number of jamming sources set is n, and the number m of DEM grid points is obtained by reading the digital elevation model DEM of the area. n points are selected from the m DEM grid points to form C m n deployment schemes;

[0010] Next, for each deployment scheme, the visible area of each jamming source is obtained using the visibility analysis method, and the number of visible jamming sources at each grid point is calculated one by one; according to the maximum or minimum number of visible jamming sources at each point, a deployment scheme that simultaneously satisfies the jamming source generation constraint conditions for the key suppression area and the key protection area is obtained, and the combined jamming area of all jamming sources under the constraint conditions of the general control area is calculated. According to the combined jamming area and the total area of the mission target area, the mission area coverage rate is calculated;

[0011] Then, the optimal deployment scheme of the jamming source is output according to the mission area coverage rate.

[0012] As the multi-GNSS jamming source deployment optimization method of the present invention, further, when the type of jamming source is fixed, the GNSS signals in the region are jammed using the minimum number of jamming sources, including:

[0013] First, for the mission target area, the fixed number of jamming sources set is n, and the number m of DEM grid points is obtained by reading the digital elevation model DEM of the area; the highest elevation point among the m DEM grid points is found, and this highest elevation point is used as the jamming source layout position point. The coverage range of the position point is obtained using the visibility analysis method, and this coverage range of the position point is used as the original coverage range;

[0014] Then, the new highest elevation point outside the original coverage range is found, and the coverage range of the new position point is obtained using the visibility analysis method. The coverage range of the new position point is added to the original coverage range. The grid points with the accumulated area greater than the area of the original position point coverage range are used as the jamming source layout position points, and the accumulated result is used as the original coverage range for the next round of iteration to find the new highest elevation point until the area of the original coverage range obtained in the iteration and the area of the mission target area meet the preset proximity threshold, and the optimal deployment of the jamming source is completed.

[0015] As the multi-GNSS jamming source deployment optimization method of the present invention, further, in obtaining the coverage range of the position points by using the visibility analysis method, it includes: for the case where the coverage range of the position points does not meet the key protection area, by looking down for the second-highest elevation point among the grid points as the position point for jamming source layout, and by re-obtaining the coverage range of the position points to meet the key protection area.

[0016] As the multi-GNSS jamming source deployment optimization method of the present invention, further, in obtaining the coverage range of the position points by using the visibility analysis method, it also includes: for the case where the coverage range of the position points does not meet the key jamming area, taking the elevation from high to low in sequence as the position points for jamming source layout, by re-obtaining the coverage range of the position points to determine the jamming source layout scheme that meets the conditions of the key jamming area, and simultaneously checking whether the coverage range of the position points meets the key protection area until all the key protection areas and key jamming areas are within the coverage range of the position points.

[0017] As the multi-GNSS jamming source deployment optimization method of the present invention, further, in obtaining the visibility area of each jamming source by using the visibility analysis method, it includes:

[0018] First, obtain relevant jamming parameters and the maximum jamming distance according to the type of jamming source;

[0019] Next, take the position of the jamming source as the center of a circle and the maximum jamming distance as the radius to construct an analysis circle, take the circumscribed square corresponding to the analysis circle as the analysis area, and read the digital elevation model (DEM) data corresponding to the analysis area;

[0020] Then, starting from the position directly in front of the jamming source and using eight direction lines evenly divided in the 360-degree direction as the dividing lines, divide the analysis area, conduct visibility analysis on the grid points on the direction lines according to the point-to-point visibility algorithm, and conduct visibility analysis on each grid point within each divided area outside the direction lines, and merge the visibility analysis results to construct the visibility area of the jamming source.

[0021] As the multi-GNSS jamming source deployment optimization method of the present invention, further, the visibility analysis of the grid points on the direction lines according to the point-to-point visibility algorithm includes:

[0022] First, connect the jamming source with the nearest grid point on the direction line to form a line of sight, and record the elevation angle of the nearest grid point as the current maximum elevation angle;

[0023] Next, calculate the elevation angles between the interference source and each grid point on the direction line outward from the interference source in sequence, and compare each calculated elevation angle with the current maximum elevation angle. If the calculated elevation angle is greater than the current maximum elevation angle, it is determined that the corresponding grid point is visible, and the current maximum elevation angle is updated using the calculated elevation angle. If the calculated elevation angle is less than the current maximum elevation angle, it is determined that the corresponding grid point is not visible;

[0024] Then, obtain all visible grid points on the direction line and record them using the visibility matrix, and perform corresponding marking in the analysis area through the visibility matrix.

[0025] As the multi-GNSS jamming source deployment optimization method of the present invention, further, perform visibility analysis on each grid point in each divided area outside the direction line, including:

[0026] First, connect the interference source and the target grid point to form a connection line, and obtain the plane coordinates of the intersection points of the connection line and the grid lines;

[0027] Next, determine the elevation values of the intersection points based on the coordinates of the two grid points adjacent to the intersection points, and obtain the elevation values of all intersection points through linear interpolation;

[0028] Then, form an equivalent direction line with the interference source, all intersection points, and the target grid point;

[0029] Finally, perform visibility analysis on the equivalent direction line according to the point-to-point visibility algorithm to obtain the visibility of each grid point in each divided area, and record it using the visibility matrix.

[0030] Further, the present invention also provides a multi-GNSS jamming source deployment optimization system, including: a constraint condition setting module and a jamming scheme optimization module, where,

[0031] The constraint condition setting module is used to set the constraint conditions for generating the jamming source according to the regional protection type, where the regional protection type includes the key jamming area type, the key protection area type, and the general control area type;

[0032] The jamming scheme optimization module is used to generate constraint conditions based on the jamming source, and use each jamming source type to jointly jam the regional GNSS signal when the number of jamming sources is fixed or use the minimum number of jamming sources to jam the regional GNSS signal when the jamming source type is fixed.

[0033] Further, the present invention also provides a multi-GNSS jamming system, implemented based on the above method, including: a site selection and planning module, a deployment generation module, and a real-time monitoring module, where,

[0034] The site selection and planning module is used to set the site selection and planning of the jamming source device in the target area according to the operation service requirements;

[0035] A deployment generation module, configured to obtain an optimal deployment path of interference source devices according to the site selection planning result and by using the above-mentioned method;

[0036] A real-time monitoring module, configured to monitor and display in real time the working parameters of the interference source devices and the interference effect in the target area, wherein the working parameters at least include the position, status and parameters of the interference source devices.

[0037] Advantages of the present invention:

[0038] By refining different purposes of regional suppression operations, the present invention distinguishes three types of areas: key suppression, key protection, and general control, so as to achieve more targeted GNSS regional signal suppression interference control; combined with actual control conditions, the deployment plan is optimized according to two categories: fixed number of interference sources and fixed type of interference sources, ensuring that the deployment plan has good practical application value; fully considering the influence of terrain conditions to give play to the collaborative advantages of multiple GNSS suppression interference sources, it can realize the automatic generation of the deployment plan of multiple GNSS interference sources, reduce manual intervention, improve the efficiency of generating the deployment plan, and has good application prospects. Description of the Drawings

[0039] Figure 1 It is a schematic diagram of the deployment optimization process of multiple GNSS suppression interference sources in the embodiment;

[0040] Figure 2 It is a schematic diagram of the principle of the automatic generation algorithm of the deployment plan under the condition of fixed number of interference sources in the embodiment;

[0041] Figure 3 It is a schematic diagram of the principle of the automatic generation algorithm of the deployment plan under the condition of fixed type of interference sources in the embodiment;

[0042] Figure 4 It is a schematic diagram of the analysis area division in the embodiment;

[0043] Figure 5 It is a schematic diagram of the visibility analysis of the interference source and the grid points on the direction line in the embodiment;

[0044] Figure 6 It is a schematic diagram of the principle of the GNSS suppression interference source visibility analysis algorithm in the embodiment. Detailed Embodiments

[0045] To make the purpose, technical solutions and advantages of the present invention clearer and more understandable, the present invention will be further described in detail below with reference to the drawings and technical solutions.

[0046] For a traditional single interference source, to achieve effective interference suppression, the required interference power increases with the increase of the interference distance, which brings the problem that it is easier to expose its own position and is easily detected. "Easy to detect" is obviously not conducive to the operation of the interference source. At the same time, considering that the operation object may be a receiver with a certain anti-interference ability, the effectiveness of a single interference source will be more restricted. For an anti-interference receiver with N array elements, it is usually required that the number of interference sources is greater than or equal to N to possibly generate effective interference. Therefore, the influence of a single interference source is relatively "easy to resist" in practice. For this reason, in the embodiments of the present invention, as shown in Figure 1 shown, a method for optimizing the deployment of multiple GNSS jamming sources is provided, including:

[0047] S101. Set the constraints for generating jamming sources according to the type of area protection, where the type of area protection includes the type of key suppression area, the type of key protection area, and the type of general control area;

[0048] S102. Based on the constraints for generating jamming sources, and when the number of jamming sources is fixed, use various jamming source types to jointly suppress the regional GNSS signals, or when the type of jamming source is fixed, use the minimum number of jamming sources to suppress the regional GNSS signals.

[0049] In the embodiments of this case, through the collaborative operation of multiple dispersed small GNSS jamming sources, the defects of "easy to detect" and "easy to resist" of traditional single jamming sources can be overcome, providing a new technical route for GNSS signal suppression control and effectively improving the security of spatio-temporal information applications.

[0050] As a preferred embodiment, further, the constraints for generating jamming sources set according to the type of area protection include: for the key suppression area, set the area greater than X-fold interference coverage according to the application requirements; for the key protection area, set the area less than Y-fold interference coverage according to the anti-interference performance of the protected target object in the area; for the general control area, set the area with single-fold jamming for the area outside the key suppression area and the key protection area, where X and Y are user-defined parameters.

[0051] To better achieve GNSS signal regional suppression and fully consider the actual application operation requirements, the entire area is divided into three types: (1) Key suppression area (can be regarded as the blue area): This area is the key suppression area. Even if it is considered that the suppression target may use an anti-jamming receiver, the suppression effect must be achieved normally. It is required to be greater than X times the interference coverage. X is an input parameter given by the user according to the application requirements. (2) Key protection area (can be regarded as the red area): This area is the key protection area, mainly considering facilities or areas that need to be key protected within the area, such as mobile communication base stations, etc. It is necessary to keep the GNSS signal continuous and effective. The user can require that the interference coverage within this area must be less than Y times according to the anti-jamming performance of the protected object. Y is an input parameter. (3) General control area (can be regarded as the green area): This area is the general control area, mainly considering other areas outside the red and blue areas, and only single-layer suppression interference needs to be achieved.

[0052] When automatically generating a multi-GNSS suppression interference source deployment plan, the above constraints are considered to achieve the optimal deployment of the interference source.

[0053] The main factors affecting the cost performance of the deployment plan: one is the number of interference sources. The more the number, the greater the cost. The other is the type of interference source. Different types have different performances. Usually, the better the performance, the greater the cost, but at the same time, the suppression distance of a single interference source is farther.

[0054] On the basis of fully considering the deployment constraints and combining the actual application situation, the deployment plan can be optimized in two cases: (1) The number of interference sources is fixed: The number of interference sources is fixed, and the type of interference source is not fixed. The actual situation is: The suppression interference equipment of the operation entity is certain (a fixed number of devices and their types), and coordinated deployment suppression needs to be achieved. (2) The type of interference source is fixed: The type of interference source is fixed, and the number is not fixed. The actual situation is: The operation entity selects a fixed type of interference source, aiming to achieve the optimal deployment and seeking the minimum number of interference sources.

[0055] When using various interference source types to suppress GNSS signals in the area collaboratively under the condition that the number of interference sources is fixed, it can be designed to include the following content:

[0056] First, for the task target area, the fixed number of interference sources set is n, and the number of DEM grid points m is obtained by reading the digital elevation model DEM of the area. C m n deployment plans are formed by selecting n points from the m DEM grid points.

[0057] Next, for each deployment plan, the visible area analysis method is used to obtain the visible areas of each interference source, and the number of visible interference sources at each grid point is calculated one by one; according to the maximum or minimum number of visible interference sources at each point, a deployment plan that simultaneously satisfies the suppression interference source generation constraint conditions of the key suppression area and the key protection area is obtained, and the combined interference area of all interference sources under the constraint conditions of the general control area is calculated. Based on this combined interference area and the total area of the mission target area, the mission area coverage rate is calculated;

[0058] Then, based on the mission area coverage rate, the optimal deployment plan of the suppression interference source is output.

[0059] As Figure 2 shown, assuming that the number of interference sources is fixed, which is n, the steps of the automatic deployment plan generation algorithm can be described as follows:

[0060] (1) Read the digital elevation model (DEM) of the mission area, calculate the number of DEM grid points in the mission area. Assuming there are m grid points in total, and arbitrarily select n (the given number of interference sources) points from them, then there are C m n possibilities.

[0061] (2) Select one from the C m n possible plans, and use the visible area analysis method to analyze the visible areas of the n interference sources at this time.

[0062] (3) According to the visible area analysis results, judge whether the red area (must be less than the Y heavy interference coverage) and the blue area (must be greater than the X heavy interference coverage) meet the conditions. The method is as follows: Divide it into the red area and the blue area. According to the visible area analysis results of each interference source, calculate the number of visible interference sources at each grid point in the area one by one. For the red area, find the maximum number N of visible interference sources at each point in the area. If N < Y, for the blue area, find the minimum number N' of visible interference sources at each point in the area. If N' > X, and both of the above two conditions are met, then continue; otherwise, discard the alternative plan.

[0063] (4) Under the condition of only considering single - layer coverage, calculate the combined interference area S' of all interference sources, compare it with the total area S of the mission area, and calculate the mission area coverage rate

[0064] ξ=(S' / S)×100% (1)

[0065] where ξ represents the mission area coverage rate, and its initial value is set to 0.

[0066] (5) Repeat steps (2)-(4) to calculate the new task area coverage rate ξ′, compare it with the existing task area coverage rate. If ξ′ > ξ, then let ξ = ξ′ and return to step (2); otherwise, directly return to step (2). And so on until: the task area coverage rate ξ reaches 100%, or all searches are completed. At this time, the green area coverage condition reaches the optimum.

[0067] (6) Output the finally retained solution, which is the optimal solution.

[0068] Furthermore, when the type of interference source is fixed, using the least number of interference sources to suppress the GNSS signal in the area can be designed to include the following content:

[0069] First, for the task target area, set the fixed number of interference sources as n, and obtain the number of DEM grid points m by reading the digital elevation model DEM of the area; find the highest elevation point among the m DEM grid points, use this highest elevation point as the location point for deploying the interference source, use the visibility analysis method to obtain the coverage range of the location point, and use this coverage range of the location point as the original coverage range;

[0070] Then, find the new highest elevation point outside the original coverage range, use the visibility analysis method to obtain the coverage range of the new location point, accumulate the coverage range of the new location point and the original coverage range, use the grid points whose accumulated area is greater than the area of the original location point coverage range as the location points for deploying the interference source, and use the accumulated result as the original coverage range for the next round of iteration to find the new highest elevation point until the area of the original coverage range obtained in the iteration is close enough to the area of the task target area, and complete the optimal deployment of the suppression interference source.

[0071] For the case where the coverage range of the location point does not meet the requirements of the key protection area, the second highest elevation point can be found downward among the grid points as the location point for deploying the interference source, and the coverage range of the location point can be re-obtained to meet the key protection area. And for the case where the coverage range of the location point does not meet the requirements of the key suppression area, the location points for deploying the interference source are arranged in descending order of elevation, the coverage range of the location point is re-obtained to determine the deployment plan of the interference source that meets the conditions of the key suppression area, and at the same time, check whether the coverage range of the location point meets the requirements of the key protection area until all key protection areas and key suppression areas are within the coverage range of the location point.

[0072] As Figure 3 shown, assuming that the type of interference source is fixed, the farthest interference distance of a single interference source is certain. Under the condition that there is no special limit on the number of interference sources, in order to achieve the optimal suppression control coverage, the deployment scheme automatic generation algorithm is as follows:

[0073] (1) Read the digital elevation model (DEM) of the mission area, calculate the number of DEM grid points in the mission area. Assume there are m grid points in total, and record the total area S of the mission area.

[0074] (2) Find the highest point P0 among the m grid points. Then, use this point as the interference source layout point, and use the visibility analysis method to analyze the coverage range of the found location point, and record the coverage area S0. At the same time, check whether the red area no longer meets the conditions. If not, find the second-highest point in the area until the red area conditions are met.

[0075] (3) Outside the above coverage range, find the highest point P1, layout the interference source, and determine its actual coverage range through the visibility analysis method. At the same time, check whether the red area no longer meets the conditions. If not, find the second-highest point in the area until the red area conditions are met. Then, add the new coverage range to the existing coverage range, record the coverage area S1, compare this area with the original coverage area. If S1 ≤ S0, discard this point. If S1 > S0, use the corresponding grid point as the interference source layout point, and let S0 = S1.

[0076] (4) Continue to find the highest point outside the coverage range according to the method in step (3) until the coverage area is close enough to the coverage analysis area, that is, |S1 - S| ≤ 1e-6, or the analysis of the last grid point is completed, and then stop the analysis.

[0077] (5) In the above process, it can be ensured that the red area meets the conditions and the green area reaches the optimal coverage. On this basis, check whether the blue area meets the conditions. If not, in the interference coverage range that affects the area where the blue area does not meet the conditions, layout the interference source from high to low according to the elevation, and determine its actual coverage range through the visibility analysis method.

[0078] At the same time, check again whether the red area meets the conditions. And so on until all red areas and blue areas meet the conditions.

[0079] As a preferred embodiment, further, the visibility analysis method can be used to obtain the following contents that can be designed in the visibility area of each interference source:

[0080] First, obtain the relevant interference parameters and the maximum interference distance according to the type of interference source;

[0081] Next, construct an analysis circle with the location of the interference source as the center and the maximum interference distance as the radius, use the circumscribed square of this analysis circle as the analysis area, and read the DEM data corresponding to the analysis area.

[0082] Then, starting from the position directly in front of the interference source and using eight direction lines evenly divided in the 360-degree direction as the dividing lines, the analysis area is divided. The visibility analysis is carried out for the grid points on the direction lines according to the point-to-point visibility algorithm, and the visibility analysis is carried out for each grid point in each divided area outside the direction lines. The visibility analysis results are combined to construct the visibility area of the interference source.

[0083] In the GNSS jamming source visibility area analysis method, as shown in Figure 6 Using the DEM model data of the analysis area, with the grid point where the interference source is located as the center, eight directions are divided. The visibility on the direction lines is solved by point-to-point visibility analysis, and the visibility analysis of other areas is normalized into the visibility analysis problem of the direction lines by mathematical methods. Finally, the visibility of each grid point in the area is analyzed clearly, and the situation in the grid can be determined by the visibility of the grid point in the lower left corner.

[0084] Taking the visibility area analysis of one interference source as an example, a specific implementation method is given. When analyzing the visibility areas of multiple interference sources, the following method can be repeated multiple times:

[0085] (1) According to the type of the interference source, obtain the relevant interference parameters, and calculate the maximum interference distance of the GNSS interference source according to Equation (2):

[0086]

[0087] In the formula, d s represents the maximum interference distance of the GNSS interference source; λ is the wavelength of the GNSS interference signal; P JT is the transmission power of the GNSS interference source; G T is the transmitting antenna gain of the GNSS interference source; G R is the receiving antenna gain; P S is the actual GNSS signal power received by the receiver; M J is the power drop between the interference signal and the real signal that the receiver can tolerate, which is determined by the anti-jamming ability of the receiver itself.

[0088] (2) Taking the position of the interference source as the center of the circle and the maximum interference distance as the radius, form an analysis circle, and use the circumscribed square corresponding to this circle as the analysis area, and read in the digital elevation model (DEM) data corresponding to the analysis area;

[0089] (3) Taking the eight direction lines where the GNSS interference source is located as the boundaries, divide the analysis area into eight areas M1, M2, …, M8, as shown in Figure 4 ;

[0090] (4) For the target points on the eight direction lines, according to the point-to-point visibility algorithm, starting from the source, conduct visibility analysis. Taking one of the direction lines as an example, as shown inFigure 5 as shown, where m1, m2, …, m p are grid points on the direction line.

[0091] As a preferred embodiment, further, visibility analysis is performed on the grid points on the direction line according to the point-to-point visibility algorithm, including:

[0092] First, connect the interference source with the nearest grid point on the direction line to form a line of sight, and record the elevation angle of this nearest grid point as the current maximum elevation angle;

[0093] Then, calculate the elevation angles between the interference source and each grid point on the direction line in turn from the interference source, and compare each calculated elevation angle with the current maximum elevation angle. If the calculated elevation angle is greater than the current maximum elevation angle, it is determined that the corresponding grid point is visible, and the current maximum elevation angle is updated with the calculated elevation angle. If the calculated elevation angle is less than the current maximum elevation angle, it is determined that the corresponding grid point is not visible;

[0094] Then, obtain all visible grid points on the direction line and record them using a visibility matrix, and make corresponding markings in the analysis area through the visibility matrix.

[0095] Figure 5 As shown, the interference source m0 can be connected to the nearest grid point m1 on the direction line to form a line of sight, record the elevation angle α1 of m1, and record it as the current maximum elevation angle β. Here, the elevation angle is defined as the angle between the line of sight and the plumb line passing through the viewing point (the position of the interference source m0). The interference source position and grid point coordinates are given by the DEM model. Assume that the coordinates of the interference source m0 and any grid point m i are (x0, y0, H0) and (x i , y i , H i ), respectively, where x i , y i are the plane coordinates and H i is the elevation, then

[0096]

[0097] where α i is the elevation angle of grid point m i relative to the interference source m0.

[0098] Then, from the source outwards, calculate the elevation angles α i between the interference source m0 and each grid point m i on the direction line in turn, and compare with the recorded maximum elevation angle. If α i > β, then the grid point m i is visible, and update the maximum elevation angle, that is, let β = α i ; if α i < β, then the grid point mi Invisible, keep the maximum elevation angle β unchanged. And so on, obtain the visibility of all grid points on the direction line and record it in the visibility matrix A. Number from the upper left corner of the analysis area, and use k and j to represent the row and column numbers of the grid point m in the analysis area respectively. i If m i is visible, then A(k,j)=1; if m i is invisible, then A(k,j)=0.

[0099] As a preferred embodiment, further, perform visibility analysis on each grid point in each divided area outside the direction line, including:

[0100] First, connect the interference source and the target grid point to form a connection line, and obtain the plane coordinates of the intersection points of the connection line and the grid lines.

[0101] Next, determine the elevation value of the intersection point based on the coordinates of the two adjacent grid points of the intersection point, and obtain the elevation values of all intersection points through linear interpolation.

[0102] Then, form an equivalent direction line with the interference source, all intersection points and the target grid point.

[0103] Finally, perform visibility analysis on the equivalent direction line according to the point-to-point visibility algorithm to obtain the visibility of each grid point in each divided area, and record it using the visibility matrix.

[0104] For the area outside the direction line, it is divided into 8 blocks by 8 direction lines. As shown, for the grid points in each area, in the order from the source outwards, perform visibility analysis on each grid point m Figure 3 as follows: i′ The method is as follows:

[0105] First, connect the interference source m0 and the grid point m i′ to form a connection line, and find the plane coordinates of the intersection points of the connection line and the grid lines. The coordinates of the interference source and the grid point are given by the DEM model. The method for determining the plane coordinates of the intersection points can be described as follows:

[0106] Assume that the coordinates of m0 and m i′ are (x0,y0,H0) and (x i′ ,y i′ ,H i′ ) respectively, then the equation of the connection line between m0 and m i′ on the plane is as follows:

[0107]

[0108] where x and y represent the coordinates of the points on the connection line on the plane.

[0109] Take one of the grid lines as an example. The other grid lines are similar. Assume that the coordinates of the two grid points are (x j ,y j ,H j )、(x k ,y k ,H k ), then the corresponding equation of the grid line on the plane is as follows:

[0110]

[0111] By solving equations (4) and (5) simultaneously, we can obtain the plane coordinates of the intersection of the connecting line and the grid line. The plane coordinates of the intersection of the connecting line and other grid lines can be obtained by the above method.

[0112] Then, the elevation values ​​of the grid points adjacent to the intersection are used for linear interpolation to obtain the elevation values ​​of all intersections. q ,y q ) represents the plane coordinates of the intersection point obtained above, then the elevation value of the intersection point H q Determined by the following formula:

[0113]

[0114] Among them, (x j ,y j ,H j )、(x k ,y k ,H k ) are the coordinates of the two grid points adjacent to the intersection.

[0115] Secondly, the interference source m0, all intersection points and target grid point m i′ A direction line equivalent to the above is formed, and the visual analysis can be performed according to the visual situation analysis method of the grid points on the direction line to obtain the target grid point m i′ The visual situation.

[0116] Finally, according to the above method, the visual condition of each grid point in the eight blocks is obtained. k′ and j′ represent the grid point m, i′ The number of rows and columns in the analysis area continues to be assigned to the visual matrix A. If m i′ Visible, then A(k′,j′)=1; if m i′ If it is not visible, then A(k′,j′)=0.

[0117] By combining the visibility of the grid points on the direction lines and the grid points in the divided area outside the direction lines, the visibility of each grid in the analysis area can be obtained and represented by the visibility matrix A. The visibility in each grid is determined by the grid point in the lower left corner.

[0118] Furthermore, based on the above method, an embodiment of the present invention further provides a multi-GNSS jamming source deployment optimization system, including: a constraint condition setting module and a jamming scheme optimization module, where,

[0119] The constraint condition setting module is used to set the jamming source generation constraint conditions according to the area protection type, where the area protection type includes key jamming area type, key protection area type, and general control area type;

[0120] The jamming scheme optimization module is used to generate jamming source constraint conditions, and use each jamming source type to jointly jam the regional GNSS signals when the number of jamming sources is fixed, or use the minimum number of jamming sources to jam the regional GNSS signals when the jamming source type is fixed.

[0121] Furthermore, based on the above method, an embodiment of the present invention further provides a multi-GNSS jamming system, including: a site selection and planning module, a deployment generation module, and a real-time monitoring module, where,

[0122] The site selection and planning module is used to set the site selection and planning of the jamming source equipment in the target area according to the operation business requirements;

[0123] The deployment generation module is used to obtain the optimal jamming source equipment deployment path according to the site selection and planning results and using the above deployment optimization method;

[0124] The real-time monitoring module is used to monitor and display the working parameters of the jamming source equipment and the interference effect in the target area in real time, where the working parameters at least include the position, status, and parameters of the jamming source equipment.

[0125] The site selection and planning of the jamming source equipment can be realized by manual, semi-automatic, full-automatic and other methods for the target area to meet the simple, fast and efficient site selection and planning business requirements; according to the results of the site selection and planning, the optimal and feasible interference equipment deployment planning route can be generated in an interactive or automatic manner based on the above deployment optimization method to guide the operator to deploy on site according to the planning route; and the working parameters such as the position, status, and parameter information of the deployed satellite jamming equipment can be viewed in real time, and the interference effect, position, status, parameters, equipment type, quantity, etc. of the satellite jamming source equipment can be intuitively displayed, and the position, status, and parameters of a single or multiple satellite jamming source equipment can be adjusted in real time to meet the safety business requirements of various daily and sudden satellite interferences, maximize the savings of funds, improve the layout efficiency of small equipment, and improve the equipment supervision efficiency.

[0126] Unless otherwise specifically stated, the relative steps, numerical expressions, and values of the components and steps set forth in these embodiments do not limit the scope of the present invention.

[0127] In the present specification, the various embodiments are described in a progressive manner. Each embodiment focuses on the differences from other embodiments. For the same or similar parts among the embodiments, reference can be made to each other. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple. For the relevant parts, reference can be made to the description in the method section.

[0128] The units and method steps of the examples described in combination with the embodiments disclosed herein can be implemented by electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described according to functions in the above description. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Those of ordinary skill in the art can use different methods to implement the described functions for each specific application, but such implementation is not considered to exceed the scope of the present invention.

[0129] Those of ordinary skill in the art can understand that all or part of the steps in the above methods can be completed by instructing relevant hardware through a program, and the program can be stored in a computer-readable storage medium, such as a read-only memory, a magnetic disk, or an optical disc, etc. Optionally, all or part of the steps of the above embodiments can also be implemented using one or more integrated circuits. Correspondingly, the various modules / units in the above embodiments can be implemented in the form of hardware or in the form of software function modules. The present invention is not limited to any specific form of the combination of hardware and software.

[0130] Finally, it should be noted that the above-described embodiments are only specific implementation manners of the present invention, used to illustrate the technical solutions of the present invention, rather than to limit it. The protection scope of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that any person skilled in the art within the technical scope disclosed by the present invention can still modify the technical solutions recorded in the foregoing embodiments, or can easily think of changes, or make equivalent replacements for some of the technical features; and these modifications, changes, or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered by the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.

Claims

1. A method for optimizing the deployment of multi-GNSS jamming sources, characterized in that Comprising: Setting suppression interference source generation constraint conditions according to the regional protection type, where the regional protection type includes the key suppression area type, the key protection area type, and the general control area type. The key suppression area is set according to the application requirements with an area greater than X - fold interference coverage. The key protection area is set according to the anti - interference performance of the protected target objects in the area with an area less than Y - fold interference coverage. The general control area is set for the area outside the key suppression area and the key protection area with single - fold suppression interference, where X and Y are user - preset parameters; Based on the suppression interference source generation constraint conditions, and when the number of interference sources is fixed, using each interference source type to synergistically suppress the regional GNSS signal, or when the interference source type is fixed, using the minimum number of interference sources to suppress the regional GNSS signal. When using the minimum number of interference sources to suppress the regional GNSS signal with the interference source type fixed, it includes: First, for the task target area, the fixed number of interference sources set is n, and the number of DEM grid points m is obtained by reading the digital elevation model (DEM) of the area; find the highest elevation point among the m DEM grid points, use this highest elevation point as the interference source layout position point, obtain the coverage range of the position point using the visibility analysis method, and take this coverage range of the position point as the original coverage range. And when obtaining the coverage range of the position point using the visibility analysis method, for the case where the coverage range of the position point does not meet the key protection area, find the second - highest elevation point in the grid points as the interference source layout position point, and meet the key protection area by re - obtaining the coverage range of the position point; then, find the new highest elevation point outside the original coverage range, obtain the new coverage range of the position point through the visibility analysis method, accumulate the new coverage range of the position point and the original coverage range, take the grid points with the accumulated area greater than the area of the original position point coverage range as the interference source layout position points, and take the accumulation result as the original coverage range for the next round of iteration to find the new highest elevation point until the area of the original coverage range obtained in the iteration and the area of the task target area meet the preset proximity threshold, and complete the optimal deployment of the suppression interference source.

2. The method for optimizing the deployment of multi-GNSS jamming sources according to claim 1, characterized in that When the number of interference sources is fixed, using each interference source type to synergistically suppress the regional GNSS signal includes the following content: First, for the task target area, the fixed number of interference sources set is n, and the number of DEM grid points m is obtained by reading the digital elevation model DEM of the area. C deployment schemes are formed by selecting n points from the m DEM grid points. n m kinds of deployment schemes; Next, for each deployment plan, use the visibility analysis method to obtain the visibility of each interference source, and calculate the number of visible interference sources for each grid point one by one; obtain the deployment plan that simultaneously meets the suppression interference source generation constraint conditions of the key suppression area and the key protection area according to the maximum or minimum number of visible interference sources for each point, and calculate the combined interference area of all interference sources under the constraint conditions of the general control area, and calculate the task area coverage rate based on this combined interference area and the total area of the task target area; Then, output the optimal deployment plan of the suppression interference source according to the task area coverage rate.

3. The method for optimizing the deployment of multi-GNSS jamming sources according to claim 1, characterized in that When obtaining the coverage range of location points using the visible area analysis method, it also includes: in the case where the coverage range of location points does not meet the key suppression area, the location points are arranged as interference sources in descending order of elevation. By re-obtaining the coverage range of location points, a layout plan of interference sources that meets the conditions of the key suppression area is determined, and at the same time, it is checked whether the coverage range of location points meets the key protection area until all key protection areas and key suppression areas are within the coverage range of location points.

4. The method for optimizing the deployment of multi-GNSS jamming sources according to claim 1 or 2, characterized in that When obtaining the visible areas of each interference source using the visible area analysis method, it includes: First, relevant interference parameters and the maximum interference distance are obtained according to the type of interference source. Next, an analysis circle is constructed with the location of the interference source as the center and the maximum interference distance as the radius. The circumscribed square corresponding to this analysis circle is used as the analysis area, and the digital elevation model (DEM) data corresponding to the analysis area is read. Then, starting from the position directly in front of the interference source and using 8 direction lines evenly divided in the 360-degree direction as the dividing lines, the analysis area is divided. The visibility analysis of the grid points on the direction lines is performed according to the point-to-point visibility algorithm, and the visibility analysis of each grid point in each divided area outside the direction lines is performed. The visibility analysis results are combined to construct the visible area of the interference source.

5. The method for optimizing the deployment of multi-GNSS jamming sources according to claim 4, characterized in that The visibility analysis of the grid points on the direction lines according to the point-to-point visibility algorithm includes: First, a line of sight is formed by connecting the interference source and the nearest grid point on the direction line, and the elevation angle of this nearest grid point is recorded as the current maximum elevation angle. Next, the elevation angles between the interference source and each grid point on the direction line are calculated sequentially outward from the interference source, and the calculated elevation angles are compared with the current maximum elevation angle. If the calculated elevation angle is greater than the current maximum elevation angle, the corresponding grid point is determined to be visible, and the current maximum elevation angle is updated using the calculated elevation angle. If the calculated elevation angle is less than the current maximum elevation angle, the corresponding grid point is determined to be invisible. Then, all visible grid points on the direction line are obtained and recorded using the visibility matrix, and corresponding markings are made in the analysis area through the visibility matrix.

6. The method for optimizing the deployment of multi-GNSS jamming sources according to claim 4, characterized in that The visibility analysis of each grid point in each divided area outside the direction line includes: First, a connection line is formed by connecting the interference source and the target grid point, and the plane coordinates of the intersection points of the connection line and the grid lines are obtained. Next, the elevation values of the intersection points are determined based on the coordinates of the two adjacent grid points of the intersection points, and all elevation values of the intersection points are obtained through linear interpolation. Then, an equivalent direction line is formed by the interference source, all intersection points, and the target grid point. Finally, visibility analysis is performed on the equivalent direction line according to the point-to-point visibility algorithm to obtain the visibility of each grid point in each divided area, and it is recorded using the visibility matrix.

7. A system for optimizing the deployment of multi-GNSS jamming sources, characterized in that Implemented based on the method described in claim 1, it includes: a constraint condition setting module and an interference scheme optimization module, where The constraint condition setting module is used to set the generation constraints of suppression interference sources according to the type of area protection. Among them, the type of area protection includes the key suppression area type, the key protection area type, and the general control area type. The interference scheme optimization module is used to generate constraint conditions based on the jamming sources, and when the number of jamming sources is fixed, use each type of jamming source to perform collaborative jamming on the regional GNSS signals, or when the type of jamming source is fixed, use the minimum number of jamming sources to jam the regional GNSS signals.

8. A multi-GNSS jamming system, characterized in that Implemented based on the method described in claim 1, including: a site selection and planning module, a deployment generation module, and a real-time monitoring module, where The site selection and planning module is used to set the site selection and planning of the jamming source equipment in the target area according to the operation service requirements; The deployment generation module is used to obtain the optimal deployment path of the jamming source equipment according to the site selection and planning results and by using the method described in claim 1; The real-time monitoring module is used to monitor and display the working parameters of the jamming source equipment and the interference effect of the target area in real time, where the working parameters at least include the position, status, and parameters of the jamming source equipment.

Images