Projection dimension reduction rapid calculation method and system for large-scale spherical column phased array scanning beam dynamic activation subarray plane

Abstract

The invention discloses a projection dimensionality reduction rapid calculation method and system for a large-scale spherical column phased array scanning beam dynamic activation subarray plane, and belongs to the field of phased array measurement and control performance modeling, and the method comprises the steps: S1, setting parameters of a scheduling problem model; s2, acquiring orbit data of the aircraft, obtaining a visible time window set and calculating a visible arc section; s3, establishing a phased array resource set, and searching for an activated array element; and S4, establishing a large-scale phased-array antenna multi-beam scheduling objective function, and obtaining a final calculation result. According to the invention, the multi-beam measurement and control performance can be improved.

Description

Technical Field

[0001] The present invention relates to the field of phased array measurement and control performance modeling, and more specifically, to a projection dimension reduction fast calculation method and system for dynamically activated sub-array surfaces of a large-scale spherical-cylindrical phased array scanning beam. Background Art

[0002] According to the development trend of aerospace at home and abroad, large-scale constellation construction is a key development area in the future. The Starlink low-orbit satellite constellation of SpaceX in the United States deployed a total of 11,943 satellites in two phases: 1,584 Ka / Ku-band satellites were launched in the first phase; 2,841 Ka / Ku-band satellites and 7,518 V-band satellites were launched in the second phase. my country has also successively announced some satellite development plans. It can be foreseen that the number of low-orbit satellites in my country will usher in rapid, long-term and sustainable growth. It is expected that by 2027, the number of low-orbit satellites in orbit in my country will increase several times. The large-scale increase in the number of satellites in orbit has led to a sharp increase in measurement and control needs, posing a huge challenge to the capabilities of existing satellite measurement and control systems. Starting with existing equipment, improving the measurement and control efficiency of equipment is a necessary means to alleviate measurement and control pressure. Due to its simultaneous multi-beam forming capability, large-scale phased arrays have obvious advantages over traditional parabolic antenna measurement and control systems in terms of beam number, beam flexibility, and equipment integration, and are the main equipment for large-scale constellation measurement and control in the future.

[0003] Taking the hemispherical + cylindrical phased array antenna as an example, it usually forms a beam by activating array elements and uses window sliding to complete continuous target tracking. When beam scheduling is performed according to the measurement and control task, the phased array antenna array is activated according to the satellite trajectory, and there is a conflict in the array resources activated at the same time for different tasks. Therefore, it is necessary to model the scheduling problem based on task requirements and array resource constraints to reasonably allocate beams to perform satellite measurement and control tasks and reduce array resource conflicts. In scenarios with large-scale phased arrays and a large number of measurement and control tasks, how to quickly and effectively calculate the activation of array elements is a difficulty in modeling due to the large number of array elements and task targets. Summary of the invention

[0004] The purpose of the present invention is to overcome the deficiencies of the prior art and provide a method and system for fast calculation of projection dimension reduction of dynamically activated sub-array planes of large-scale spherical-cylindrical phased array scanning beams, thereby improving the multi-beam measurement and control performance of the system.

[0005] The object of the present invention is achieved through the following solutions:

[0006] A fast calculation method for projection dimension reduction of dynamically activated sub-array planes of large-scale spherical cylindrical phased array scanning beams comprises the following steps:

[0007] S1. Set the parameters of the scheduling problem model;

[0008] S2, obtain the spacecraft orbit data, obtain the visible time window set and calculate the visible arc segment;

[0009] S3, establishing a phased array resource set and searching for activated array elements;

[0010] S4. Establish the multi-beam scheduling objective function of large-scale phased array antennas and obtain the final calculation results.

[0011] Furthermore, in step S1, the parameters of the scheduling problem model specifically include:

[0012] Scheduling cycle [S T ,E T ,d T ]:S T is the starting time, E T is the end time, d T is the total duration of the time period;

[0013] Satellite set S: the set of all satellites participating in the scheduling period. The visible window of each satellite corresponds to a subset S. i ;

[0014] Visible time window set W: the visible time window subset S between the ground station and different satellites during the scheduling period i , the subset contains the visible time window identifier id wl , the visible window corresponds to the satellite identifier, st wl 、et wl Respectively represent the start time and end time of the visible time window, dt wl Indicates the duration of the visible time window;

[0015] Visible arc segment set H: corresponding to the visible time window, the visible arc segment subset hl between the ground station and different satellites, the subset contains the visible time window identifier id wl , the azimuth, elevation and altitude of the visible arc segment, which is used in conjunction with the visible time window set;

