CN109870676A - Radar-communication integration system performance boundary measurement method based on location estimation rate - Google Patents
Radar-communication integration system performance boundary measurement method based on location estimation rate Download PDFInfo
- Publication number
- CN109870676A CN109870676A CN201910187415.7A CN201910187415A CN109870676A CN 109870676 A CN109870676 A CN 109870676A CN 201910187415 A CN201910187415 A CN 201910187415A CN 109870676 A CN109870676 A CN 109870676A
- Authority
- CN
- China
- Prior art keywords
- radar
- communication
- bandwidth
- rate
- signal
- 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 199
- 230000010354 integration Effects 0.000 title claims abstract description 9
- 238000000691 measurement method Methods 0.000 title abstract 2
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- 239000002131 composite material Substances 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 34
- 230000014509 gene expression Effects 0.000 claims description 21
- 230000005540 biological transmission Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 abstract description 6
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 9
- 230000006872 improvement Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 240000005523 Peganum harmala Species 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention discloses a kind of radar-communication integration system performance boundary measurement method based on location estimation rate, applied to Radar Technology field, for the prior art when measuring radar-communication integration system performance boundary, only target range information, there is no the problem of angle information, the present invention has comprehensively considered target range information and angle information, the radar signal that receiver is received first, signal of communication and its composite signal are modeled, further according to the mode of radar and communications frequency spectrum share in radar-communication integration system, the case where frequency spectrum share is divided into: individually subband is wide for tradition, it is wide wide with radar individual subbands to communicate individual subbands, then location estimation rate is respectively adopted, communication rate is to radar function, communication functionality can be carried out characterization, it is finally logical to radar under different frequency bands sharing mode using location estimation rate and communication rate Believe that the performance limit of system integration system carries out theory deduction, realizes the measurement to the performance of radar-communication integration system.
Description
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a performance limit measuring technology of a radar communication integrated system.
Background
Modern radar devices have a variety of functions, such as target detection, device deployment, and sharing of acquired information with other devices. How to combine radar equipment and communication equipment to form a complete system and solve the problem of reasonable utilization of resources is an increasing research subject in recent years. The radar equipment and the communication equipment are combined to form an integrated system, and how to solve the problem of reasonable utilization of resources becomes one of the research hotspots of experts at home and abroad.
The sharing of radar and communication frequency spectrum is one of key researches for solving the shortage of frequency spectrum resources and improving the utilization rate of frequency spectrum. The method is a research subject with practical significance, but the mutual interference problem caused by the sharing of the radar and the communication frequency spectrum is still a difficult problem of the research, and the development of a characterization method of the performance of the radar communication integrated system and the theoretical research of the performance limit of the radar communication integrated system are still difficult problems of the research. In the literature "Inner Bounds on Performance of Radar and communications Co-Existence, IEEE Transactions on Signal Processing, vol.64, No.2, pp.464-474,2015", the concept of Radar estimated information rate is proposed, and the Performance limit achievable by the Radar communication integrated system is derived by estimating the information rate by the Radar system and the communication rate of the communication system, but only the Time of Arrival (TOA) and not the beam Arrival angle (DOA), i.e. only the target distance information and no angle information, are considered. From the published article at present, the measurement of performance limit of radar communication integration system by using positioning estimation rate has not been studied.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method for measuring the performance limit of a radar communication integrated system based on a positioning estimation rate, which considers the joint information of the arrival time and the arrival angle of a wave beam and realizes the performance limit measurement of the radar communication integrated system in different frequency band sharing modes.
The method for measuring the performance limit of the radar communication integrated system based on the positioning estimation rate comprises the following steps:
s1, receiving the signals by the receiving window at the same time: radar signals, communication signals, composite signals;
s2, dividing a spectrum sharing mode: a traditional independent sub-bandwidth mode, a communication independent sub-bandwidth mode and a radar independent sub-bandwidth mode;
under the traditional independent sub-bandwidth mode, the total bandwidth is divided into two parts, one part is used for communication, and the other part is used for radar; the two parts work in respective sub-frequency bands without any interference;
under the communication independent sub-bandwidth mode, the total bandwidth is divided into two parts, one part is only used for communication and is called communication independent sub-bandwidth, and the other part is simultaneously used for radar and communication and is called mixed sub-bandwidth;
under the radar independent sub-bandwidth mode, the total bandwidth is divided into two parts, one part is only used for radar and is called radar independent sub-bandwidth, and the other part is simultaneously used for radar and communication and is called hybrid sub-bandwidth;
s3, obtaining a positioning estimation rate boundary expression representing radar function performance by considering the beam arrival time and the beam arrival angle joint information;
s4, representing the performance of the communication function by adopting the communication rate;
s5, calculating the communication rate and the estimated location rate boundary under different spectrum sharing modes in step S2.
Further, in step S3, the process of obtaining the expression of the positioning estimation value characterizing the performance of the radar function is:
a1, obtaining probability density functions of the radar signals received in the step S1 with respect to tau and theta; where τ represents the beam arrival time and θ represents the beam arrival angle;
a2, obtaining a corresponding snow information matrix according to the probability density function of the step A1;
a3, obtaining a Cramer-Rao boundary of beam arrival time and a Cramer-Rao boundary of beam arrival angle with the Fisher-snow information matrix of the step A2;
a4, Cramer-Role bound according to beam arrival time, and Cramer-Role bound according to beam arrival angle; and obtaining an upper bound expression of the positioning estimation rate. The step a4 specifically includes:
a41, the positioning estimation rate is:
wherein h isrsRepresenting the entropy of the received signal, hestRepresenting the estimated entropy, TpriRepresenting a pulse repetition interval;
a42, obtaining the minimum mean square error of the positioning estimator according to the Cramer-Rao bound of the arrival time of the wave beam and the Cramer-Rao bound of the arrival angle of the wave beam;
a43, respectively obtaining h of the received radar signals according to the minimum mean square error of the positioning estimatorrsAnd hest：
Wherein, JminRepresenting the minimum mean square error, n, of the positioning estimatorproRepresenting the remaining noise except the estimated noise;
a44, h obtained in the step A43rsAnd hestSubstituting into the expression of the step A41, obtaining an upper bound expression of the positioning estimation rate as follows:
where T represents the pulse width, γ represents the signal-to-noise ratio,
further, in step S4, the communication rate expression is:
wherein, BcDenotes the communication independent sub-bandwidth, h denotes the communication signal propagation gain, k is the boltzmann constant, T0Indicating the absolute temperature, p, of the communication systemcRepresenting the communication power in the total bandwidth.
Further, the communication rate and the positioning estimation rate boundary in the conventional individual sub-bandwidth manner in step S5 are respectively:
wherein α is a bandwidth adjustment factor, 0 is equal to or more than α is equal to or more than 1, B represents a bandwidth, n is equal to or more than 1proRepresenting the remaining noise, except the estimated noise, gamma representing the signal-to-noise ratio, BrRepresenting the radar independent sub-bandwidth, delta representing the duty cycle, T representing the pulse width,n denotes the number of antenna elements, d denotes the antenna element spacing, c denotes the speed of light, λ is the plane wave wavelength, prDenotes the radar transmission power, and g denotes the propagation gain of the radar signal.
Further, in step S5, the communication rate of the communication independent sub-bandwidth in the communication independent sub-bandwidth mode is:
wherein p isc,oAnd the communication power distributed in the communication independent sub-bandwidth mode.
Further, the calculation process of the communication rate and the position estimation rate boundary of the hybrid sub-bandwidth in the communication independent sub-bandwidth manner in step S5 is as follows:
b1, processing the received composite signal in the mixed frequency band by using SIC;
b2, the joint receiver firstly decodes the received communication signal;
b3, when the communication signal is successfully decoded, subtracting the communication signal from the composite signal to obtain a radar signal without communication signal interference;
b4, using the radar signal in the step B3 as interference, and obtaining the communication rate in the mixed sub-bandwidth as:
wherein p isc,mCommunication power allocated in the mixed sub-bandwidth in the communication independent sub-bandwidth mode;
b5, since the communication signal has been removed, according to the positioning estimation rate boundary expression of step S3, the positioning estimation rate boundary in the mixed sub-bandwidth is obtained as:
further, in step S5, the positioning estimation rate of the radar independent sub-bandwidth in the radar independent sub-bandwidth mode is:
wherein, γr,oFor radar in the radar independent sub-bandwidth modeCommunication signal-to-noise ratio, p, of independent sub-bandwidthsr,oAnd allocating the radar power in the radar independent sub-bandwidth mode.
Further, the boundary calculation process of the communication rate and the positioning estimation rate of the hybrid sub-bandwidth in the radar-independent sub-bandwidth manner in step S5 is as follows:
c1, processing the received composite signal in the mixed frequency band by using SIC;
c2, the joint receiver firstly decodes the received radar signal;
c3, when the radar signal is decoded successfully, subtracting the radar signal from the composite signal to obtain a communication signal without radar signal interference;
c4, using the communication signal in the step C3 as interference, obtaining the positioning estimation rate in the mixed sub-bandwidth:
wherein, γr,mCommunication signal-to-noise ratio, p, for mixed sub-bandwidths in radar independent sub-bandwidth moder,mAllocating radar power in the mixed sub-bandwidth in a radar independent sub-bandwidth mode;
c5, since the radar signal has been removed, according to the communication rate expression of step S4, the communication rate in the mixed sub-bandwidth is found to be:
the invention has the beneficial effects that: the method comprises the steps of firstly modeling radar signals, communication signals and composite signals thereof received by a receiver, dividing frequency spectrum sharing conditions into three types according to a radar and communication frequency spectrum sharing mode in the radar communication integrated system, namely traditional single sub-bandwidth, communication independent sub-bandwidth and radar independent sub-bandwidth, then formulating a performance characterization method of the radar communication integrated system, respectively characterizing radar functions and communication function performances by adopting positioning estimation rate and communication rate, and finally carrying out theoretical derivation on performance limits of the radar communication integrated system in different frequency spectrum sharing modes by utilizing the positioning estimation rate and the communication rate to realize the measurement of the performance of the radar communication integrated system; the method has the advantages that the combined information of the time of arrival (TOA) of the wave beam and the angle of arrival (DOA) of the wave beam is considered, and the performance limit measurement of the radar communication integrated system under different frequency band sharing modes is realized; the method is suitable for the fields of civil military and the like.
Drawings
FIG. 1 is a block diagram showing the general structure of the method of the present invention.
Fig. 2 is a block diagram of a conventional ISB spectrum sharing method according to the present invention.
Fig. 3 is a block diagram of a structure of the present invention adopting a CIB spectrum sharing method.
Fig. 4 is a block diagram of the RIB spectrum sharing method according to the present invention.
Fig. 5 is a simulation result of performance limits of the radar communication integrated system in three spectrum sharing modes of ISB, CIB, and RIB in the embodiment of the present invention.
Fig. 6 shows a simulation result of the improvement of the communication rate performance of the CIB and RIB spectrum sharing modes compared with the ISB spectrum sharing mode in the embodiment of the present invention.
Fig. 7 shows a simulation result of the improved positioning estimation rate performance of the CIB and RIB spectrum sharing modes compared with the ISB spectrum sharing mode in the embodiment of the present invention.
Detailed Description
In order to facilitate the understanding of the technical contents of the present invention by those skilled in the art, the present invention will be further explained with reference to the accompanying drawings.
The invention mainly adopts a simulation experiment method for verification, and all the steps and conclusions are verified to be correct on Matlab 2014. The present invention will be described in further detail with reference to specific embodiments.
Step 1: modeling a received signal:
receiving window simultaneously receiving radar signal zr(t) and communication signal zc(t), the composite signal is z (t) is:
z(t)＝zr(t)+zc(t) (1)
defining the eigenvector a (θ) of the antenna:
a(θ)＝[1,e-j2πdcosθ/λ,···,e-j2π(N-1)dcosθ/λ]T(2)
where d is the antenna element spacing, λ is the plane wave wavelength, θ is DOA, and the received radar signal zr(t) is given by:
wherein p isrRepresenting radar transmission power, frRepresenting the carrier frequency of the radar signal, g representing the propagation gain of the radar signal, τrRepresenting the propagation delay time, phi, of the radar signalrRepresents the carrier phase caused by the propagation delay time and satisfies phir＝-2πfrτr。nr(t) is a mean of 0 and a variance of σ2＝κTsysWhite Gaussian Noise (AWGN) of B, where κ is Boltzmann constant, TsysRepresenting the radar system noise temperature and B the bandwidth.
Mathematical expression of a received communication signal:
wherein f iscRepresenting the carrier frequency, s, of the communication signalc(t) represents an equivalent communication baseband signal, h represents a communication signal propagation gain, τcRepresenting propagation delay time, phi, of the communication signalrRepresents the carrier phase caused by propagation delay and satisfies phic＝-2πfcτc。
Step 2: dividing a spectrum sharing mode:
according to the method of sharing the radar and the communication frequency band in the radar communication integrated system, the frequency spectrum sharing is divided into three types, namely, a traditional independent sub-band (ISB), a communication independent sub-band (CIB) and a radar independent sub-band (RIB).
Step 2.1: traditional individual sub-bandwidths:
in the ISB mode, the total bandwidth is divided into two parts, one part is used for communication and the other part is used for radar, and the two parts of communication and radar operate without any interference in the respective sub-bands.
Step 2.2: communication independent sub-bandwidths:
in the CIB mode, the total bandwidth is divided into two parts, one part is used for communication only, called communication independent sub-bandwidth, and the other part is used for radar and communication simultaneously, called hybrid sub-bandwidth.
Step 2.3: radar independent sub-bandwidth:
in the RIB mode, the total bandwidth is divided into two parts, one part is only used for radar and is called radar independent sub-bandwidth, and the other part is used for radar and communication simultaneously and is called mixed sub-bandwidth.
And step 3: characterization of radar communication integrated system:
step 3.1: performance characterization of radar function:
in view of AWGN in equation (3), radar reception signal zr(t) the probability density function (pdf) with respect to τ and θ is:
assuming that p (z; τ, θ) satisfies the "normal" condition, the corresponding Fisher Information Matrix (FIM) can be obtained:
after some algebraic derivations, the clarmero bound (CRB) for TOA and DOA are:
where γ represents the Signal to noise ratio (SNR),represents the mean square effective bandwidth and satisfies:
defining Δ τ and Δ θ as the estimation errors of τ and θ, respectively, the location estimate can be expressed as:
minimum Mean Squared Error (MMSE) J of a position estimatorminComprises the following steps:
according to the excitation of the estimated entropy and the random process entropy to the communication rate, the positioning estimation rate is defined as:
wherein, TpriWhere T/δ denotes pulse repetition interval, T denotes pulse width, δ denotes duty cycle, hrsAnd hestRepresenting the received signal entropy and the estimated entropy, respectively:
wherein n isproIs the remaining noise in addition to the estimated noise. The upper bound of the positioning estimation rate can be obtained from equations (13) and (14):
wherein,
step 3.2: performance characterization of communication functions:
using the communication rate RcomTo measure the performance of the communication function, the larger the communication rate is, the larger the maximum number of bits that can be transmitted on the channel in a unit time is, the stronger the information transmission capability is. According to Shannon's theorem:
wherein, T0Indicating the absolute temperature of the communication system, BcRepresenting the communication channel bandwidth.
And 4, step 4: performance limit analysis of radar communication integrated systems in different spectrum sharing modes:
step 4.1: analyzing performance limit of the radar communication integrated system in an ISB mode;
in ISB mode, let the radar bandwidth be BrA communication bandwidth of Bc。
Bc＝αB,Br＝(1-α)B (17)
Wherein α (0 is equal to or more than α is equal to or less than 1) is a bandwidth adjustment factor, BcFor communication of independent sub-bandwidths, BrIs a radar independent sub-bandwidth.
According to the formulas (15) and (16), the corresponding communication rate RcomAnd the positioning estimation rate is respectively:
step 4.2: performance limit analysis of the radar communication integrated system in the CIB mode:
in CIB mode, let communication independent sub-bandwidth be BcThe hybrid sub-bandwidth is BmixAnd simultaneously, the communication power is distributed by adopting a water injection method, the bandwidth regulating factor α is given (0 is more than or equal to α is more than or equal to 1), and the mathematical expressions of the two sub-bandwidths are as follows:
Bc＝αB,Bmix＝(1-α)B (20)
step 4.2.1: communication rate of communication independent sub-bandwidths
According to equation (16), the communication rate of the available communication independent sub-bandwidth is:
wherein p isc,oCommunication power allocated in the communication independent sub-bandwidth.
Step 4.2.2: communication rate and positioning estimation rate of hybrid sub-bandwidth:
the SIC is adopted to process the received composite signal in the mixed frequency band, and the joint receiver firstly decodes the received communication signal. Once the communication signal is successfully decoded, it is subtracted from the composite signal to obtain a radar signal free of communication signal interference.
Regarding radar signals as interference, the communication signal-to-interference ratio (SIR) of the hybrid sub-bandwidth is:
wherein p isc,mThe allocated communication power in the mixed sub-bandwidth.
Communication rates in the mixed sub-bandwidths according to equation (16)The mathematical expression of (a) is:
since the communication signal has been removed, the positioning estimation rate R in the mixed sub-bandwidth according to equation (15)estComprises the following steps:
step 4.2.3: communication power allocation:
communication power is distributed by water injection method to optimize pc,o、pc,mTo achieve maximum communication rates. The corresponding lagrange function is:
where λ is the lagrange multiplier. By solving the KKT point of the Lagrangian function, p can be obtainedc,oThe value range is as follows:
because of pc,o+pc,m＝pc，pc,mNot less than 0, further obtaining pcAnd prThe constraint relationship is as follows:
step 4.3: performance limit analysis of the radar communication integrated system in the RIB mode:
in RIB mode, let thunderUp to independent sub-bandwidth of BrThe hybrid sub-bandwidth is BmixGiven a bandwidth adjustment factor α (0 ≦ α ≦ 1), the mathematical expression for the two sub-bandwidths is:
Br＝αB,Bmix＝(1-α)B (28)
step 4.3.1: positioning estimation rate of radar independent sub-bandwidth
For radar signals in the radar independent sub-bandwidth, the radar signals are not interfered by communication signals, and the signal-to-noise ratio is as follows:
wherein p isr,oThe allocated radar power in the radar-independent sub-bandwidth. According to equation (15), the corresponding location estimation rate is:
step 4.3.2: communication rate and positioning estimation rate of hybrid sub-bandwidth:
and processing the received composite signal in the mixed frequency band by using the SIC, and decoding the received radar signal by using the combined receiver firstly. Once the radar signal is successfully decoded, it is subtracted from the composite signal to yield a communication signal free of radar signal interference.
Considering the communication signal as interference of a radar signal, the radar signal-to-interference ratio (SIR) in the mixed sub-bandwidth is:
wherein p isr,mIs prepared by mixingThe radar power allocated in the combined bandwidth, according to equation (15), the corresponding positioning estimation rate is:
since the radar signal has been removed, the communication rate R in the mixed sub-bandwidth is according to equation (16)comThe mathematical expression of (a) is:
step 4.3.3: radar power distribution:
like the principle of communication power distribution, the radar power is distributed by adopting a water injection method to optimize pr,o、pr,mThereby maximizing the positioning estimation rate.
The effect of the invention is further illustrated by the following simulation test:
simulation scene: assuming that the target cross-sectional area is known, the received power of the receiver follows a typical propagation loss model, i.e., the received signal power and r-nProportional, where r is the radar detection range or communication transmission range and n is the path loss exponent. The path loss exponents for radar and communication are set to 4 and 2, respectively. The working parameters of the radar communication integrated system are shown in table 1.
The simulation results of the performance limit of the radar communication integrated system in the three frequency spectrum sharing modes of ISB, CIB and RIB are shown in fig. 5.
Fig. 6 shows simulation results of the improvement of communication rate performance of the CIB and RIB spectrum sharing methods compared with the ISB spectrum sharing method.
The simulation result of the improved positioning estimation rate performance of the CIB and RIB spectrum sharing modes compared with the ISB spectrum sharing mode is shown in FIG. 7.
TABLE 1 Radar communication Integrated System working parameters
| Parameter(s) | Numerical value |
| Bandwidth (B) | 8MHz |
| Wavelength (lambda) | 0.3m |
| Absolute temperature (T)0) | 290K |
| Radar detection power (p)r) | 10KW |
| Radar detection range (r)r) | 10km |
| Radar detection antenna gain (g) | 1000 |
| Antenna element spacing (d) | 0.3m |
| Communication transmission power (p)c) | 50W |
| Communication transmission distance (r)c) | 10km |
| Number of antenna elements (N) | 5 |
| Pulse width (T) | 20μs |
| Target area of area (sigma) | 2m2 |
| Pulse duty cycle (delta) | 0.05 |
| DOA | π/2 |
Fig. 5 plots the performance limit of the radar communication integrated system in three frequency spectrum sharing modes of ISB, CIB and RIB. The external constraint indicates that under three frequency spectrum sharing modes, the maximum positioning estimation rate of 7890bit/s and the maximum communication data rate of 1.445 multiplied by 10 can be realized by allocating the total bandwidth to communication or radar respectively7bit/s as shown at points A and D, respectively, in FIG. 4. It can be seen that the performance limit in the conventional ISB mode is lower than that in the RIB mode and the CIB mode. Wherein, at point B, because all the bandwidth is allocated to the radar, the communication rate in the ISB and RIB modes is 0bit/s, but because the radar and the communication share all the bandwidth in the CIB mode, the communication rate can still reach 4.301 multiplied by 106bit/s. Similarly, at point C, since all the bandwidth is allocated to the communication, the positioning estimation rate in the ISB and CIB modes is 0bit/s, but since the radar and the communication share all the bandwidth in the RIB mode, the positioning estimation rate can still reach 2591 bit/s.
Fig. 6 is a graph showing the improvement of the communication rate performance in the CIB and RIB spectrum sharing scheme as compared with the ISB spectrum sharing scheme. It can be seen that compared with the ISB mode, the communication rates in the CIB and RIB modes are bothThe performance of the composite material is obviously improved. Wherein, the performance improvement effect is most obvious under the CIB mode, when the communication rate is 4.301 multiplied by 106When the signal is at bit/s, the positioning estimation rate can reach the maximum value.
Fig. 7 plots the improvement of the positioning estimation rate performance of the CIB and RIB spectrum sharing modes compared with the ISB spectrum sharing mode. Compared with the ISB mode, the performance of the positioning estimation rate under the CIB mode and the RIB mode is obviously improved. The performance improvement effect is most obvious in the RIB mode, and when the positioning estimation rate is 2591bit/s, the communication rate can reach the maximum value.
In conclusion, the performance limit of the radar communication integrated system can be well measured.
It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Various modifications and alterations to this invention will become apparent to those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the claims of the present invention.
Claims (9)
1. The method for measuring the performance limit of the radar communication integrated system based on the positioning estimation rate is characterized by comprising the following steps:
s1, receiving the signals by the receiving window at the same time: radar signals, communication signals, composite signals;
s2, the spectrum sharing mode is divided into: a traditional independent sub-bandwidth mode, a communication independent sub-bandwidth mode and a radar independent sub-bandwidth mode;
under the traditional independent sub-bandwidth mode, the total bandwidth is divided into two parts, one part is used for communication, and the other part is used for radar; the two parts work in respective sub-frequency bands without any interference;
under the communication independent sub-bandwidth mode, the total bandwidth is divided into two parts, one part is only used for communication and is called communication independent sub-bandwidth, and the other part is simultaneously used for radar and communication and is called mixed sub-bandwidth;
under the radar independent sub-bandwidth mode, the total bandwidth is divided into two parts, one part is only used for radar and is called radar independent sub-bandwidth, and the other part is simultaneously used for radar and communication and is called hybrid sub-bandwidth;
s3, obtaining a positioning estimation rate boundary expression representing radar function performance by considering the beam arrival time and the beam arrival angle joint information;
s4, representing the performance of the communication function by adopting the communication rate;
s5, calculating the communication rate and the estimated location rate boundary under different spectrum sharing modes in step S2.
2. The method for measuring performance limit of integrated radar communication system based on positioning estimation rate as claimed in claim 1, wherein the step S3 of obtaining the positioning estimation rate expression characterizing the radar function performance comprises:
a1, obtaining probability density functions of the radar signals received in the step S1 with respect to tau and theta; where τ represents the beam arrival time and θ represents the beam arrival angle;
a2, obtaining a corresponding snow information matrix according to the probability density function of the step A1;
a3, obtaining a Cramer-Rao bound of the arrival time of the wave beam and a Cramer-Rao bound of the arrival angle of the wave beam according to the Fisher-snow information matrix of the step A2;
a4, Cramer-Role bound according to beam arrival time, and Cramer-Role bound according to beam arrival angle; and obtaining an upper bound expression of the positioning estimation rate.
3. The method for measuring performance limit of radar communication integration system based on positioning estimation rate as claimed in claim 2, wherein step a4 specifically comprises:
a41, the positioning estimation rate is:
wherein h isrsRepresenting the entropy of the received signal, hestRepresenting the estimated entropy, TpriRepresenting a pulse repetition interval;
a42, obtaining the minimum mean square error of the positioning estimator according to the Cramer-Rao bound of the arrival time of the wave beam and the Cramer-Rao bound of the arrival angle of the wave beam;
a43, respectively obtaining h of the received radar signals according to the minimum mean square error of the positioning estimatorrsAnd hest；
A44, h obtained in the step A43rsAnd hestAnd substituting the expression of the step A41 to update the expression of the positioning estimation rate.
4. The method of claim 3, wherein the communication rate expression of step S4 is as follows:
wherein, BcDenotes the communication independent sub-bandwidth, h denotes the communication signal propagation gain, k is the boltzmann constant, T0Indicating the absolute temperature, p, of the communication systemcRepresenting the communication power in the total bandwidth.
5. The method of claim 4, wherein the communication rate and the position estimation rate boundary in the conventional single sub-bandwidth manner in step S5 are respectively:
wherein α is a bandwidth adjustment factor, 0 is equal to or more than α is equal to or more than 1, B represents a bandwidth, n is equal to or more than 1proRepresenting the remaining noise, except the estimated noise, gamma representing the signal-to-noise ratio, BrRepresenting the radar independent sub-bandwidth, delta representing the duty cycle, T representing the pulse width,n denotes the number of antenna elements, d denotes the antenna element spacing, c denotes the speed of light, λ is the plane wave wavelength, prDenotes the radar transmission power, and g denotes the propagation gain of the radar signal.
6. The method for measuring performance limit of integrated radar communication system according to claim 5, wherein the communication rate of the communication independent sub-bandwidth in the communication independent sub-bandwidth manner in step S5 is:
wherein p isc,oAnd the communication power distributed in the communication independent sub-bandwidth mode.
7. The method of claim 6, wherein the boundary calculation process of the communication rate and the position estimation rate of the mixed sub-bandwidth in the communication independent sub-bandwidth mode in step S5 is as follows:
b1, processing the received composite signal in the mixed frequency band by using SIC;
b2, the joint receiver firstly decodes the received communication signal;
b3, when the communication signal is successfully decoded, subtracting the communication signal from the composite signal to obtain a radar signal without communication signal interference;
b4, using the radar signal in the step B3 as interference, and obtaining the communication rate in the mixed sub-bandwidth as:
wherein p isc,mCommunication power allocated in the mixed sub-bandwidth in the communication independent sub-bandwidth mode;
b5, since the communication signal has been removed, according to the positioning estimation rate boundary expression of step S3, the positioning estimation rate boundary in the mixed sub-bandwidth is obtained as:
8. the method for measuring performance boundary of radar-communication integrated system according to claim 4, wherein the positioning estimation rate of radar independent sub-bandwidth in radar independent sub-bandwidth mode in step S5 is:
wherein, γr,oFor the communication signal-to-noise ratio, p, of the radar independent sub-bandwidth in the radar independent sub-bandwidth moder,oAnd allocating the radar power in the radar independent sub-bandwidth mode.
9. The method for measuring performance boundary of radar communication integration system based on positioning estimation rate as claimed in claim 8, wherein the boundary calculation procedure of communication rate and positioning estimation rate of mixed sub-bandwidth in radar independent sub-bandwidth mode in step S5 is as follows:
c1, processing the received composite signal in the mixed frequency band by using SIC;
c2, the joint receiver firstly decodes the received radar signal;
c3, when the radar signal is decoded successfully, subtracting the radar signal from the composite signal to obtain a communication signal without radar signal interference;
c4, using the communication signal in the step C3 as interference, obtaining the positioning estimation rate in the mixed sub-bandwidth:
wherein, γr,mCommunication signal-to-noise ratio, p, for mixed sub-bandwidths in radar independent sub-bandwidth moder,mAllocating radar power in the mixed sub-bandwidth in a radar independent sub-bandwidth mode;
c5, since the radar signal has been removed, according to the communication rate expression of step S4, the communication rate in the mixed sub-bandwidth is found to be:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910187415.7A CN109870676B (en) | 2019-03-13 | 2019-03-13 | Radar communication integrated system performance limit measuring method based on positioning estimation rate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910187415.7A CN109870676B (en) | 2019-03-13 | 2019-03-13 | Radar communication integrated system performance limit measuring method based on positioning estimation rate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109870676A true CN109870676A (en) | 2019-06-11 |
| CN109870676B CN109870676B (en) | 2021-01-12 |
Family
ID=66920328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910187415.7A Active CN109870676B (en) | 2019-03-13 | 2019-03-13 | Radar communication integrated system performance limit measuring method based on positioning estimation rate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109870676B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111025276A (en) * | 2019-11-21 | 2020-04-17 | 南京航空航天大学 | Bistatic radar optimal radio frequency stealth power distribution method under frequency spectrum coexistence environment |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107167785A (en) * | 2017-07-11 | 2017-09-15 | 吉林大学 | A kind of sane big array MIMO radar target transmitting-receiving angle combined estimation method |
| US20180188365A1 (en) * | 2013-10-13 | 2018-07-05 | Oculii Corp | Systems and methods for 4-dimensional radar tracking |
| CN108983226A (en) * | 2018-07-20 | 2018-12-11 | 北京航空航天大学 | It is a kind of to be structured the formation the MIMO radar communicating integral method of modulation based on antenna |
-
2019
- 2019-03-13 CN CN201910187415.7A patent/CN109870676B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20180188365A1 (en) * | 2013-10-13 | 2018-07-05 | Oculii Corp | Systems and methods for 4-dimensional radar tracking |
| CN107167785A (en) * | 2017-07-11 | 2017-09-15 | 吉林大学 | A kind of sane big array MIMO radar target transmitting-receiving angle combined estimation method |
| CN108983226A (en) * | 2018-07-20 | 2018-12-11 | 北京航空航天大学 | It is a kind of to be structured the formation the MIMO radar communicating integral method of modulation based on antenna |
Non-Patent Citations (2)
| Title |
|---|
| ALEX R. CHIRIYATH等: ""Inner Bounds on Performance of Radar and Communications Co-Existence"", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
| 陈伯孝: "《现代雷达系统分析与设计》", 30 September 2012, 西安电子科技大学出版社 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111025276A (en) * | 2019-11-21 | 2020-04-17 | 南京航空航天大学 | Bistatic radar optimal radio frequency stealth power distribution method under frequency spectrum coexistence environment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109870676B (en) | 2021-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105978641B (en) | The signal-noise ratio estimation method of time-frequency overlapped signal in a kind of cognitive radio | |
| Wu et al. | FILA: Fine-grained indoor localization | |
| CN105939188B (en) | The carrier frequency estimation method of time-frequency overlapping MASK signal under a kind of underlay frequency spectrum share | |
| US9071359B2 (en) | Method for noise floor and interference estimation | |
| Paul et al. | Extending joint radar-communications bounds for FMCW radar with Doppler estimation | |
| CN100518012C (en) | Authorization user signal detecting method for cognitive radio system | |
| Wyne et al. | A statistical model for indoor office wireless sensor channels | |
| WO2015096074A1 (en) | Method, base station, positioning server, and system for uplink reference signal configuration | |
| Angjelicinoski et al. | Comparative analysis of spatial interpolation methods for creating radio environment maps | |
| US9800279B2 (en) | Method and apparatus for automatic gain control in wireless receiver | |
| CN109861768B (en) | Radar communication integrated system performance analysis method based on mutual information | |
| CN109870676B (en) | Radar communication integrated system performance limit measuring method based on positioning estimation rate | |
| CN111398946A (en) | Shared signal design method based on L FM signal modulation frequency modulation | |
| CN1235429C (en) | Method for estimating position | |
| Saito et al. | A computational study of indoor-to-outdoor propagation in office environment at 2.4 GHz and 5.2 GHz bands | |
| Hu et al. | Secure crowdsourced radio environment map construction | |
| WO2018035757A1 (en) | Pilot frequency allocation method and pilot frequency allocation device | |
| Keunecke et al. | Deriving 2D TOA/TDOA IEEE 802.11 g/n/ac location accuracy from an experimentally verified fading channel model | |
| CN105915299B (en) | Spectrum prediction method based on time-frequency two-dimensional LMBP neural networks in ISM band | |
| CN109041215B (en) | CSI indoor AP positioning algorithm based on probability density | |
| CN103491550A (en) | Cognitive radio dynamic spectrum distribution method integrated with system overall transmission speed and distribution equity | |
| Vanhoy et al. | A spatial interpolation method for radio frequency maps based on the discrete cosine transform | |
| Shan et al. | The method of electromagnetic environment map construction based on Kriging spatial interpolation | |
| Nnebe et al. | Empirical analysis of signal-to-interference ratio variations in IEEE 802.11 b WLAN | |
| KR20080005685A (en) | Apparatus and method for bemaforming in a wireless communication network |
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 |