A high-throughput satellite frequency planning method based on backtracking search

Abstract

The present invention discloses a high-throughput satellite frequency plan arrangement method based on backtracking search, including inputting the total available frequency band, the number of regions, the number of beams and beam width in each region, the required bandwidth of a single beam in each region and the corresponding number combination, the beam topology relationship, the total bandwidth of the system, whether the gateway uses dual polarization, and whether each region uses a multi-feed synthetic antenna; in order from large beams to small beams, from more frequency band taboos to fewer frequency band taboos, frequency bands are allocated to each beam in turn, and when a beam has no optional frequency band, it traces back to the previous beam adjacent to it that has been allocated and re-allocates other frequency bands to the beam, and finally outputs the frequency bands allocated to each beam and the corresponding gateway. In the present invention, a frequency plan arrangement scheme that meets the feasibility constraints of multiple transponders and antennas can be obtained, and it gets rid of the dependence on manual arrangement, and the time period is less than 10 minutes, which is conducive to the optimization of the scheme and expansion to large-scale non-equal width beam scenarios.

Description

Technical Field

[0001] The present invention relates to the technical field of high-throughput satellite frequency plan arrangement methods, and in particular to a high-throughput satellite frequency plan arrangement method based on backtracking search. Background Art

[0002] After decades of development, the service capabilities of high-throughput satellites have been significantly improved. However, improper allocation of on-board resources will prevent us from fully utilizing their various advantages. As one of the most valuable resources on board, the planning and management of frequency resources is a key issue in the design and optimization of high-throughput satellites. If the frequency is not well managed during the overall design and planning of high-throughput satellites, it will cause the satellite system's carrier-to-interference ratio, effective capacity and other performance indicators to deteriorate, thus affecting user experience and failing to fully utilize satellite resources.

[0003] The frequency plan arrangement of traditional high-throughput communication satellite transponders mostly adopts a uniform multi-color multiplexing allocation scheme. However, once the user demand between different beams is unbalanced, it will lead to "uneven busy and idle" resources, affecting the satellite resource utilization and user service quality. Therefore, it is necessary to study the flexible arrangement method of on-board frequency plan to effectively alleviate the "uneven busy and idle" phenomenon between beams and improve the system's matching rate with user capacity requirements.

[0004] When it comes to flexibly arranging frequency plans for high-throughput satellites, existing methods have two important drawbacks:

[0005] (1) The feasibility constraints of transponders and antennas are not fully considered, such as: co-frequency and co-polarization beams are not adjacent; all gateway stations must have full two-polarization frequency resources, and the upper sideband and lower sideband must be aligned in segments as much as possible; the uplink and downlink frequencies must be frequency-changed as much as possible, and the frequency-changing relationships must not exceed 3; the minimum protection bandwidth between beams is 10% of the bandwidth of adjacent wide beams; if the antenna is a multi-feed synthesis method, the beams synthesized by multi-feed must be distributed "co-polarized by column"; the co-frequency interval of small beams adjacent to large beams must be greater than the width of a large beam; the frequencies of adjacent large beams must be different and as far apart as possible. Therefore, the generated frequency plan cannot be used in the actual model.

[0006] (2) The practical frequency plan for the model relies too much on manual arrangement, which takes a long time and cannot form a solution optimization. In addition, it is difficult to obtain a feasible solution when large-scale non-uniform width beams are used. Summary of the invention

[0007] The purpose of the present invention is to solve the above problems and propose a high-throughput satellite frequency plan arrangement method based on backtracking search for scenarios with large-scale non-equal-width beams and uneven bandwidth requirements of each beam, so as to improve satellite design capabilities and achieve better satellite resource utilization and user service quality.

[0008] In order to achieve the above object, the present invention adopts the following technical solutions:

[0009] A high-throughput satellite frequency plan arrangement method based on back-tracking search includes inputting the total available frequency band, the number of regions, the number of beams and beam widths in each region, the required bandwidth of a single beam in each region and the corresponding number combination, the beam topology relationship, the total bandwidth of the system, whether the gateway uses dual polarization, and whether each region uses a multi-feed synthetic antenna;

[0010] Frequency bands are allocated to each beam in order from large beams to small beams and from more frequency band taboos to fewer frequency band taboos. When a beam has no optional frequency band, it traces back to the previous adjacent beam that has completed allocation and reallocates other frequency bands to the beam. Finally, the frequency bands allocated to each beam and the corresponding gateway stations are output.

[0011] Preferably, the steps include:

[0012] S1. Initialize input parameters:

[0013] Manually input the user's downlink available frequency band [f1, f2], where f1 and f2 are positive, in GHz, and f2 is greater than f1;

[0014] The number of regions m is manually input, where m is a positive integer, the beam width in a single region is the same, and the set of all region numbers M is {1, 2, …, m};

[0015] Manually input the number of beams in each area vector N, the i-th vector element N i represents the number of beams in region i, N i is a positive integer and i∈M; the set of all beam numbers Q is {1, 2, …, q},

[0016] Manually input the beam width vector B of each region, the i-th vector element B i represents the beam width in region i, B i is positive, in degrees, and i∈M;

[0017] Manually input the required bandwidth of a single beam in each area and the corresponding number of cell arrays S, which contains a total of m cells; the i-th cell is a matrix S (i) ,i∈M,S (i) It contains u rows and 2 columns, where u is a positive integer and different i can correspond to different u, which means that area i contains u types of single-beam required bandwidths; the matrix element in the vth row and the first column is Represents the required bandwidth value of the vth single beam in region i, the matrix element in the vth row and the second column represents the number of beams corresponding to the vth single-beam required bandwidth value in region i, i is a positive integer and v≤u;

[0018] S2, divide the total frequency band resources of the system into a set of sub-bands Δ according to the beam bandwidth requirements; the number of available frequency bands [f1, f2] δ is the total system bandwidth divided by the available frequency band bandwidth, that is, If it is dual polarization, the number of frequency bands for left and right polarization is δ / 2 respectively;

[0019] S3, initializing the polarization definition library, the frequency band taboo library, the allocated beam number library and the frequency band total resource library;

[0020] S4. Sort the regions according to beam width from large to small. The region numbers are λ1, λ2, Lλ. m ;

[0021] S5. Use the most central beam in area λ1 as the beam to be allocated;

[0022] S6. From the matrix Randomly select a row, the row number is θ, and the element value of the θth row and the first column is used as the bandwidth requirement of the beam to be allocated;

[0023] S7. Randomly select a sub-band that meets the following conditions from the total frequency band resource library and allocate it to the beam to be allocated: the bandwidth of the sub-band is equal to the bandwidth requirement of the beam to be allocated; if the element value in the polarization definition library corresponding to the beam to be allocated is not 0, the polarization value of the sub-band must be the same as the element value, otherwise there is no requirement;

[0024] S8, adding the beam to be allocated to the allocated beam number library;

[0025] S9, taking the beams in region λ1 that are adjacent to the beam to be allocated and do not belong to the allocated beam number library as a beam set to be allocated;

[0026] S10, re-determining the first beam in the set of beams to be allocated as the beam to be allocated;

[0027] S11, according to the method of the above steps S6-S8, allocate sub-bands for the beams to be allocated and update the allocated beam number library and matrix If there is no sub-band that meets the conditions in the frequency band total resource library, the frequency band taboo library, and the polarization definition library, then trace back to the beam that is adjacent to the beam to be allocated and belongs to the allocated beam number library, determine it as the new beam to be allocated, and restart the allocation according to the method of steps S6-S8;

[0028] S12, adding beams in region λ1 that are adjacent to the beam to be allocated and do not belong to the allocated beam number library or the beam set to be allocated to the beam set to be allocated, and deleting the beam to be allocated from the beam set to be allocated, to form a new beam set to be allocated;

[0029] S13, repeating the method of steps S10-S12 until sub-band allocation is completed for all beams in region λ1;

[0030] S14, according to the method of steps S5-S13, the area λ2, Lλ m The allocation is completed, wherein step S5 needs to be adjusted as follows: the area λ2, ..., or λ m The beam with the most taboos in the inner frequency band is used as the beam to be allocated;

[0031] S15. Output the sub-frequency band corresponding to each beam as the final frequency plan arrangement scheme.

[0032] Preferably, in step S1, the topological relationship matrix T between the beams is manually input, and the matrix element T xy Indicates whether the beam numbered x is adjacent to the beam numbered y. If so, T xy =1, if not adjacent, then T xy =0; if the beam numbered x and the beam numbered y have different widths, they are considered adjacent if the interval between them is less than the maximum beam width;

[0033] Manually input the total bandwidth of the system G, where G is a positive number in GHz;

[0034] Manual input of whether the gateway uses dual polarization α, if yes, then α=1, if no, then α=0;

[0035] Manual input whether the user uses dual polarization β, if yes, β = 1, if no, β = 0;

[0036] Manually input the vector Y of whether each area uses a multi-feed synthetic antenna, the i-th vector element Y i Indicates whether region i uses a multi-feed synthetic antenna. If so, Y i =1, if not used, then Y i =0.

[0037] Preferably, in step S2, sub-frequency bands with bandwidths equal to the maximum value of the beam bandwidth requirement are first divided in sequence from the minimum frequency points of the δ available frequency bands. If the number ε of beams with bandwidth requirements equal to the maximum value is less than δ, two values ​​whose sum is closest to the maximum value of the beam bandwidth requirement are selected from the beam bandwidth requirement values, and then the δ-ε available frequency bands are divided in sequence from the minimum frequency points into sub-frequency bands with bandwidths equal to the two values, completing the first round of division of the δ available frequency bands; the δ available frequency bands and the beam bandwidth requirement matrix are adjusted, and multiple rounds of division are completed in sequence according to the above method until the total frequency band resources of the system are divided into a set of sub-frequency bands that meet all beam bandwidth requirements;

[0038] Δ contains c rows and 3 columns, where c is a natural number. Each row corresponds to an optional sub-band. The matrix element Δ in the dth row and the first column isd1 Represents the minimum frequency of the dth optional sub-band, in GHz, and the matrix element Δ in the dth row and second column d2 Represents the maximum frequency of the dth optional sub-band, in GHz, and the matrix element Δ in the dth row and third column d3 Represents the polarization value of the dth optional sub-band. If it is left polarization, Δ d3 =1, if it is right polarization, Δ d3 =-1, d is a natural number and d≤c.

[0039] Preferably, the polarization definition library in step S3 is a vector P, and the jth vector element P j Represents the polarization value corresponding to the beam numbered j, j∈Q, if it is left polarization, P j =1, if it is right polarization, P j =-1, when initialized, P j =0;

[0040] The frequency band taboo library is a cell array F, which contains q cells in total; the jth cell is a matrix F (j) ,j∈Q,F (j) Contains a rows and 2 columns, where a is a natural number and different j can correspond to different a, which means that the beam numbered j contains a taboo frequency bands, and the matrix element in the bth row and the first column is Represents the minimum frequency of the b-th taboo frequency band of the beam numbered j, the matrix element in the b-th row and the second column represents the maximum frequency of the bth taboo frequency band of the beam numbered j, where b is a natural number and b≤a; during initialization, F (j) =Φ,j∈Q;

[0041] The allocated beam number library is a vector D, whose elements are the beam numbers that have been allocated in sequence. When initialized, D = Φ;

[0042] The total frequency band resource library is a matrix R, and when initialized, R=Δ.

[0043] Preferably, in step S8, the matrix The value of the element in the θth row and the 2nd column is reduced by 1 and stored again. If the result is equal to 0, the matrix Delete the θth row; delete the allocated sub-band from the total frequency band resource library; expand the allocated sub-band by 10% of the protection bandwidth and add it to the frequency band taboo library of all other beams adjacent to the beam to be allocated; if The elements in the polarization definition library corresponding to the beams in the same column as the beam to be allocated in region λ1 are reassigned to be equal to the polarization value of the allocated sub-frequency band, and the elements in the polarization definition library corresponding to the beams in different columns are reassigned to be opposite to the polarization value of the allocated sub-frequency band.

[0044] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0045] 1. This application proposes a sub-band division method that matches the beam bandwidth requirements. According to the idea of ​​beam bandwidth requirements from large to small, multiple rounds of division, and strict alignment of each round, the division of the total frequency band resources of the system is matched with the beam bandwidth requirements, while the two polarization frequency resources are filled as much as possible, and the upper sideband and the lower sideband are aligned in segments as much as possible.

[0046] 2. The present invention proposes a frequency plan arrangement method based on a backtracking search algorithm, which determines two beams of different widths whose interval is less than the maximum beam width as adjacent, sets a "polarization definition library" and a "frequency band taboo library", considers protection bandwidth and other methods, and solves a frequency plan arrangement scheme that meets the feasibility constraints of multiple repeaters and antennas.

[0047] 3. The present invention proposes a frequency plan arrangement method based on a backtracking search algorithm. When a sub-frequency band cannot be allocated for a certain beam, a backtracking algorithm is used to comprehensively search the solution space and automatically generate a feasible solution. At the same time, the beam that is adjacent to the beam that cannot be allocated and belongs to the "allocated beam number library" is backtracked to be reallocated instead of the previous beam, which can narrow the feasible solution space and improve the efficiency of the algorithm. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] Figure 1 A flowchart of a high-throughput satellite frequency plan arrangement method based on backtracking search according to an embodiment of the present invention is shown;

[0049] Figure 2 A schematic diagram of the topological relationship between beams provided according to an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0050] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0051] See also Figure 1-2 , the present invention provides a technical solution:

[0052] A high-throughput satellite frequency plan arrangement method based on back-tracking search includes inputting the total available frequency band, the number of regions, the number of beams and beam widths in each region, the required bandwidth of a single beam in each region and the corresponding number combination, the beam topology relationship, the total bandwidth of the system, whether the gateway uses dual polarization, and whether each region uses a multi-feed synthetic antenna;

[0053] Frequency bands are allocated to each beam in order from large beams to small beams and from more frequency band taboos to fewer frequency band taboos. When a beam has no optional frequency band, it traces back to the previous adjacent beam that has completed allocation and reallocates other frequency bands to the beam. Finally, the frequency bands allocated to each beam and the corresponding gateway stations are output.

[0054] Specifically, Figure 1 As shown, a high-throughput satellite frequency plan arrangement method based on back-tracking search includes the following steps:

[0055] S1. Initialize input parameters:

[0056] Manually input the user's downlink available frequency band [f1, f2], where f1 and f2 are positive, in GHz, and f2 is greater than f1;

[0057] The number of regions m is manually input, where m is a positive integer, the beam width in a single region is the same, and the set of all region numbers M is {1, 2, …, m};

[0058] Manually input the number of beams in each area vector N, the i-th vector element N i represents the number of beams in region i, N i is a positive integer and i∈M; the set of all beam numbers Q is {1, 2, …, q},

[0059] Manually input the beam width vector B of each region, the i-th vector element B i represents the beam width in region i, B i is positive, in degrees, and i∈M;

[0060] Manually input the required bandwidth of a single beam in each area and the corresponding number of cell arrays S, which contains a total of m cells; the i-th cell is a matrix S (i) ,i∈M,S (i) It contains u rows and 2 columns, where u is a positive integer and different i can correspond to different u, which means that area i contains u types of single-beam required bandwidths; the matrix element in the vth row and the first column is Represents the required bandwidth value of the vth single beam in region i, the matrix element in the vth row and the second column represents the number of beams corresponding to the vth single-beam required bandwidth value in region i, i is a positive integer and v≤u;

[0061] In step S1, the topological relationship matrix T between the beams is manually input, and the matrix element T xy Indicates whether the beam numbered x is adjacent to the beam numbered y. If so, T xy =1, if not adjacent, then T xy=0; if the beam numbered x and the beam numbered y have different widths, they are considered adjacent if the interval between them is less than the maximum beam width;

[0062] Manually input the total bandwidth of the system G, where G is a positive number in GHz;

[0063] Manual input of whether the gateway uses dual polarization α, if yes, then α=1, if no, then α=0;

[0064] Manual input whether the user uses dual polarization β, if yes, β = 1, if no, β = 0;