[0016] Tracking and control task set F: Tracking and control task set F contains the execution status of all tasks of all satellites in the scheduling time period. 1 means that the task is executed in the visible time window, and 0 means that the task is not executed in the visible time window. The format is as follows:

[0017]

[0018] in, n is the number of given tasks, j k It represents the number of visible windows of the kth satellite in the scheduling period; the matrix F corresponds to the visible window set W, and once the corresponding visible window is activated, a tracking and control beam can be determined one by one.

[0019] Furthermore, in step S1, the setting of the parameters of the scheduling problem model further includes the following sub-steps:

[0020] The angle between the beam pointing direction and the normal direction of the activated array element is defined as the activation angle. When the angle between the normal direction of the array element and the beam pointing direction is smaller than the activation angle, the array element is activated.

[0021] The set of activated array elements corresponding to each measurement and control task or each measurement and control window is represented as Z i , all array element sets are represented as follows:

[0022] Z={z1,z2,…,z i ,…,z m}.

[0023] Furthermore, in step S2, the acquisition of the aircraft orbit data, obtaining the visible time window set and calculating the visible arc segment specifically includes the following sub-steps:

[0024] Obtain aircraft trajectory data to obtain the aircraft position, determine the visible time window set and calculate the corresponding visible arc segment through the aircraft position at different times; wherein the specific calculation method of the visible arc segment is as follows:

[0025] Step 1: Convert the ground-fixed coordinate system of the aircraft into the station-centered coordinate system: Establish a coordinate system with a known measuring station (X r ,Y r ,Z r ) is the origin of the station-centered rectangular coordinate system, then the coordinate X of the satellite in this coordinate system L (E,N,U) is:X L =R(X s -X r );where X s is the coordinate vector of the satellite in the Earth-fixed coordinate system; X r is the coordinate vector of the measuring station in the ground-fixed coordinate system; X L is the coordinate vector of the satellite in the station center coordinate system; R is the rotation matrix, and:

[0026] R=R x (90-B)·R z (L+90);

[0027]

[0028] In the above formula, B is the geodetic latitude of the measuring station, and L is the geodetic accuracy of the measuring station;

[0029] Step 2: Calculate the pitch angle, azimuth angle and distance of the aircraft relative to the ground station: The formula for converting the satellite from the station-centered rectangular coordinate system to the station-centered polar coordinate system is:

[0030]

[0031] The visible arc segments are obtained through the visible positions of multiple time points in the same time window, and the design code obtains the visible arc segments by inputting the visible window.

[0032] Furthermore, in step S3, the establishing of the phased array resource set specifically includes the following sub-steps:

[0033] All array elements are numbered according to coordinates, array element coordinate information on the phased array is stored, and numbered; the established array element resource set includes array element numbers, coordinate information, and activation status.

[0034] Furthermore, in step S3, the searching for activated array elements specifically includes the sub-steps of searching for activated array elements by plane projection method, and determining by mapping and using binary search method.

[0035] Furthermore, the search for the activated array element by plane projection method is determined by mapping and using a binary search method, which specifically includes the following sub-steps:

[0036] For hemispherical and cylindrical phased arrays, the following situations exist during the search:

[0037] Case 1: All the activated array surfaces are on the hemisphere: the coordinates of the satellite and the center of the sphere are known, and the coordinates of the activation center point are obtained according to the equation of the line connecting the two points and the spherical equation; the mapping surface of the activated array element on the spherical surface mapped to the bottom surface of the cylinder is an ellipse, and the X and Y axis coordinates x1, y1 of the activation center mapping point on this ellipse are the same as the coordinates of the corresponding axes of the activation center point on the activation surface. The X and Y axis coordinates of the activation center mapping point can be used to obtain the angle A between the mapping line connecting the satellite and the center of the sphere and the X axis; (x1, y1), (x3, y1), (x1, y3), (x4, y4) are located on the mapped ellipse, and the three points on the ellipse and the center point determine an ellipse. The center point is the activation center mapping point, and then the three points x3, y3, x4, y4 are obtained by searching x3, y3, x4, y4 by binary search.

[0038] Furthermore, the search for the activated array element by plane projection method is determined by mapping and using a binary search method, which specifically includes the following sub-steps:

[0039] For hemispherical and cylindrical phased arrays, the following situations exist during the search:

[0040] Case 2: The active array is partially on the hemisphere and partially on the cylinder: The hemisphere is the same as case 1; the mapping surface on the hemisphere is the intersection of the cylinder base and the ellipse, that is, the X and Y axis coordinates do not exceed the bottom circle; (X3, Y3, H), (X2, Y2, H) are the edge points of the active surface and are located on the intersection of the hemisphere and the cylinder;

[0041] The satellite coordinates are known, and the angle h between the line from the satellite to the center of the sphere and its plane projection on the cylinder can be calculated; the angle between the line from (X2, Y2, H) to the center of the sphere and the line from the satellite to the center of the sphere is equal to θ degrees, so the angle a between the plane projection of the line from the satellite to the center of the sphere on the cylinder and the line from (X2, Y2, H) to the center of the sphere is calculated by the following formula:

[0042] cosθ=cosa*cosh;

[0043] Then use the equation of the plane circle on the cylinder to find (X2, Y2, H), and then find (X3, Y3, H);

[0044] For the cylindrical part, obtain X2, Y2, X3, and Y3, and then use the cylinder height H and 0 to binary search Z1. Z1 is the point with the smallest Z-axis coordinate on the active surface. Find the equation of the line connecting (X3, Y3, H) and (X2, Y2, H), and search for points on the cylinder whose X and Y axis coordinates are greater than this line and whose Z axis coordinate is greater than Z1 and less than the cylinder height, and whose angle between the line connecting the array element to the center of the sphere and the line connecting the satellite to the center of the sphere is less than or equal to θ degrees.

[0045] Further, in step S4, the establishment of a large-scale phased array antenna multi-beam scheduling objective function specifically includes the following sub-steps:

[0046] During the scheduling process, the comprehensive index η4 = w1η1 + w2η2 - w3η3 is evaluated based on the measurement and control task satisfaction rate η1, array resource utilization rate η2, and array resource conflict rate η3. The measurement and control task satisfaction rate η1 is calculated as:

[0047]

[0048] in, represents the number of elements that are 1 in the F matrix; the array element resource utilization rate η2 represents the ratio of all used array elements in the Z matrix to the total number of array elements, expressed as:

[0049]

[0050] in, represents the total number of array elements used, |Z| represents the total number of array elements; the array element resource conflict rate η3 represents the ratio of the number of conflicting array elements to the number of array elements used, which is expressed as:

[0051]

[0052] in, Indicates the number of conflicting array elements;

[0053] The objective function is set according to the comprehensive indicator η4=w1η1+w2η2-w3η3, which is used to comprehensively evaluate the measurement and control task satisfaction rate η1, array element resource utilization rate η2 and array element resource conflict rate η3.

[0054] A fast calculation method for projection dimension reduction of dynamically activated sub-array surfaces of a large-scale spherical-cylindrical phased array scanning beam includes a computer device, a computer program is stored in the memory of the computer device, and when the computer program is loaded by the processor of the computer device, any of the methods described above is executed.

[0055] The beneficial effects of the present invention include:

[0056] Aiming at the rapid search for activated array elements of large-scale phased array antennas, a fast calculation method of projection dimensionality reduction of dynamically activated sub-array planes of large-scale spherical-cylindrical phased array scanning beams is proposed. The activated array elements are quickly searched through plane projection and combined with multi-objective optimization functions to complete the large-scale phased array beam scheduling modeling, which provides a basis for improving the multi-beam measurement and control performance of the system. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0058] Figure 1 is a flowchart of the implementation steps of an embodiment of the present invention;

[0059] Figure 2 The active surface is located in the hemisphere schematic diagram and mapping diagram;

[0060] Figure 3 This is a schematic diagram showing that the active surface is partially located on the hemisphere and partially on the cylinder;

[0061] Figure 4 This is a large-scale phased array antenna model. DETAILED DESCRIPTION

[0062] All features disclosed in all embodiments in this specification, or steps in all methods or processes implicitly disclosed, except for mutually exclusive features and / or steps, can be combined and / or expanded or replaced in any manner.

[0063] In view of the problems in the background, the present invention proposes a fast calculation method for projection dimensionality reduction of dynamically activated sub-array surfaces of large-scale spherical-cylindrical phased array scanning beams. Under the conditions of target dynamic motion, array element channel power and weighted resource constraints, the multi-beam scheduling problem modeling is completed under the constraints of task requirements and array resources, which provides a basis for achieving optimal scheduling under the goals of task satisfaction rate, resource utilization rate and resource conflict rate, and ultimately improves the multi-beam measurement and control performance of the system.

