CN108333583B - Resource allocation method based on dual-target optimization of phased array radar search and tracking - Google Patents
Resource allocation method based on dual-target optimization of phased array radar search and tracking Download PDFInfo
- Publication number
- CN108333583B CN108333583B CN201810057487.5A CN201810057487A CN108333583B CN 108333583 B CN108333583 B CN 108333583B CN 201810057487 A CN201810057487 A CN 201810057487A CN 108333583 B CN108333583 B CN 108333583B
- Authority
- CN
- China
- Prior art keywords
- allocation
- search
- during
- sub
- array radar
- 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.)
- Active
Links
- 238000013468 resource allocation Methods 0.000 title claims abstract description 104
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005457 optimization Methods 0.000 title claims abstract description 37
- 239000013598 vector Substances 0.000 claims abstract description 64
- 230000009977 dual effect Effects 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 2
- 239000002245 particle Substances 0.000 claims description 2
- 239000000243 solution Substances 0.000 description 60
- 238000004088 simulation Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- -1 period Chemical class 0.000 description 2
- 239000002994 raw material Substances 0.000 description 2
- 239000012088 reference solution Substances 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241000764238 Isis Species 0.000 description 1
- 241000022852 Letis Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- XDDAORKBJWWYJS-UHFFFAOYSA-N glyphosate Chemical compound OC(=O)CNCP(O)(O)=O XDDAORKBJWWYJS-UHFFFAOYSA-N 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a resource allocation method based on phased array radar searching and tracking dual-target optimization, which belongs to the technical field of radar and mainly comprises the following steps: setting search area existence Q of time-controlled array radar distributed for the kth timekThe search area of the phased array radar is divided into N when the k-th distribution is carried outkA plurality of non-overlapping search sectors; let K be the kth distribution, K is more than or equal to 1 and less than or equal to K, the initial value of K is 1, and K is an even number greater than 0; obtaining the phased array radar allocation N during the 1 st allocation period1Optimal search time resources for non-overlapping search sectorsPhased array radar allocation to N during allocation up to KthKOptimal search time resources for non-overlapping search sectorsAnd phased array radar allocation to Q during the 1 st allocation1Tracking time resource column vector optimal solution of each targetPhased array radar allocation to Q during the Kth allocationKTracking time resource column vector optimal solution of each targetAnd recording as a resource allocation result based on phased array radar searching and tracking dual-target optimization.
Description
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a resource allocation method based on phased array radar searching and tracking dual-target optimization, which is suitable for resource allocation under the condition of limited time resource budget and improves the searching capability of a phased array radar and the tracking precision of a target to the maximum extent.
Background
Recently, advances in technology have enabled the development of agile, multi-tasking phased array radar systems; generally, a phased array radar uses an electronically steered array antenna and thus has extremely high beam flexibility, a characteristic that enables the phased array radar to perform multiple tasks; in fact, different radar functions may have competing demands on system resources, and therefore require the use of automated techniques to allocate resources depending on the capabilities of the radars and their goals; for radar search and tracking (SAT) applications, if one phased array radar detects targets with insufficient time resources, multiple low detectable targets may still not be found; at the same time, if a phased array radar does not have enough time resources to illuminate the previously tracked target, a discontinuous trajectory may result.
Heretofore, many approaches have been applied to the problem of resource allocation, whether a search function or a tracking function or both functions act together; for radar search applications, the challenge is to expand the radar surveillance area and improve the probability of target detection while using as few resources as possible; for target tracking, one may seek to minimize the total radar resources needed to track a target by optimizing the tracking-revisit interval, target signal strength, and detection threshold, and existing work, i.e., resource allocation schemes designed for SAT applications, can be broadly divided into two categories, rule-based and optimization-based; in rule-based schemes, a series of rules are formulated according to some operation requirements or radar features, and although the methods are very effective, the methods are not optimal in Bayesian theory and unpredictable behaviors can occur; another approach is to compute the SAT task by using a single cost function; the radar resource allocation problem may be formulated as a mathematical optimization scheme, and typically the cost function is a weighted sum of metrics corresponding to search capability and tracking accuracy, e.g., probability of detecting a target, tracking bayesian cralmelo lower bound (BCRLB), and expected measured signal-to-noise ratio (SNR); however, the disadvantage of these methods is the selection of weights and the meaningless aggregation of non-corresponding metrics.
Disclosure of Invention
In view of the problems in the prior art, the present invention is directed to a resource allocation method based on dual target search and tracking optimization for a phased array radar, which can solve the problem of resource allocation in a limited illumination time budget and simultaneously improve the search capability of the phased array radar and the tracking accuracy of a target to the maximum extent.
In view of the above problems with the prior art, it is an object of the present invention to design the SAT task resource allocation scheme as a dual target constrained optimization problem and to determine its well-known pareto subset (BK-PS) using the pareto theory. The resource allocation scheme employs two cost functions: (i) emphasizing that the target search capability is maximized in terms of multi-search sector search minimum signal-to-noise ratio (worst case search signal-to-noise ratio (WCS-SNR)); (ii) emphasizing multi-target tracking mean square error minimization in terms of a worst-case tracking Bayesian Cramer lower bound (WCT-BCRLB); in order to discuss the multiple balance between the two targets, a pareto optimal solution set of the two-target problem needs to be found; however, for many dual-objective problems, because the number of solutions in a solution set is large, determining the entire pareto optimal set is practically impossible; thus, a practical approach to dual target optimization is to calculate the BK-PS and represent the pareto optimal set with BK-PS as much as possible, with which one can find a suitable compromise between SAT tasks and select a resource allocation scheme accordingly to meet the specific application requirements.
In order to achieve the technical purpose, the invention is realized by adopting the following technical scheme.
A resource allocation method based on phased array radar searching and tracking dual-target optimization comprises the following steps:
Step 7, letAdds 1 to the value of (1), returns to step 2 until the phased array radar assignment during the 1 st assignment is obtainedOptimal search time resources for non-overlapping search sectorsTo the firstPhased array radar allocation during sub-allocationOptimal search time resources for non-overlapping search sectorsAnd phased array radar assignment during the 1 st assignmentTracking time resource column vector optimal solution of each targetTo the firstPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each targetAnd recording the resource allocation result as a resource allocation result based on dual target searching and tracking of the phased array radar.
The invention has the beneficial effects that:
firstly, the method utilizes the unique structure of a dual-target constraint optimization problem, and can obtain a pareto subset BK-PS by solving a salient pole small maximum optimization (CMO) problem in parallel; in previous studies, in order to obtain the pareto subset BK-PS, researchers developed many methods such as the weighted sum method, the ant colony optimization algorithm, and the genetic algorithm; the method of the invention develops a parallel minimization scheme to research the pareto subset BK-PS by utilizing the unique structure of the dual-objective optimization problem, so that the method of the invention is used as an alternative method for obtaining BK-PS, different pareto solutions can be obtained by solving a plurality of dual-objective problems with different SAT requirements in parallel, each complex dual-objective problem can be divided into two CMO problems, one for a search task and the other for a tracking task, and the solutions of search resource allocation (S-RA) problems corresponding to different SAT requirements are proportional.
Secondly, for different SAT parameters, the minimum search resource allocation problem only needs to be solved once, and as the target function of target tracking resource allocation is nonlinear, the method can obtain a pareto subset BK-PS with a base number of M by solving M +1 CMO problems in parallel; the result shows that the infinitesimal maximum radar search resource allocation problem generates a linear programming model, so that the problem can be solved easily by a famous linear programming method; for the target tracking resource allocation scheme, the generated minimum maximum problem consists of a set of separable monotonically decreasing convex functions, and a minimum maximum solution algorithm can be used for solving the T-RA problem.
Drawings
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
FIG. 1 is a flow chart of a resource allocation method based on dual target optimization for phased array radar search and tracking according to the present invention;
FIG. 2 is a schematic diagram of quantizing the search area of a phased array radar to non-overlapping sectors;
FIG. 3 is a schematic diagram of target deployment within a detection range of a phased array radar;
FIG. 4(a) is a schematic diagram comparing the performance of the resource-averaging allocation scheme with the pareto-based dual-target-optimized resource allocation scheme during the 5 th allocation;
fig. 4(b) is a schematic diagram comparing the performance of the resource average allocation scheme with the pareto-based dual target optimized resource allocation scheme during the 15 th allocation.
Detailed Description
Referring to fig. 1, it is a flowchart of a resource allocation method based on dual target optimization of phased array radar search and tracking according to the present invention; the resource allocation method based on the phased array radar searching and tracking dual-target optimization comprises the following steps:
Specifically, a phased array radar is determined, and the south-plus-south direction 93km and the west-plus-45 km of the phased array radar are taken as the origin pointsAnd establishing a plane rectangular coordinate system by taking the north direction as the Y axis and the east direction as the X axis.
As shown in FIG. 2, letIs as followsThe number of the sub-distribution is equal to the number of the sub-distribution,is set to an initial value of 1,is an even number greater than 0, in the present embodimentA value of 20; is set toSearch area existence of time-controlled array radar in sub-distributionObject, in the present embodimentThe value is 5; and a firstThe search area of the sub-distributed time-controlled array radar is divided intoNon-overlapping search sectorsFront sideThe search regions of the sub-distributed time phased array radar are respectively divided intoA number of non-overlapping search sectors that are,rear endThe search regions of the sub-distributed time phased array radar are respectively divided intoA plurality of non-overlapping search sectors; in this example
Specifically, in order to find an undetected target, the phased array radar needs to allocate its time resources to different search sectors; to further obtainTo the firstTime phased array radar in sub-distributionSearch sectorThe search model of (2) is:
wherein,a represents an intermediate variable which is,is shown asTime-controlled array radar allocation to the secondSearch sectorThe time resources of the search of (2),is shown asTime phased array radar in sub-distributionSearch sectorThe target search signal-to-noise ratio of (c),represents the average transmit power of the phased array radar,indicating the set effective receive aperture of the phased array radar antenna,is shown asTime phased array radar in sub-distributionSearch sectorThe cross-sectional area of scattering of the target,which represents the boltzmann constant, represents,indicating the set phased array radar temperature,the loss of the phased array radar is shown,is shown asTime-controlled array radar scanned by sub-distributionSearch sectorAngle of (2), in the present embodimentTaking the empirical value of the raw material to be tested,the empirical values of (A) are 6 °, 8 °, 12 °, 16 °;is shown asTime phased array radar in sub-distributionSearch sectorTarget distance search value of, in the present embodimentTaking the empirical value of the raw material to be tested,the empirical values of (A) are 200km, 240km, 295km, 300km and 310 km;can be directly measured according to prior information.
Is set toNumber of targets during sub-allocationIs known, in this embodimentWill be firstDuring the sub-distribution periodThe state vector of each object is recorded asWhich represents the transpose of the row vector,is shown asDuring the sub-distribution periodThe X-axis directional position of the individual target,is shown asDuring the sub-distribution periodThe Y-axis directional position of the individual target,is shown asDuring the sub-distribution periodThe velocity of the individual target in the direction of the X-axis,is shown asDuring the sub-distribution periodThe speed of the object along the Y-axis direction with the constraint of the firstDuring the sub-distribution periodState vector of individual targetDimension of (2)In this exampleIs shown asAnd distributing the number of targets in the search area of the phased array radar in time.
wherein,is shown asDuring the sub-distribution periodThe process noise of the individual target is,is shown asDuring the sub-distribution periodThe state vector of the individual objects is,is as followsDuring the sub-distribution periodThe transformation matrix of the individual objects is,the expression of the kronecker operator,representing a 2 x 2 dimensional identity matrix, assuming the process is noisySubject to a mean-zero Gaussian process that is noisyHas a covariance matrix ofIndicating the duration of each dispensing period, in this embodimentWhen in useThe 0 th allocation period when the value is 1The state vector of the individual target is noted asInitial state vector of individual targetFirst, theProcess noise of initial state vector of individual targetIs as followsThe process of each target is noisy with respect to the random number during initial assignment.
First, theDuring the sub-distribution periodThe measured value of each target isThe expression is as follows:
in the formula (3)
Wherein,is shown asDuring the sub-distribution periodState vector of individual targetIs/are as followsThe dimensional non-linear distance and orientation measurement functions,representing the coordinates of the phased array radar in a planar rectangular coordinate system,the position of the phased array radar in the X-axis direction in a plane rectangular coordinate system is represented,the position of the phased array radar in the Y-axis direction in a plane rectangular coordinate system is represented,is shown asDuring the sub-distribution periodThe position information of the individual objects is determined,is shown asDuring the sub-distribution periodThe radial distance of an individual target from the phased array radar,is shown asDuring the sub-distribution periodThe X-axis directional position of the individual target,is shown asDuring the sub-distribution periodY-axis position of individual target, superscriptIndicating transpose and arctan for arctan.
Will be firstDuring the sub-distribution periodThe error of each target is recorded asSet errorIs mean zero uncoupledResultant measurement error ofDuring the sub-distribution periodError of individual targetIs a diagonal covariance matrix of
Wherein,is shown asDuring the sub-distribution periodThe range of each target estimates the lower boundary of the cramer-circle of mean square error,is shown asDuring the sub-distribution periodThe lower boundary of the Cramer-Rao bound of the mean square error of the azimuth information estimation of each target is respectively as follows:
wherein,the speed of light is indicated and is,it is shown that the set constant is,is shown asDuring the sub-distribution periodThe expected measured echo signal to noise ratio (signal-to-noise ratio) for an individual target,is shown asDuring the sub-distribution periodThe reflectivity of the individual target is such that,is shown asDuring the second allocation period phased array radar is allocated to the secondThe tracking time resources of the individual targets,is shown asDuring the sub-distribution periodThe radial distance of an individual target to the phased array radar,is shown asThe-3 dB bandwidth of the electromagnetic wave signals transmitted by the phased array radar during the sub-distribution period is marked with the-1 to represent inversion,is shown as-3dB beamwidth of the phased array radar antenna during the secondary allocation; in this example
Due to the fact thatDuring the sub-distribution periodCramer-Rao bound lower bound of mean squared error of range estimates for individual targetsFirst, theDuring the sub-distribution periodCramer-Rao bound lower bound of mean square error of azimuth information estimation of individual targetsElectromagnetic wave signal-3 dB bandwidth-3dB beamwidthSum signal to noise ratioAll with tracking time resourcesIs inversely proportional, therefore will beDuring the sub-distribution periodError of individual targetThe diagonal covariance matrix of (A) extracts common factorsThe post rewrite is:
wherein,is shown asDuring the sub-distribution periodThe remaining matrix of the individual objects is,
Specifically, for the secondDuring the sub-distribution periodThe search sectors are not overlapped, and the search resource allocation of the phased array radar aims to optimally allocate search time resources to a plurality of areas and maximize the search signal-to-noise ratio under the worst condition; by usingIs shown asPhased array radar allocation during sub-allocationSearch time column vectors of non-overlapping search sectors, aThe objective function of searching the resource allocation scheme during the secondary allocation is:
wherein,is shown asDuring the sub-distribution periodA set of non-overlapping search sector numbers, a constraint condition is expressed in terms of the number of the elements,is shown asPhased array radar allocation during sub-allocationThe total search time resources of the non-overlapping search sectors,is shown asDuring the second allocation period phased array radar is allocated to the secondSearch time resources of the respective search sectors; the first constraint in equation (8) indicates thatPhased array radar allocation during sub-allocationThe total search time resources of the non-overlapping search sectors areThe second constraint states thatThe search time resources allocated to each search sector by phased array radar during the sub-allocation are limited by a minimum value, i.e. the secondThe search time resources allocated to each search sector by the phased array radar during the sub-allocation are greater than or equal to 0.
As is readily known, a maximized object type can be converted into a minimized object type by inversion; therefore, the search resource allocation problem of equation (8) can be newly formulated as the secondSearch for the transfer objective function of the resource allocation scheme during the secondary allocation:
wherein,is shown asThe convex function during the sub-allocation period,is shown as Search sector 1 to search sector 1 during sub-allocationThe search sectors search the accumulated sum of time resources, alpha represents an intermediate variable,representation calculationThe minimum value of (a) is determined,representing a set of computationsRatio of each search sector in the searchThen obtainA ratio is then comparedThe ratio value is then used to select the maximum value operation.
Specifically, for multi-target tracking, the time resources can be optimized according to the previous tracking informationThe tracking performance of a plurality of targets under the worst condition is improved; here, WCT-BCRLB is used as a standard function, and a target function of a target tracking resource allocation problem is made to be the secondTarget criteria function for tracking resource allocation scheme during sub-allocation:
wherein,is shown asDuring the sub-distribution periodA set of object numbers, each object number being,is shown asPhased array radar allocation during sub-allocationA tracking time resource column vector for each target,is shown asDuring the second allocation period phased array radar is allocated to the secondThe tracking time resources of the individual targets,is shown asPhased array radar allocation during sub-allocationTotal tracking time resources of the individual targets;expression solutionThe minimum value of (a) is determined,representing normalized worst caseDuring the sub-distribution periodTracking a Bayesian Claritrol lower bound convex function of each target; the first constraint indicatesPhased array radar allocation during sub-allocationThe total tracking time resource of each target isAnd the second constraint denotesThe tracking time resources allocated to each target by the phased array radar during the secondary allocation are limited by a minimum value, i.e. the secondThe tracking time resource allocated to each target by the phased array radar during the secondary allocation period is greater than or equal to 0; the normalized worst case firstDuring the sub-distribution periodTracking Bayesian Claritrol lower bound convex function of each targetThe expression is as follows:
wherein,representation matrixAnd Λ represents a standardized matrix, and shows that the elements of the bayesian cramer lower bound matrix are on different scales, and the expression is as follows:
indicating the duration of each dispensing period, in this embodimentRepresentation collectionEach of the other eyeThe marks correspond toA medium maximum value;is shown asDuring the second allocation period phased array radar is allocated to the secondTracking time resource of individual targetThe Bayesian Clarithrome lower bound matrix of (2) has the expression:
wherein,is shown asDuring the sub-distribution periodA predicted Bayesian information matrix of the observed state of each target,obtained by the following formula:
wherein,is to show toDuring the sub-distribution periodState vector of individual targetIs/are as followsDimensional non-linear distance and orientation measurement functionIs transferred toWith respect to state vectorsΔ represents the amount of change of the measured value,is shown asDuring the sub-distribution periodThe state vector of the individual objects is,is shown asDuring the sub-distribution periodOf a single objectA matrix of the Wijacobi is formed,it is shown that the set constant is,show thatValue ofIn (1),to representIs given a value ofCalculating the value of the point; when in useFirst 0 time during the dispensingState vector of individual targetIs shown asDuring the sub-distribution periodProcess noise of individual targetThe covariance matrix of (a) is determined,is as followsDuring the sub-distribution periodThe transformation matrix of the individual objects is,is shown asPhased array radar allocation during sub-allocationA tracking time resource column vector for each target,is shown asDuring the second allocation period phased array radar is allocated to the secondThe tracking time resources of the individual targets,is shown asDuring the sub-distribution periodThe remaining matrix of each object, superscript-1, represents the inversion,representing a row vector transpose.
In particular, the capabilities of phased array radars present a significant challenge to radar resource managers, which must determine during each assignment whether the radar should search for new targets or track existing targets; under the ideal condition, the maximization of the searching capability and the multi-target tracking precision are two mutually conflicting targets and must be considered simultaneously; thus, firstThe mathematical model of a dual target resource allocation scheme for phased array radar integrated SAT application during sub-allocation can be written as:
wherein,showing obtainedAndat the same time makeAndthe size of the particles is minimized and,is shown asPhased array radar allocation during sub-allocationThe total search time resources of the non-overlapping search sectors,is shown asPhased array radar allocation during sub-allocationThe total tracking time resources of the individual targets,is shown asIntegrating the total time resources of phased array radar search and tracking applications during the secondary allocation; the last constraint is indicated inThe total resource of the integrated phased array radar search and tracking application during the sub-allocation isIs shown asThe duty cycle during the sub-dispensing period,is shown asDuring the second allocation period phased array radar is allocated to the secondThe search time resources of each search sector,is shown asDuring the second allocation period phased array radar is allocated to the secondThe tracking time resources of the individual targets,is shown asThe convex function during the sub-allocation period,representing normalized worst caseDuring the sub-distribution periodThe tracked Bayesian Clarithrome lower bound convex function of each target,is shown asPhased array radar allocation during sub-allocationA search time column vector of non-overlapping search sectors,is shown asPhased array radar allocation during sub-allocationA tracking time resource column vector for each target.
Will be firstPhased array radar allocation during sub-allocationSearch time column vector of non-overlapping search sectorsAnd a firstPhased array radar allocation during sub-allocationTracking time resource column vector of individual targetsIntegrated into a single vector, denotedDimension vectorRepresenting a row vector transpose; further obtain the firstThe mathematical optimization model of the dual-target resource allocation scheme during the secondary allocation period is as follows:
wherein,is shown asThe convex function during the sub-allocation period,representing normalized worst caseDuring the sub-distribution periodThe tracked Bayesian Clarithrome lower bound convex function of each target,is expressed as length ofAnd beforeA row vector with 1 element and zero elements,is expressed as length ofAnd from the secondElement to elementA row vector with 1 element and 0 elements,representing row vectorsTo middleAn element; the first constraint indicatesPhased array radar allocation during sub-allocationThe total search time resources of the non-overlapping search sectors areFirst, thePhased array radar allocation during sub-allocationThe total tracking time resource of each target isThe second constraint indicatesThe searching time of each searching sector and the tracking time of each target in the secondary distribution period are both more than or equal to 0; the last constraint indicatesPhased array radar allocation during sub-allocationTotal search time resources for non-overlapping search sectorsAnd a firstPhased array radar allocation during sub-allocationTotal tracking time resource of each targetIs a sum ofThe total time resource of the integrated phased array radar search and tracking application during the sub-allocation is
The substep of step 6 is:
6.1 in order to obtain pareto subsets (the solutions obtained by different overall search budget solutions (9) are mutually independent, and the solutions obtained by different overall tracking budget solution problems (10) are mutually independent); is set toDuring the sub-distribution periodTotal search budget andan overall tracking budget, in this embodimentFrom 1 st total search budget toThe total search budget satisfies:whereinIs shown asDuring the sub-distribution periodTotal search budget, willDuring the sub-distribution periodTotal tracking budget asAnd the firstDuring the sub-distribution periodTotal search budgetAnd a firstDuring the sub-distribution periodTotal tracking budgetSatisfies the following conditions:
whereinIs shown asThe total time resources of the phased array radar search and tracking application are integrated during the secondary allocation.
6.2 according to the firstDuring the sub-distribution periodTotal search budgetAnd linear programming to solve equation (9)During the sub-distribution periodTotal search budgetSubstituting the right side of the first constraint condition in the formula (9) to obtain the second constraint condition according to a linear programming methodPhased array radar allocation during sub-allocationSearch time vector pair of non-overlapping search sectorsOptimal solution for individual search budget resource allocationAccording to the firstDuring the sub-distribution periodTotal tracking budgetAnd maximum and minimum solution algorithm solving equation (10), i.e. the firstDuring the sub-distribution periodTotal tracking budgetSubstituting the right side of the first constraint condition in the step (10), and obtaining the second constraint condition according to a maximum and minimum solution algorithmPhased array radar allocation during sub-allocationTracking time resource column vector pair of individual targetsOptimal solution for individual total tracking budget resource allocationSubscriptIs shown asSub-distribution period, subscriptDenotes a search, subscriptDenotes trace, subscriptRepresenting an optimal solution; will be the firstTotal searchOptimal solution for allocation of cable budget resourcesAnd the said firstOptimal solution for individual total tracking budget resource allocationForm two target resourcesThe second of the sub-distribution (equation (16))Pareto optimal solution (parallel minimization scheme) in which superscripts are appliedSecond to indicate transposed, dual target resource allocationPareto optimal solutionIncludedAnd (4) each element.
6.3 orderAre respectively 1 toRepeatedly executing 6.2 until obtaining the second target resourceFirst pareto optimal solution of sub-distributionTo the two target resourcesSecond of the sub distributionPareto optimal solutionIs expressed as a base numberPareto subsets of
Due to the fact that for the firstAny two different total search budgets during sub-allocationAndupper labelIs shown asIndividual total search budget, superscriptIs shown asThe total search budget for the search is,is shown asDuring the sub-distribution periodThe total search budget for the search is,is shown asDuring the sub-distribution periodTotal search budget, ofDuring the sub-distribution periodOptimal solution for individual search budget resource allocationAnd a firstDuring the sub-distribution periodOptimal solution for individual search budget resource allocationHas the following relationshipTherefore, it is caused byOne of the total search budgets during the sub-allocation and the optimal solution sum of the search budget resource allocationThe proportional relationship between the total search budgets can be obtainedThe total search budget is related to the solution of equation (9) respectively.
According to the cardinality ofPareto subsets ofThe calculation base isPareto subsets ofThe function value of each optimal solution in the system is defined as the base numberPareto subsets ofThe function value of each optimal solution is respectively marked as a pareto point, and then the second solution is obtainedPareto point set of time-controlled array radar for sub-distribution (just before calculation)One element), a represents an intermediate variable,is shown asTime phased array radar in sub-distributionSearch sectorIs calculated from the target distance of (a) to the target distance search value,is shown asTime-controlled array radar scanned by sub-distributionSearch sectorThe angle of (a) is determined,is shown asDuring the sub-distribution periodA set of non-overlapping search sector numbers,is shown asSub-distribution phase-controlled array radar 1 st pareto optimal solutionThe function value of (a) is determined,is shown asTime phased array radar with sub-distributionPareto optimal solutionThe function value of (a) is determined,is shown asTime phased array radar with sub-distributionPareto optimal solutionThe function value of (1).
By using the firstDuring the sub-distribution periodConvex function ofMonotonicity of, i.e. Indicating two target resourcesSub-assigned [ gamma ] pareto optimal solutionThe function value of (a) is determined,indicating two target resourcesPareto optimal solution of the beta of the sub-distributionThe function value of (a) is determined, is shown asThe optimal solution for the beta total search budget resource allocation during the secondary allocation,represents the optimal solution for the beta total tracking budget resource allocation,is shown asThe optimal solution for the gamma th overall search budget resource allocation during the sub-allocation,an optimal solution representing the gamma total tracking budget resource allocation; the upper level t represents the transpose,is shown asDuring the sub-distribution periodThe worst-case search signal-to-noise ratio in the non-overlapping search sectors,andare the function values of two adjacent pareto optimal solutions.
From two pareto optimal solutionsAndtwo total search budgets can be derivedAnd is shown asThe beta total search budget during the secondary allocation period is the second of the dual target resourcesPareto optimal solution of the beta of the sub-distributionMiddle frontA cumulative sum of the elements;is shown asThe gamma total search budget during the sub-allocation period is the second of the dual target resourcesSub-assigned [ gamma ] pareto optimal solutionMiddle frontA cumulative sum of the elements; budget for two total searchesAndusing a dichotomy to obtainPhased array radar allocation during sub-allocationMaximum of non-overlapping search sectorsOptimal total search time resourcesAnd a firstPhased array radar allocation during sub-allocationOptimal total tracking time resource of each target
The optimal total search time resourceSubstitute intoSearching the conversion objective function of the resource allocation scheme during the sub-allocation period, and solving by using a linear programming method to obtain the secondPhased array radar allocation during sub-allocationOptimal search time resources for non-overlapping search sectorsThe optimal total tracking time resourceSubstitute intoTracking the target standard function of the resource allocation scheme during the sub-allocation period, and obtaining the first time by utilizing a minimum maximum solution algorithmPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each target
Wherein,is the firstThe optimal search resource allocation result during the sub-allocation,is the firstOptimal target tracking resource allocation results during sub-allocation, the firstPhased array radar allocation during sub-allocationOptimal total search time resources for non-overlapping search sectorsAnd the said firstPhased array radar allocation during sub-allocationOptimal total tracking time resource of each targetAs a sum ofTotal time resources for integrated phased array radar search and tracking applications during sub-allocation
The method comprises the following specific steps:
6a) separately setting iteration indexesAnd a stop threshold ε of 10-3And is provided withSearching for lower bounds of resources during secondary allocationAnd a firstSearching for upper bounds of resources during secondary allocation Is shown asThe beta total search budget during the secondary allocation,is shown asA gamma total search budget during the secondary allocation; wherein, gamma is beta +1,is shown asThe total search budget number or the total tracking budget number set in the secondary allocation period.
6b) According to the firstSearching for lower bounds of resources during secondary allocationAnd a firstSearching for upper bounds of resources during secondary allocationCalculate the firstAfter the second iterationTotal search resources for phased array radar during secondary allocationThen use the firstAfter the second iterationTotal search resources for phased array radar during secondary allocationAnd linear programming method to solveSearching resource allocation during secondary allocationConverting the objective function of the scheme to obtainAfter the second iterationPhased array radar allocation during sub-allocationSearch time column vector optimal solution for non-overlapping search sectors
6e) If it is notGetSearch for the optimal solution for resource allocation asIs shown asPhased array radar allocation during sub-allocationThe optimal search time resources of the non-overlapping search sectors,is shown asAfter the second iterationPhased array radar allocation during sub-allocationThe search time vector optimal solution of the non-overlapping search sectors.
6f) Calculate the firstPhased array radar allocation during sub-allocationOptimal total tracking time resource of each target Wherein,is shown asThe total resources of the integrated phased array radar search and tracking application during the sub-allocation,is shown asThe duty cycle during the sub-dispensing period,indicating the duration of each dispensing period, in this embodiment
6g) The optimal total tracking time resourceSubstitute intoTracking the target standard function of the resource allocation scheme during the sub-allocation period, and obtaining the first time by utilizing a minimum maximum solution algorithmPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each target
Step 7, letAdds 1 to the value of (1), returns to step 2 until the phased array radar assignment during the 1 st assignment is obtainedOf non-overlapping search sectorsOptimal search time resourceTo the firstPhased array radar allocation during sub-allocationOptimal search time resources for non-overlapping search sectorsAnd phased array radar assignment during the 1 st assignmentTracking time resource column vector optimal solution of each targetTo the firstPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each targetRecording as a resource allocation result based on phased array radar searching and tracking dual-target optimization; this example is to get
The invention makes a resource allocation scheme into a dual-target optimization framework, and obtains a cardinal number of the resource allocation scheme according to a parallel minimization schemePareto subsets of (a); then using linear programmingThe method and the minimum and maximum solution algorithm effectively solve the problem of dual-target resource allocation.
The effects of the present invention are further verified and explained by the following simulations.
1. Simulation parameters:
the fixed position of the phased array radar is (93,45) km, and the effective signal bandwidth and the half-power beam width are respectively set to beAndthe total time resource isThe simulation was performed using a 20 dispense period, with the duration of each dispense being set toIn practice, phased array radars may have different search requirements. Thus, two types of search models are consideredAndto describe two different search requirements, see tables I and II for details, Table I is a modelTable II is a model of the search parameters for each sector inOf the search parameter per sector.
The first 10 dispensing periodsPhased array radar adopts a first search modelWherein the number of sectors is set toDuring the last 10 dispensing sessionsThe phased array radar adopts a second search modelThe number of sectors is
TABLE I
|
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
Ri,k(km) | 240 | 240 | 240 | 300 | 200 | 300 | 310 | 295 |
θi,k(○) | 8 | 8 | 8 | 8 | 8 | 6 | 8 | 12 |
TABLE II
|
1 | 2 | 3 | 4 | 5 | 6 |
Ri,k(km) | 240 | 240 | 240 | 300 | 200 | 300 |
θi,k(○) | 8 | 8 | 16 | 8 | 8 | 8 |
In the simulation, the number of targets to be tracked is set toEach target state parameter is given in table III.
TABLE III
|
1 | 2 | 3 | 4 | 5 |
Location (Km) | (3,55) | (-23,85) | (233,60) | (300,30) | (52,80) |
Speed (m/s) | (300,0) | (100,-150) | (200,-200) | (10,-200) | (200,100) |
Initial Range (m) | 127.28 | 126.32 | 169.19 | 86.45 | 122.09 |
|
2 | 0.8 | 1.5 | 1 | 1 |
Referring to fig. 3, a schematic diagram of target deployment in a detection range of a phased array radar; figure 3 shows the angular distribution of these targets with respect to the radar system.
To obtain a base number ofThe pareto sub-set of (a),individual search time budget settingFor a given Denotes a reference solution (solution of the resource equal allocation scheme) in which Is shown as havingA vector of all 1 columns of the elements,is shown asDuring the sub-distribution periodThe total resource budget is solved by a resource equal allocation scheme,is shown asDuring the sub-distribution periodThe total search resource budget is a solution to the search resource even allocation scheme,is shown asDuring the sub-distribution periodThe total tracking resource budget adopts the solution of the tracking resource average allocation scheme, and the reference set isThe objective function value of each reference solution is called the reference result, the curve formed by the reference results is called the reference curve, and the same pareto subsetThe objective function value of each pareto solution is called a pareto result, and a curve formed by the pareto results is called a pareto curve.
2. Simulation content:
the invention aims at the resource average allocation scheme and the allocation result of the dual-target optimization resource allocation scheme based on the pareto theory to compare simulation experiments.
3. And (3) simulation result analysis:
the results of FIG. 4(a) show that the demand model was searchedFor the worst case search signal-to-noise ratio and the worst case tracking bayesian cramer lower bound, the pareto curve is significantly better than the reference curve; the results of FIG. 4(b) show that the demand model was searchedFor the worst case search signal-to-noise ratio and the worst case lower bound of the Bayesian-Cramer-Role boundary, PaThe cumulative-over curve is obviously superior to the reference curve; with the pareto curve, the best search and tracking performance can be easily obtained using the dichotomy for any mission requirements.
In conclusion, the simulation experiment verifies the correctness, the effectiveness and the reliability of the method.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (7)
1. A resource allocation method based on phased array radar searching and tracking dual-target optimization is characterized by comprising the following steps:
step 1, initialization: order toIs as followsThe number of the sub-distribution is equal to the number of the sub-distribution,is set to an initial value of 1,an even number greater than 0; is set toSearch area existence of time-controlled array radar in sub-distributionAn object, andtime control of sub-distributionThe search area of the array radar is divided intoA plurality of non-overlapping search sectors;
step 2, determiningTime phased array radar in sub-distributionSearch model and method for searching sectorDuring the sub-distribution periodA tracking model of the individual target; wherein,is shown asDistributing the number of targets in a search area of the time-controlled array radar in a secondary mode;is shown asThe time-controlled array radar searches the total number of sectors in the time distribution;
step 3, according toTime phased array radar in sub-distributionA search model of the search sector, getSearching for objective function and second order of resource allocation scheme during secondary allocationSearching a conversion objective function of the resource allocation scheme during the secondary allocation period;
step 4, according toDuring the sub-distribution periodA tracking model of the object, determiningTracking a target standard function of a resource allocation scheme during the secondary allocation;
step 5, according toSearch for transfer objective function and the second of resource allocation scheme during sub-allocationTracking a target criteria function of the resource allocation scheme during the sub-allocation period to obtain a secondA mathematical optimization model of a dual target resource allocation scheme during secondary allocation;
step 6, solvingThe mathematical optimization model of the double-target resource allocation scheme in the secondary allocation period respectively obtainsPhased array radar allocation during sub-allocationOptimal search time resources for non-overlapping search sectorsAnd a firstPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each target
Step 7, letAdds 1 to the value of (1), returns to step 2 until the phased array radar assignment during the 1 st assignment is obtainedOptimal search time resources for non-overlapping search sectorsTo the firstPhased array radar allocation during sub-allocationOptimal search time for non-overlapping search sectors(Resource)And phased array radar assignment during the 1 st assignmentTracking time resource column vector optimal solution of each targetTo the firstPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each targetAnd recording the resource allocation result as a resource allocation result based on dual target searching and tracking of the phased array radar.
2. The method for resource allocation based on dual target optimization for phased array radar search and tracking as claimed in claim 1, wherein in step 2, the first stepTime phased array radar in sub-distributionThe search model for each search sector is:
wherein, alpha represents an intermediate variable,is shown asTime-controlled array radar allocation to the secondSearch sectorThe time resources of the search of (2),is shown asTime phased array radar in sub-distributionSearch sectorThe target search signal-to-noise ratio of (c),represents the average transmit power of the phased array radar,indicating the set effective receive aperture of the phased array radar antenna,is shown asTime phased array radar in sub-distributionSearch sectorThe cross-sectional area of scattering of the target,which represents the boltzmann constant, represents,indicating the set phased array radar temperature,the loss of the phased array radar is shown,is shown asTime-controlled array radar scanned by sub-distributionSearch sectorThe angle of (a) is determined,is shown asTime phased array radar in sub-distributionSearch sectorThe target distance search value of (1);
the first mentionedDuring the sub-distribution periodThe tracking model of each target is as follows:
wherein,is shown asDuring the sub-distribution periodThe process noise of the individual target is,is shown asDuring the sub-distribution periodThe state vector of the individual objects is,is as followsDuring the sub-distribution periodThe transformation matrix of the individual objects is,the expression of the kronecker operator,representing a 2 x 2 dimensional identity matrix,indicating the duration of each dispense.
3. The method for resource allocation based on dual target optimization for phased array radar search and tracking as claimed in claim 1, wherein in step 3, the first stepThe objective function of searching the resource allocation scheme during the secondary allocation is:
wherein,is shown asDuring the sub-distribution periodA set of non-overlapping search sector numbers, a constraint condition is expressed in terms of the number of the elements,is shown asPhased array radar allocation during sub-allocationThe total search time resources of the non-overlapping search sectors,is shown asDuring the second allocation period phased array radar is allocated to the secondThe search time resources of each search sector,is shown asTime phased array radar in sub-distributionSearch sectorThe target search signal-to-noise ratio of (c),is shown asPhased array radar allocation during sub-allocationA search time column vector of non-overlapping search sectors,upper labelThe transpose is represented by,
the first mentionedThe transfer objective function for searching the resource allocation scheme during the secondary allocation is:
wherein,is shown asThe convex function during the sub-allocation period,is shown asSearch sector 1 to search sector 1 during sub-allocationThe accumulated sum of the search time resources of each search sector, wherein alpha represents an intermediate variable;is shown asTime phased array radar in sub-distributionSearch sectorIs calculated from the target distance of (a) to the target distance search value,is shown asTime-controlled array radar scanned by sub-distributionSearch sectorThe angle of (d);representation calculationThe minimum value of (a) is determined,representing a set of computationsRatio of each search sector in the searchThen obtainA ratio is then comparedThe ratio value is then used to select the maximum value operation.
4. The method for resource allocation based on dual target optimization for phased array radar search and tracking as claimed in claim 1, wherein in step 4, the first stepThe objective criteria function for tracking the resource allocation scheme during the secondary allocation is:
wherein,is shown asDuring the sub-distribution periodA set of object numbers, each object number being,is shown asPhased array radar allocation during sub-allocationA tracking time resource column vector for each target,upper labelThe transpose is represented by,is shown asDuring the second allocation period phased array radar is allocated to the secondThe tracking time resources of the individual targets,is shown asPhased array radar allocation during sub-allocationTotal tracking time resources of the individual targets;expression solutionThe minimum value of (a) is determined,representing normalized worst caseDuring the sub-distribution periodAnd tracking a Bayesian Claritrol Laurve lower bound convex function of each target.
5. The method for resource allocation based on dual target optimization for phased array radar search and tracking as claimed in claim 1, wherein in step 5, the first stepThe mathematical optimization model of the double-target resource allocation scheme during the secondary allocation period comprises the following processes:
firstly, first, theThe mathematical model of the dual target resource allocation scheme during the secondary allocation is represented as:
wherein,showing obtainedAndat the same time makeAndthe size of the particles is minimized and,is shown asPhased array radar allocation during sub-allocationThe total search time resources of the non-overlapping search sectors,is shown asPhased array radar allocation during sub-allocationThe total tracking time resources of the individual targets,is shown asThe total time resources of the integrated phased array radar search and tracking application during the sub-allocation,is shown asThe duty cycle during the sub-dispensing period,is shown asDuring the second allocation period phased array radar is allocated to the secondThe search time resources of each search sector,is shown asDuring the second allocation period phased array radar is allocated to the secondThe tracking time resources of the individual targets,is shown asThe convex function during the sub-allocation period,representing normalized worst caseDuring the sub-distribution periodThe tracked Bayesian Clarithrome lower bound convex function of each target,is shown asIs divided into sub-divisionsAllocation of phased array radar to allocation periodsA search time column vector of non-overlapping search sectors,is shown asPhased array radar allocation during sub-allocationA tracking time resource column vector for each target;
then it will bePhased array radar allocation during sub-allocationSearch time vector of non-overlapping search sectorsAnd a firstPhased array radar allocation during sub-allocationTracking time resource column vector of individual targetsIntegrated into a single vector, denotedDimension vectorRepresenting a row vector transpose; thereby obtaining the firstThe mathematical optimization model of the dual-target resource allocation scheme during the secondary allocation period is as follows:
wherein,is shown asThe convex function during the sub-allocation period,representing normalized worst caseDuring the sub-distribution periodThe tracked Bayesian Clarithrome lower bound convex function of each target,is expressed as length ofAnd beforeA row vector with 1 element and zero elements,is expressed as length ofAnd from the secondElement to elementA row vector with 1 element and 0 elements,representing row vectorsTo middleThe number of the elements is one,is shown asThe total time resources of the phased array radar search and tracking application are integrated during the secondary allocation.
6. The method for resource allocation based on dual target optimization for phased array radar search and tracking as claimed in claim 1, wherein in step 6, the first stepPhased array radar allocation during sub-allocationOptimal search time resources for non-overlapping search sectorsAnd a firstPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each targetThe obtaining process comprises the following steps:
6a) separately setting iteration indexesAnd a stop threshold epsilon, and settingSearching for lower bounds of resources during secondary allocationAnd a firstSearching for upper bounds of resources during secondary allocationIs shown asThe beta total search budget during the secondary allocation,is shown asA gamma total search budget during the secondary allocation; wherein, gamma is beta +1,is shown asThe total search budget number or the total tracking budget number set in the secondary allocation period;
6b) according to the firstSearching for lower bounds of resources during secondary allocationAnd a firstSearching for upper bounds of resources during secondary allocationCalculate the firstAfter the second iterationTotal search resources for phased array radar during secondary allocationThen use the firstAfter the second iterationTotal search resources for phased array radar during secondary allocationAnd linear programming method to solveSearching a conversion objective function of the resource allocation scheme during the sub-allocation period to obtain a first stepAfter the second iterationPhased array radar allocation during sub-allocationSearch time column vector optimal solution for non-overlapping search sectors
6c) If it is notUpdatingOrder toAdd 1 to the value of (6 b); otherwise 6d) is executed,is shown asAfter the second iterationPhased array radar allocation during sub-allocationFunction values of the search time vector optimal solution of the non-overlapping search sectors;is shown asDuring the sub-distribution periodWorst case search signal-to-noise ratios in non-overlapping search sectors; a represents an intermediate variable which is,is shown asTime phased array radar in sub-distributionSearch sectorIs calculated from the target distance of (a) to the target distance search value,is shown asTime-controlled array radar scanned by sub-distributionSearch sectorThe angle of (d);
6e) if it is notGetSearch for the optimal solution for resource allocation asIs shown asPhased array radar allocation during sub-allocationThe optimal search time resources of the non-overlapping search sectors,is shown asAfter the second iterationPhased array radar allocation during sub-allocationSearch time vector optimal solutions of non-overlapping search sectors;
6f) calculate the firstPhased array radar allocation during sub-allocationOptimal total tracking time resource of each target Wherein,is shown asThe total resources of the integrated phased array radar search and tracking application during the sub-allocation,is shown asThe duty cycle during the sub-dispensing period,indicating the duration of each dispensing period;
6g) tracking time resources according to the optimal total tracking time resourcesAnd minimum maximum solution methodTracking a target criteria function of the resource allocation scheme during the sub-allocation period to obtain a secondPhased array radar allocation during sub-allocationTracking time resource column vector optimal solution of each target
7. The method for resource allocation based on dual target optimization for phased array radar search and tracking according to claim 6, wherein in 6(a), the method comprisesIs shown asBeta total search budget during secondary allocation and theIs shown asA γ th total search budget during the secondary allocation, further comprising:
is shown asThe beta total search budget during the secondary allocation period is the second of the dual target resourcesPareto optimal solution of the beta of the sub-distributionMiddle frontA cumulative sum of the elements;represents the gamma total search budget during the kth allocation, is the second target resourceSub-assigned [ gamma ] pareto optimal solutionMiddle frontA cumulative sum of the elements;
the two target resources are treated asSecond of the sub distributionThe pareto optimal solution is recorded asThe obtaining process comprises the following steps:
6.1 setting ofSub-distribution periodThere is aTotal search budget andtotal tracking budget from 1 st total search budget toThe total search budget satisfies:whereinIs shown asDuring the sub-distribution periodTotal search budget, willDuring the sub-distribution periodTotal tracking budget asAnd the firstDuring the sub-distribution periodTotal search budgetAnd a firstDuring the sub-distribution periodTotal tracking budgetSatisfies the following conditions:
whereinIs shown asIntegrating the total time resources of phased array radar search and tracking applications during the secondary allocation;
6.2 according to the firstDuring the sub-distribution periodTotal search budgetAnd linear programming method to solveSearching a conversion objective function of the resource allocation scheme during the sub-allocation period to obtain a first stepPhased array radar allocation during sub-allocationSearch time vector pair of non-overlapping search sectorsOptimal solution for individual search budget resource allocationAccording to the firstDuring the sub-distribution periodTotal tracking budgetAnd maximum and minimum solution algorithm solvingTracking a target criteria function of the resource allocation scheme during the sub-allocation period to obtain a secondPhased array radar allocation during sub-allocationTracking time resource column vector pair of individual targetsOptimal solution for individual total tracking budget resource allocationWill be the firstOptimal solution for individual search budget resource allocationAnd the said firstOptimal solution for individual total tracking budget resource allocationForm two target resourcesSecond of the sub distributionPareto optimal solutionUpper labelSecond to indicate transposed, dual target resource allocationPareto optimal solutionIncludedAnd (4) each element.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810057487.5A CN108333583B (en) | 2018-01-22 | 2018-01-22 | Resource allocation method based on dual-target optimization of phased array radar search and tracking |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810057487.5A CN108333583B (en) | 2018-01-22 | 2018-01-22 | Resource allocation method based on dual-target optimization of phased array radar search and tracking |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108333583A CN108333583A (en) | 2018-07-27 |
CN108333583B true CN108333583B (en) | 2021-08-03 |
Family
ID=62926491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810057487.5A Active CN108333583B (en) | 2018-01-22 | 2018-01-22 | Resource allocation method based on dual-target optimization of phased array radar search and tracking |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108333583B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112034448B (en) * | 2020-08-10 | 2022-11-04 | 西安电子科技大学 | Networking radar resource allocation optimization method based on tracking precision and resource limitation |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7555383B2 (en) * | 2003-05-28 | 2009-06-30 | Northrop Grumman Corporation | Target acquisition and tracking system |
CN104077488A (en) * | 2014-07-05 | 2014-10-01 | 中国船舶重工集团公司第七二四研究所 | Rotary phased array radar sliding window resource scheduling technique based on sectors |
CN105046412A (en) * | 2015-06-29 | 2015-11-11 | 中国船舶重工集团公司第七二四研究所 | Passive phased array radar multistation joint resource scheduling and distribution method |
CN106199579A (en) * | 2016-06-22 | 2016-12-07 | 中国人民解放军信息工程大学 | Distributed MIMO radar target tracking precision method for joint optimization of resources |
CN106990399A (en) * | 2017-05-11 | 2017-07-28 | 西安电子科技大学 | Radar network system power and bandwidth combined distributing method for target following |
CN107167798A (en) * | 2017-05-05 | 2017-09-15 | 电子科技大学 | Based on the cognitive tracking of the controllable many radars of tracking accuracy |
CN107450070A (en) * | 2017-04-14 | 2017-12-08 | 电子科技大学 | Phased-array radar wave beam and residence time combined distributing method based on target following |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8634982B2 (en) * | 2009-08-19 | 2014-01-21 | Raytheon Company | System and method for resource allocation and management |
-
2018
- 2018-01-22 CN CN201810057487.5A patent/CN108333583B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7555383B2 (en) * | 2003-05-28 | 2009-06-30 | Northrop Grumman Corporation | Target acquisition and tracking system |
CN104077488A (en) * | 2014-07-05 | 2014-10-01 | 中国船舶重工集团公司第七二四研究所 | Rotary phased array radar sliding window resource scheduling technique based on sectors |
CN105046412A (en) * | 2015-06-29 | 2015-11-11 | 中国船舶重工集团公司第七二四研究所 | Passive phased array radar multistation joint resource scheduling and distribution method |
CN106199579A (en) * | 2016-06-22 | 2016-12-07 | 中国人民解放军信息工程大学 | Distributed MIMO radar target tracking precision method for joint optimization of resources |
CN107450070A (en) * | 2017-04-14 | 2017-12-08 | 电子科技大学 | Phased-array radar wave beam and residence time combined distributing method based on target following |
CN107167798A (en) * | 2017-05-05 | 2017-09-15 | 电子科技大学 | Based on the cognitive tracking of the controllable many radars of tracking accuracy |
CN106990399A (en) * | 2017-05-11 | 2017-07-28 | 西安电子科技大学 | Radar network system power and bandwidth combined distributing method for target following |
Non-Patent Citations (2)
Title |
---|
Simultaneous Multibeam Resource Allocation Scheme for Multiple Target Tracking;Junkun Yan等;《IEEE Transactions on Signal Processing》;20150326;第3110-3122页 * |
相控阵雷达多目标跟踪模式波束资源管理与优化研究;牛俊翔;《中国优秀硕士学位论文全文数据库信息科技辑》;20170315;第43-60页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108333583A (en) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yan et al. | Joint beam selection and power allocation for multiple target tracking in netted colocated MIMO radar system | |
Daun et al. | Tracking in multistatic passive radar systems using DAB/DVB-T illumination | |
Chai et al. | Reference tag supported RFID tracking using robust support vector regression and Kalman filter | |
US7671789B1 (en) | Method and system for target detection and angle estimation based on a radar signal | |
Chen et al. | A modified probabilistic data association filter in a real clutter environment | |
CN109752710B (en) | Rapid target angle estimation method based on sparse Bayesian learning | |
CN109581354B (en) | Multi-target tracking resource management method for simultaneous multi-beam co-location MIMO radar | |
CN108761381B (en) | Target information source detection method based on array antenna | |
CN108168564A (en) | A kind of Data Association based on LHD grey relational grades | |
CN106842184A (en) | A kind of multiple target detection and tracking based on beam dispath | |
CN106021697B (en) | A kind of rapid phase-control battle array radar Time-energy resource joint management method | |
CN109581355B (en) | Centralized MIMO radar self-adaptive resource management method for target tracking | |
CN105137418A (en) | Multi-object tracking and data interconnection method based on whole neighborhood fuzzy clustering | |
KR20220023749A (en) | Co-optimization method of integrated network radar dwell time and radiated power for low detection | |
CN110244273A (en) | It is a kind of based on the target angle estimation method for being uniformly distributed formula array | |
CN107229032A (en) | A kind of method and apparatus for building four array element solid arrays | |
CN110376580A (en) | A kind of isomery radar fence resource allocation methods of asynchronous multiple target tracking performance driving | |
CN107037398A (en) | A kind of two-dimentional MUSIC algorithms estimate the parallel calculating method of direction of arrival | |
Chen et al. | Iterative subspace alternating projection method for GNSS multipath DOA estimation | |
CN111274529B (en) | Robust Gao Sini Weisal PHD multi-expansion target tracking algorithm | |
CN108333583B (en) | Resource allocation method based on dual-target optimization of phased array radar search and tracking | |
CN116148837A (en) | Robust multi-sensor multi-target tracking method under heterogeneous clutter background | |
CN113325410A (en) | Radar antenna signal processing method and device, control equipment and storage medium | |
CN108535704A (en) | A kind of signal Pre-sorting method based on self-adaption two-dimensional cluster | |
CN110488276B (en) | Heterogeneous radar network optimal resource on-demand distribution method oriented to multi-target tracking task |
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 |