[0065] Manually input the vector Y of whether each area uses a multi-feed synthetic antenna, the i-th vector element Y i Indicates whether region i uses a multi-feed synthetic antenna. If so, Y i =1, if not used, then Y i =0.

[0066] S2, divide the total frequency band resources of the system into a set of sub-bands Δ according to the beam bandwidth requirements; the number of available frequency bands [f1, f2] δ is the total system bandwidth divided by the available frequency band bandwidth, that is, If it is dual polarization, the number of frequency bands for left and right polarization is δ / 2 respectively;

[0067] In step S2, first, sub-frequency bands with bandwidths equal to the maximum value of the beam bandwidth requirement are divided in sequence from the minimum frequency points of the δ available frequency bands. If the number of beams ε with bandwidth requirements equal to the maximum value is less than δ, two values ​​whose sum is closest to the maximum value of the beam bandwidth requirement are selected from the beam bandwidth requirement values, and then the δ-ε available frequency bands are divided in sequence from the minimum frequency points into sub-frequency bands with bandwidths equal to the two values, completing the first round of division of the δ available frequency bands; the δ available frequency bands and the beam bandwidth requirement matrix are adjusted, and multiple rounds of division are completed in sequence according to the above method until the total frequency band resources of the system are divided into a set of sub-frequency bands that meet all beam bandwidth requirements;

[0068] Δ contains c rows and 3 columns, where c is a natural number. Each row corresponds to an optional sub-band. The matrix element Δ in the dth row and the first column is d1 Represents the minimum frequency of the dth optional sub-band, in GHz, and the matrix element Δ in the dth row and second column d2 Represents the maximum frequency of the dth optional sub-band, in GHz, and the matrix element Δ in the dth row and third column d3 Represents the polarization value of the dth optional sub-band. If it is left polarization, Δ d3 =1, if it is right polarization, Δ d3 =-1, d is a natural number and d≤c.

[0069] S3, initializing the polarization definition library, the frequency band taboo library, the allocated beam number library and the frequency band total resource library;

[0070] The polarization definition library in step S3 is a vector P, and the jth vector element P j Represents the polarization value corresponding to the beam numbered j, j∈Q, if it is left polarization, P j =1, if it is right polarization, P j =-1, when initialized, P j =0;

[0071] The frequency band taboo library is a cell array F, which contains q cells in total; the jth cell is a matrix F (j) ,j∈Q,F (j) Contains a rows and 2 columns, where a is a natural number and different j can correspond to different a, which means that the beam numbered j contains a taboo frequency bands, and the matrix element in the bth row and the first column is Represents the minimum frequency of the b-th taboo frequency band of the beam numbered j, the matrix element in the b-th row and the second column represents the maximum frequency of the bth taboo frequency band of the beam numbered j, where b is a natural number and b≤a; during initialization, F (j) =Φ,j∈Q;

[0072] The allocated beam number library is a vector D, whose elements are the beam numbers that have been allocated in sequence. When initialized, D = Φ;

[0073] The total frequency band resource library is a matrix R, and when initialized, R=Δ.

[0074] S4. Sort the regions according to beam width from large to small. The region numbers are λ1, λ2, Lλ. m ;

[0075] S5. Use the most central beam in area λ1 as the beam to be allocated;

[0076] S6. From the matrix Randomly select a row, the row number is θ, and the element value of the θth row and the first column is used as the bandwidth requirement of the beam to be allocated;

[0077] S7. Randomly select a sub-band that meets the following conditions from the total frequency band resource library and allocate it to the beam to be allocated: the bandwidth of the sub-band is equal to the bandwidth requirement of the beam to be allocated; if the element value in the polarization definition library corresponding to the beam to be allocated is not 0, the polarization value of the sub-band must be the same as the element value, otherwise there is no requirement;

[0078] S8, adding the beam to be allocated to the allocated beam number library;

[0079] In step S8, the matrix The value of the element in the θth row and the 2nd column is reduced by 1 and stored again. If the result is equal to 0, the matrix Delete the θth row; delete the allocated sub-band from the total frequency band resource library; expand the allocated sub-band by 10% of the protection bandwidth and add it to the frequency band taboo library of all other beams adjacent to the beam to be allocated; if The elements in the polarization definition library corresponding to the beams in the same column as the beam to be allocated in region λ1 are reassigned to be equal to the polarization value of the allocated sub-frequency band, and the elements in the polarization definition library corresponding to the beams in different columns are reassigned to be opposite to the polarization value of the allocated sub-frequency band.