[0064] The present invention is aimed at large-scale phased arrays. When performing measurement and control of a large number of tasks, when performing beam scheduling under the task requirements and array resource constraints, it is necessary to calculate a large number of array elements, determine which array elements need to be activated and whether there are conflicts, and consider the dynamic motion of the target, the array element channel power and the weighted resource constraints. It is necessary to model the beam scheduling problem under multiple constraints and propose the following technical solutions, such as Figure 1 As shown, the following steps are included:

[0065] Step A, initializing the parameters of the scheduling problem model;

[0066] Step B, obtaining the spacecraft orbit data, and calculating the visible time window set and the visible arc segment;

[0067] Step C, establishing a phased array resource set, and searching for activated array elements by plane projection method;

[0068] Step D: Establish a large-scale phased array antenna multi-beam scheduling objective function.

[0069] Furthermore, step A defines the parameters of the scheduling problem model. The specific parameters are defined as follows:

[0070] Scheduling cycle [S T ,E T ,d T ]:S T is the starting time, E T is the end time, d T is the total duration of the time period;

[0071] Satellite set S: the set of all satellites that can participate in the scheduling period. The visible window of each satellite corresponds to a subset S. i ;

[0072] Visible time window set W: the visible time window subset A between the ground station and different satellites during the scheduling period i , the subset contains the visible time window identifier id wl , the visible window corresponds to the satellite identifier, st wl 、et wl Respectively represent the start time and end time of the visible time window, dt wl Indicates the duration of the visible time window;

[0073] Visible arc segment set H: corresponding to the visible time window, the visible arc segment subset hl between the ground station and different satellites, the subset contains the visible time window identifier id wl , the azimuth (from the ground station), pitch (from the ground station), and altitude of the visible arc. This set is used in conjunction with the visible time window set;

[0074] Tracking and control task set F: Tracking and control task set F contains the execution status of all tasks of all satellites in the scheduling time period. 1 means that the task is executed in the visible time window, and 0 means that the task is not executed in the visible time window. The format is as follows:

[0075]

[0076] in, n is the number of given tasks, j k Represents the number of visible windows of the kth satellite in the scheduling period. The matrix F corresponds to the visible window set W. Once the corresponding visible window is activated, a telemetry beam can be determined one by one. In order to form this beam, a set of array elements that can synthesize the beam is obtained. Generally, the angle between the beam pointing and the normal direction of the activated array element is defined as the activation angle, and the activation angle is set to 60 degrees. When the angle between the normal direction of the array element and the beam pointing is less than the activation angle, the array element is activated. Therefore, the set of activated array elements corresponding to each telemetry task or each telemetry window is expressed as Z i , all array element sets can be expressed as the following set:

[0077] Z={z1,z2,…,z i ,…,z m}.

[0078] Furthermore, step B obtains the aircraft trajectory data in order to obtain the aircraft position, and calculates the visible time window set and the corresponding visible arc segment through the aircraft position at different times. For different aircraft, the trajectory data obtained are different. For example, for satellites, TLE data is obtained for calculation to obtain the satellite position at different times; if it is a drone, its trajectory data can be obtained by setting parameters such as heading and speed. Through the aircraft trajectory data, combined with the ground station position information, the visible time window and the corresponding visible arc segment can be calculated. The specific calculation method of the visible arc segment information is as follows.

[0079] 1. Convert the ground-fixed coordinate system of the aircraft into the station-centered coordinate system: Establish a coordinate system with a known measuring station (X r ,Y r ,Z r ) is the origin of the station-centered rectangular coordinate system, then the coordinate X of the satellite in this coordinate system L (E,N,U) is:

[0080] X L =R(X s -X r );

[0081] In the above formula, X s is the coordinate vector of the satellite in the Earth-fixed coordinate system; X r is the coordinate vector of the measuring station in the ground-fixed coordinate system; X Lis the coordinate vector of the satellite in the station center coordinate system; R is the rotation matrix, that is:

[0082] R=R x (90-B)·R z (L+90);

[0083]

[0084] In the above formula, B is the geodetic latitude of the measuring station, and L is the geodetic accuracy of the measuring station;

[0085] 2. Calculate the pitch angle, azimuth angle and distance of the aircraft relative to the ground station: The formula for converting the satellite from the station-centered rectangular coordinate system to the station-centered polar coordinate system is:

[0086]

