A high-throughput satellite frequency allocation method based on carrier-to-interference-to-noise ratio calculation and optimization

Abstract

The present invention provides a high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization, comprising randomly generating a plurality of feasible solutions for frequency allocation under constraints such as different frequencies of adjacent carriers and bandwidth requirements of each carrier as an initial population, performing frequency optimization allocation based on a differential evolution algorithm, iteratively selecting individuals with the best system capacity matching rate as a system frequency allocation plan through operations such as mutation and crossover, and in the iterative adjustment process of the frequency allocation plan, considering the influence of the frequency plan adjustment on the system carrier-to-interference-noise ratio, performing dynamic quantitative analysis on the system carrier-to-interference-noise ratio based on statistics, and then accurately predicting the real-time system capacity matching rate, finding a nearly optimal frequency allocation plan, and improving the effective capacity of the satellite and resource utilization.

Description

Technical Field

[0001] The present invention belongs to the technical field of high-throughput satellite frequency control, and in particular to a high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization. Background Art

[0002] The frequency allocation of traditional high-throughput communication satellite transponders is mostly fixed throughout their life cycle, such as using a static, uniform multi-color multiplexing allocation scheme. However, once large-scale adjustments in user demand between different beams occur, or user demand within the same beam is seriously unbalanced, it will lead to "uneven busy and idle" resources, affecting satellite resource utilization and user service quality. Therefore, it is necessary to study the on-board frequency optimization allocation method to effectively alleviate the "uneven busy and idle" phenomenon between beams and improve the system's matching rate with user capacity requirements.

[0003] With the improvement of flexible payload capacity on satellites, many researchers have begun to study flexible and dynamic frequency plans. When flexibly allocating frequency resources of high-throughput satellites, the existing methods have two important defects:

[0004] (1) The carrier-to-interference-to-noise ratio is fixedly selected as a typical value, without considering the impact of frequency plan adjustment on the system carrier-to-interference-to-noise ratio;

[0005] (2) The method only focuses on finding feasible solutions for frequency arrangement under constraints such as different frequencies of adjacent carriers and bandwidth requirements of each carrier, but does not accurately predict the capacity matching rate after the frequency plan is adjusted, making it impossible to optimize the capacity matching rate. Summary of the invention

[0006] The technical problem solved by the present invention is: to overcome the shortcomings of the prior art and provide a high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization, thereby improving the effective capacity and resource utilization of satellites through dynamic and quantitative analysis of the carrier-to-interference-noise ratio and accurate prediction of the capacity matching rate.

[0007] The technical solution of the present invention is:

[0008] A high-throughput satellite frequency allocation method based on carrier-to-interference-to-noise ratio calculation and optimization comprises the following steps:

[0009] 1) Obtain the number of users to be allocated n and the required bandwidth B of each user i ,i=1,2...,n;get the frequency range available for allocation [f min , f max ] and evenly divided into N sub-bands, obtaining the frequency range corresponding to each sub-band, and each sub-band is numbered in ascending order according to the center frequency point; the required bandwidth B i Constructed by at least one subband, N is a positive integer;

[0010] 2) Randomly sort all users, and assign subbands to each user in order from the front to the back according to the sorting. Repeat this process W times, and obtain the starting subband number set assigned to each user each time as the current population, where W is a positive integer;

[0011] 3) Obtaining the set of termination subband numbers allocated to each user each time according to the current population, and forming W kinds of initial frequency allocation schemes;

[0012] 4) adjusting the W initial frequency allocation schemes respectively so that the subbands allocated to any two users in adjacent beams or belonging to the same beam have no overlapping parts, thereby forming W optimized frequency allocation schemes;

[0013] 5) Calculating the system carrier-to-interference-to-noise ratio of each optimized frequency allocation scheme in step 4) respectively, and calculating the system capacity matching rate of each optimized frequency allocation scheme based on the above ratio;

[0014] 6) The current population obtained in step 2) is subjected to mutation and crossover operations in sequence using a differential evolution algorithm to obtain a candidate population, and the candidate population is subjected to processing in steps 3) to 5) to obtain a system capacity matching rate corresponding to each optimized frequency allocation scheme formed by the candidate population;

[0015] 7) Compare and screen the current population and the candidate population, and replace the new population composed of the screened individuals as the new current population. The screening basis is the size of the system capacity matching rate corresponding to each optimized frequency allocation scheme formed by the current population and the candidate population;

