CN112213718B - Networking radar node selection and radiation resource joint optimization method under multi-target tracking - Google Patents

Networking radar node selection and radiation resource joint optimization method under multi-target tracking Download PDF

Info

Publication number
CN112213718B
CN112213718B CN202011022261.5A CN202011022261A CN112213718B CN 112213718 B CN112213718 B CN 112213718B CN 202011022261 A CN202011022261 A CN 202011022261A CN 112213718 B CN112213718 B CN 112213718B
Authority
CN
China
Prior art keywords
radar
target
networking
matrix
radiation
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
Application number
CN202011022261.5A
Other languages
Chinese (zh)
Other versions
CN112213718A (en
Inventor
时晨光
丁琳涛
王奕杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202011022261.5A priority Critical patent/CN112213718B/en
Publication of CN112213718A publication Critical patent/CN112213718A/en
Application granted granted Critical
Publication of CN112213718B publication Critical patent/CN112213718B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/66Radar-tracking systems; Analogous systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/87Combinations of radar systems, e.g. primary radar and secondary radar

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 networking radar node selection and radiation resource joint optimization method under multi-target tracking, which comprises the following steps: s1, determining the composition of a networking radar system and a multi-target tracking task; s2, constructing a predicted Bayesian-Luo lower bound matrix of the target motion state estimation error by taking radar node selection, each radar dwell time and radiation power as independent variables, and taking the trace sum of the predicted Bayesian-Luo lower bound matrix of each target as a measurement index of multi-target tracking accuracy; s3, establishing a networking radar node selection and radiation resource joint optimization model under multi-target tracking; and S4, solving the networking radar node selection and radiation resource joint optimization model under multi-target tracking by using a two-step optimization algorithm. The method reduces the estimation error of the networking radar system on the motion state of each target, and effectively improves the multi-target tracking precision.

Description

Networking radar node selection and radiation resource joint optimization method under multi-target tracking
Technical Field
The invention relates to a radar signal processing technology, in particular to a networking radar node selection and radiation resource joint optimization method under multi-target tracking.
Background
Radio frequency radiation resources are very important resources in radar systems. Therefore, the radio frequency radiation resource optimal allocation in a multi-target tracking scene or the radio frequency radiation resource optimal allocation in a networking radar tracking scene is a research hotspot problem at present.
The networking radar is a radar system of a new system developed on the basis of the traditional phased array radar, can realize the independent control of each radar transmitter and receiver, and synthesizes expected beams according to the total amount of radar system resources, battlefield environment, task types, quantity and the like so as to realize the organic allocation of the radio frequency radiation resources of the networking radar system.
In order to improve the utilization rate of radar radio frequency radiation resources, research institutions and colleges at home and abroad have already done a great deal of work, and a solid theoretical foundation is laid for subsequent research. However, no method for joint optimization of networking radar node selection and radiation resource under multi-target tracking exists in the prior art.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a networking radar node selection and radiation resource joint optimization method under multi-target tracking, which takes the radiation resources of a networking radar system as constraint conditions, takes the sum of traces of Bayesian-Luo lower bound matrixes predicted by each target at each moment as an optimization target, and establishes a networking radar node selection and radiation resource joint optimization model under multi-target tracking, thereby reducing the motion state estimation error of each target by the networking radar system and effectively improving the multi-target tracking precision.
The technical scheme is as follows: the invention discloses a networking radar node selection and radiation resource joint optimization method under multi-target tracking, which comprises the following steps:
s1, determining the composition of a networking radar system and a multi-target tracking task;
s2, constructing a prediction Bayes-Lame lower bound matrix of the target motion state estimation error by taking radar node selection, radar residence time and radiation power as independent variables, and taking the track sum of the prediction Bayes-Lame lower bound matrix of each target as a measurement index of multi-target tracking accuracy;
s3, establishing a networking radar node selection and radiation resource joint optimization model under multi-target tracking;
and S4, solving the networking radar node selection and radiation resource joint optimization model under multi-target tracking by using a two-step optimization algorithm.
Furthermore, in the step S1, a networking radar system formed by M phased array radars is considered to track Q targets, the M phased array radars can keep accurate time, space and frequency synchronization, and each phased array radar can only receive and process target echoes from its own transmitted signal.
Further, step S2 specifically includes:
the Bayesian information matrix calculation expression of the qth target is as follows:
Figure BDA0002701012070000021
wherein, (. Cndot.) -1 The inverse operation of the matrix is represented,
Figure BDA0002701012070000022
a Bayesian information matrix for a qth target;
Figure BDA0002701012070000023
predicting the motion state vector of the qth object at time k for time k-1, wherein T Representing a transpose operation of a matrix or vector,
Figure BDA0002701012070000024
indicating the location of the qth object at time k-1 at which time k is predicted,
Figure BDA0002701012070000025
representing the motion speed of the qth target at the predicted k moment at the k-1 moment;
Figure BDA0002701012070000026
selecting variables for radar nodes when
Figure BDA0002701012070000027
When the target is shot and tracked by the nth radar at the time point k, the shot and the tracked target are shot and tracked
Figure BDA0002701012070000028
Then, the nth radar does not perform irradiation tracking on the qth target at the time k; q is a process noise covariance matrix, and the mathematical expression is as follows:
Figure BDA0002701012070000029
wherein, T 0 In order to be the sampling interval of the sample,
Figure BDA00027010120700000210
representing a matrix direct product operation, I 2 Is an identity matrix of order 2, r q Process noise strength for the qth target; f is a target state transition matrix, and the mathematical expression of the target state transition matrix is as follows:
Figure BDA00027010120700000211
Figure BDA00027010120700000212
as a function of the nth radar measurement
Figure BDA00027010120700000213
The jacobian matrix of (a), wherein,
Figure BDA00027010120700000214
representing the q-th target state vector
Figure BDA00027010120700000215
First order partial derivative, nth radar measurement function is obtained
Figure BDA00027010120700000216
The mathematical expression of (a) is:
Figure BDA0002701012070000031
wherein (x) n ,y n ) Position coordinates of the nth part of radar in the space;
Figure BDA0002701012070000032
the measured noise covariance matrix of the nth radar to the qth target is expressed by the following mathematical expression:
Figure BDA0002701012070000033
wherein, c =3 × 10 8 m/s, B is effective bandwidth of radar emission signal, lambda is radar working wavelength, and gamma isThe aperture of the antenna is set to be,
Figure BDA0002701012070000034
predicting the echo signal-to-noise ratio of the k moment to the q target for the nth radar k-1 moment, wherein the mathematical expression is as follows:
Figure BDA0002701012070000035
wherein,
Figure BDA0002701012070000036
and
Figure BDA0002701012070000037
respectively the dwell time and the radiation power of the nth radar irradiating the qth target at the time k, T r For each radar pulse repetition period, G t And G r For each radar transmit antenna gain and receive antenna gain,
Figure BDA0002701012070000038
radar cross section of the qth target relative to the nth radar, G RP Processing gain, k, for radar receivers 0 And T 0 Respectively Boltzmann constant and noise temperature of each radar receiver, B r Matching the filter bandwidth, F, for each radar receiver r Noise figure for each radar receiver;
Figure BDA0002701012070000039
predicting the distance between the k moment and the q target for the nth radar k-1 moment;
the formula (1) is inverted, namely a Bayesian Classmei-Rou lower bound matrix for predicting the estimation error of the q-th target motion state is obtained, and the mathematical expression is as follows:
Figure BDA00027010120700000310
here, bayesian clar is predicted using each targetMei-Luo lower bound matrix
Figure BDA0002701012070000041
The multi-target tracking precision is characterized by the sum of the traces of (1), namely:
Figure BDA0002701012070000042
wherein, tr (-) represents the trace operation of the matrix.
Further, step S3 specifically includes:
the method comprises the following steps of taking radiation resources of a networking radar system as constraint conditions, taking the sum of traces of a Bayesian Classman-Luo lower bound matrix for minimizing target prediction at each moment as an optimization target, and establishing a networking radar node selection and radiation resource combined optimization model under multi-target tracking as follows:
Figure BDA0002701012070000043
wherein, T d,min And T d,max Respectively the lower limit and the upper limit of the residence time of each radar; p min And P max Respectively the lower limit and the upper limit of the radiation power of each radar; t is a unit of total Is the total dwell time of the networked radar; p total Is the total radiated power of the networking radar; l is max And the number of radar nodes for carrying out irradiation tracking on the q-th target at the time k.
Further, step S4 specifically includes:
first, select
Figure BDA0002701012070000044
L with the largest value max The radar irradiates and tracks the q-th target and makes the corresponding target
Figure BDA0002701012070000045
Predicting the distance between the k moment and the q target for the nth radar k-1 moment;
and secondly, combining an interior point method and a cycle minimum method to obtain the residence time and radiation power optimal distribution result of each radar under the selection of the given radar node.
Furthermore, the specific steps of obtaining the residence time and the radiation power optimal distribution result of each radar under the selection of the given radar node by combining the interior point method and the cyclic minimum method are as follows:
(a) For the obtained radar node selection result, the networking radar node selection and radiation resource joint optimization model under multi-target tracking is simplified as follows:
Figure BDA0002701012070000051
(b) Order to
Figure BDA0002701012070000052
Adopting an interior point method in MATLAB software to call a function fmincon to calculate and solve the model simplified formula (10), namely obtaining an optimal solution of radiation power distribution when the residence time distribution of the networking radar is fixed
Figure BDA0002701012070000053
(c) Order to
Figure BDA0002701012070000054
Adopting an interior point method in MATLAB software to call a function fmincon to calculate and solve the model simplified formula (10), namely obtaining an optimal solution of residence time distribution when the distribution of the networking radar radiation power is fixed
Figure BDA0002701012070000055
(d) An optimal solution for assigning residence time
Figure BDA0002701012070000056
Substituting into a model simplified formula (10), jumping to the step (a), and circulating in sequence until the difference of the multi-target tracking precision obtained by two times of iterative calculation is less than a fixed value epsilon 0 At this time can be obtainedAnd determining the residence time and the radiation power optimization distribution result of each radar under the selection of the radar nodes.
Has the advantages that: compared with the prior art, the invention has the following advantages:
(1) The invention provides a method for joint optimization of node selection and radiation resources of networking radar under multi-target tracking, which completes the main task of tracking multiple targets by considering a networking radar system consisting of a plurality of phased array radars, wherein the radars can keep accurate time, space and frequency synchronization, and each radar can only receive and process target echoes of self-transmitted signals. Firstly, constructing a prediction Bayesian-Lame lower bound matrix of a target motion state estimation error by taking radar node selection, radar residence time and radiation power as independent variables, and taking the sum of trails of the prediction Bayesian-Lame lower bound matrix of each target as a measurement index of multi-target tracking accuracy; secondly, by taking the radiation resources of the networking radar system as constraint conditions, and by taking the sum of the traces of the Bayesian-Luo lower bound matrix for minimizing the prediction of each target at each moment as an optimization target, a networking radar node selection and radiation resource combined optimization model under multi-target tracking is established, so that the motion state estimation error of each target by the networking radar system is reduced, and the multi-target tracking precision is effectively improved.
The method has the advantages that the radiation resource of the networking radar system can be met, the motion state estimation error of the networking radar system on each target can be reduced, and the multi-target tracking precision of the networking radar system is effectively improved. The method takes the radiation resources of a networking radar system as a constraint condition, takes the sum of the traces of a Bayesian-Ro lower bound matrix for each target prediction at each moment as an optimization target, and establishes a networking radar node selection and radiation resource combined optimization model under multi-target tracking. By adopting a two-step optimization algorithm to solve the optimization model, a networking radar node selection, residence time and radiation power combined optimization result meeting constraint conditions can be obtained, and therefore the multi-target tracking precision of the networking radar system is effectively improved.
(2) Compared with the prior art, the networking radar node selection and radiation resource joint optimization method under multi-target tracking provided by the invention not only can meet the radiation resource requirement of a networking radar system, but also can reduce the motion state estimation error of each target of the networking radar system, thereby effectively improving the multi-target tracking precision.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention is further described with reference to the following figures and specific examples.
The invention provides a networking radar node selection and radiation resource combined optimization method under multi-target tracking based on the current military application requirements, and establishes a networking radar node selection and radiation resource combined optimization model under multi-target tracking by taking the radiation resources of a networking radar system as constraint conditions and minimizing the sum of the traces of Bayesian-Rale lower bound matrixes predicted by each target at each moment as an optimization target, so that the motion state estimation error of each target by the networking radar system is reduced, and the multi-target tracking precision of the networking radar system is effectively improved.
As shown in fig. 1, the method for jointly optimizing networking radar node selection and radiation resources under multi-target tracking of the present invention includes the following steps:
s1, determining the composition of a networking radar system and a multi-target tracking task;
a networking radar system consisting of M phased array radars is considered to track Q targets, the radars can keep accurate time, space and frequency synchronization, and each radar can only receive and process target echoes from own transmitted signals.
S2, constructing a predicted Bayesian-Luo lower bound matrix of the target motion state estimation error by taking radar node selection, radar dwell time and radiation power as independent variables, and taking the track sum of the predicted Bayesian-Luo lower bound matrix of each target as a measurement index of the multi-target tracking accuracy, wherein the measurement index is as follows:
the Bayesian information matrix calculation expression of the qth target is as follows:
Figure BDA0002701012070000071
wherein, (.) -1 The inverse operation of the matrix is represented,
Figure BDA0002701012070000072
a Bayesian information matrix for a qth target;
Figure BDA0002701012070000073
predicting the motion state vector of the qth target at time k for time k-1, wherein · T Representing a transpose operation of a matrix or vector,
Figure BDA0002701012070000074
indicating the position of the qth object at time k-1,
Figure BDA0002701012070000075
representing the motion speed of the qth target at the predicted k moment at the k-1 moment;
Figure BDA0002701012070000076
selecting variables for radar nodes when
Figure BDA0002701012070000077
When the target is located at the k time, the nth radar performs irradiation tracking on the qth target
Figure BDA0002701012070000078
Then, the nth radar does not perform irradiation tracking on the qth target at the time k; q is a process noise covariance matrix, and the mathematical expression of the process noise covariance matrix is as follows:
Figure BDA0002701012070000079
wherein, T 0 In order to be the sampling interval of the sample,
Figure BDA00027010120700000710
representing a matrix direct product operation, I 2 Is an identity matrix of order 2, r q Process noise strength for the qth target; f is a target state transition matrix, and the mathematical expression of the target state transition matrix is as follows:
Figure BDA00027010120700000711
Figure BDA00027010120700000712
as a function of nth radar measurement
Figure BDA00027010120700000713
The jacobian matrix of (a), wherein,
Figure BDA00027010120700000714
representing the q-th target state vector
Figure BDA00027010120700000715
First order partial derivative, nth radar measurement function is obtained
Figure BDA00027010120700000716
The mathematical expression of (a) is:
Figure BDA0002701012070000081
wherein (x) n ,y n ) Position coordinates of the nth part of radar in the space;
Figure BDA0002701012070000082
the measured noise covariance matrix of the nth radar to the qth target is expressed by the following mathematical expression:
Figure BDA0002701012070000083
wherein, c =3 × 10 8 m/s, B is the effective bandwidth of the radar emission signal, lambda is the radar operating wavelength, gamma is the antenna aperture,
Figure BDA0002701012070000084
predicting the echo signal-to-noise ratio of the k moment to the q target at the k-1 moment of the nth radar, wherein the mathematical expression is as follows:
Figure BDA0002701012070000085
wherein,
Figure BDA0002701012070000086
and
Figure BDA0002701012070000087
respectively the dwell time and the radiation power of the nth radar irradiating the qth target at the time k, T r For each radar pulse repetition period, G t And G r For each radar transmit antenna gain and receive antenna gain,
Figure BDA0002701012070000088
radar cross section of the qth target relative to the nth radar, G RP Processing gain, k, for radar receivers 0 And T 0 Respectively Boltzmann constant and noise temperature of each radar receiver, B r Matching the filter bandwidth, F, for each radar receiver r For each of the radar receiver noise figure(s),
Figure BDA0002701012070000089
predicting the distance between the k moment and the q target for the nth radar k-1 moment; .
The formula (1) is inverted, so that a Bayesian Classmen-Rou lower bound matrix for predicting the estimation error of the q-th target motion state can be obtained, and the mathematical expression is as follows:
Figure BDA00027010120700000810
here, a Bayesian-Lame lower bound matrix is predicted using each target
Figure BDA0002701012070000091
The multi-target tracking precision is characterized by the sum of the traces of (1), namely:
Figure BDA0002701012070000092
wherein, tr (·) represents a trace operation of matrix calculation.
S3, establishing a networking radar node selection and radiation resource joint optimization model under multi-target tracking, as follows:
by taking the radiation resources of the networking radar system as constraint conditions and taking the sum of the traces of the Bayesian-Ro lower bound matrix for minimizing the prediction of each target at each moment as an optimization target, a networking radar node selection and radiation resource combined optimization model under multi-target tracking is established, and is as follows:
Figure BDA0002701012070000093
wherein, T d,min And T d,max Respectively the lower limit and the upper limit of the residence time of each radar; p min And P max Respectively the lower limit and the upper limit of the radiation power of each radar; t is total Is the total dwell time of the networked radar; p total Is the total radiated power of the networking radar; l is max And the number of radar nodes for carrying out irradiation tracking on the q-th target at the time k.
S4, solving a networking radar node selection and radiation resource combined optimization model under multi-target tracking by using a two-step optimization algorithm, namely a formula (9):
first, select
Figure BDA0002701012070000094
L with the largest value max The radar irradiates the q-th targetShoot-follow and order the corresponding
Figure BDA0002701012070000095
And secondly, combining an interior point method and a cycle minimum method to obtain the residence time and radiation power optimal distribution result of each radar under the selection of the given radar node.
(a) For the obtained radar node selection result, the model for jointly optimizing networking radar node selection and radiation resources under multi-target tracking, namely the formula (9), can be simplified as follows:
Figure BDA0002701012070000101
(b) Order to
Figure BDA0002701012070000102
Adopting an interior point method in MATLAB software to call a function fmincon to calculate and solve the model simplified formula (10), and obtaining an optimal solution of radiation power distribution when the residence time distribution of the networking radar is fixed
Figure BDA0002701012070000103
(c) Order to
Figure BDA0002701012070000104
Adopting an interior point method in MATLAB software to call a function fmincon to calculate and solve the model simplified formula (10), namely obtaining an optimal solution of residence time distribution when the distribution of the networking radar radiation power is fixed
Figure BDA0002701012070000105
(d) An optimal solution for assigning residence time
Figure BDA0002701012070000106
Substituting into a model simplified formula (10), jumping to the step (a), and circulating in sequence until the difference of the multi-target tracking precision obtained by two times of iterative calculationLess than a fixed value epsilon 0 And at the moment, the residence time and radiation power optimization distribution result of each radar under the selection of the given radar node can be obtained.
The working principle and the working process of the invention are as follows:
the invention considers that a networking radar system formed by a plurality of phased array radars tracks multiple targets, the radars can keep accurate time, space and frequency synchronization, and each radar can only receive and process target echoes of self-transmitted signals. Firstly, constructing a predicted Bayesian-Luo lower bound matrix of a target motion state estimation error by taking radar node selection, each radar dwell time and radiation power as independent variables, and taking the trace sum of the predicted Bayesian-Luo lower bound matrix of each target as a measurement index of multi-target tracking accuracy; secondly, by taking the radiation resources of the networking radar system as constraint conditions and taking the sum of the traces of the Bayesian-Rauwolf lower bound matrix for predicting each target at each moment as an optimization target, a networking radar node selection and radiation resource combined optimization model under multi-target tracking is established, so that the motion state estimation error of each target by the networking radar system is reduced, and the multi-target tracking precision of the networking radar system is effectively improved; and finally, solving the optimization model by adopting a two-step optimization algorithm to obtain a networking radar node selection, residence time and radiation power combined optimization result which meets the constraint condition.
The invention is characterized in that:
1. tracking multiple targets by a networking radar system consisting of multiple phased array radars, wherein the radars can keep accurate time, space and frequency synchronization, and each radar can only receive and process target echoes of a signal transmitted by the radar; constructing a prediction Bayesian-Lame lower bound matrix of a target motion state estimation error by taking radar node selection, radar residence time and radiation power as independent variables, and taking the sum of traces of the prediction Bayesian-Lame lower bound matrix of each target as a measurement index of multi-target tracking accuracy;
2. the method comprises the steps of taking radiation resources of a networking radar system as constraint conditions, taking the sum of the traces of a Bayesian-Luo lower bound matrix of target prediction at each moment as an optimization target, establishing a networking radar node selection and radiation resource combined optimization model under multi-target tracking, solving the optimization model by adopting a two-step optimization algorithm, and determining the networking radar node selection, residence time and radiation power which enable the sum of the traces of the Bayesian-Luo lower bound matrix of target prediction to be minimized as an optimal solution.

Claims (2)

1. The method for jointly optimizing networking radar node selection and radiation resources under multi-target tracking is characterized by comprising the following steps:
s1, determining the composition of a networking radar system and a multi-target tracking task, and tracking Q targets by considering the networking radar system consisting of M phased array radars;
s2, constructing a predicted Bayesian-Luo lower bound matrix of the target motion state estimation error by taking radar node selection, each radar dwell time and radiation power as independent variables, and taking the trace sum of the predicted Bayesian-Luo lower bound matrix of each target as a measurement index of multi-target tracking accuracy; the method specifically comprises the following steps:
the Bayesian information matrix calculation expression of the qth target is as follows:
Figure FDA0003735113530000011
wherein, (.) -1 The inverse operation of the matrix is represented,
Figure FDA0003735113530000012
a Bayesian information matrix for a qth target;
Figure FDA0003735113530000013
predicting the motion state vector of the qth object at time k for time k-1, wherein T Representing a transpose operation of a matrix or vector,
Figure FDA0003735113530000014
indicating the location of the qth object at time k-1 at which time k is predicted,
Figure FDA0003735113530000015
representing the motion speed of the qth target at the predicted k moment at the k-1 moment;
Figure FDA0003735113530000016
selecting variables for radar nodes when
Figure FDA0003735113530000017
When the target is located at the k time, the nth radar performs irradiation tracking on the qth target
Figure FDA0003735113530000018
When the target is shot, the nth radar does not carry out irradiation tracking on the qth target at the moment k; q is a process noise covariance matrix, and the mathematical expression is as follows:
Figure FDA0003735113530000019
wherein, T 0 In order to be the sampling interval of the sample,
Figure FDA00037351135300000110
representing a matrix direct product operation, I 2 Is an identity matrix of order 2, r q Process noise strength for the qth target; f is a target state transition matrix, and the mathematical expression of the target state transition matrix is as follows:
Figure FDA00037351135300000111
Figure FDA0003735113530000021
as a function of the nth radar measurement
Figure FDA0003735113530000022
The jacobian matrix of (a), wherein,
Figure FDA0003735113530000023
representing the q-th target state vector
Figure FDA0003735113530000024
To find the first order partial derivative, the nth radar measurement function
Figure FDA0003735113530000025
The mathematical expression of (a) is:
Figure FDA0003735113530000026
wherein (x) n ,y n ) Position coordinates of the nth part of radar in the space;
Figure FDA0003735113530000027
the measured noise covariance matrix of the nth radar to the qth target is expressed by the following mathematical expression:
Figure FDA0003735113530000028
wherein, c =3 × 10 8 m/s, B is the effective bandwidth of the radar emission signal, lambda is the radar operating wavelength, gamma is the antenna aperture,
Figure FDA0003735113530000029
predicting the echo signal-to-noise ratio of the k moment to the q target at the k-1 moment of the nth radar, wherein the mathematical expression is as follows:
Figure FDA00037351135300000210
wherein,
Figure FDA00037351135300000211
and
Figure FDA00037351135300000212
respectively the dwell time and the radiation power of the nth radar irradiating the qth target at the time k, T r For each radar pulse repetition period, G t And G r For each radar transmit antenna gain and receive antenna gain,
Figure FDA00037351135300000213
radar cross section of the qth target relative to the nth radar, G RP Processing gain, k, for radar receivers 0 And T 0 Respectively Boltzmann constant and noise temperature of each radar receiver, B r Matching the filter bandwidth, F, for each radar receiver r Noise figure for each radar receiver;
Figure FDA00037351135300000214
predicting the distance between the k moment and the q target for the nth radar k-1 moment;
the formula (1) is inverted, namely a Bayesian Classmei-Rou lower bound matrix for predicting the estimation error of the q-th target motion state is obtained, and the mathematical expression is as follows:
Figure FDA0003735113530000031
here, a Bayesian-Lame lower bound matrix is predicted using each target
Figure FDA0003735113530000032
The multi-target tracking precision is characterized by the sum of the traces of (1), namely:
Figure FDA0003735113530000033
wherein, tr (-) represents the trace operation of matrix calculation;
s3, establishing a networking radar node selection and radiation resource joint optimization model under multi-target tracking; the method specifically comprises the following steps:
the method comprises the following steps of taking radiation resources of a networking radar system as constraint conditions, taking the sum of traces of a Bayesian Classman-Luo lower bound matrix for minimizing target prediction at each moment as an optimization target, and establishing a networking radar node selection and radiation resource combined optimization model under multi-target tracking as follows:
Figure FDA0003735113530000034
wherein, T d,min And T d,max Respectively the lower limit and the upper limit of the residence time of each radar; p min And P max Respectively the lower limit and the upper limit of the radiation power of each radar; t is a unit of total Total residence time for the networked radar; p total Is the total radiated power of the networking radar; l is max The number of radar nodes for performing irradiation tracking on the q-th target at the time k;
s4, solving a networking radar node selection and radiation resource joint optimization model under multi-target tracking by using a two-step optimization algorithm; the method comprises the following specific steps:
first, select
Figure FDA0003735113530000041
L with the largest value max The radar irradiates and tracks the q-th target and makes the corresponding target
Figure FDA0003735113530000042
Figure FDA0003735113530000043
Predicting the distance between the k moment and the q target for the nth radar k-1 moment;
secondly, combining an interior point method and a cycle minimum method to obtain the residence time and radiation power optimal distribution result of each radar under the selection of the given radar node;
the method comprises the following specific steps of combining an interior point method and a cycle minimum method to obtain the residence time and radiation power optimal distribution result of each radar under the condition that a given radar node is selected:
(a) For the obtained radar node selection result, the networking radar node selection and radiation resource joint optimization model under multi-target tracking is simplified as follows:
Figure FDA0003735113530000044
(b) Order to
Figure FDA0003735113530000045
Adopting an interior point method in MATLAB software to call a function fmincon to calculate and solve the model simplified formula (10), namely obtaining an optimal solution of radiation power distribution when the residence time distribution of the networking radar is fixed
Figure FDA0003735113530000046
(c) Order to
Figure FDA0003735113530000047
Adopting an interior point method in MATLAB software to call a function fmincon to calculate and solve the model simplified formula (10), namely obtaining an optimal solution of residence time distribution when the distribution of the networking radar radiation power is fixed
Figure FDA0003735113530000048
(d) An optimal solution for assigning residence time
Figure FDA0003735113530000049
Substituting the model into a simplified formula (10), skipping to the step (a), and circulating in sequence until the difference of the multi-target tracking precision obtained by two times of iterative computation is less than a fixed value epsilon 0 At this time, can be obtainedAnd determining the residence time and the radiation power optimization distribution result of each radar under the selection of the radar nodes.
2. The method for joint optimization of node selection and radiation resource of networking radar under multi-target tracking according to claim 1, wherein in the step S1, the M phased array radars can keep accurate time, space and frequency synchronization, and each phased array radar can only receive and process target echoes from its own transmitted signal.
CN202011022261.5A 2020-09-25 2020-09-25 Networking radar node selection and radiation resource joint optimization method under multi-target tracking Active CN112213718B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011022261.5A CN112213718B (en) 2020-09-25 2020-09-25 Networking radar node selection and radiation resource joint optimization method under multi-target tracking

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011022261.5A CN112213718B (en) 2020-09-25 2020-09-25 Networking radar node selection and radiation resource joint optimization method under multi-target tracking

Publications (2)

Publication Number Publication Date
CN112213718A CN112213718A (en) 2021-01-12
CN112213718B true CN112213718B (en) 2022-10-11

Family

ID=74052275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011022261.5A Active CN112213718B (en) 2020-09-25 2020-09-25 Networking radar node selection and radiation resource joint optimization method under multi-target tracking

Country Status (1)

Country Link
CN (1) CN112213718B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113093172B (en) * 2021-03-15 2023-12-05 南京航空航天大学 Multi-platform radar path and waveform selection combined optimization method for target tracking
CN113253239B (en) * 2021-05-26 2024-06-07 中国人民解放军空军工程大学 Node scheduling and transmitting resource allocation method for centralized MIMO radar network
CN114666219B (en) * 2022-03-15 2023-04-28 南京航空航天大学 Multi-radar network power and bandwidth joint optimization allocation method and system under non-ideal detection
CN118244258A (en) * 2024-03-14 2024-06-25 南京航空航天大学 Networked radar detection threshold and radio frequency resource joint optimization method under multi-target tracking

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110084871A1 (en) * 2009-10-13 2011-04-14 Mcmaster University Cognitive tracking radar
DE102013209708A1 (en) * 2013-05-24 2014-11-27 Robert Bosch Gmbh Method for operating a MIMO radar
CN110412534B (en) * 2019-08-08 2023-04-07 南京航空航天大学 Networking radar multi-target tracking residence time optimization method based on radio frequency stealth
CN110412515B (en) * 2019-08-19 2023-02-03 南京航空航天大学 Multi-target tracking transmission power distribution method for networking radar based on radio frequency stealth
CN111025275B (en) * 2019-11-21 2021-10-08 南京航空航天大学 Multi-base radar radiation parameter multi-target joint optimization method based on radio frequency stealth
CN111323773A (en) * 2020-02-20 2020-06-23 南京航空航天大学 Networking radar power and bandwidth joint optimization distribution method based on radio frequency stealth

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Joint Subcarrier Assignment and Power Allocation Strategy for Integrated Radar and Communications System Based on Power Minimization;Chenguang Shi 等;《IEEE SENSORS JOURNAL》;20191201;第19卷(第23期);11167-11179 *
基于协同噪声干扰的机载雷达射频隐身性能优化算法;时晨光 等;《数据采集与处理》;20170331;第32卷(第2期);321-329 *

Also Published As

Publication number Publication date
CN112213718A (en) 2021-01-12

Similar Documents

Publication Publication Date Title
CN112213718B (en) Networking radar node selection and radiation resource joint optimization method under multi-target tracking
Shi et al. Joint transmit resource management and waveform selection strategy for target tracking in distributed phased array radar network
CN111323773A (en) Networking radar power and bandwidth joint optimization distribution method based on radio frequency stealth
CN112068124B (en) Networking radar residence time and radiation power combined optimization method for low interception
CN106199579B (en) Distributed MIMO radar target tracking precision method for joint optimization of resources
CN106405536B (en) A kind of MIMO radar multiple target tracking method for managing resource
Yan et al. Power allocation algorithm for target tracking in unmodulated continuous wave radar network
CN111025275B (en) Multi-base radar radiation parameter multi-target joint optimization method based on radio frequency stealth
CN108896985B (en) Networking radar multi-target tracking sampling interval control method based on radio frequency stealth
CN111090078B (en) Networking radar residence time optimal control method based on radio frequency stealth
CN110412515B (en) Multi-target tracking transmission power distribution method for networking radar based on radio frequency stealth
CN109283522B (en) Co-location MIMO radar target tracking method based on joint time-space resource management
CN111190176A (en) Adaptive resource management method of co-location MIMO radar networking system
CN108562894B (en) Method for distributing radar beam pointing and transmitting power
CN106021697B (en) A kind of rapid phase-control battle array radar Time-energy resource joint management method
Narykov et al. Algorithm for resource management of multiple phased array radars for target tracking
CN112147600B (en) Multi-base radar transmission parameter optimization method facing radio frequency stealth and target tracking
CN110794395B (en) Networking radar multi-target tracking time resource and signal bandwidth combined optimization method
CN108398678A (en) A kind of quick self-adapted method for tracking target of centralization MIMO radar
CN114706045A (en) Networking radar power time joint optimization method for multi-target tracking under space frequency perception
Junkun et al. Joint power and bandwidth allocation for centeralized target tracking in multiple radar system
CN115561748A (en) Networked radar target search tracking resource allocation method based on radio frequency stealth
CN110488276A (en) The optimal resource allocation method based on demand of isomery radar fence towards multiple target tracking task
CN116577730A (en) Phased array radar combined space-time resource allocation method based on variable false alarm detection
Du et al. Joint power and bandwidth allocation algorithm for radio frequency stealth radar

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