[0087] The visible arc is obtained through the visible positions of multiple time points in the same time window, and the design code obtains the visible arc (pitch angle, azimuth angle, distance) by inputting the visible window.

[0088] Furthermore, one method of establishing a phased array resource set in step C is: numbering all array elements according to coordinates, storing array element coordinate information on the phased array, and numbering the array element resource set established includes array element numbers, coordinate information, and activation status.

[0089] One method of searching for activated array elements in step C is to determine them by mapping using a binary search method, as follows: for a hemisphere + cylindrical configuration phased array, the following situations exist during the search: 1. The activated array surface is entirely on the hemisphere; 2. The activated array surface is partially on the hemisphere and partially on the cylinder.

[0090] 1. All activation arrays are on the hemisphere: Figure 2 The activation surfaces shown are all on the hemisphere. The coordinates of the satellite and the center of the sphere are known. The coordinates of the activation center point can be obtained according to the equation of the line connecting the two points and the spherical equation. The mapping surface of the activation array element on the spherical surface mapped to the bottom surface of the cylinder is an ellipse. The X and Y axis coordinates x1, y1 of the activation center mapping point on this ellipse are the same as the coordinates of the corresponding axes of the activation center point on the activation surface. The X and Y axis coordinates of the activation center mapping point can be used to obtain the angle A between the mapping line connecting the satellite and the center of the sphere and the X axis. (x1, y1)(x3, y1)(x1, y3)(x4, y4) are located on the mapped ellipse. The three points on the ellipse and the center point determine an ellipse. The center point is the activation center mapping point, and then the three points x3, y3, x4, y4 are obtained by searching x3, y3, x4, y4 by binary search.

[0091] Take y3 as an example: It is known that the value of y3 is between y1 and y2. y1 can be obtained based on the activation center point, and y2 can be obtained based on the equation of the outer circle of the cylinder and x1. Take y3 = (y2 + y1) / 2. If the angle between the line connecting the corresponding activated array element and the satellite and the center of the phased array is exactly equal to the activation angle θ, then (x1, y3) is the mapping point of the activated array element outside the activation surface on the phased array; if it is less than the activation angle θ, y1 ​​is covered with (y2 + y1) / 2; if it is greater than the activation angle θ, y2 is covered with (y2 + y1) / 2. Repeat this cycle until y3 is found. After determining the ellipse equation according to the above method, the X and Y axis coordinates of the activated array element can be determined, and then the corresponding Z axis coordinate can be obtained according to the spherical equation.

[0092] 2. The activation array is partially on the hemisphere and partially on the cylinder: Figure 3 As shown, when the active array surface is partially located on the hemisphere and partially on the cylinder, the hemisphere part is the same as the above method, but because part of the active surface is located on the cylinder, the mapping surface on the hemisphere is the intersection of the cylinder base and the ellipse, that is, the X and Y axis coordinates do not exceed the bottom circle. (X3, Y3, H), (X2, Y2, H) are the edge points of the active surface and are located on the intersection of the hemisphere and the cylinder. From the known satellite coordinates, the angle h between the satellite-sphere center line and its plane projection on the cylinder can be obtained. The angle between the (X2, Y2, H)-sphere center line and the satellite-sphere center line is equal to θ degrees, then the angle a between the satellite-sphere center line plane projection on the cylinder and the (X2, Y2, H)-sphere center line can be calculated by the following formula:

[0093] cosθ=cosa*cosh;

[0094] Combined with the equation of the plane circle on the cylinder, we can find (X2, Y2, H). (X3, Y3, H) can be found in the same way;

[0095] For the cylindrical part, X2, Y2, X3, and Y3 are obtained according to the above process. Then, the cylinder height H is used to divide 0 by binary search (binary search in z = [0, H] to determine the value of Z1) to find Z1, which is the point with the smallest Z-axis coordinate on the active surface; the equation of the line connecting (X3, Y3, H) and (X2, Y2, H) is obtained, and the point whose angle between the line connecting the array element and the center of the sphere and the line connecting the satellite and the center of the sphere is less than or equal to θ degrees is searched on the cylinder whose X and Y axis coordinates are greater than this line and whose Z axis coordinates are greater than Z1 and less than the cylinder height.

[0096] Furthermore, a method for establishing the scheduling objective function in step D is as follows: During the scheduling process, the comprehensive index η4 = w1η1 + w2η2 - w3η3 can be evaluated based on the measurement and control task satisfaction rate η1, the array resource utilization rate η2, and the array resource conflict rate η3. The measurement and control task satisfaction rate η1 is calculated as:

[0097]

[0098] in Indicates the number of elements in the F matrix that are 1. The array element resource utilization rate η2 represents the ratio of all used array elements in the Z matrix to the total number of array elements, expressed as:

[0099]

[0100] in represents the total number of array elements used, and |Z| represents the total number of array elements. The array element resource conflict rate η3 represents the ratio of the number of conflicting array elements to the number of array elements used, which is expressed as:

[0101]

[0102] in Indicates the number of conflicting elements.

[0103] The above analysis of the target function data can be set according to the comprehensive index η4 = w1η1 + w2η2 - w3η3, which can comprehensively evaluate the measurement and control task satisfaction rate η1, array element resource utilization rate η2, and array element resource conflict rate η3.

[0104] In other preferred embodiments of the present invention, the specific embodiment includes the following steps:

[0105] Step 1: Set the parameters of the scheduling problem model:

[0106] The number of measurement and control tasks is 400;

[0107] The satellite collection uses satellites in the Starlink Phase I constellation;

[0108] Problem scheduling time period start time S_T: 2024,3,3,23:59:59;

[0109] Problem scheduling time period end time E_T: 2024,3,4,23:59:59;

[0110] Total duration of the period d_T: 1 day;

[0111] Set the activation angle θ to 60 degrees.

[0112] Step 2: Get the spacecraft orbit parameters and calculate the visible time window set and visible arc segment:

[0113] Find the Starlink satellite TLE file, set the ground station position (latitude 31°N, longitude 103°E), and take one of the satellites (numbered sa101 in the simulation) as an example. Its visible time window is shown in Table 1:

[0114] Table 1

[0115]

[0116] Calculate the visible arc (A is the pitch angle, h is the azimuth angle, and r is the distance) as shown in Table 2:

[0117] Table 2

[0118]

[0119] Step 3: Create a phased array resource set and search for activated array elements:

[0120] The large-scale phased array antenna array elements are set to 26469 in total, and all the elements are evenly distributed as follows: Figure 4 , and number them accordingly.

[0121] Combined with the phased array antenna elements, search for its activated elements. Taking the first row in the visible arc of the sa101 satellite as an example, the activated elements are obtained, and some are shown below.

[0122] Table 3

[0123]

[0124] Step 4: Establish the multi-beam scheduling objective function of large-scale phased array antennas:

[0125] Here, the weights w1=0.2, w2=0.4, w3=0.2 are set, and the task performance is evaluated according to the comprehensive indicator η4=w1η1+w2η2-w3η3.

[0126] Based on the first three steps, the final evaluation results are as follows.

[0127] Table 4

[0128] Number of tasks <![CDATA[η1]]> <![CDATA[η2]]> <![CDATA[η3]]> <![CDATA[η4]]> Detection time / s 20 1 0.9112 0.0188 0.5607 69.8281 100 0.73 1 0.0129 0.5434 216.5313 400 0.9875 1 0.0649 0.5845 1427.8

[0129] The units involved in the embodiments of the present invention may be implemented by software or hardware, and the units described may also be arranged in a processor. The names of these units do not, in some cases, limit the units themselves.

[0130] According to one aspect of an embodiment of the present invention, a computer program product or a computer program is provided, the computer program product or the computer program includes a computer instruction, and the computer instruction is stored in a computer-readable storage medium. A processor of a computer device reads the computer instruction from the computer-readable storage medium, and the processor executes the computer instruction, so that the computer device executes the method provided in the above various optional implementations.

[0131] As another aspect, an embodiment of the present invention further provides a computer-readable medium, which may be included in the electronic device described in the above embodiment; or may exist independently without being assembled into the electronic device. The above computer-readable medium carries one or more programs, and when the above one or more programs are executed by an electronic device, the electronic device implements the method described in the above embodiment.

Claims

1. A fast calculation method for projection dimension reduction of dynamically activated sub-arrays of large-scale spherical-cylindrical phased array scanning beams, characterized in that: The steps include: S1. Set the parameters of the scheduling problem model; S2, obtain the spacecraft orbit data, obtain the visible time window set and calculate the visible arc segment; S3, establishing a phased array resource set and searching for activated array elements; S4. Establish the multi-beam scheduling objective function of large-scale phased array antennas and obtain the final calculation results.