[0016] 8) Repeat steps 6) to 7) R times, and select the scheme with the largest system capacity matching rate from the W optimized frequency allocation schemes formed by the final current population as the frequency allocation scheme of the system, where R is a positive integer.

[0017] Preferably, the allocation of subbands to each user in step 2) is specifically as follows:

[0018] The subbands occupied by the users that are ranked before the user and belong to the same beam or adjacent beam are removed from the frequency range available for allocation [f min , f max ] to obtain the remaining available sub-bands, and determine whether there is a frequency interval in the remaining available sub-bands whose continuous frequency is greater than or equal to the user's required bandwidth. If so, select the frequency interval with the smallest continuous frequency, and allocate the sub-bands to the user according to the required bandwidth starting from the smallest sub-band number. Otherwise, allocate the sub-band corresponding to the frequency interval with the largest remaining available continuous frequency to the user.

[0019] Preferably, the set of termination subband numbers allocated to each user each time is obtained in step 3), satisfying the following expression:

[0020]

[0021] In the formula, The starting subband number assigned to the user numbered i for the uth time, The number of the terminating subband allocated for the uth time to the user numbered i, u = 1, 2..., W, B uni is the sub-band granularity, N = (f max -f min ) / B uni .

[0022] Preferably, in step 4), the subbands allocated to any two users in adjacent beams or belonging to the same beam are free of overlap, specifically:

[0023] Two users with user numbers a and b are randomly selected to form a user group. If the two users are located in adjacent beams or belong to the same beam, and there is an overlap in the subband number ranges of the two users, the subband numbers of the two users are adjusted according to the following formula. Otherwise, no adjustment is required until all user groups are traversed and adjusted:

[0024]

[0025] in, Indicates the starting subband number of the u-th allocation to user a before adjustment, represents the adjusted starting subband number of the u-th allocation to user a, Indicates the number of the terminating subband allocated to user a for the uth time before adjustment. represents the adjusted terminating subband number of the u-th allocation to user a. represents the starting subband number of the u-th allocation to user b before adjustment, represents the adjusted starting subband number of the u-th allocation to user b, Indicates the ending subband number of the u-th allocation to user b before adjustment. represents the adjusted end subband number of user b u's allocation, B a represents the bandwidth requirement of user a, B b represents the bandwidth requirement of user b, Indicates the rounding operation. Indicates floor operation.

[0026] Preferably, the system capacity matching ratio corresponding to each optimized frequency allocation scheme is calculated in step 5), specifically:

[0027]

[0028] In the formula, θ xis the system capacity matching rate of the xth frequency allocation scheme, x=1,2...,W,C′ i is the capacity available to the ith user, C i is the required capacity of the i-th user.

[0029] Preferably, in the step 7), comparison and screening are performed according to the system capacity matching rate, specifically: the system capacity matching rate of each optimized frequency allocation scheme formed by the current population is compared with the optimized frequency allocation scheme formed by the corresponding candidate population, and the starting subband number set of each user corresponding to the optimized frequency allocation scheme with a higher system capacity matching rate is screened as an individual of the new population.

[0030] Preferably, in step 6), operations are performed according to the mutation operator and the crossover operator of the differential evolution algorithm, and the boundary constraint is: the subband number is greater than or equal to 1 and less than or equal to N.

[0031] Preferably, the value range of W is 5n≤W≤20n.

[0032] Preferably, the value range of R is 100-500.

[0033] Preferably, the sub-band granularity B uni For the frequency range available for allocation [f min , f max ] is used as the scale for averaging, and its value range is 1 to 5 MHz.

[0034] The advantages of the present invention compared with the prior art are:

[0035] (1) The present invention provides a high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization. During the iterative adjustment of the frequency plan, the system carrier-to-interference-noise ratio is dynamically and statistically quantitatively analyzed, thereby accurately predicting the real-time capacity matching rate, making the evaluation of the frequency allocation effect more realistic.

[0036] (2) The high-throughput satellite frequency allocation method adopted in the present invention performs frequency optimization allocation based on the differential evolution algorithm, randomly generates multiple feasible solutions for frequency allocation under constraints such as different frequencies of adjacent carriers and bandwidth requirements of each carrier as the initial population, and iteratively selects individuals with the best capacity matching rate through operations such as mutation and crossover. This method can find a nearly optimal frequency allocation solution and improve the effective capacity and resource utilization of the satellite. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 It is a flow chart of the present invention;

[0038] Figure 2 A schematic diagram of beam topology relationship in an embodiment of the present invention;