[0080] S9, taking the beams in region λ1 that are adjacent to the beam to be allocated and do not belong to the allocated beam number library as a beam set to be allocated;

[0081] S10, re-determining the first beam in the set of beams to be allocated as the beam to be allocated;

[0082] S11, according to the method of the above steps S6-S8, allocate sub-bands for the beams to be allocated and update the allocated beam number library and matrix If there is no sub-band that meets the conditions in the frequency band total resource library, the frequency band taboo library, and the polarization definition library, then trace back to the beam that is adjacent to the beam to be allocated and belongs to the allocated beam number library, determine it as the new beam to be allocated, and restart the allocation according to the method of steps S6-S8;

[0083] S12, adding beams in region λ1 that are adjacent to the beam to be allocated and do not belong to the allocated beam number library or the beam set to be allocated to the beam set to be allocated, and deleting the beam to be allocated from the beam set to be allocated, to form a new beam set to be allocated;

[0084] S13, repeating the method of steps S10-S12 until sub-band allocation is completed for all beams in region λ1;

[0085] S14, according to the method of steps S5-S13, the area λ2, Lλ m The allocation is completed, wherein step S5 needs to be adjusted as follows: the area λ2, ..., or λ m The beam with the most taboos in the inner frequency band is used as the beam to be allocated;

[0086] S15. Output the sub-frequency band corresponding to each beam as the final frequency plan arrangement scheme.

[0087] Embodiment 1:

[0088] The input parameters initialized in step S1 are as follows:

[0089] The available frequency band for user downlink is [17.7, 20.2], in GHz;

[0090] The number of regions m is 4;

[0091] The vector of the number of beams in each region N = [65 8 9 7];

[0092] The beam width vector of each region B = [0.6 1.0 1.3 2.0];

[0093] The required bandwidth of a single beam in each area and the corresponding number cell array S are:

[0094] S (2) =[200 8], S (4) = [100 7];

[0095] The topological relationship between beams is as follows: Figure 2 As shown by Figure 2 The matrix T can be obtained;

[0096] The total system bandwidth G is 10GHz;

[0097] The gateway uses dual polarization, then α = 1;

[0098] If the user uses dual polarization, β = 1;

[0099] The vector Y=[0 0 0 0] indicating whether each area uses a multi-feed synthetic antenna.

[0100] In step S2, the total frequency band resources of the system are divided into a sub-frequency band set Δ according to the beam bandwidth requirement:

[0101]

[0102]

[0103]

[0104] In steps S3-S14, a computer program is written according to the algorithm flow to complete the automatic calculation.

[0105] In step S15, the sub-frequency bands corresponding to each beam are output as the final frequency plan arrangement scheme, as shown in Table 1.

[0106] Table 1 Frequency allocation plan for users

[0107]

[0108]

[0109] The method of the present invention can obtain a frequency plan arrangement scheme that meets the feasibility constraints of various repeaters and antennas, and gets rid of the dependence on manual arrangement. The time cycle is less than 10 minutes, which is conducive to scheme optimization and expansion to large-scale non-uniform width beam scenarios.

[0110] The present invention proposes a frequency plan arrangement method based on a backtracking search algorithm. When a sub-frequency band cannot be allocated for a certain beam, a backtracking algorithm is used to comprehensively search the solution space and automatically generate a feasible solution. At the same time, the beam that is adjacent to the beam that cannot be allocated and belongs to the "allocated beam number library" is backtracked to be reallocated instead of the previous beam, which can narrow the feasible solution space and improve the efficiency of the algorithm.