2. The fast calculation method for projection dimension reduction of dynamic activation sub-array planes of large-scale spherical cylindrical phased array scanning beam according to claim 1 is characterized in that: In step S1, the parameters of the scheduling problem model specifically include: Scheduling cycle [S T ,E T ,d T ]:S T is the starting time, E T is the end time, d T is the total duration of the time period; Satellite set S: the set of all satellites participating in the scheduling period. The visible window of each satellite corresponds to a subset S. i ; Visible time window set W: the visible time window subset S between the ground station and different satellites during the scheduling period i , the subset contains the visible time window identifier id wl , the visible window corresponds to the satellite identifier, st wl 、et wl Respectively represent the start time and end time of the visible time window, dt wl Indicates the duration of the visible time window; Visible arc segment set H: corresponding to the visible time window, the visible arc segment subset hl between the ground station and different satellites, the subset contains the visible time window identifier id wl , the azimuth, elevation and altitude of the visible arc segment, which is used in conjunction with the visible time window set; Tracking and control task set F: Tracking and control task set F contains the execution status of all tasks of all satellites in the scheduling time period. 1 means that the task is executed in the visible time window, and 0 means that the task is not executed in the visible time window. The format is as follows: Among them, f ijk ∈{0,1},i=1,2,…,n,j k =1,2,…,j i , n is the number of given tasks, j k It represents the number of visible windows of the kth satellite in the scheduling period; the matrix F corresponds to the visible window set W, and once the corresponding visible window is activated, a tracking and control beam can be determined one by one.

3. The fast calculation method for projection dimension reduction of large-scale spherical cylindrical phased array scanning beam dynamically activated sub-array surface according to claim 1 is characterized in that: In step S1, the parameters of the scheduling problem model are set, and the following sub-steps are also included: The angle between the beam pointing direction and the normal direction of the activated array element is defined as the activation angle. When the angle between the normal direction of the array element and the beam pointing direction is smaller than the activation angle, the array element is activated. The set of activated array elements corresponding to each measurement and control task or each measurement and control window is represented as Z i , all array element sets are represented as follows: Z={z1,z2,…,z i ,…,With m }。 4. The method for fast calculation of projection dimension reduction of dynamic activation sub-array planes of large-scale spherical cylindrical phased array scanning beam according to claim 1 is characterized in that: In step S2, the acquisition of the aircraft orbit data, obtaining the visible time window set and calculating the visible arc segment specifically includes the following sub-steps: Obtain aircraft trajectory data to obtain the aircraft position, determine the visible time window set and calculate the corresponding visible arc segment through the aircraft position at different times; wherein the specific calculation method of the visible arc segment is as follows: Step 1: Convert the ground-fixed coordinate system of the aircraft into the station-centered coordinate system: Establish a coordinate system with a known measuring station (X r ,Y r ,Z r ) is the origin of the station-centered rectangular coordinate system, then the coordinate X of the satellite in this coordinate system L (E,N,U) is:X L =R(X s -X r );where X s is the coordinate vector of the satellite in the Earth-fixed coordinate system; X r is the coordinate vector of the measuring station in the ground-fixed coordinate system; X L is the coordinate vector of the satellite in the station center coordinate system; R is the rotation matrix, and: R=R x (90-B)·R z (L+90); In the above formula, B is the geodetic latitude of the measuring station, and L is the geodetic accuracy of the measuring station; Step 2: Calculate the pitch angle, azimuth angle and distance of the aircraft relative to the ground station: The formula for converting the satellite from the station-centered rectangular coordinate system to the station-centered polar coordinate system is: The visible arc segments are obtained through the visible positions of multiple time points in the same time window, and the design code obtains the visible arc segments by inputting the visible window.

5. The fast calculation method for projection dimension reduction of large-scale spherical cylindrical phased array scanning beam dynamically activated sub-array surface according to claim 1 is characterized in that: In step S3, the phased array resource set is established, specifically including the sub-steps of: All array elements are numbered according to coordinates, array element coordinate information on the phased array is stored, and numbered; the established array element resource set includes array element numbers, coordinate information, and activation status.

6. The fast calculation method for projection dimension reduction of large-scale spherical cylindrical phased array scanning beam dynamically activated sub-array surface according to claim 1 is characterized in that: In step S3, the searching for activated array elements specifically includes the sub-steps of searching for activated array elements by plane projection method and determining by mapping and using binary search method.