[0039] Figure 3 This is a diagram comparing the capacity matching of each user in the embodiment of the present invention with the capacity matching of each user that can be obtained by four-color multiplexing under the same input conditions. DETAILED DESCRIPTION

[0040] A high-throughput satellite frequency allocation method based on carrier-to-interference-to-noise ratio calculation and optimization takes into account the system carrier-to-interference-to-noise ratio changes caused by frequency allocation, and then adjusts the beam frequency band with the accurately predicted system capacity matching rate as the optimization target, such as Figure 1 As shown, the specific steps are:

[0041] 1) Obtain input parameters for frequency allocation, including:

[0042] Obtain m beams and n users, where m and n are positive integers and m and n ≥ 4;

[0043] Get the required bandwidth B of each user i and the required capacity C i , i is the user number, i=1,2...,n;

[0044] Get the frequency range available for allocation [f min , f max ];

[0045] Get sub-band granularity B uni , B uni Determined by the satellite's flexible payload capacity, the value range is 1 to 5 MHz;

[0046] Get the number of subbands available for allocation N = (f max -f min ) / B uni , all sub-bands are sorted from small to large according to the center frequency, and the sub-band numbers are obtained, which are 1, 2..., N, where N is a positive integer;

[0047] The topological relationship between beams and the topological relationship between users and beams are obtained.

[0048] 2) Randomly sort all users and allocate subbands to each user in order from the first to the last order. Specifically, the subbands occupied by users that are ranked before the user and belong to the same beam or adjacent beams are removed from the available frequency range [f min , f max ] to obtain the remaining available sub-bands, and determine whether there is a frequency interval with a continuous frequency greater than or equal to the user's required bandwidth. If so, select the frequency interval with the smallest continuous frequency, and allocate the sub-band to the user according to the required bandwidth of the user to be allocated starting from the minimum frequency of the minimum frequency interval. Otherwise, allocate the sub-band corresponding to the largest continuous frequency interval in the remaining available frequency bands to the user;

[0049] After obtaining the allocation, the frequency band starting subband number set Q of each user is 1 ,…,Q i ,…,Q n}, Q i The starting subband number of the user numbered i.

[0050] 3) Repeat step 2) W times, where W is a positive integer and 5n≤W≤20n, to obtain the frequency band starting subband number set Q allocated to each user each time (1) ,Q (2) ,…,Q (W) , which constitutes the current population, Q (u) is the set of starting sub-band numbers of the frequency band allocated to each user for the uth time, The starting sub-band number of the frequency band for the user numbered i in the u-th allocation, u = 1, 2, ..., W.

[0051] 4) According to the frequency band starting subband number set of each user in the current population, obtain the corresponding frequency band ending subband number set of each user Specifically by Q (u) get The expression is:

[0052]

[0053] In the formula, The starting sub-band number of the frequency band allocated for the uth time to the user with user number i, The number of the frequency band termination subband allocated for the uth time to the user numbered i, u = 1, 2, ..., W.

[0054] 5) Adjust Q (1) ,Q (2) ,…,Q (W) and Make sure that there is no overlap between any two user frequency bands in adjacent beams or in the same beam. Specifically, for the set Q (1) ,Q (2) ,…,Q (W) as well as Group in pairs, Q (u) and Form a group, divided into W groups in total, and complete the adjustment in groups one by one. (u) and When two users with user numbers a and b are selected to form a user group, if the two users are located in adjacent beams or belong to the same beam, and the subband number range of user a is Subband number range with user b If there is an overlap, the frequency sub-band numbers of users a and b are adjusted according to the following formula. Otherwise, no adjustment is required until all user groups are traversed and adjusted:

[0055]

[0056] in, Indicates the starting subband number of the u-th allocation of user a before adjustment, Indicates the adjusted starting subband number of the u-th allocation to user a, Indicates the end subband number of the u-th allocation of user a before adjustment, Indicates the adjusted end subband number of user a for the uth allocation. Indicates the starting subband number of the u-th allocation to user b before adjustment, represents the adjusted starting subband number of the u-th allocation to user b, Indicates the ending subband number of the u-th allocation to user b before adjustment, represents the adjusted termination subband number of user b for the uth allocation, u = 1, 2, ..., W, B a represents the bandwidth requirement of user a, B b represents the bandwidth requirement of user b, Indicates the rounding operation. Indicates floor operation.

[0057] The adjusted sets are Q′ (1) ,Q′ (2) ,…,Q′ (W) and Group the adjusted sets into pairs, Q′ (u) and Form a group, which is divided into W groups in total, each group corresponds to a frequency allocation scheme, forming W frequency allocation schemes.

[0058] 6) Calculate the system capacity matching rate θ corresponding to each frequency allocation scheme in step 5) respectively x , specifically:

[0059] First, the system carrier-to-interference-to-noise ratio under this frequency allocation scheme is calculated, and the link budget is performed, specifically:

[0060] Calculate the system carrier-to-interference-to-noise ratio Φ under this frequency allocation scheme. The calculation method is:

[0061]

[0062] In the formula, χ u is the uplink carrier-to-noise ratio, χ d is the downlink carrier-to-noise ratio, δu is the uplink carrier-to-interference ratio, δ d is the downlink carrier-to-interference ratio. u or d The method of obtaining is: the equivalent isotropic radiated power of the transmitting antenna minus the free space propagation loss and other additional losses, plus the quality factor of the receiving antenna, minus 10lg(kB), where k is the Boltzmann constant and B is the receiver bandwidth; δ u or δ d The method for obtaining is the ratio of the transmitted signal power to the sum of the co-channel interference signal powers.

[0063] Then the link budget is performed to obtain the capacity C′ available to each user i The method is: select a modulation coding method from the DVB modulation coding table, and look up the table to obtain the spectrum efficiency S corresponding to the modulation coding method E , roll-off coefficient α, signal-to-noise ratio Eb / N0 threshold H u , Eb is the signal energy per bit of the satellite signal, N0 is the noise power spectrum density; calculate the carrier-to-noise ratio threshold T Φ :T Φ =10lg(S E )+H u ; Select T Φ Under the premise that S is less than Φ E The maximum modulation coding method, and the corresponding spectrum efficiency is recorded as The roll-off factor is denoted as The capacity C′ available to each user i for: Where B′ i The bandwidth allocated to user numbered i under this frequency allocation plan.

[0064] Then the system capacity matching rate is calculated as follows:

[0065]

[0066] In the formula, θ x is the system capacity matching ratio of the xth frequency allocation scheme, x=1,2...,W.

[0067] 7) The current population Q obtained in step 3) (1) ,Q (2) ,…,Q (W) The mutation operator and crossover operator of the differential evolution algorithm are operated in sequence, and the boundary constraints of the subband number being greater than or equal to 1 and less than or equal to N are considered to obtain the candidate population. Specifically:

[0068] The mutation operator is: randomly select three numbers from 1, 2..., W, denoted as ε, η, μ, and perform the following calculation: Q (ε) +F×(Q (η) -Q (μ) ), and a mutant individual is obtained, where F is the scaling factor, ranging from 0 to 2; repeat W times to obtain the mutant population.

[0069] The crossover operator is: (1) Taking the crossover of mutant individual 1 as an example, a row vector containing W elements is generated, and the element value is randomly selected from 0 to 1; the crossover probability ξ is set, and the value range is 0 to 1; each element of the row vector is compared in turn to see if it is greater than ξ. If so, the element is updated to Q (1) If not, the element is updated and assigned to the element at the corresponding position of the variant individual 1, and finally the crossover individual 1 is generated; Q is completed in sequence (2) Crossover with mutant individual 2, ..., Q (W) Crossover with the mutant individual W to generate the post-crossover population.

[0070] 8) For candidate population Perform steps 4) to 6) to calculate the system capacity matching rate corresponding to each frequency allocation scheme formed by the candidate population.

[0071] 9) Compare the system capacity matching rates of the frequency allocation schemes formed by the current population and the frequency allocation schemes formed by the corresponding candidate populations, and select the frequency band starting sub-band number sets of each user corresponding to the frequency allocation scheme with a higher system capacity matching rate as individuals of the new population, specifically:

[0072] x takes values ​​from 1 to W, and compares θ x and like Then filter out As the xth individual; otherwise, Q (x) As the xth individual, the selected individuals form a new population and replace the current population Q (1) ,Q (2) ,…,Q (W) , the system capacity matching rate corresponding to the current population is reassigned to θ in turn 1 ,θ 2 ,…,θ W .

[0073] 10) Repeat steps 7) to 9) R times, where the value range of R is 100 to 500, and select the frequency allocation scheme with the largest system capacity matching rate from the final current population as the frequency allocation scheme of the system, specifically:

[0074] The system capacity matching rate set {θ 1 ,θ 2 ,…,θ W} select the maximum value θ q ,q∈{1,2,…,W}, then the optimal individual in the population is Q (q) , the corresponding frequency band termination subband number set of each user is Q (q) and The combination forms the optimal frequency allocation plan for the system.

[0075] In a specific embodiment, the satellite downlink has 19 beams and a total of 37 users; the frequency range available for downlink allocation is 17.7 GHz to 20.2 GHz, and the sub-band granularity is 3.125 MHz; the topological relationship between users and beams and the bandwidth requirements and capacity requirements of each user are shown in Table 1. Figure 2 It is a schematic diagram of the topological relationship between beams.

[0076] Table 1 Topological relationship between users and beams and bandwidth and capacity requirements of each user

[0077]

[0078] According to the method proposed in the present invention, the frequency allocation plan of the final user is calculated as shown in Table 2. The capacity matching of each user obtained by the method of the present invention is compared with the capacity matching of each user obtained by four-color multiplexing under the same input conditions. Figure 3 The total system capacity matching rate obtained by the method proposed in the present invention is 92%, which is about 7% higher than that of four-color multiplexing.

[0079] Table 2 User frequency allocation plan

[0080]

[0081]

[0082] Although the present invention has been disclosed as above in the form of a preferred embodiment, it is not intended to limit the present invention. Any person skilled in the art may make possible changes and modifications to the technical solution of the present invention by using the methods and technical contents disclosed above without departing from the spirit and scope of the present invention. Therefore, any simple modifications, equivalent changes and modifications made to the above embodiments based on the technical essence of the present invention without departing from the content of the technical solution of the present invention shall fall within the protection scope of the technical solution of the present invention.

[0083] The contents not described in detail in the specification of the present invention belong to the common knowledge of those skilled in the art.

Claims

1. A high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization, characterized in that: The following steps are involved: 1) Obtain the number of users to be allocated n and the required bandwidth B of each user i ,i=1,2...,n;get the frequency range available for allocation [f min , f max ] and evenly divided into N sub-bands, obtaining the frequency range corresponding to each sub-band, and each sub-band is numbered in ascending order according to the center frequency point; the required bandwidth B i Constructed by at least one subband, N is a positive integer; 2) Randomly sort all users, and assign subbands to each user in order from the front to the back according to the sorting. Repeat this process W times, and obtain the starting subband number set assigned to each user each time as the current population, where W is a positive integer; 3) Obtaining the set of termination subband numbers allocated to each user each time according to the current population, and forming W kinds of initial frequency allocation schemes; 4) adjusting the W initial frequency allocation schemes respectively so that the subbands allocated to any two users in adjacent beams or belonging to the same beam have no overlapping parts, thereby forming W optimized frequency allocation schemes; 5) Calculating the system carrier-to-interference-to-noise ratio of each optimized frequency allocation scheme in step 4) respectively, and calculating the system capacity matching rate of each optimized frequency allocation scheme based on the above ratio; The system carrier-to-interference-to-noise ratio Φ is calculated as: In the formula, χ u is the uplink carrier-to-noise ratio, χ d is the downlink carrier-to-noise ratio, δ u is the uplink carrier-to-interference ratio, δ d is the downlink carrier-to-interference ratio; u or d The method of obtaining is: the equivalent isotropic radiated power of the transmitting antenna minus the free space propagation loss and other additional losses, plus the quality factor of the receiving antenna, minus 10lg(kB), where k is the Boltzmann constant and B is the receiver bandwidth; δ u or δ d The method for obtaining is the ratio of the transmitted signal power to the sum of the co-frequency interference signal power; Then the link budget is performed to obtain the capacity C available to each user i ', the method is: polling a modulation coding method from the DVB modulation coding table, and looking up the table to obtain the spectrum efficiency S corresponding to the modulation coding method E , roll-off coefficient α, signal-to-noise ratio Eb / N0 threshold H u , Eb is the signal energy per bit of satellite signal, N0 is the noise power spectral density; Calculate the carrier-to-interference-to-noise ratio threshold T Φ :T Φ =10lg(S E )+H u ; Select T Φ Under the premise that S is less than Φ E The maximum modulation coding method, and the corresponding spectrum efficiency is recorded as The roll-off factor is denoted as The capacity C available to each user i 'for: Among them B i ' is the bandwidth allocated to user numbered i under this frequency allocation plan; Then the system capacity matching rate is calculated as follows: In the formula, θ x is the system capacity matching rate of the xth frequency allocation scheme, x=1,2...,W; 6) The current population obtained in step 2) is subjected to mutation and crossover operations in sequence using a differential evolution algorithm to obtain a candidate population, and the candidate population is subjected to processing in steps 3) to 5) to obtain a system capacity matching rate corresponding to each optimized frequency allocation scheme formed by the candidate population; 7) Compare and screen the current population and the candidate population, and replace the new population composed of the screened individuals as the new current population. The screening basis is the size of the system capacity matching rate corresponding to each optimized frequency allocation scheme formed by the current population and the candidate population; 8) Repeat steps 6) to 7) R times, and select the scheme with the largest system capacity matching rate from the W optimized frequency allocation schemes formed by the final current population as the frequency allocation scheme of the system, where R is a positive integer.