[0111] The above description of the embodiments enables those skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A high-throughput satellite frequency plan arrangement method based on backtracking search, characterized in that: Including inputting the total available frequency band, the number of regions, the number and width of beams in each region, the required bandwidth of a single beam in each region and the corresponding number combination, the beam topology relationship, the total bandwidth of the system, whether the gateway uses dual polarization, and whether each region uses a multi-feed synthetic antenna; Allocate frequency bands to each beam in order from large beams to small beams and from those with more frequency band taboos to those with fewer frequency band taboos. When a beam has no optional frequency bands, trace back to the previous adjacent beam that has been allocated and re-allocate other frequency bands to the beam. Finally, output the frequency bands allocated to each beam and the corresponding gateways. The following steps are involved: S1. Initialize input parameters: Manually input the user's downlink available frequency band [f1, f2], where f1 and f2 are positive, in GHz, and f2 is greater than f1; The number of regions m is manually input, where m is a positive integer, the beam width in a single region is the same, and the set of all region numbers M is {1, 2, …, m}; Manually input the number of beams in each area vector N, the i-th vector element N i represents the number of beams in region i, N i is a positive integer and i∈M; the set of all beam numbers Q is {1, 2, …, q}, Manually input the beam width vector B of each region, the i-th vector element B i represents the beam width in region i, B i is positive, in degrees, and i∈M; Manually input the required bandwidth of a single beam in each area and the corresponding number of cell arrays S, which contains a total of m cells; the i-th cell is a matrix S (i) ,i∈M,S (i) It contains u rows and 2 columns, where u is a positive integer and different i can correspond to different u, which means that area i contains u types of single-beam required bandwidths; the matrix element in the vth row and the first column is Represents the required bandwidth value of the vth single beam in region i, the matrix element in the vth row and the second column represents the number of beams corresponding to the vth single-beam required bandwidth value in region i, i is a positive integer and v≤u; S2, divide the total frequency band resources of the system into a set of sub-bands Δ according to the beam bandwidth requirements; the number of available frequency bands [f1, f2] δ is the total system bandwidth divided by the available frequency band bandwidth, that is, If it is dual polarization, the number of frequency bands for left and right polarization is δ / 2 respectively; S3, initializing the polarization definition library, the frequency band taboo library, the allocated beam number library and the frequency band total resource library; S4. Sort the regions according to the beam width from large to small, and the region numbers are λ1, λ2, …λ m ; S5. Use the most central beam in area λ1 as the beam to be allocated; S6. From the matrix Randomly select a row, the row number is θ, and the element value of the θth row and the first column is used as the bandwidth requirement of the beam to be allocated; S7. Randomly select a sub-band that meets the following conditions from the total frequency band resource library and allocate it to the beam to be allocated: the bandwidth of the sub-band is equal to the bandwidth requirement of the beam to be allocated; if the element value in the polarization definition library corresponding to the beam to be allocated is not 0, the polarization value of the sub-band must be the same as the element value, otherwise there is no requirement; S8, adding the beam to be allocated to the allocated beam number library; S9, taking the beams in region λ1 that are adjacent to the beam to be allocated and do not belong to the allocated beam number library as a beam set to be allocated; S10, re-determining the first beam in the set of beams to be allocated as the beam to be allocated; S11, according to the method of the above steps S6-S8, allocate sub-bands for the beams to be allocated and update the allocated beam number library and matrix If there is no sub-band that meets the conditions in the frequency band total resource library, the frequency band taboo library, and the polarization definition library, then trace back to the beam that is adjacent to the beam to be allocated and belongs to the allocated beam number library, determine it as the new beam to be allocated, and restart the allocation according to the method of steps S6-S8; S12, adding beams in region λ1 that are adjacent to the beam to be allocated and do not belong to the allocated beam number library or the beam set to be allocated to the beam set to be allocated, and deleting the beam to be allocated from the beam set to be allocated, to form a new beam set to be allocated; S13, repeating the method of steps S10-S12 until sub-band allocation is completed for all beams in region λ1; S14, according to the method of steps S5-S13, the regions λ2, ...λ m The allocation is completed, wherein step S5 needs to be adjusted as follows: the area λ2, ..., or λ m The beam with the most taboos in the inner frequency band is used as the beam to be allocated; S15. Output the sub-frequency band corresponding to each beam as the final frequency plan arrangement scheme.