7. The fast calculation method for projection dimension reduction of dynamic activation sub-array planes of large-scale spherical cylindrical phased array scanning beam according to claim 6 is characterized in that: The method of searching for the activated array element by plane projection method and determining it by mapping and using binary search method specifically includes the following sub-steps: For hemispherical and cylindrical phased arrays, the following situations exist during the search: Case 1: All the activated array surfaces are on the hemisphere: the coordinates of the satellite and the center of the sphere are known, and the coordinates of the activation center point are obtained according to the equation of the line connecting the two points and the spherical equation; the mapping surface of the activated array element on the spherical surface mapped to the bottom surface of the cylinder is an ellipse, and the X and Y axis coordinates x1, y1 of the activation center mapping point on this ellipse are the same as the coordinates of the corresponding axes of the activation center point on the activation surface. The X and Y axis coordinates of the activation center mapping point can be used to obtain the angle A between the mapping line connecting the satellite and the center of the sphere and the X axis; (x1, y1), (x3, y1), (x1, y3), (x4, y4) are located on the mapped ellipse, and the three points on the ellipse and the center point determine an ellipse. The center point is the activation center mapping point, and then the three points x3, y3, x4, y4 are obtained by searching x3, y3, x4, y4 by binary search.

8. The fast calculation method for projection dimension reduction of dynamic activation sub-array planes of large-scale spherical cylindrical phased array scanning beam according to claim 7 is characterized in that: The method of searching for the activated array element by plane projection method and determining it by mapping and using binary search method specifically includes the following sub-steps: For hemispherical and cylindrical phased arrays, the following situations exist during the search: Case 2: The active array is partially on the hemisphere and partially on the cylinder: The hemisphere is the same as case 1; the mapping surface on the hemisphere is the intersection of the cylinder base and the ellipse, that is, the X and Y axis coordinates do not exceed the bottom circle; (X3, Y3, H), (X2, Y2, H) are the edge points of the active surface and are located on the intersection of the hemisphere and the cylinder; The satellite coordinates are known, and the angle h between the line from the satellite to the center of the sphere and its plane projection on the cylinder can be calculated; the angle between the line from (X2, Y2, H) to the center of the sphere and the line from the satellite to the center of the sphere is equal to θ degrees, so the angle a between the plane projection of the line from the satellite to the center of the sphere on the cylinder and the line from (X2, Y2, H) to the center of the sphere is calculated by the following formula: cosθ=cosa*cosh; Then use the equation of the plane circle on the cylinder to find (X2, Y2, H), and then find (X3, Y3, H); For the cylindrical part, obtain X2, Y2, X3, and Y3, and then use the cylinder height H and 0 to binary search Z1. Z1 is the point with the smallest Z-axis coordinate on the active surface. Find the equation of the line connecting (X3, Y3, H) and (X2, Y2, H), and search for points on the cylinder whose X and Y axis coordinates are greater than this line and whose Z axis coordinate is greater than Z1 and less than the cylinder height, and whose angle between the line connecting the array element to the center of the sphere and the line connecting the satellite to the center of the sphere is less than or equal to θ degrees.

9. The fast calculation method for projection dimension reduction of dynamic activation sub-array planes of large-scale spherical cylindrical phased array scanning beam according to claim 1 is characterized in that: In step S4, the establishment of a large-scale phased array antenna multi-beam scheduling objective function specifically includes the following sub-steps: During the scheduling process, the comprehensive index η4 = w1η1 + w2η2 - w3η3 is evaluated based on the measurement and control task satisfaction rate η1, array resource utilization rate η2, and array resource conflict rate η3. The measurement and control task satisfaction rate η1 is calculated as: in, represents the number of elements that are 1 in the F matrix; the array element resource utilization rate η2 represents the ratio of all used array elements in the Z matrix to the total number of array elements, expressed as: in, represents the total number of array elements used, |Z| represents the total number of array elements; the array element resource conflict rate η3 represents the ratio of the number of conflicting array elements to the number of array elements used, which is expressed as: in, Indicates the number of conflicting array elements; The objective function is set according to the comprehensive indicator η4=w1η1+w2η2-w3η3, which is used to comprehensively evaluate the measurement and control task satisfaction rate η1, array element resource utilization rate η2 and array element resource conflict rate η3.

10. A fast calculation method for projection dimension reduction of dynamically activated sub-arrays of large-scale spherical cylindrical phased array scanning beams, comprising a computer device, characterized in that: A computer program is stored in the memory of the computer device, and when the computer program is loaded by the processor of the computer device, the method according to any one of claims 1 to 9 is executed.

Images