2. The high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to claim 1, characterized in that: The specific steps of allocating subbands to each user in step 2) are as follows: The subbands occupied by the users that are ranked before the user and belong to the same beam or adjacent beam are removed from the frequency range available for allocation [f min , f max ] to obtain the remaining available sub-bands, and determine whether there is a frequency interval in the remaining available sub-bands whose continuous frequency is greater than or equal to the user's required bandwidth. If so, select the frequency interval with the smallest continuous frequency, and allocate the sub-bands to the user according to the required bandwidth starting from the smallest sub-band number. Otherwise, allocate the sub-band corresponding to the frequency interval with the largest remaining available continuous frequency to the user.

3. The high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to claim 2, characterized in that: In step 3), the set of termination subband numbers allocated to each user each time is obtained, which satisfies the following expression: In the formula, The starting subband number assigned to the user numbered i for the uth time, The number of the terminating subband allocated for the uth time to the user numbered i, u = 1, 2..., W, B uni is the sub-band granularity, N = (f max -f min ) / B uni .

4. The high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to claim 3, characterized in that: In step 4), the subbands allocated to any two users in adjacent beams or belonging to the same beam are made to have no overlapping parts, specifically: Two users with user numbers a and b are randomly selected to form a user group. If the two users are located in adjacent beams or belong to the same beam, and there is an overlap in the subband number ranges of the two users, the subband numbers of the two users are adjusted according to the following formula. Otherwise, no adjustment is required until all user groups are traversed and adjusted: in, Indicates the starting subband number of the u-th allocation to user a before adjustment, represents the adjusted starting subband number of the u-th allocation to user a, Indicates the number of the terminating subband allocated to user a for the uth time before adjustment. represents the adjusted terminating subband number of the u-th allocation to user a. represents the starting subband number of the u-th allocation to user b before adjustment, represents the adjusted starting subband number of the u-th allocation to user b, Indicates the ending subband number of the u-th allocation to user b before adjustment. represents the adjusted end subband number of user b u's allocation, B a represents the bandwidth requirement of user a, B b represents the bandwidth requirement of user b, Indicates the rounding operation. Indicates floor operation.

5. The high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to claim 4, characterized in that: The system capacity matching ratio corresponding to each optimized frequency allocation scheme is calculated in step 5), specifically: In the formula, θ x is the system capacity matching rate of the xth frequency allocation scheme, x=1,2...,W,C i ′ is the capacity available to the i-th user, C i is the required capacity of the i-th user.

6. The high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to claim 5, characterized in that: In the step 7), comparison and screening are performed according to the system capacity matching rate, specifically: the system capacity matching rate of each optimized frequency allocation scheme formed by the current population is compared with the system capacity matching rate of the optimized frequency allocation scheme formed by the corresponding candidate population, and the starting subband number set of each user corresponding to the optimized frequency allocation scheme with a higher system capacity matching rate is selected as the individual of the new population.

7. A high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to any one of claims 1 to 6, characterized in that: In the step 6), operations are performed according to the mutation operator and the crossover operator of the differential evolution algorithm, and the boundary constraint is: the subband number is greater than or equal to 1 and less than or equal to N.

8. A high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to any one of claims 1 to 6, characterized in that: The value range of W is 5n≤W≤20n.

9. A high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to any one of claims 1 to 6, characterized in that: The value range of R is 100-500.

10. A high-throughput satellite frequency allocation method based on carrier-to-interference-noise ratio calculation and optimization according to any one of claims 3 to 6, characterized in that: The sub-band granularity B uni For the frequency range available for allocation [f min , f max ] is used as the scale for averaging, and its value range is 1 to 5 MHz.

Images