2. The high-throughput satellite frequency plan arrangement method based on backtracking search according to claim 1 is characterized in that: In step S1, the topological relationship matrix T between the beams is manually input, and the matrix element T xy Indicates whether the beam numbered x is adjacent to the beam numbered y. If so, T xy =1, if not adjacent, then T xy = 0; if the beam numbered x and the beam numbered y have different widths, they are considered adjacent if the interval between them is less than the maximum beam width; Manually input the total bandwidth of the system G, where G is a positive number in GHz; Manually input whether the gateway uses dual polarization α, if yes, then α=1, if no, then α=0; Manual input of whether the user uses dual polarization β, if yes, β = 1, if no, β = 0; Manually input the vector Y of whether each area uses a multi-feed synthetic antenna, the i-th vector element Y i Indicates whether region i uses a multi-feed synthetic antenna. If so, Y i =1, if not used, then Y i =0.

3. The high-throughput satellite frequency plan arrangement method based on backtracking search according to claim 1 is characterized in that: In step S2, first, sub-frequency bands with bandwidths equal to the maximum value of the beam bandwidth requirement are divided in sequence from the minimum frequency points of the δ available frequency bands. If the number of beams ε with bandwidth requirements equal to the maximum value is less than δ, two values ​​whose sum is closest to the maximum value of the beam bandwidth requirement are selected from the beam bandwidth requirement values, and then the δ-ε available frequency bands are divided in sequence from the minimum frequency points into sub-frequency bands with bandwidths equal to the two values, completing the first round of division of the δ available frequency bands; the δ available frequency bands and the beam bandwidth requirement matrix are adjusted, and multiple rounds of division are completed in sequence according to the above method until the total frequency band resources of the system are divided into a set of sub-frequency bands that meet all beam bandwidth requirements; Δ contains c rows and 3 columns, where c is a natural number. Each row corresponds to an optional sub-band. The matrix element Δ in the dth row and the first column is d1 Represents the minimum frequency of the dth optional sub-band, in GHz, and the matrix element Δ in the dth row and second column d2 Represents the maximum frequency of the dth optional sub-band, in GHz, and the matrix element Δ in the dth row and third column d3 Represents the polarization value of the dth optional sub-band. If it is left polarization, Δ d3 =1, if it is right polarization, Δ d3 =-1, d is a natural number and d≤c.

4. The high-throughput satellite frequency plan arrangement method based on backtracking search according to claim 1, characterized in that: The polarization definition library in step S3 is a vector P, and the jth vector element P j Represents the polarization value corresponding to the beam numbered j, j∈Q, if it is left polarization, P j =1, if it is right polarization, P j =-1, when initialized, P j =0; The frequency band taboo library is a cell array F, which contains q cells in total; the jth cell is a matrix F (j) ,j∈Q,F (j) Contains a rows and 2 columns, where a is a natural number and different j can correspond to different a, which means that the beam numbered j contains a taboo frequency bands, and the matrix element in the bth row and the first column is Represents the minimum frequency of the b-th taboo frequency band of the beam numbered j, the matrix element in the b-th row and the second column represents the maximum frequency of the bth taboo frequency band of the beam numbered j, where b is a natural number and b≤a; during initialization, F (j) =Φ,j∈Q; The allocated beam number library is a vector D, whose elements are the beam numbers that have been allocated in sequence. When initialized, D = Φ; The total frequency band resource library is a matrix R, and when initialized, R=Δ.

5. The high-throughput satellite frequency plan arrangement method based on backtracking search according to claim 1, characterized in that: In step S8, the matrix The value of the element in the θth row and the 2nd column is reduced by 1 and stored again. If the result is equal to 0, the matrix Delete the θth row; delete the allocated sub-band from the total frequency band resource library; expand the allocated sub-band by 10% of the protection bandwidth and add it to the frequency band taboo library of all other beams adjacent to the beam to be allocated; if Y λ1 =1, the elements in the polarization definition library corresponding to the beams in the same column as the beam to be allocated in region λ1 are reassigned to be equal to the polarization value of the allocated sub-frequency band, and the elements in the polarization definition library corresponding to the beams not in the same column are reassigned to be opposite to the polarization value of the allocated sub-frequency band.

Images