CN111008446A - Speed optimization system based on grid ship position layout calculation - Google Patents
Speed optimization system based on grid ship position layout calculation Download PDFInfo
- Publication number
- CN111008446A CN111008446A CN201911242385.1A CN201911242385A CN111008446A CN 111008446 A CN111008446 A CN 111008446A CN 201911242385 A CN201911242385 A CN 201911242385A CN 111008446 A CN111008446 A CN 111008446A
- Authority
- CN
- China
- Prior art keywords
- grid
- area
- ship
- day
- fine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to the technical field of distribution calculation of ship positions of measuring ships in the technical field of aerospace measurement and control, and particularly relates to a speed optimization system based on grid ship position distribution calculation. The working process of the optimization system comprises the following steps: calculating ship distribution areas of the ballistic offshore measurement and control arc section from head to tail to the lower sea surface, respectively taking intersections, and roughly dividing grids; carrying out mesh subdivision on the coarse mesh meeting the conditions according to the precision requirement to obtain a fine mesh; recording fine grids meeting the constraint for N consecutive days to obtain fine grid points; calculating a circumscribed convex polygon formed by the deployment region according to the fine grid points meeting the requirements; the area in the polygon is an area where the measuring vessel needs to measure and control multiple lanes for N consecutive days. The system reduces the calculated amount to a great extent and improves the calculation efficiency on the premise of ensuring the comprehensiveness and accuracy of the ship position deployment calculation.
Description
Technical Field
The invention belongs to the technical field of distribution calculation of ship positions of measuring ships in the technical field of aerospace measurement and control, and particularly relates to a speed optimization system based on grid ship position distribution calculation.
Background
The spacecraft has a plurality of key measurement and control arc sections from launching to entering the mission orbit, wherein the key measurement and control arc sections such as the orbit entering section, the near-location orbital transfer sections and the like only can provide measurement and control support for the offshore measuring vessel. The traditional survey vessel position arrangement is mostly manually calculated and drawn by survey crew, and the method has the advantages of large workload, strong experience and difficulty in adapting to the requirements of large-scale and high-frequency offshore survey tasks. The existing ship position layout calculation method based on the fine grid solves the problem of measuring the ship position layout area through manual calculation by using a computer. However, the calculation magnitude of the complete cycle traversal based on the fine grid grows exponentially with the increase of the number of the measuring ships and the increase of the number of the task days, and the calculation efficiency is low. In order to meet the requirement of a space measurement and control task and improve the ship position layout calculation efficiency, an implementation method based on fine grid ship position layout calculation is improved, the calculation amount is reduced, and the calculation speed is improved.
Disclosure of Invention
Technical problem to be solved
The technical problem to be solved by the invention is as follows: in the arrangement calculation of the ship positions of the measuring ships based on the fine grids, the problems of large calculation amount and long calculation time are solved along with the increase of the number of the measuring ships and the increase of the number of task days.
(II) technical scheme
In order to solve the above technical problem, the present invention provides a speed optimization system based on grid berth layout calculation, wherein the optimization system comprises: the system comprises a grid rough division module, a first operation module, a grid subdivision module, a second operation module, a fine grid point acquisition module and a measurement and control required area acquisition module;
the grid rough-dividing module is used for calculating ship distribution areas of a plurality of ballistic offshore measurement and control arc sections from head to tail to the lower sea surface in one day, respectively taking intersections, and roughly dividing the intersections of the ship distribution areas into grids;
the first operation module is used for circularly traversing the visibility of each coarse division grid in the calculation area to the measurement and control arc section and the link communication condition, and recording the coarse grids meeting the requirements;
the grid subdivision module is used for carrying out grid subdivision on the coarse grid meeting the conditions according to the precision requirement to obtain a fine grid;
the second operation module is used for circularly traversing and calculating the visibility of each fine grid to the measurement and control arc section and the link communication condition in the coarse grid meeting the conditions, and recording the fine grids meeting the requirements;
repeating the work of the grid subdivision module and the second operation module until the grid precision of the fine grid meets the task requirement, and stopping calculation;
the fine grid point acquisition module is used for recording the fine grids meeting the constraint for N consecutive days to obtain fine grid points;
the measurement and control requirement region acquisition module is used for calculating a circumscribed convex polygon formed by the deployment region according to the fine grid points meeting the requirements; the area in the polygon is an area where the measuring vessel needs to measure and control multiple lanes for N consecutive days.
In the working process of the grid rough division module, a plurality of sea surface ship distribution areas corresponding to the starting points of a plurality of ballistic maritime flight sections in the first day are calculated, and the intersection of the areas is taken as an area A; and simultaneously calculating a plurality of sea surface ship distribution areas corresponding to the terminal points of the plurality of ballistic maritime flight sections on the first day, and taking the intersection of the areas as an area B.
In the working process of the grid rough division module, two conditions exist: in the first situation, if the area A and the area B have an intersection area C, the measurement and control requirements of the first day can be met by arranging one measuring ship in the area C; in the second case, if the area a and the area B do not intersect, it indicates that the task needs two or more measurement ships to perform measurement and control.
Under the first condition in the working process of the grid rough division module, the area C can be divided into warps and wefts according to the calculation precision requirement only by one measuring ship due to the fact that the area is small, a plurality of fine grids are directly formed, visibility of each fine grid in the area C to the measurement and control arc section and the link communication condition are calculated according to the circular traversal of the first operation module, and the fine grids capable of meeting the full arc section visibility and the link communication are recorded.
Under the second condition in the working process of the grid rough dividing module, carrying out preliminary warp and weft grid rough dividing on the area A and the area B, and circularly traversing and calculating the measuring ship S in the area A according to the first operation moduleAFor visibility and link communication condition of corresponding measurement and control arc sections, simultaneously calculating a measurement ship S in the area BBThe visibility and the link communication condition of the corresponding measurement and control arc section; finding a grid M in an area A that can satisfy full coverage for rocket trajectory measurement and control1AiAnd a grid M in the region B1BjForming coarse mesh pairs M1AiAnd M1Bj。
Wherein, in the working process of the grid subdivision module and the second operation module, all the found coarse grid pairs M1AiAnd M1BjCarrying out grid subdivision, and circularly traversing and calculating the measuring ship S in each subdivided gridAVisibility and link communication conditions to the corresponding measurement and control arc segments, and a survey vessel SBAnd recording the visibility of the corresponding measurement and control arc section and the link communication condition in the coarse grid pair of the area A and the area B, wherein the fine grid pair can meet the visibility of the full arc section and the link communication.
In the working process of the repeated grid subdivision module and the second operation module, the fine grid pairs are further divided until the grid precision meets the task requirement and simultaneously meets the visibility of corresponding measurement and control arc sections and the full coverage of link communication, the calculation is stopped, and all the found fine grid pairs M are recorded1SAiAnd M1SBj。
Wherein, theIn the working process of the fine grid point acquisition module, according to the working processes of the grid rough division module, the first operation module, the grid subdivision module and the second operation module, the measurement and control coverage performance of the launching trajectory on the second day, the third day … and the Nth day is calculated respectively; recording day two, day three …, day N fine mesh pairs M satisfying visibility and link communication constraints2SAiAnd M2SBj，M3SAiAnd M3SBj…，MNSAiAnd MNSBj。
In the working process of the fine grid point acquisition module, according to the distance D of the measurement ship moving every day as the radius, a certain fine grid M in the ship distribution area of the first day is used1SAiMaking a circle as the circle center, wherein the intersection of the area covered by the circle and the ship distribution area on the second day is the layout-available area of the measuring ship corresponding to a certain grid on the second day, traversing and calculating all grids in the ship distribution area on the first day to obtain the ship distribution area on the second day corresponding to each fine grid in the ship distribution area on the first day of the measuring ship, and obtaining the fine grid points of the ship distribution of the measuring ship corresponding to the second day;
and then calculating the fine grid points of the measuring ship on the third day, … and the Nth day according to the method.
In the working process of the measurement and control required area acquisition module, calculating an external convex polygon formed in a ship position layout area of the measurement ship for the fine grid points recorded by the fine grid point acquisition module according to the following rule;
enveloping the areas constructed by all the discrete points to form a final ship position layout area; defining a longitude axis of a coordinate system as an X axis and a latitude axis as a Y axis;
the calculation process is as follows:
1) for the fine grid points in the discrete state recorded by the fine grid point acquisition module, finding one point in the discrete points, and recording the minimum point of the x coordinate as A under the condition of ensuring the maximum y coordinate1Point;
2) with A1Point is the origin, X-axis positive ray A1x clockwise scanningFinding the point scanned when the rotation angle is minimum, and marking as B1Point;
3) with point B1 as the origin, A1B1Direction ray A1B1Scanning clockwise, finding the point scanned when the rotation angle is minimum, and recording as C1Point;
4) with C1Point is the origin, B1C1Directional ray B1C1Scanning clockwise, finding the scanned point with the minimum rotation angle, and recording as D1Point;
and so on until finding the starting point A1(ii) a Thereby forming a starting point A1Point, the most important point is A1The points are circumscribed with convex polygons.
(III) advantageous effects
Compared with the prior art, the invention adopts a grid-based ship position layout calculation optimization scheme, greatly reduces the calculation amount and improves the calculation efficiency on the premise of ensuring the comprehensiveness and the accuracy of ship position layout calculation.
Drawings
Fig. 1 is a schematic diagram of coverage area subdivision optimization.
FIG. 2 is a schematic diagram of a convex polygon calculation method.
Detailed Description
In order to make the objects, contents, and advantages of the present invention clearer, the following detailed description of the embodiments of the present invention will be made in conjunction with the accompanying drawings and examples.
In order to solve the above technical problem, the present invention provides a speed optimization system based on grid berth layout calculation, wherein the optimization system comprises: the system comprises a grid rough division module, a first operation module, a grid subdivision module, a second operation module, a fine grid point acquisition module and a measurement and control required area acquisition module;
the grid rough-dividing module is used for calculating ship distribution areas of a plurality of ballistic offshore measurement and control arc sections from head to tail to the lower sea surface in one day, respectively taking intersections, and roughly dividing the intersections of the ship distribution areas into grids;
the first operation module is used for circularly traversing the visibility of each coarse division grid in the calculation area to the measurement and control arc section and the link communication condition, and recording the coarse grids meeting the requirements;
the grid subdivision module is used for carrying out grid subdivision on the coarse grid meeting the conditions according to the precision requirement to obtain a fine grid;
the second operation module is used for circularly traversing and calculating the visibility of each fine grid to the measurement and control arc section and the link communication condition in the coarse grid meeting the conditions, and recording the fine grids meeting the requirements;
repeating the work of the grid subdivision module and the second operation module until the grid precision of the fine grid meets the task requirement, and stopping calculation;
the fine grid point acquisition module is used for recording the fine grids meeting the constraint for N consecutive days to obtain fine grid points;
the measurement and control requirement region acquisition module is used for calculating a circumscribed convex polygon formed by the deployment region according to the fine grid points meeting the requirements; the area in the polygon is an area where the measuring vessel needs to measure and control multiple lanes for N consecutive days.
In the working process of the grid rough division module, a plurality of sea surface ship distribution areas corresponding to the starting points of a plurality of ballistic maritime flight sections in the first day are calculated, and the intersection of the areas is taken as an area A; and simultaneously calculating a plurality of sea surface ship distribution areas corresponding to the terminal points of the plurality of ballistic maritime flight sections on the first day, and taking the intersection of the areas as an area B.
In the working process of the grid rough division module, two conditions exist: in the first situation, if the area A and the area B have an intersection area C, the measurement and control requirements of the first day can be met by arranging one measuring ship in the area C; in the second case, if the area a and the area B do not intersect, it indicates that the task needs two or more measurement ships to perform measurement and control.
Under the first condition in the working process of the grid rough division module, the area C can be divided into warps and wefts according to the calculation precision requirement only by one measuring ship due to the fact that the area is small, a plurality of fine grids are directly formed, visibility of each fine grid in the area C to the measurement and control arc section and the link communication condition are calculated according to the circular traversal of the first operation module, and the fine grids capable of meeting the full arc section visibility and the link communication are recorded.
Under the second condition in the working process of the grid rough dividing module, carrying out preliminary warp and weft grid rough dividing on the area A and the area B, and circularly traversing and calculating the measuring ship S in the area A according to the first operation moduleAFor visibility and link communication condition of corresponding measurement and control arc sections, simultaneously calculating a measurement ship S in the area BBThe visibility and the link communication condition of the corresponding measurement and control arc section; finding a grid M in an area A that can satisfy full coverage for rocket trajectory measurement and control1AiAnd a grid M in the region B1BjForming coarse mesh pairs M1AiAnd M1Bj。
Wherein, in the working process of the grid subdivision module and the second operation module, all the found coarse grid pairs M1AiAnd M1BjCarrying out grid subdivision, and circularly traversing and calculating the measuring ship S in each subdivided gridAVisibility and link communication conditions to the corresponding measurement and control arc segments, and a survey vessel SBAnd recording the visibility of the corresponding measurement and control arc section and the link communication condition in the coarse grid pair of the area A and the area B, wherein the fine grid pair can meet the visibility of the full arc section and the link communication.
In the working process of the repeated grid subdivision module and the second operation module, the fine grid pairs are further divided until the grid precision meets the task requirement and simultaneously meets the visibility of corresponding measurement and control arc sections and the full coverage of link communication, the calculation is stopped, and all the found fine grid pairs M are recorded1SAiAnd M1SBj。
In the working process of the fine grid point acquisition module, according to the working processes of the grid rough division module, the first operation module, the grid subdivision module and the second operation module, the measurement and control coverage performance of the launching trajectory on the second day, the third day … and the Nth day is calculated respectively; recording day two, day three …, day N fine mesh pairs M satisfying visibility and link communication constraints2SAiAnd M2SBj，M3SAiAnd M3SBj…，MNSAiAnd MNSBj。
In the working process of the fine grid point acquisition module, according to the distance D of the measurement ship moving every day as the radius, a certain fine grid M in the ship distribution area of the first day is used1SAiMaking a circle as the circle center, wherein the intersection of the area covered by the circle and the ship distribution area on the second day is the layout-available area of the measuring ship corresponding to a certain grid on the second day, traversing and calculating all grids in the ship distribution area on the first day to obtain the ship distribution area on the second day corresponding to each fine grid in the ship distribution area on the first day of the measuring ship, and obtaining the fine grid points of the ship distribution of the measuring ship corresponding to the second day;
and then calculating the fine grid points of the measuring ship on the third day, … and the Nth day according to the method.
In the working process of the measurement and control required area acquisition module, calculating an external convex polygon formed in a ship position layout area of the measurement ship for the fine grid points recorded by the fine grid point acquisition module according to the following rule;
enveloping the areas constructed by all the discrete points to form a final ship position layout area; defining a longitude axis of a coordinate system as an X axis and a latitude axis as a Y axis;
the calculation process is as follows:
1) for the fine grid points in the discrete state recorded by the fine grid point acquisition module, finding one point in the discrete points, and recording the minimum point of the x coordinate as A under the condition of ensuring the maximum y coordinate1Point;
2) with A1Point is the origin, X-axis positive ray A1x scanning clockwise, finding the scanned point with the smallest rotation angle, and recording as B1Point;
3) with point B1 as the origin, A1B1Direction ray A1B1Scanning clockwise, finding the point scanned when the rotation angle is minimum, and recording as C1Point;
4) with C1The point is taken as the origin point,B1C1directional ray B1C1Scanning clockwise, finding the scanned point with the minimum rotation angle, and recording as D1Point;
and so on until finding the starting point A1(ii) a Thereby forming a starting point A1Point, the most important point is A1The points are circumscribed with convex polygons.
In addition, the invention also provides a speed optimization method based on grid ship position layout calculation, and the optimization method can solve the problems of repeated calculation and large calculation amount in the ship position based on grid calculation and effectively improve the calculation efficiency;
the optimization method comprises the following steps:
step 1: calculating ship distribution areas of a plurality of ballistic offshore measurement and control arc sections from head to tail to the lower sea surface in one day, respectively taking intersections, and roughly dividing the intersections of the ship distribution areas into grids;
step 2: circularly traversing the visibility of each coarse grid in the calculation area to the measurement and control arc section and the link communication condition, and recording the coarse grids meeting the requirements;
and step 3: carrying out mesh subdivision on the coarse mesh meeting the conditions according to the precision requirement to obtain a fine mesh;
and 4, step 4: in the coarse grids meeting the conditions, circularly traversing and calculating the visibility of each fine grid to the measurement and control arc section and the link communication condition, and recording the fine grids meeting the requirements;
and 5: repeating the step 3 and the step 4 until the grid precision of the fine grid meets the task requirement, and stopping calculation;
step 6: recording fine grids meeting the constraint for N consecutive days to obtain fine grid points, and eliminating ship position points which do not meet the requirement;
and 7: calculating a circumscribed convex polygon formed by the deployment region according to the fine grid points meeting the requirements;
the area in the polygon is an area where the measuring vessel needs to measure and control multiple lanes for N consecutive days.
In the step 1, a plurality of sea surface ship distribution areas corresponding to the starting points of the plurality of ballistic offshore flight sections on the first day are calculated, and the intersection of the areas is taken as an area A; and simultaneously calculating a plurality of sea surface ship distribution areas corresponding to the terminal points of the plurality of ballistic maritime flight sections on the first day, and taking the intersection of the areas as an area B.
In step 1, there are two situations: in the first situation, if the area A and the area B have an intersection area C, the measurement and control requirements of the first day can be met by arranging one measuring ship in the area C; in the second case, if the area a and the area B do not intersect, it indicates that the task needs two or more measurement ships to perform measurement and control.
In the first case of step 1, the area C can be divided into warp and weft according to the calculation accuracy requirement because the area is small and only one measuring vessel is needed, so as to directly form a plurality of fine grids, and the fine grids capable of meeting the visibility of the full arc section and the link communication are recorded according to the visibility of each fine grid in the step 2 to the measurement and control arc section and the link communication condition.
In the second case of the step 1, the area A and the area B are subjected to preliminary rough division of warp and weft grids, and the survey ship S in the area A is calculated according to the loop traversal of the step 2AFor visibility and link communication condition of corresponding measurement and control arc sections, simultaneously calculating a measurement ship S in the area BBThe visibility and the link communication condition of the corresponding measurement and control arc section; finding a grid M in an area A that can satisfy full coverage for rocket trajectory measurement and control1AiAnd a grid M in the region B1BjForming coarse mesh pairs M1AiAnd M1Bj。
Wherein, in the step 3 and the step 4, all the coarse grid pairs M to be found1AiAnd M1BjCarrying out grid subdivision, and circularly traversing and calculating the measuring ship S in each subdivided gridAVisibility and link communication conditions to the corresponding measurement and control arc segments, and a survey vessel SBAnd recording the visibility of the corresponding measurement and control arc section and the link communication condition in the coarse grid pair of the area A and the area B, wherein the fine grid pair can meet the visibility of the full arc section and the link communication.
Wherein, in the step 5, the fine grid pair is processedFurther dividing, repeating the step 3 and the step 4 until the grid precision meets the task requirement, and stopping the calculation when the visibility of the corresponding measurement and control arc section and the full coverage of the link communication are met, and recording all the found fine grid pairs M1SAiAnd M1SBj。
In the step 6, according to the methods in the steps 1 to 5, the measurement and control coverage performance of the launching trajectory on the second day, the third day … and the Nth day is calculated respectively; recording day two, day three …, day N fine mesh pairs M satisfying visibility and link communication constraints2SAiAnd M2SBj，M3SAiAnd M3SBj…，MNSAiAnd MNSBj。
In step 6, a certain fine grid M in the first-day ship distribution area is used as a radius according to the distance D of the measurement ship moving every day1SAiMaking a circle as the circle center, wherein the intersection of the area covered by the circle and the ship distribution area on the second day is the layout-available area of the measuring ship corresponding to a certain grid on the second day, traversing and calculating all grids in the ship distribution area on the first day to obtain the ship distribution area on the second day corresponding to each fine grid in the ship distribution area on the first day of the measuring ship, and obtaining the fine grid points of the ship distribution of the measuring ship corresponding to the second day;
and then calculating the fine grid points of the measuring ship on the third day, … and the Nth day according to the method.
In the step 7, for the fine grid points recorded in the step 6, calculating an external convex polygon formed in the ship position layout area of the measurement ship according to the following rule;
enveloping the areas constructed by all the discrete points to form a final ship position layout area; defining a longitude axis of a coordinate system as an X axis and a latitude axis as a Y axis;
the calculation process is as follows:
1) for the fine grid points recorded in the step 6 and in the discrete state, finding one point in the discrete points, and recording the minimum point of the x coordinate as A under the condition that the y coordinate is maximum1Point;
2) with A1Point is the origin, X-axis positive ray A1x scanning clockwise, finding the scanned point with the smallest rotation angle, and recording as B1Point;
3) with point B1 as the origin, A1B1Direction ray A1B1Scanning clockwise, finding the point scanned when the rotation angle is minimum, and recording as C1Point;
4) with C1Point is the origin, B1C1Directional ray B1C1Scanning clockwise, finding the scanned point with the minimum rotation angle, and recording as D1Point;
and so on until finding the starting point A1(ii) a Thereby forming a starting point A1Point, the most important point is A1The points are circumscribed with convex polygons.
Example 1
The specific implementation of the invention is explained by taking the example that the carrier rocket is launched to the front and back of the separation of the satellite and the rocket, and the marine flight section of the carrier rocket needs more than 1 measuring ship to complete the continuous tracking measurement and control of a plurality of launching trajectories in a certain interval for a plurality of days, and combining the attached drawings. The method comprises the following specific steps:
step one, calculating a plurality of sea surface ship distribution areas corresponding to starting points of a plurality of ballistic maritime flight sections on the first day, and taking the intersection of the areas as an area A; and simultaneously calculating a plurality of sea surface ship distribution areas corresponding to the terminal points of the plurality of ballistic maritime flight sections on the first day, and taking the intersection of the areas as an area B.
There are two cases at this time: in the first situation, if the area A and the area B have an intersection area C, the measurement and control requirements of the first day can be met by arranging one measuring ship in the area C; in the second case, if the area a and the area B do not intersect, it indicates that the task needs two or more measurement ships to perform measurement and control.
And step two, under the first condition, because the area is smaller and only one measuring ship is needed, dividing the area C into warps and wefts according to the calculation precision requirement to form a plurality of fine grids, circularly traversing and calculating the visibility of each grid in the area C to the key arc section and the link communication condition, and recording the grids capable of meeting the full arc section visibility and the link communication.
Step three, under the second condition, roughly dividing the area A and the area B into a primary warp grid and a primary weft grid, and circularly traversing and calculating the measuring ship S in the area AAFor visibility and link communication condition of corresponding measurement and control trajectory, simultaneously calculating measurement ship S in region BBAnd correspondingly measuring and controlling the visibility of the trajectory and the link communication condition. Finding a grid M in an area A that can satisfy full coverage for rocket trajectory measurement and control1AiAnd a grid M in the region B1BjForming coarse mesh pairs M1AiAnd M1Bj；
The grid rough division method is as follows, taking a sea surface ship distribution circular area corresponding to the starting point of the offshore flight segment as an example:
and calculating the diameter D of the circular area, wherein the value of D is between 1000 km and 3000 km after a plurality of task analyses. In order to ensure the calculation efficiency of the first rough division, the calculation test is carried out for a plurality of times, and when D belongs to (1000,2000) kilometers, the step length D of the rough division grid is 100 kilometers; when D belongs to (2000,3000) kilometers, the step length D of the rough division grid is 200 kilometers;
step four, all the coarse grid pairs M found in the step three are used1AiAnd M1BjCarrying out grid subdivision, and circularly traversing and calculating the measuring ship S in each subdivided gridAVisibility and link communication conditions of the corresponding measurement and control trajectory, and a measurement vessel SBAnd recording the visibility of the corresponding measurement and control trajectory and the link communication condition in the coarse grid pair of the area A and the area B, wherein the fine grid pair can meet the full arc visibility and the link communication.
The mesh subdivision method is as follows.
As shown in FIG. 1, A, B is a starting point and an end point of an offshore flight segment, a circle C and a circle D are respectively an offshore coverage area of the starting point and the end point, ∠α is a maximum elevation angle value of an antenna, the height of a sub-satellite point of the starting point is h, and an O point is a midpoint of a connecting line of circle centers of the circle C and the circle D. for a rocket launching offshore measurement and control task, relevant constraints include the number (N) of measuring ships, the distance (h) of the sub-satellite points, the elevation angle range of the antenna, the length of a task arc segment and the like, and the rocket launching offshore:
the subdivision grid step length d' is:
d'＝d·Δ
step five, further dividing the fine grids, repeating the step four and the step five until the grid precision meets the task requirement, and stopping calculation when the visibility of the corresponding measurement and control trajectory and the full coverage of link communication are met, and recording all the found fine grid pairs M1SAiAnd M1SBj。
And step six, calculating the measurement and control coverage performance of the launching trajectory on the second day, the third day … and the Nth day according to the method of the step one to the step five. Recording day two, day three …, day N fine mesh pairs M satisfying visibility and link communication constraints2SAiAnd M2SBj，M3SAiAnd M3SBj…，MNSAiAnd MNSBj；
Step seven, taking a certain grid M in the first ship distribution area as a radius according to the distance D of the measurement ship moving every day1SAiAnd making a circle as the center of the circle, wherein the intersection of the area covered by the circle and the ship distribution area on the second day is the layout-available area of the measuring ship corresponding to a certain grid on the second day, and traversing and calculating all grids in the ship distribution area on the first day to obtain the ship distribution area on the second day corresponding to each grid in the ship distribution area on the first day of the measuring ship.
And step eight, calculating fine grid points of the measuring ship on the third day, … and the Nth day according to the method in the step seven.
And step nine, calculating the circumscribed convex polygon formed by the ship position layout area of the measuring ship according to the following rules for the points recorded in the step eight.
And enveloping the areas constructed by all the discrete points to form a final berth layout area. The calculation method is shown in fig. 2, in which the X axis is the longitude axis and the Y axis is the latitude axis.
The calculation process is as follows:
1) finding out the minimum point of the x coordinate in the discrete points and recording the minimum point as the point A under the condition of ensuring the maximum y coordinate;
2) taking the point A as an origin, scanning the X-axis positive ray Ax clockwise, finding a scanned point with the minimum rotation angle, and marking as a point B;
3) taking the point B as an original point, scanning clockwise by using rays AB in the direction AB, finding a scanned point when the rotation angle is minimum, and recording the point as a point C;
4) taking the point C as an origin, scanning the ray BC in the BC direction clockwise, finding a scanned point with the minimum rotation angle, and recording the scanned point as a point D;
and so on until a starting point a is found.
Step ten, if the grid pair M is not found in the step five until the grid precision meets the task requirement1AiAnd M1BjThen, it means that at least three measuring vessels are needed to perform measurement and control tasks. At the moment, the intersection D of the head-tail middle points of the ballistic offshore flight measurement and control sections in the measurement and control area on the sea surface is taken as the deployable area of the third measuring vessel. Dividing the area A, B, D into coarse grids and fine grids according to the method of the fourth step to the tenth step, and calculating the measuring ship S at the same timeASurvey vessel SBAnd a survey vessel SDAnd the visibility of the corresponding measurement and control trajectory and the link communication condition. And recording grid pairs which can fully cover the rocket trajectory for N consecutive days, and finding a proper ship position deployment area according to the daily moving distance constraint of the measuring ship.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (10)
1. A grid berth layout calculation-based speed optimization system, the optimization system comprising: the system comprises a grid rough division module, a first operation module, a grid subdivision module, a second operation module, a fine grid point acquisition module and a measurement and control required area acquisition module;
the grid rough-dividing module is used for calculating ship distribution areas of a plurality of ballistic offshore measurement and control arc sections from head to tail to the lower sea surface in one day, respectively taking intersections, and roughly dividing the intersections of the ship distribution areas into grids;
the first operation module is used for circularly traversing the visibility of each coarse division grid in the calculation area to the measurement and control arc section and the link communication condition, and recording the coarse grids meeting the requirements;
the grid subdivision module is used for carrying out grid subdivision on the coarse grid meeting the conditions according to the precision requirement to obtain a fine grid;
the second operation module is used for circularly traversing and calculating the visibility of each fine grid to the measurement and control arc section and the link communication condition in the coarse grid meeting the conditions, and recording the fine grids meeting the requirements;
repeating the work of the grid subdivision module and the second operation module until the grid precision of the fine grid meets the task requirement, and stopping calculation;
the fine grid point acquisition module is used for recording the fine grids meeting the constraint for N consecutive days to obtain fine grid points;
the measurement and control requirement region acquisition module is used for calculating a circumscribed convex polygon formed by the deployment region according to the fine grid points meeting the requirements; the area in the polygon is an area where the measuring vessel needs to measure and control multiple lanes for N consecutive days.
2. The grid berth layout calculation-based speed optimization method as claimed in claim 1, wherein in the working process of the grid rough separation module, a plurality of sea surface ship distribution areas corresponding to the starting points of a plurality of ballistic maritime flight segments on the first day are calculated, and the intersection of the areas is taken as an area A; and simultaneously calculating a plurality of sea surface ship distribution areas corresponding to the terminal points of the plurality of ballistic maritime flight sections on the first day, and taking the intersection of the areas as an area B.
3. The grid berth layout calculation-based speed optimization method as claimed in claim 2, wherein during the operation of the grid rough separation module, two conditions exist: in the first situation, if the area A and the area B have an intersection area C, the measurement and control requirements of the first day can be met by arranging one measuring ship in the area C; in the second case, if the area a and the area B do not intersect, it indicates that the task needs two or more measurement ships to perform measurement and control.
4. The grid ship position layout calculation-based speed optimization method according to claim 3, wherein in a first situation during the operation of the grid rough separation module, because the area is small and only one measuring ship is needed, the area C can be divided into warps and wefts according to the calculation accuracy requirement to directly form a plurality of fine grids, and the fine grids capable of meeting the full arc visibility and the link communication condition of each fine grid in the calculation area C are recorded according to the cyclic traversal of the first operation module on the visibility and the link communication condition of the measurement and control arc.
5. The grid berth layout calculation-based speed optimization method according to claim 4, wherein in the second situation in the working process of the grid rough-dividing module, the area A and the area B are subjected to preliminary warp and weft grid rough division, and a measuring ship S in the area A is circularly traversed and calculated according to the first operation moduleAFor visibility and link communication condition of corresponding measurement and control arc sections, simultaneously calculating a measurement ship S in the area BBThe visibility and the link communication condition of the corresponding measurement and control arc section; finding a grid M in an area A that can satisfy full coverage for rocket trajectory measurement and control1AiAnd a grid M in the region B1BjForming coarse mesh pairs M1AiAnd M1Bj。
6. The grid berth layout calculation-based speed optimization method of claim 5, wherein all the coarse grid pairs M to be found in the working process of the grid subdivision module and the second operation module1AiAnd M1BjCarrying out grid subdivision, and circularly traversing and calculating the measuring ship S in each subdivided gridATo corresponding measuring and controlling arc sectionVisibility and link communication conditions, and survey vessel SBAnd recording the visibility of the corresponding measurement and control arc section and the link communication condition in the coarse grid pair of the area A and the area B, wherein the fine grid pair can meet the visibility of the full arc section and the link communication.
7. The grid berth layout calculation-based speed optimization method of claim 6, wherein in the working process of repeating the grid subdivision module and the second operation module, the fine grid pairs are further divided until the grid precision meets the task requirement and simultaneously meets the visibility of corresponding measurement and control arc sections and the full coverage of link communication, the calculation is stopped, and all the found fine grid pairs M are recorded1SAiAnd M1SBj。
8. The grid berth layout calculation-based speed optimization method according to claim 5, wherein in the working process of the fine grid point acquisition module, the measurement and control coverage performance of the launching trajectory on the second day, the third day … and the Nth day is calculated respectively according to the working processes of the grid rough division module, the first operation module, the grid subdivision module and the second operation module; recording day two, day three …, day N fine mesh pairs M satisfying visibility and link communication constraints2SAiAnd M2SBj，M3SAiAnd M3SBj…，MNSAiAnd MNSBj。
9. The method as claimed in claim 5, wherein the fine grid point obtaining module is used for obtaining a fine grid M in the first-day distribution area according to the distance D of the ship moving each day as the radius1SAiMaking a circle as the center of the circle, wherein the intersection of the area covered by the circle and the ship distribution area on the second day is the layout-available area of the second-day measuring ship corresponding to a certain grid, and traversing and calculating all grids in the ship distribution area on the first day to obtain each fine grid in the ship distribution area on the first day of the measuring shipObtaining a corresponding fine grid point of the ship distribution of the survey ship in the corresponding ship distribution area of the second day;
and then calculating the fine grid points of the measuring ship on the third day, … and the Nth day according to the method.
10. The speed optimization method based on grid berth layout calculation of claim 9, characterized in that, in the working process of the measurement and control required region acquisition module, for the fine grid points recorded by the fine grid point acquisition module, the circumscribed convex polygon formed by the ship berth layout region of the measurement ship is calculated according to the following rules;
enveloping the areas constructed by all the discrete points to form a final ship position layout area; defining a longitude axis of a coordinate system as an X axis and a latitude axis as a Y axis;
the calculation process is as follows:
1) for the fine grid points in the discrete state recorded by the fine grid point acquisition module, finding one point in the discrete points, and recording the minimum point of the x coordinate as A under the condition of ensuring the maximum y coordinate1Point;
2) with A1Point is the origin, X-axis positive ray A1x scanning clockwise, finding the scanned point with the smallest rotation angle, and recording as B1Point;
3) with point B1 as the origin, A1B1Direction ray A1B1Scanning clockwise, finding the point scanned when the rotation angle is minimum, and recording as C1Point;
4) with C1Point is the origin, B1C1Directional ray B1C1Scanning clockwise, finding the scanned point with the minimum rotation angle, and recording as D1Point;
and so on until finding the starting point A1(ii) a Thereby forming a starting point A1Point, the most important point is A1The points are circumscribed with convex polygons.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911242385.1A CN111008446B (en) | 2019-12-06 | 2019-12-06 | Speed optimization system based on grid ship position layout calculation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911242385.1A CN111008446B (en) | 2019-12-06 | 2019-12-06 | Speed optimization system based on grid ship position layout calculation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111008446A true CN111008446A (en) | 2020-04-14 |
| CN111008446B CN111008446B (en) | 2023-03-10 |
Family
ID=70114969
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911242385.1A Active CN111008446B (en) | 2019-12-06 | 2019-12-06 | Speed optimization system based on grid ship position layout calculation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111008446B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111125948A (en) * | 2019-12-06 | 2020-05-08 | 北京京航计算通讯研究所 | Speed optimization method based on grid ship position layout calculation |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07253329A (en) * | 1994-03-16 | 1995-10-03 | Tsurumi Seiki:Kk | Apparatus for controlling position of underwater object with respect to ground |
| US20100118644A1 (en) * | 2008-11-07 | 2010-05-13 | Seale Daniel B | Method and system for controlling streamers |
| CN102520688A (en) * | 2011-12-12 | 2012-06-27 | 中国人民解放军92859部队 | Marine comprehensive survey ship measurement control method |
| CN105571405A (en) * | 2015-12-16 | 2016-05-11 | 北京航天控制仪器研究所 | Optical fiber underwater sound collision part and target-missing amount measuring system |
| CN105929709A (en) * | 2016-04-13 | 2016-09-07 | 中国人民解放军63680部队 | Multi-model rocket remote measurement data simulation system |
| CN105956286A (en) * | 2016-05-06 | 2016-09-21 | 北京航空航天大学 | Prediction method of forecabin thermal protection system whole trajectory temperature boundary of hypersonic velocity aircraft |
| US20170090053A1 (en) * | 2015-09-28 | 2017-03-30 | Conocophillips Company | 3d seismic acquisition |
| CN106649992A (en) * | 2016-11-17 | 2017-05-10 | 复旦大学 | Ship and wake mesh model merging and optimization method |
| CN107942743A (en) * | 2017-11-20 | 2018-04-20 | 中国运载火箭技术研究院 | A kind of measurement and control system of aircraft and method based on steel tower base station |
| CN109061583A (en) * | 2018-06-26 | 2018-12-21 | 中国卫星海上测控部 | A kind of boat-carrying Large Radar pedestal vertical section deformation measurement method |
| CN110414343A (en) * | 2019-06-24 | 2019-11-05 | 西安电子科技大学 | A kind of satellite in orbit remote sensing images Ship Detection |
| CN111125948A (en) * | 2019-12-06 | 2020-05-08 | 北京京航计算通讯研究所 | Speed optimization method based on grid ship position layout calculation |
-
2019
- 2019-12-06 CN CN201911242385.1A patent/CN111008446B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07253329A (en) * | 1994-03-16 | 1995-10-03 | Tsurumi Seiki:Kk | Apparatus for controlling position of underwater object with respect to ground |
| US20100118644A1 (en) * | 2008-11-07 | 2010-05-13 | Seale Daniel B | Method and system for controlling streamers |
| CN102520688A (en) * | 2011-12-12 | 2012-06-27 | 中国人民解放军92859部队 | Marine comprehensive survey ship measurement control method |
| US20170090053A1 (en) * | 2015-09-28 | 2017-03-30 | Conocophillips Company | 3d seismic acquisition |
| CN105571405A (en) * | 2015-12-16 | 2016-05-11 | 北京航天控制仪器研究所 | Optical fiber underwater sound collision part and target-missing amount measuring system |
| CN105929709A (en) * | 2016-04-13 | 2016-09-07 | 中国人民解放军63680部队 | Multi-model rocket remote measurement data simulation system |
| CN105956286A (en) * | 2016-05-06 | 2016-09-21 | 北京航空航天大学 | Prediction method of forecabin thermal protection system whole trajectory temperature boundary of hypersonic velocity aircraft |
| CN106649992A (en) * | 2016-11-17 | 2017-05-10 | 复旦大学 | Ship and wake mesh model merging and optimization method |
| CN107942743A (en) * | 2017-11-20 | 2018-04-20 | 中国运载火箭技术研究院 | A kind of measurement and control system of aircraft and method based on steel tower base station |
| CN109061583A (en) * | 2018-06-26 | 2018-12-21 | 中国卫星海上测控部 | A kind of boat-carrying Large Radar pedestal vertical section deformation measurement method |
| CN110414343A (en) * | 2019-06-24 | 2019-11-05 | 西安电子科技大学 | A kind of satellite in orbit remote sensing images Ship Detection |
| CN111125948A (en) * | 2019-12-06 | 2020-05-08 | 北京京航计算通讯研究所 | Speed optimization method based on grid ship position layout calculation |
Non-Patent Citations (3)
| Title |
|---|
| K.L. FANCHER等: "Spatial relationship surveys of navigation and surveying instrumentation of NOAA hydrographic ships and launches", 《PROCEEDINGS OF OCEANS 2005 MTS/IEEE》 * |
| 向颉等: "一种航天测量船船位数据替代方法", 《测控技术》 * |
| 周大尉等: "测量船校飞试验的视景仿真系统实现", 《系统仿真学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111125948A (en) * | 2019-12-06 | 2020-05-08 | 北京京航计算通讯研究所 | Speed optimization method based on grid ship position layout calculation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111008446B (en) | 2023-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101758934B (en) | Star sensor setting angle determining method based on mission planning | |
| CN103700286A (en) | Automatic carrier-landing guiding method of carrier-borne unmanned aircraft | |
| CN103674029B (en) | A kind of method of many ships collaborative navigation formation configuration based on underwater sound communication | |
| CN106681342B (en) | Airplane search and rescue method | |
| CN110866015B (en) | Moving target moving range recording method based on local grid | |
| CN103727937B (en) | Star sensor based naval ship attitude determination method | |
| CN103985952B (en) | Boat-carrying A-E-C three-axis satellite communication antenna polarization angle of deviation real-time correcting method | |
| CN111008446B (en) | Speed optimization system based on grid ship position layout calculation | |
| CN104348539A (en) | Satellite communication point wave beam covering region calculation method | |
| CN114326774A (en) | Method and system for generating collision avoidance strategy of spacecraft | |
| CN111125948A (en) | Speed optimization method based on grid ship position layout calculation | |
| CN108489497B (en) | Safe navigation aiding method for preventing reef touch by using map | |
| Hayes et al. | Sonic‐boom propagation through a stratified atmosphere | |
| CN113359737A (en) | Ship formation self-adaptive event trigger control method considering formation expansion | |
| CN114200929B (en) | Rapid comb-type path planning method for maximum detection coverage rate of multi-underwater robot | |
| Zhang et al. | Optimising feeder routing for container ships through an electronic chart display and information system | |
| CN111239777B (en) | Satellite cluster hierarchical positioning method based on position fingerprint | |
| CN112330091B (en) | Satellite-borne SAR imaging autonomous task planning method | |
| CN110988789A (en) | Space survey ship position layout method based on grids | |
| Fogarty et al. | Computer mechanization of six-degree of freedom flight equations | |
| CN105353767A (en) | Satellite formation keeping control method capable of improving payload observation efficiency | |
| CN113655786B (en) | Unmanned boat group control method based on African bee intelligent algorithm | |
| CN209112419U (en) | One kind being used for Naval Architecture and Ocean Engineering dynamic positioning system | |
| Zhao et al. | A Star Pattern Recognition Algorithm Based on the Radial Companion‐Circumferential Feature | |
| Nebylov et al. | Control aspects of aerospace plane docking with ekranoplane for water landing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |