CN115441906A - Cooperative MIMO radar communication integrated system power distribution method based on cooperative game - Google Patents
Cooperative MIMO radar communication integrated system power distribution method based on cooperative game Download PDFInfo
- Publication number
- CN115441906A CN115441906A CN202210873425.8A CN202210873425A CN115441906A CN 115441906 A CN115441906 A CN 115441906A CN 202210873425 A CN202210873425 A CN 202210873425A CN 115441906 A CN115441906 A CN 115441906A
- Authority
- CN
- China
- Prior art keywords
- communication
- radar
- target
- representing
- cooperative
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 127
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000009826 distribution Methods 0.000 title claims abstract description 21
- 238000001514 detection method Methods 0.000 claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims description 32
- 239000013598 vector Substances 0.000 claims description 25
- 230000005540 biological transmission Effects 0.000 claims description 18
- 238000007781 pre-processing Methods 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 5
- 230000000996 additive effect Effects 0.000 claims description 5
- 238000005562 fading Methods 0.000 claims description 5
- 230000000295 complement effect Effects 0.000 claims description 4
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 101100497221 Bacillus thuringiensis subsp. alesti cry1Ae gene Proteins 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 5
- 238000009827 uniform distribution Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000005336 cracking Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0426—Power distribution
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0619—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal using feedback from receiving side
- H04B7/0621—Feedback content
- H04B7/0634—Antenna weights or vector/matrix coefficients
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/021—Estimation of channel covariance
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a cooperative MIMO radar communication integrated system power distribution method based on cooperative game, and belongs to the field of communication. The power distribution method of the cooperative MIMO radar communication integrated system is characterized in that under the Neyman-Pearson criterion, the target detection probability is used as a radar performance evaluation index, and the communication mutual information quantity is used as a communication performance evaluation index. In the method, a radar subsystem and a communication subsystem are regarded as two cooperative game parties, the total power limitation of the system is considered, and an optimization problem based on Nash bargaining solution is established. According to the cooperative game theory, the Nash bargaining solution obtained by utilizing the iterative NBS algorithm meets the properties of pareto optimality, fairness and the like. The iterative NBS algorithm can adjust whether system performance is more prone to radar or communication by negotiating the point of failure.
Description
Technical Field
The invention belongs to the field of signal processing, relates to the power distribution problem of a cooperative MIMO radar communication integrated system, and is suitable for designing the antenna transmitting power of an MIMO radar subsystem and an MIMO communication subsystem in the cooperative integrated system.
Background
In the past research, radar and communication generally work independently, and in recent years, many experts and scholars focus on radar and communication integrated systems. Radar systems and communication systems have many similarities, for example, in hardware, radio frequency parts of the two systems are similar, signals are radiated to the space through an antenna and then processed at a receiving end, and the antenna, a transmitter, a receiver and the like have common possibility. The difference of the two is reflected in signal processing, information of radar comes from a target, and information of communication comes from a transmitter. Therefore, radar and communication systems have a basis for achieving integration.
In recent years, the commercial 5G technology and the research of 6G technology and the development of millimeter wave radar provide fertile soil for the development of radar communication integration. Through extensive research and verification, the existing research divides radar communication integrated systems into the following three categories: the system comprises (1) a radar communication coexistence system, (2) a cooperative radar communication coexistence system and (3) a dual-function radar communication integrated system. The radar communication in the radar communication coexistence system shares resources, but the radar and communication systems are designed separately. The present invention is based on a cooperative coexistence system that utilizes paths that would otherwise be considered interference in the coexistence system to help improve system performance. The radar and the communication system in the dual-function integrated system share the same hardware platform, and the same platform can simultaneously realize radar and communication functions, so that the dual-function integrated system has high integration.
In the field of communications, as early as the 4G era, MIMO (Multiple Input Multiple Out) technology is widely used, and the key technology of 5G is massive MIMO antenna array and beam forming technology. At present, due to the maturity of millimeter wave technology, millimeter wave radars are rapidly developed in the fields of automobile auxiliary driving, security inspection, medical detection and the like. Similar to communication, millimeter wave technology brings higher bandwidth, but transmission loss is high, and combining millimeter wave and massive MIMO technology can complement advantages. The MIMO technology can improve the precision of target parameter estimation and the target detection capability. The MIMO technology applied in the communication field can improve the channel capacity, improve the reliability of the channel and reduce the error rate.
The power allocation problem has been an important issue in the research of various types of systems. The MIMO system that this patent was considered is the split antenna MIMO system, and present MIMO radar divide into common antenna MIMO radar and split antenna MIMO radar. The multi-angle target detection can be realized for the MIMO radar with the split antennas, and the target detection performance is improved. The method considers the power distribution problem under the condition that the total power of the radar and the communication is limited, designs by taking the overall performance of the optimized radar and the communication as a target, and provides an algorithm for realizing the optimal overall performance in the pareto meaning.
Current applications to gaming theory focus on cooperative gaming and non-cooperative gaming. In the field of radar, game theory is often used for researching the mutual influence process between different antennas of the MIMO radar and different base radars, the mutual influence process between the radar and jammers, the cooperation and confrontation between the radar and targets, and the like. In the field of communication, cooperative game is often adopted to study how to reasonably and fairly distribute resources of each system to finally realize the optimal overall performance. The game theory provides a brand-new visual angle and a model, and provides an efficient solution for solving complex problems such as a non-convex multi-target optimization problem and the like.
Therefore, aiming at the problem of optimizing the power distribution of the communication integrated transmitting antenna of the cooperative MIMO radar, the invention provides a cooperative game model, and provides an iterative Nash locking Solution (NBS) algorithm to realize the optimization problem of the overall performance of the radar and the communication under the condition that the total power of the system is limited.
Disclosure of Invention
The invention provides a power distribution method of a cooperative MIMO radar communication integrated system, which combines cooperative game knowledge and takes target detection probability as a radar performance evaluation index and communication mutual information quantity as a communication performance evaluation index under a Neyman-Pearson criterion. In the method, a radar subsystem and a communication subsystem are regarded as two cooperative game parties, the total power limitation of the system is considered, and an optimization problem based on Nash bargaining solution is established. According to the cooperative game theory, the Nash bargaining solution obtained by utilizing the iterative NBS algorithm meets the properties of pareto optimality, fairness and the like. The iterative NBS algorithm can adjust whether system performance is more prone to radar or communication by negotiating the point of failure.
The technical scheme of the invention is that a cooperative MIMO radar communication integrated system power distribution method based on cooperative game comprises the following steps:
step 1: let MIMO radar have N R A receiving antenna, M R One transmitting antenna, MIMO communication has N C A receiving antenna, M C Transmitting antennas, each antenna position being known for both radar and communication; the total transmitting power of the integrated system isThe total transmitting power of the radar isThe transmission power of the mth radar transmission antenna is denoted as E R,m Total transmission power of the communication isThe m' th communication transmitting antenna has the transmitting power of E C,m′ The radar transmission power distribution weight is defined as eta R ,
Step 2: defining the observed value of the MIMO radar end receiving signal as a vector r R The observed value of the MIMO communication terminal receiving signal is a vector r C :
wherein URt and UR Representing a radar-radar channel matrix, comprising a block diagonal matrix of target reflection coefficients, fading coefficients, respectively, U Ct and UC A block diagonal matrix s representing a communication-radar channel matrix and including a target reflection coefficient and a fading coefficient respectively Rt and sR Representing the warp of a targetRadar signal vector of the direct and direct paths, s Ct and sC Representing vectors of communication signals, U, reflected and direct paths through the object Rt and UR Representing a radar-communication channel matrix, comprising a block diagonal matrix of target reflection channel gains, direct path channel gains, respectively, U Ct and UC A block diagonal matrix representing a communication-communication channel matrix containing channel gains for target reflections, and direct path channel gains, respectively,andrepresenting the radar signal vectors received by the receiving end of the communication via the target reflection and direct path,andvector of communication signals, w, representing reflected and direct paths through the target received at the receiving end of the communication R Represents an additive white Gaussian noise vector with a mean of zero and a covariance matrix of Q R ,w C Representing an additive white Gaussian noise vector, obeying a zero mean and having a covariance matrix of Q C ;
And step 3: the assumption defining the existence of the target isThe assumption that the target does not exist isEstablishing a hypothesis test problem to obtain a log-likelihood ratio, and writing the detection statistics as:
and 4, step 4: in general, NP criteria specify false alarm profilesRate P FA Below a certain value alpha f In the case of (2), the maximum detection probability is calculated, giving a false alarm probability α f And introducing the detection statistic to obtain the radar target detection probability P D :
The function Q (-) represents a complementary distribution function,σ is the detection statistic T R Standard deviation of (2), operatorExpressing the mathematical expectation of solving random variables; and because ofOc indicates that the two are in direct proportion, indicated byTo represent the detection probability P D In the case of a change in the state of the art,
wherein the parameter alpha m,m′ ,β m,m′ ,γ m,m′ Relative to target reflection coefficient, radar and communication emission signals, time delay;
and 5: according to the mutual information quantity definition, the communication mutual information quantity MI is calculated as:
wherein Comprising radar hairThe method comprises the steps of transmitting power, radar transmitting signals, reflection coefficients, target time delay and target estimation error items, assuming that communication receiving signals are Gaussian signals with the mean value of zero, at different observation moments, the receiving signals of different communication receivers are independent and uncorrelated, assuming that the time delay estimation error of a communication receiving end is in a reasonable range, and meeting the requirement that the time delay estimation error of the communication receiving end is in the reasonable rangeIs much smaller than Q C Each element of (1) willApproximated as a diagonal matrix, using the diagonal matrix algorithm, the mutual information MI can be written as:
can see chi m,m′,n′,k Communication signal vector with communication receiving endThe direct path channel gain is related to the target reflection channel gain; channel gain via target reflection and direct path gain via pre-processing means, psi, due to radar communication cooperation n′,m,k The channel gain reflected by the target in the system is obtained in a preprocessing mode and is influenced by time delay estimation errors;
step 6: based on cooperative game knowledge in game theory, there is a unique and fair Nash bargained solutionCan be obtained by maximizing the nash product:
wherein the point of failure is negotiated asAnd (5) solving by adopting an iterative NBS algorithm.
The iterative NBS algorithm provided by the invention is suitable for mutual cooperation between the radar and the communication subsystem, communication signals can be accurately decoded at a radar receiving end, and the target reflection coefficient and channel gain of the radar and the communication can be obtained through preprocessing. Under the condition that the total power of the system is limited, the radar and communication overall performance is optimal, and whether the overall performance emphasizes radar or communication can be realized by adjusting negotiation breaking points. The iterative NBS algorithm emphasizes the overall performance of the system, can realize pareto optimality, and has higher precision compared with a genetic algorithm. Therefore, the method provided by the invention is an efficient and effective method for solving the problem of power distribution of the cooperative MIMO radar communication integrated system.
Drawings
Fig. 1 shows a distribution diagram of the positions of the transmitting and receiving antennas of the cooperative integrated system.
Fig. 2 is a graph of iteration number versus transmit power using an iterative NBS algorithm.
FIG. 3 is a graph showing the relationship between the target detection probability and the mutual information amount under the NSGA-II algorithm and the iterative NBS algorithm.
Fig. 4 is a graph of the variation of the target detection probability and the communication mutual information amount with the signal-to-noise ratio SCNR under the iterative NBS algorithm, the trisection search nash equalization algorithm, the uniform distribution method and the random distribution method.
FIG. 5 is a graph showing the overall performance of the system varying with SCNR under the iterative NBS algorithm, the trisection search Nash equilibrium algorithm, the uniform distribution method and the random distribution method.
Detailed Description
Defining the signal transmitted by the m-th radar transmitting antenna asBy s R,m (T) denotes a radar transmission signal, T s Denotes the time sampling interval, with K (K =1, 2.., K) denoting the different sampling samples, and s for the communication end C,m′ (t) represents a communication transmission signal, and a transmission signal of a communication transmission antenna is represented asConsidering that both radar and communication signals have been normalized,according to a model of the radar received signal, at kT s Time, nth (N =1,., N) R ) The received signal of each radar receiving antenna is represented as
wherein ζRt,nm and ζCt,nm′ Representing the reflection coefficient of the target for different paths. Similarly, ζ R,nm and ζC,nm′ Representing the radar fading coefficients, tau, under different paths Rt,nm and τCt,nm′ Representing the time delay, τ, of reflections off the target under different paths R,nm and τC,nm′ Representing the time delay of the direct path under the different paths. w is a R,n [k]Gaussian white noise signals, noise of different paths is zero-mean and independently and uniformly distributed.
The antenna position, the transmitting signal, the transmitting power and other information of the cooperative integrated system are mutually shared between the radar and the communication. And the communication signal can be decoded and reconstructed more accurately by using the communication information at the radar receiving end. The target reflection coefficient of the radar can be known by means of preprocessing.
The observation vector of the nth receiving antenna at the MIMO radar end at different time can represent r R,n =(r R,n [1],...,r R,n [K]) T The observation vector of the MIMO radar end can be written as
wherein
Diag {. Is } representation BlockDiagonal matrix, w R Is an additive white Gaussian noise vector with a mean of zero and a covariance matrix of Q R ,
At the radar receiving end, if the radar task is to detect the presence of a target, then the hypothesis testing problem is expressed as
Accordingly, the log-likelihood ratio is expressed as
The sum r is found in formula (7) R Related terms, resulting in detection statistics expressed as
In a general radar detection problem, judgment can be performed only through a received signal, the prior probability and the cost cannot be known in advance, and a Neyman-Pearson criterion is adopted as a main judgment criterion at this time. In general, the Neyman-Pearson criterion specifies that the false alarm probability is below a certain value α f And calculating the maximum detection probability. Because the false alarm probability and the detection probability are usually increased at the same time, and the actual requirement is that the lower the false alarm probability is, the higher the detection probability is, the maximum detection probability under the premise can be calculated on the premise of giving the maximum acceptable false alarm probability. The false alarm probability is expressed as
wherein Indicating that there is actually a target present. The detection probability is expressed as
Bringing the detection statistic into (9)
Where beta represents the detection threshold and where,
β=σQ -1 (α f )+μ 0 (12)
thus, the target detection probability of a radar is expressed as
wherein
It is known that the standard Gaussian complement distribution function Q (-) is a monotonically decreasing function when the argument is greater than zero, and thus P D And (Q) -1 (α f )+(μ 0 -μ 1 ) /σ) is inversely proportional to the region greater than zero. Suppose Q -1 (α f ) A value of greater than (mu) 0 -μ 1 ) σ (in actual simulation signal selection, this assumption generally holds).
wherein
At kT s Time, N '(N' =1, 2., N) C ) The received signal of a communication can be expressed as
wherein ζCt,n′m′ and ζRt,n′m Representing the channel gain reflected by the target for different channels. ζ represents a unit R,n′m and ζC,n′m′ Representing the direct path channel gain at different channels,andrepresenting the time delay of the reflection off the target under different channels,andrepresenting the time delay of the direct path under different channels. w is a C,n′ [k]The Gaussian white noise signals received by the communication receiving end are zero-mean and independently distributed.
In the cooperative integrated system, the antenna positions, the transmission signals and the transmission power of the radar and the communication are known through cooperation. The communication receiving end can estimate the target position by using the target echo signal of the radar, thereby eliminating the signal interference of the radar to the communication direct path and the target reflected path.
The observation vectors of the n' th receiving antenna at different time at the MIMO communication end are represented as r C,n′ =(r C,n′ [1],...,r C,n′ [K]) T The observation vector of the MIMO communication terminal can be written as
wherein UCt 、U C 、U Rt 、U R Definition of (2) and U Ct 、U C 、U Rt 、U R In a similar manner to the above-described embodiments,definition of (1) and s Ct 、s C 、s Rt 、s R Similarly. w is a C Subject to mean of zero and covariance matrix of Q C The white gaussian noise signal of (a) is,
at MIMO communication receiving end, the direct wave interference of the communication receiving end radar can be eliminated by cooperatively sharing the radar positionAnd further eliminating the interference of the radar target echo according to the estimation result of the target position. Assuming that the target position to be estimated is θ = (x, y), the communication transmission signal is a gaussian signal, and the communication reception signal r C Can be expressed as
the patent adopts ML estimation method, ML estimation of the target position,
wherein Is an estimate of the target location. The estimated value of the target position influences the time delay obtained by radar and communication through a target reflection path, and the estimation error n of the target parameter estimation is assumed Ct,n′m′ and nRt,n′m Obeying a Gaussian distribution, the estimated delay is
The parameter estimation of the target at the communication end can be used for eliminating the target echo interference from radar transmission signals at the communication end, and communication signals reflected by the target can also be utilized. Echo signalTime delay inSubstitution by estimated time delayTo obtainIs estimated value ofTo simplify the analysis, the estimation error of the reflection of the communication signal by the target is assumedCan be omitted. Due to estimation errors, the communication reception signal can be written as
wherein
wherein Denotes s R,m (t) calculating a partial derivative of t. The amount of mutual communication information can be written as
Since the communication reception signal is a gaussian signal with a mean value of zero, the reception signals of different communication receivers are independent and uncorrelated with each other at different times. For different N 1 and N2 ,k 1 and k2 ,Further assume that the time delay estimation error n arriving at the receiving end of the communication via the target reflection Ct,n′m′ and nRt,n′m Within a reasonable range, so thatAre much smaller than Q C Each element of (a). Thus, can be connectedApproximated as a diagonal matrix. Further using the diagonal matrix algorithm, the amount of mutual information can be expressed as
wherein
According to the cooperative game theory, a cooperative game model is established, game participants comprise two parties of radar and communication,the strategy set of the radar isThe strategy set of communication isSymbol x represents a cartesian product. The utility function u of the radar is expressed by the simplified radar detection probability (see equation (15)) R (E R ,E C ) The utility function of the communication is represented as u by the simplified mutual traffic information quantity (see equation (25)) C (E R ,E C )。
The power allocation problem is how to allocate the transmit power of the radar and communication systems to achieve fairness under total power constraints. Defining negotiation failure point asNash product isThe optimization problem can be described as
The steps for the iterative NBS algorithm are as follows:
initialization: initial valueLagrange multiplierStep length s t Step size k of iteration 1 The number of iterations n =0, the convergence factor epsilon > 0.
And 3, step 3: if obtained in steps 1 and 2Beyond the feasible domain, i.e.Updating with orthogonal projection operatorsThe orthogonal projection operator P is
P=I n -A T (AA T ) -1 A
And (4) finishing the algorithm: and obtaining the optimal power distribution.
Three simulation examples and comparison curves are given for the cooperative game based power distribution method.
The simulation parameters are set as follows: the system antennas are assumed to be in a two-dimensional cartesian coordinate system and are each 70 kilometers from the origin of coordinates. Consider that MIMO radar has M R =2 transmitting antennas located at (70, 0) km and (-70, 0) km, with N R =3 receiving antennas, positions are (66, 24) km, (-54, 45) km and (-12, -69) km; consider that MIMO communication has M C =2 transmitting antennas, located at (0, 70) km and (0, -70) km, having N C =3 receive antennas, positions (-24,66) km, (-45,64) km and (69, -12) km; the coordinates of the target position to be estimated are (50, 30) meters, and the antenna position and the target position are shown in fig. 1. The MIMO radar transmits signals which are single Gaussian pulse signals,f Δ =125hz, t =0.01s. The transmission signal of MIMO communication employs an orthogonal frequency division multiplexing signal,T′=0.01s,Δf=125Hz,N f and (6). Covariance matrixSignal to clutter noise ratio ofTotal power of transmitting antenna of cooperative integrated systemKilowatts, SCNR = -6dB.
In simulation 1, the convergence condition is set to be 10 as the difference between the powers of two adjacent iterations -5 The number of the tiles is such that,to negotiate the break point. As shown in fig. 2, it can be concluded that the transmit power tends to stabilize after about 25 iterations.
In simulation 2, the optimal individual coefficient pareto frame using the NSGA-II algorithm is 0.8, the population size popup size is 100, the maximum genetic Generations are 200, the stop algebra stalgenlimit is 200, the fitness function deviation TolFun is 1e to 100, and the resulting pareto boundary is as shown in fig. 3 (a). Modifying the pareto fraction to 0.7, results in fig. 3 (b), which shows that the theoretical pareto optimum has been achieved for agreement at the negotiation breaking point (0, 0), and that the results obtained using the NSGA-II algorithm are still inferior to the results obtained using the NBS algorithm at the negotiation breaking point (1, 0).
In the simulation 3, under different SCNRs, changes of radar detection probability and communication mutual information amount are analyzed by adopting an NE algorithm based on three-section searching, an iterative NBS algorithm (negotiation cracking points (1, 0)), uniform distribution and random distribution under four power distribution methods. In fig. 4, when the SCNR is less than-2 dB, the random allocation method can obtain better detection performance than the uniform allocation and iterative NBS algorithm, but the communication performance is poor. The uniform distribution and iterative NBS algorithm can obtain similar detection probability, but the iterative NBS algorithm obtains much higher mutual information amount than other methods. Comparing the three-division search NE algorithm with the iterative NBS algorithm, the radar detection performance obtained by the three-division search NE algorithm is far better than that of the iterative NBS algorithm when the SCNR is less than 0dB, but the iterative NBS algorithm is better in a communication mutual information quantity graph. With the increasing of the SCNR, the interference signal strength of the radar is weakened, and at the moment, good target detection performance can be achieved without too much radar power. The increase of SCNR improves the cooperative gain brought by the iterative NBS algorithm, and the system distributes more power to the communication end, thereby ensuring that the performance of the communication end is optimal on the premise of optimal target detection probability. Figure 5 shows the effect of different allocation algorithms on the overall performance of the system. Fig. 5 shows that the overall performance of the system obtained by the iterative NBS algorithm is better than that of the rest of the algorithms under different SCNRs. The iterative NBS algorithm is mainly based on overall performance, and under the simulation scene of the invention, the algorithm sacrifices partial radar detection probability to obtain the optimal overall performance of the system. The algorithm proposed by the patent can provide a reference for a system designer, and an iterative NBS algorithm is adopted when the overall performance is considered.
Claims (1)
1. A cooperative MIMO radar communication integrated system power distribution method based on cooperative game comprises the following steps:
step 1: let MIMO radar have N R A receiving antenna, M R A plurality of transmitting antennas, MIMO communication having N C A receiving antenna, M C A plurality of transmit antennas, each antenna location being known for radar and communications; the total transmitting power of the integrated system isThe total transmitting power of the radar isThe transmission power of the mth radar transmission antenna is denoted as E R,m Total transmission power of communication ofTransmission power of m' th communication transmitting antennaA rate of E C,m′ The radar transmission power distribution weight is defined as eta R ,
Step 2: defining the observed value of the MIMO radar end receiving signal as a vector r R The observed value of the received signal of the MIMO communication end is a vector r C :
wherein URt and UR Representing a radar-radar channel matrix, comprising a block diagonal matrix of target reflection coefficients, fading coefficients, respectively, U Ct and UC A block diagonal matrix s representing a communication-radar channel matrix and including a target reflection coefficient and a fading coefficient respectively Rt and sR Representing the radar signal vector, s, reflected by the target and directed path Ct and sC Representing vectors of communication signals, U, reflected and direct paths through the object Rt and UR Representing a radar-communication channel matrix comprising a block diagonal matrix of target reflection channel gains, direct path channel gains, respectively, U Ct and UC A block diagonal matrix representing a communication-communication channel matrix containing channel gains for target reflections, direct path channel gains, respectively,andrepresenting a radar signal vector received by a communication receiving end via a target reflection and a direct path,andvector of communication signals, w, representing reflected and direct paths through a target received at a receiving end of the communication R Represents an additive white Gaussian noise vector with a mean of zero and a covariance matrix of Q R ,w C Representing an additive white Gaussian noise vector obeying a zero mean and a covariance matrix of Q C ;
And 3, step 3: the assumption defining the existence of the target isThe assumption that the target does not exist isEstablishing a hypothesis test problem to obtain a log-likelihood ratio, and writing the detection statistics as:
and 4, step 4: in general, the NP criterion specifies a false alarm probability P FA Below a certain value alpha f In the case of (2), the maximum detection probability is calculated, giving a false alarm probability α f And substituting detection statistic to obtain radar target detection probability P D :
The function Q (-) represents the complementary distribution function,σ is the detection statistic T R Standard deviation of (1), operatorExpressing the mathematical expectation of solving random variables; and becauseA value of-To represent the detection probability P D In the case of a change in the state of (c),
wherein the parameter alpha m,m′ ,β m,m′ ,γ m,m′ Related to target reflection coefficient, radar and communication emission signals, time delay;
and 5: according to the mutual information quantity definition, the communication mutual information quantity MI is calculated as:
wherein The method comprises radar transmitting power, radar transmitting signals, reflection coefficients, target time delay and target estimation error items, and assumes that communication receiving signals are Gaussian signals with the mean value of zero, so that the receiving signals of different communication receivers are independent and uncorrelated at different observation moments, and the time delay estimation error of a communication receiving end is assumed to be in a reasonable range, thereby meeting the requirement of the communication receiving endIs much smaller than Q C Will be each element ofApproximately as a diagonal momentArray, using the diagonal matrix algorithm, the mutual information MI can be written as:
can see chi m,m′,n′,k Communication signal vector with communication receiving endThe direct path channel gain is related to the target reflection channel gain; channel gain via target reflection and direct path gain via pre-processing means, psi, due to radar communication cooperation n′,m,k The channel gain reflected by the target is obtained in a preprocessing mode and is influenced by a time delay estimation error;
step 6: based on cooperative game knowledge in game theory, there is a unique and fair Nash bargained solutionCan be obtained by maximizing the nash product:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210873425.8A CN115441906B (en) | 2022-07-21 | 2022-07-21 | Cooperative game-based cooperative MIMO radar communication integrated system power distribution method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210873425.8A CN115441906B (en) | 2022-07-21 | 2022-07-21 | Cooperative game-based cooperative MIMO radar communication integrated system power distribution method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115441906A true CN115441906A (en) | 2022-12-06 |
CN115441906B CN115441906B (en) | 2023-08-15 |
Family
ID=84241294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210873425.8A Active CN115441906B (en) | 2022-07-21 | 2022-07-21 | Cooperative game-based cooperative MIMO radar communication integrated system power distribution method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115441906B (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106685441A (en) * | 2016-12-09 | 2017-05-17 | 南京航空航天大学 | Radio frequency stealth radar networking power distribution method based on cooperative game |
CN108717180A (en) * | 2018-07-05 | 2018-10-30 | 南京航空航天大学 | A kind of radar network power distribution method based on Stackelberg game |
CN108834209A (en) * | 2018-05-22 | 2018-11-16 | 南京航空航天大学 | Under frequency spectrum coexisted environment based on receive it is assorted negotiate a price solution radar network power distribution method |
CN109507641A (en) * | 2018-11-07 | 2019-03-22 | 电子科技大学 | A kind of MIMO radar and MIMO communication system Calculation Methods for Performance of cooperation |
CN110441740A (en) * | 2019-08-20 | 2019-11-12 | 南京航空航天大学 | Distributed MIMO radar robust power distribution method based on layering game |
US20200169458A1 (en) * | 2018-11-26 | 2020-05-28 | Eagle Technology, Llc | Radio frequency (rf) communication system providing enhanced rf equipment configuration updates for mobile vehicles based upon reward matrices and related methods |
CN112272064A (en) * | 2020-09-29 | 2021-01-26 | 电子科技大学 | Detection probability and mutual information calculation method of cooperative MIMO radar |
WO2022001148A1 (en) * | 2020-06-28 | 2022-01-06 | 南京航空航天大学 | Non-cooperative game-based radar communication integrated networking system power control method |
US20220057483A1 (en) * | 2019-07-25 | 2022-02-24 | Beijing University Of Posts And Telecommunications | Radar and communication integrated cooperative detection method and apparatus based on beam power distribution |
-
2022
- 2022-07-21 CN CN202210873425.8A patent/CN115441906B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106685441A (en) * | 2016-12-09 | 2017-05-17 | 南京航空航天大学 | Radio frequency stealth radar networking power distribution method based on cooperative game |
CN108834209A (en) * | 2018-05-22 | 2018-11-16 | 南京航空航天大学 | Under frequency spectrum coexisted environment based on receive it is assorted negotiate a price solution radar network power distribution method |
CN108717180A (en) * | 2018-07-05 | 2018-10-30 | 南京航空航天大学 | A kind of radar network power distribution method based on Stackelberg game |
CN109507641A (en) * | 2018-11-07 | 2019-03-22 | 电子科技大学 | A kind of MIMO radar and MIMO communication system Calculation Methods for Performance of cooperation |
US20200169458A1 (en) * | 2018-11-26 | 2020-05-28 | Eagle Technology, Llc | Radio frequency (rf) communication system providing enhanced rf equipment configuration updates for mobile vehicles based upon reward matrices and related methods |
US20220057483A1 (en) * | 2019-07-25 | 2022-02-24 | Beijing University Of Posts And Telecommunications | Radar and communication integrated cooperative detection method and apparatus based on beam power distribution |
CN110441740A (en) * | 2019-08-20 | 2019-11-12 | 南京航空航天大学 | Distributed MIMO radar robust power distribution method based on layering game |
WO2022001148A1 (en) * | 2020-06-28 | 2022-01-06 | 南京航空航天大学 | Non-cooperative game-based radar communication integrated networking system power control method |
CN112272064A (en) * | 2020-09-29 | 2021-01-26 | 电子科技大学 | Detection probability and mutual information calculation method of cooperative MIMO radar |
Non-Patent Citations (4)
Title |
---|
LIMING WANG ET AL.: "Transmitter selection and receiver placement for target parameter estimation in cooperative radar-communications system", 《IET SIGNAL PROCESSING》 * |
T. HARIKALA ET AL.: "Power Efficient Technique for MIMO radar using Co-operative and Non co-operative game theory in Wireless Applications", 《INTERNATIONAL JOURNAL OF RECENT TECHNOLOGY AND ENGINEERING(IJRTE)》, vol. 7, no. 5 * |
王珍等: "合作式MIMO雷达通信一体化***的检测与互信息增益", 《信号处理》, vol. 36, no. 10 * |
靳标等: "基于合作博弈的组网雷达分布式功率分配方法", 《航空学报》, vol. 43, no. 1 * |
Also Published As
Publication number | Publication date |
---|---|
CN115441906B (en) | 2023-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8380150B2 (en) | Received electric field intensity estimating device and received electric field intensity estimating program | |
US8553797B2 (en) | Channel information prediction system and channel information prediction method | |
CN111751800B (en) | Frequency control array radar angle-distance parameter decoupling method | |
CN109738856B (en) | Noise suppression Capon active target DOA estimation method based on time reversal | |
CN108469601B (en) | Transceiving joint robust optimization method for improving worst detection performance of MIMO-STAP | |
CN104794264B (en) | A kind of radar communication waveform design method based on sparse frequency | |
CN113466848B (en) | Co-location MIMO radar multi-target tracking resource optimization allocation method for angular flicker noise scene | |
CN111446999A (en) | Position-assisted beam alignment method and system based on multi-arm forced theft | |
CN112147600B (en) | Multi-base radar transmission parameter optimization method facing radio frequency stealth and target tracking | |
CN112068124A (en) | Networking radar residence time and radiation power combined optimization method for low interception | |
CN104678362B (en) | MIMO sky-wave OTH radar waveform optimization method | |
Kim et al. | RIS-aided radar sensing and object detection with single and double bounce multipath | |
WO2018006355A1 (en) | Transmission weight value selection method and base station | |
CN114662271A (en) | Wave beam-power-bandwidth joint allocation method in centralized MIMO radar multi-target tracking | |
CN112272064B (en) | Detection probability and mutual information calculation method of cooperative MIMO radar | |
Peng et al. | Semi-Passive Intelligent Reflecting Surface Enabled Sensing Systems | |
CN110161471B (en) | Method for calculating sampling rate and quantization bit of cloud MIMO radar | |
CN110146854B (en) | Robust anti-interference method for FDA-MIMO radar | |
CN115441906A (en) | Cooperative MIMO radar communication integrated system power distribution method based on cooperative game | |
CN109669167B (en) | Airborne radar emission waveform selection method based on radio frequency stealth | |
CN109669165B (en) | Multi-base radar transmitter and waveform joint selection optimization method based on radio frequency stealth | |
CN110488276B (en) | Heterogeneous radar network optimal resource on-demand distribution method oriented to multi-target tracking task | |
CN116743222A (en) | Multiplexing beam forming method, communication perception calculation integrated system and related device | |
CN116953688A (en) | Array element selection and power distribution combined optimization method under tracking and detection background | |
CN113126076B (en) | Cooperative distance and angle parameter estimation method and device based on dual-mode 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 |