Capacity Analysis for Joint Radar-Communication Capable Coherent MIMO Radars
Muharrem Arik, Ozgur B. Akan

TL;DR
This paper analyzes the communication channel capacity of joint radar-communication coherent MIMO radars under various fading conditions, revealing optimal waveform strategies for maximizing capacity.
Contribution
It provides the first detailed analysis of how unintentional amplitude modulation affects communication capacity in JRC MIMO radars under different channel conditions.
Findings
Single-side limited null direction waveform yields highest capacity.
Capacity varies with Rayleigh and Rician fading conditions.
Unintentional amplitude modulation impacts communication performance.
Abstract
Recently, huge attention is attracted to the concept of integrating communication and radar missions within the same platform. Joint radar-communications (JRC) system gives an important opportunity to reduce spectrum usage and product cost while doing concurrent operation, as target sensing via radar processing and establishing communication links. A JRC-capable coherent MIMO radar system have been proposed recently in the literature. Several methods are introduced to reach dual goal as a notable null level towards the direction of interest of the radar and MIMO radar waveform orthogonality. Due to the limitations originated form the JRC operation, communication channel may encounter unwanted amplitude variations. This unwanted modulation normally affects the communication performance by its nature, due to the fades on radiated signal amplitude towards the direction of communication.…
| Number of MIMO radar antenna element | |
| Number of sub-pulse in a radar pulse | |
| s | transmit signal vector |
| Total transmit power budget | |
| Average transmit power towards communication direction | |
| Carrier signal frequency | |
| transmit beamforming vector of the sub-pulse | |
| space-time phase coding matrix | |
| Communication receiver direction | |
| Null direction towards target or clutter | |
| steering vector of the transmit array towards | |
| Radiated signal towards the direction of | |
| Magnitude of the radiated signal towards communication direction | |
| complex channel gain with amplitude and phase at time | |
| transmitted signal towards , | |
| when communication direction fixed method in [4] is used, | |
| when null direction fixed method in [4] is used, | |
| Positive number that defines communication pre-distortion level | |
| Additive White Gaussian Noise (AWGN) variance | |
| variance of the both real/imaginary part of the complex r.v. | |
| Complex conjugate operator | |
| Complex conjugate transpose operator |
| Difference | |||
|---|---|---|---|
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Capacity Analysis for Joint Radar-Communication Capable Coherent MIMO Radars
Muharrem Arik
Ozgur B. Akan
[email protected], [email protected]
Next-generation and Wireless Communications Laboratory (NWCL)
Department of Electrical and Electronics Engineering
Koc University, Istanbul, 34450, Turkey
Internet of Everything (IoE) Group
Electrical Engineering Division, Department of Engineering
University of Cambridge, Cambridge, CB3 0FA, UK
Abstract
Recently, huge attention is attracted to the concept of integrating communication and radar missions within the same platform. Joint Radar-Communications (JRC) system gives an important opportunity to reduce spectrum usage and product cost while doing concurrent operation, as target sensing via radar processing and establishing communication links. A JRC-capable coherent MIMO radar system have been proposed recently in the literature. Several methods are introduced to reach dual goal as a notable null level towards the direction of interest of the radar and MIMO radar waveform orthogonality. Due to the limitations originated form the JRC operation, communication channel may encounter unwanted amplitude variations. This unwanted modulation normally affects the communication performance by its nature, due to the fades on radiated signal amplitude towards the direction of communication. However, the effect of this unintentional modulation on communication channel is yet to be investigated. In this paper, the communication channel for JRC capable phase-coded coherent MIMO radars is analyzed and investigated under additive white Gaussian noise and Rayleigh/Rician fading conditions. Communication capacity is evaluated for each channel condition. The results reveal that, using the single-side limited null direction fixed waveform generation method displays the best capacity performance under all channel conditions.
keywords:
MIMO radar, joint communication and radar sensing, channel capacity, fading channels.
††journal: Physical Communication Journal
\WithSuffix
[7]G^ #1,#2_#3,#4( #7—#5 #6)
url]https://nwcl.ku.edu.tr/marik.html
url]https://ioe.eng.cam.ac.uk/directory/akan
1 Introduction
Several techniques are proposed over the years to enable JRC concept with different types of radar and communication systems and different topologies [1]. Almost all types of radars and all types of modulation are becoming the subject to this research domain. Mostly, efforts are concentrated on enhancements on data rates, modulation types and transparency over radar sensing operation.
Up to now, JRC systems are investigated mostly on radar perspective and the majority of the researchers have spent more efforts on finding out how much performance the radar system loses when communication is enabled. However, due to the limited resources and being a secondary mission on radar system, communication operation has often remained behind the scenes for a JRC system. To assess the success of a communication channel, the main issue to consider is channel capacity. Basically, the channel capacity is the information theoretic performance metric for a communication channel. There are several researches on communication channel capacity or performance limits of a JRC system. In [2], radar estimation rate concept was proposed to investigate JRC information limits. This rate is a measure of the amount of information collected from target in terms of information which might be exchanged uncooperatively between target and radar as a function of time. This type of metrics are used to optimize collaborative operation between radar and communication systems. There are many researches which consider communication capable radar waveform design. Each work tries to find the best waveform that provides invisible communication when flawless radar sensing operation is conducted. In [3], authors investigate high-SNR capacity of a communications channel based on some coding alphabet in the presence of AWGN. This alphabet presents a hybrid approach which constitutes a dual optimize JRC operation.
Recently, a JRC-capable monostatic coherent MIMO radar system have been proposed in [4]. Coherent MIMO radars offer fully coherent signal processing and coherent transmit beamforming which presents radiation patterns to reduce the probability of intercept. In [4], MIMO radar waveform coding technique in [5] is modified for enabling communication without disturbing coherent MIMO radar structure. JRC ability is accomplished for coherent MIMO radar without disturbing the orthogonality and transmit beamforming requirements. Three different novel communication methods are proposed in [4]. First method uses each sub-pulse in the radar pulse for communication. This method is directional and it may provide thousands of bits per radar pulse. Second method exploits coding mechanism through the radar pulse which is only visible at communication direction. The phase of the information vector is modulated for each radar pulse. Last communication method utilizes the progressive phase difference between orthogonal waveforms radiated from the first and the last antenna element caused by small degrees of steering on transmit beampattern. Specifically, null direction fixed phase-coded coherent MIMO radar waveform generation method in [4] produces an unintentional amplitude modulation during the transmission of the radar pulse towards the communication direction. In [4], in order to provide amplitude stable channel towards receiver, some constraints are introduced during the waveform generation phase. Fig.1 displays the fading like amplitude variations during the radar pulse with and constraints. Blue curve in Fig.1 shows the amplitude modulation without any constraint applied during waveform generation. This modulation affects the communication performance by its nature due to possible fades on signal amplitude through the radar pulse.
In this paper, first, we give channel models for different cases, then, we analyze the distribution characteristics of the unintentional modulation on signal amplitude. After that, using this distribution, we try to reach the capacity expressions for different communication methods under additive white Gaussian noise (AWGN) and Rayleigh/Rician fading conditions. Lastly, we evaluate the probability density function (pdf) of the distributions and capacity expressions for various communication methods proposed for JRC capable coherent MIMO radar. While analyzing the unintentional modulation under Rayleigh/Rician fading conditions, pdf of the distribution of product of a single or double truncated Rayleigh r.v. and a Rayleigh/Rician r.v. is firstly given in the literature.
The remainder is organized as follows. Firstly, the channel model is given in Sec. II. Then, amplitude distribution of the unintentional modulation is investigated in Sec. III and capacity expressions for JRC capable coherent MIMO radars are detailed in Sec. IV. Sec. V provides the numerical evaluation results. Lastly, conclusions are discussed in Sec. VI. For ease of exposition, the major notations used throughout this paper are listed in Table 1.
2 Channel Model
We utilize the same channel model as proposed in our previous work in [4]. Corresponding JRC system is equipped with element MIMO radar antenna aligned as a uniform linear array (ULA) as in Fig. 2. In the array, each antenna element has an isotropic radiation pattern with spacing distance of in terms of a wavelength. Moreover, each antenna radiates distinct phase coded quasi-orthogonal waveforms at the same carrier frequency. transmit signal vector during the -th sub-pulse of the -th radar pulse is given as,
[TABLE]
where and are the radar pulse and sub-pulse index, respectively and is the fast time index, stands for the complex conjugate transpose and carrier signal can be given as , where is the total transmit power budget and is the radar carrier frequency. Then, is the transmit beamforming weight vector of the sub-pulse as, , , where is the phase of the sub-pulse from the antenna element . For each radar pulse, angle of the space-time phase coding matrix for antennas and sub-pulses can be given as,
[TABLE]
where is the phasor angle function.
At the receiver side, single omni-directional antenna element is connected to a communication receiver located in direction , which is assumed to be known from the receiver. Besides, JRC enabled communication receiver is equipped with a matched filter. The received signal at the communication receiver for sub-pulse can be expressed as,
[TABLE]
where is the steering vector of the transmit array toward the spatial angle ,
[TABLE]
and is the complex channel coefficient which reflects the propagation gain between the transmit coherent MIMO array and the communication receiver. Then, is the additive white Gaussian noise with zero mean and variance . By combining (1) and (3), the baseband received signal for sub-pulse can be written as,
[TABLE]
After this point, we have dropped the index in the equations for the simplicity. Using information embedding strategies via manipulating transmit signal vector , communication symbols can be buried into the radar signal.
For the coherent MIMO radar case, the channel model for JRC capable coherent MIMO radar is expressed by,
[TABLE]
where , is the channel input and can be given as information modulated carrier signal, , where is the complex information symbol. can be defined as transmitted signal towards communication direction from the MIMO array. When communication fixed waveform generation method in [4] is applied, channel model can be simplified as,
[TABLE]
where is defined as transmitted signal towards communication direction from the MIMO array when communication direction fixed waveform generation method in [4] is used. is defined as the radiated signal towards communication direction for the first sub-pulse.
Null direction fixed method in [4] aims to fix the radiation signal at null direction using the technique in [5, 6], while satisfying higher signal levels for the communication direction. This method slightly changes the channel model for the communication due to the undesired variations on the amplitude level expressed in Eq.(8). This unintentional modulation shows random behavior and for the sake of simplicity, we have assumed that this unintentional modulation is the part of the channel effect. Therefore, we consider the channel model as in Fig.3 and it is given by,
[TABLE]
where is the channel input when null direction fixed waveform generation method in [4] is used and is defined as amplitude of the radiated signal towards communication direction . Although null direction fixed method introduces unintentional modulation, a phase correction procedure is proposed to fix the phase toward receiver. In [4], the correction factor is given as the phase angle differences between radiated signal at the first and the th sub-pulse, for . After applying this procedure, the phase component is fixed to along the radar pulse. This phase correction leads phase stable radiated signal. However, radiated signal amplitude towards communication direction still contains unintentional variations which may cause a degradation on communication performance.
3 Amplitude Distribution of the Unintentional Modulation
In order to reach the channel capacity of the null direction fixed method, it is necessary to characterize the distribution of the unintentional modulation. First, we must rewrite the radiated signal . When null direction fixed method is applied, the radiated signal becomes,
[TABLE]
where are the indexes of a random permutation sequence of and is the new phase component of the -th antenna element for the -th sub-pulse [4]. Then, this unintentional modulation caused by can be rewritten as a sum of individual exponentials as,
[TABLE]
For different values, is a phase variable in the range and is the radiated signal amplitude variable in the range .
Theorem 1
Amplitude of the radiated signal towards communication direction for JRC enabled coherent MIMO radar when null direction fixed method is applied presents a Rayleigh distribution with the following probability density function,
[TABLE]
where is the number of transmit antenna elements in the MIMO array, is the sub-pulse index and is the total number of sub-pulse in a radar pulse.
Proof 1
See the A.
4 Capacity Analysis
In this part, channel capacity is analyzed with the assumption that channel distribution information is known by the receiver. Then, we assumed that, the receiver perfectly knows the Channel State Information (CSI). For two waveform generation methods in [4], the channel capacity is expressed in detail. Then, capacity expressions are given for AWGN and RayleighRician channels.
4.1 Communication Direction Fixed Method - Stable Channel
This method guarantees the radiated signal is kept constant during the radar pulse. Therefore, the channel becomes an AWGN channel since and is constant during the radar pulse. Then, the channel model in Eq.(7) can be rewritten as,
[TABLE]
The channel capacity for the communication direction fixed method with AWGN channel can be defined in nats as, , where is the channel bandwidth and is the channel Signal-Noise Ratio (SNR). SNR can be given as, , where is the average transmit power towards communication direction, is the variance of the zero mean AWGN. For the communication direction fixed method, is equal to , where is the radiated signal amplitude through communication direction of the reference transmit beampattern, . Then, the capacity expression becomes,
[TABLE]
where . can be selected as, the probability of is greater than is equal to , . Since has Rayleigh distribution with the pdf in Theorem 1, can be given as,
[TABLE]
Then, becomes . Hence, the capacity formula can be rewritten as,
[TABLE]
4.2 Null Direction Fixed Method - Unintentionally Modulated Channel
Since the amplitude distribution of , which is equal to , shows Rayleigh distribution for the null direction fixed method, the unintentional modulation affects the channel as the Rayleigh fading does. Then, using Eq.(8) the channel can be represented as,
[TABLE]
for . The channel has the Rayleigh distribution by the following pdf, if there is no constraint applied in the null direction method, i.e. ,
[TABLE]
where . In [4], null direction method requires applying amplitude constraint as , where and . Since, the distribution is modified to a double truncated Rayleigh distribution, the probability density distribution is changed to a double truncated Rayleigh pdf as,
[TABLE]
where , and . By the transformation theorem for single random variables, the channel power gain has an exponential distribution with the mean . The pdf of is
[TABLE]
Using the (19), pdf of the SNR, becomes,
[TABLE]
where .
4.2.1 Capacity for AWGN channels
For this case, is selected as in Eq.(16). Then, fast fading (ergodic) conditions can be applied, since the channel has encountered all possible fades during the radar pulse when the radar pulse is long enough. The ergodic channel capacity can be expressed as,
[TABLE]
For the numerator integral, let and . By using integral by parts we get,
[TABLE]
where is the exponential integral for . The delta constraint in [4] is given as at the null direction method. Hence, the limits, and , are given as and . Eq.(22) can be rewritten for and ,
[TABLE]
If there is no given constraint at the null direction method, the limits, and , can be given as [math] and . Eq.(22) can be rewritten for and ,
[TABLE]
and if , the capacity can be approximated to the ergodic channel capacity of the Rayleigh fading channel with known CSI as,
[TABLE]
when , where is the average SNR which is and . If the constraint is given as at the null direction method, the limits, and , becomes and . Eq.(22) can be rewritten for and ,
[TABLE]
and if , the capacity can be approximated as,
[TABLE]
where is the average SNR and given as .
4.3 Null Direction Fixed Method Under Fading Channel
In this part, we consider the capacity of the JRC capable coherent MIMO radar when the null direction fixed method is used for waveform generation, where a single-user transmission is occurred over Rician and Rayleigh fading channel. For directional communications as sidelobe based JRC systems, most of the cases, receiver has a major LoS component along with multi-path reflections. For these cases, the fading channel is assumed as Rician fading channel. However, for some applications, receiver may lost its LoS component due to the mobility. This cases is considered as Rayleigh fading channel in our derivations. Therefore, the fading channel is selected for two different conditions as Rayleigh and Rician fading channels. Moreover, slow and fast fading cases are investigated. Slow fading case occurs when fading channel coherence time is bigger than communication symbol duration , i.e. . Then, fast fading (ergodic) case is also investigated. We assumed that CSI is known by the receiver only and the receiver is able to get the information of the fading states. The main channel model in Eq.(8) can be rewritten as,
[TABLE]
where is the channel input, is the channel output, and independent AWGN random variable with mean zero and variance and at time . is a complex channel gain with amplitude and phase at time t. The phase is uniformly distributed in , and the signal amplitude is a random variable with Rician and Rayleigh pdf. The null direction fixed method also fixes the phase towards communication direction applying a correction vector. Hence, phase component of the becomes constant over the radar pulse as . The channel model in Eq.(28) becomes,
[TABLE]
where is the product of two random variable and and the phase is again uniformly distributed in . Eq.(29) becomes similar to a fading channel model. Before attempting to investigate capacity expressions for null direction fixed method under fading channel, pdf of the new random variable needs to be introduced. , is a random variable with truncated Rayleigh distribution which is given in Eq.(18).
For Rayleigh fading cases, the pdf of can be given as,
[TABLE]
where is the variance of the both real and imaginary part of the complex r.v. .
Lemma 1
Any random variable , which is a product of two random variable, i.e. , which follows double truncated Rayleigh () and Rayleigh () distribution, presents a distribution with the following probability density function (pdf),
[TABLE]
where , and is equal to and . and are the variances of the both real and the imaginary part of the complex random variables and , respectively. and are the lower and the higher limits of the double truncated Rayleigh distributed r.v. , i.e. . is the generalized incomplete gamma function [12] as .
Proof 2
See the B.
Using the Lemma.1, pdf of the in Eq.(29) can be rewritten as,
[TABLE]
where , . By the transformation theorem for single random variables, the pdf of the channel power gain becomes,
[TABLE]
Using the (33), pdf of the SNR, becomes,
[TABLE]
where . We have also investigated the pdf of the product of single side truncated Rayleigh and Rayleigh random variables.
Corollary 1
Any random variable , which is a product of two r.v., , which follows single side truncated Rayleigh (), , and Rayleigh () distribution, presents a distribution with the following probability density function (pdf),
[TABLE]
Proof 3
By replacing with in Eq.(31), and becomes and [math], respectively. Then, related distribution pdf can be simply defined as in Eq.(35).
Using the Eq.(35), pdf of the for single side truncated case in Eq.(29) can be rewritten as,
[TABLE]
by doing the same transformation step in Eq.(33,34) we get,
[TABLE]
Without constraint in null direction fixed method, random variable becomes product of two Rayleigh variable and channel becomes double Rayleigh fading channel. Integral in Eq.(76) with limits and , can be expressed from [16, p. 370] as,
[TABLE]
where is the zeroth order modified Bessel function of the second kind. Then, pdf of the SNR becomes,
[TABLE]
The same derivations must be done for Rician fading cases, the pdf of can be given as,
[TABLE]
where is the modified Bessel function of the first kind with order zero and is the non-centrality parameter. The shape parameter can be given as, and the scale parameter is written as , which is defined as the total power received from multi-paths.
Lemma 2
*Any random variable , which is a product of two random variable, i.e. , which follows double truncated Rayleigh () and Rician () distribution, presents a distribution with the following probability density function (pdf), *
[TABLE]
where , and is equal to and . and are the variances of the both real and the imaginary part of the complex random variables and , respectively. and are the lower and the higher limits of the double truncated Rayleigh distributed r.v. , i.e. . is the generalized incomplete gamma function [12] as .
Proof 4
See the C.
Using the Lemma.2, pdf of the in Eq.(29) can be rewritten as,
[TABLE]
where . By the transformation theorem for single random variables, the pdf of the channel power gain becomes,
[TABLE]
Using the (43), pdf of the SNR, becomes,
[TABLE]
where and . We have also investigated the pdf of the product of single side truncated Rayleigh and Rayleigh random variables.
Corollary 2
*Any random variable , which is a product of two r.v., , which follows single side truncated Rayleigh (), , and Rayleigh () distribution, presents a distribution with the following probability density function (pdf), *
[TABLE]
Proof 5
By replacing with in Eq.(41), and becomes and [math], respectively. Then, related distribution pdf can be simply defined as in Eq.(45).
Using the Eq.(45), pdf of the for single side truncated case in Eq.(29) can be rewritten as,
[TABLE]
by doing the same transformation step in Eq.(33,34) we get,
[TABLE]
Without constraint in null direction fixed method, random variable becomes product of Rayleigh and Rician variable. Integral in Eq.(76) with limits and , can be expressed from [17, p. 58] as,
[TABLE]
where is the i-th order modified Bessel function of the second kind. Then, pdf of the SNR becomes,
[TABLE]
4.3.1 Capacity Under Slow fading,
For the slow fading case, communication symbol duration () is smaller than the fading channel coherence time (). For slowly varying channels where the instantaneous SNR is assumed to be constant for a large number of symbols, outage capacity needs to be presented. Over a big time period, channel data rate can reach up to capacity without fading case with negligible error. However, since the transmitter does not know the instantaneous SNR, a constant transmission data rate required, and it is independent from the instantaneous received SNR.
Particularly, a design parameter has to be selected which expresses a probability that the system can be in outage. In other words, is the probability that the system cannot correctly decrypt the transmitted communication symbols. The cause of the outage is the slow fading conditions. Hence, it is independent from unintentional modulation parameter . Then, can be defined as , where is the minimum SNR level that the received symbols cannot be correctly decoded with probability . At this level, the system declares an outage. The pdf of the is,
[TABLE]
where is , and shows the expected value of a random variable and given as . is the instantaneous channel gain resulting from fading and average channel power gain for Rayleigh fading is and for Rician fading, , where is the non-centrality parameter. Then, can be defined as,
[TABLE]
Taking out from Eq.(51),
[TABLE]
Then, the capacity with outage expression can be written as,
[TABLE]
where the pdf of the as,
[TABLE]
where . Also, average outage rate [9] can be defined as,
[TABLE]
Then, the ergodic capacity for coherent JRC system under slow fading can be expressed by replacing with in Eq.(21),
[TABLE]
Capacity expressions in Eq.(23) and Eq.(24) for coherent JRC methods with different constraints can be easily modified by replacing with .
4.3.2 Capacity Under Fast fading,
For the fast fading case, communication symbol duration () is larger than the channel coherence time (). The ergodic channel capacity can be reached for Rayleigh and Rician fading cases and different limit conditions using the related pdf expressions, after the channel has encountered all possible fades and it can be expressed as,
[TABLE]
Capacity expressions for the different communication methods can be easily obtained via modifying Eq.(57) for and constraints. Moreover, the expression in Eq.(57) can be solved numerically for all cases. With only exception, capacity expression for null direction fixed method without constraint, i.e. double Rayleigh fading case, can be given using the pdf in Eq.(39) in a closed form equation as,
[TABLE]
where is the Meijer G-function defined using a contour integral as,
[TABLE]
where and is the gamma function [14].
5 Numerical Analysis
In this section, first we investigate the amplitude distribution of the unintentional modulation given in Sec.3. Then, we display and compare the results of the channel capacity for the JRC capable coherent MIMO radars when the actual channel distribution is AWGN or RayleighRician.
5.1 Unintentional Modulation Amplitude Distribution Analysis
In Sec.(3), we have proved that the probability distribution of the amplitude distribution of the unintentional modulation follows a Rayleigh distribution. In this proof, the main requirement to reach this distribution indicates that the phase distribution along the radar pulse must converge to uniform distribution when the number of antenna elements has relatively large values.
We have analyzed the phase distribution, in Eq.(10), with two different number of antennas and , which are spaced with of a wavelength. In Fig.4, phase distribution for each case spans homogeneously over interval and the values are scattered around . Moreover, the expression in Eq.(68) for the pdf of is evaluated in Fig.4. The histogram in Fig.4 and the curve in Fig.4 show similar behavior, so we can make sure that the Eq.(68) almost represents the phase distribution. Calculations are also evaluated for and illustrated in Fig.4. Phase distribution of shows uniform distributions when takes large number. The curve in Fig.4 converges to line.
After making sure about the phase distribution which converges to uniform distribution by examining the figures. We have also presented the histograms for amplitude distribution of the unintentional modulation for two different number of antenna element cases, is equal to and . The histogram in Fig.4 almost follows the Rayleigh distribution pdf curve for . However, for the case , the histogram shows exactly same behavior as the curve of the pdf in Fig.4.
For different delta constraints, unintentional amplitude distributions are calculated and displayed in Fig.5. Each one of the figure includes some part of the Rayleigh distribution. Null fixed waveform generation algorithm mainly contains random permutations. Increase on the strictness of the constraints means longer required waveform generation times. Let the radar waveform has sub-pulses, and is selected as . Then, the required iterations to reach sub-pulses becomes approximated cycles. In Fig.6, total required random permutation iteration for null direction fixed waveform generation method is displayed when and . Approximation to a curve can be easily seen from the simulation results in the Fig.6.
5.2 Amplitude Distribution under Rayleigh Fading
Previously, we have proved that the unintentional modulation presents a Rayleigh distribution. Besides, numerical analysis in Sec.5.1 is also proved that JRC systems configured with higher number of MIMO antennas shows very similar behavior with truncated Rayleigh distribution. Further, we have also investigated the channel distribution under Rayleigh and Rician fading conditions. Calculations are done for three cases.
While evaluating the pdf of amplitude distribution of the null direction fixed method under Rician fading channel, after summation steps, the resulting sum converges to a curve rapidly. Therefore, evaluations are terminated after summation iterations.
5.2.1 Case 1: With Delta constraint
The analysis is evaluated for and , with the delta constraint for different values, i.e. for , and is and . The histogram of the product distribution , and all the pdf curves in Fig.7 show exactly the same behavior, hence we can make sure that the Eq.(1) and Eq.(2) totally represent the corresponding product distribution. For Rayleigh fading conditions, selecting smaller delta values results the product distribution more Rayleigh like shape, only with larger mean. This feature is also repeated for Rician case and resulting distribution follows Rician like behavior.
5.2.2 Case 2: With Single Side constraint
The analysis is evaluated for and , with the minimum constraint as when . Then, the variance of the both real and the imaginary part of the complex channel gain for Rayleigh and Rician channel, , is selected as . For Rician cases, shape parameter is selected as . The histogram of the product distribution , and the pdf curves in Fig.8 show exactly same behavior, hence we can make sure that the Eq.(35) and Eq.(45) totally represent the product distribution. Also Eq.(35) and Eq.(45) require a numerical calculation of generalized incomplete gamma function, . In D, a method for the numerical calculation which contains infinite summation is given. For our case, in , after step summation, the resulting sum converges to a curve rapidly. Determining a lower constraint results bigger mean value for the product distribution comparing with no constraint case in Fig.8.
5.2.3 Case 3: Without any constraint
The problem becomes easier without any constraint cases. For Rayleigh fading cases, the distribution can be calculated as cascaded double Rayleigh channel. Similary, for Rician fading cases, the distribution can be calculated as cascaded Rician and Rayleigh channel. The analysis is evaluated for and . Then, mean channel gain for Rayleigh channel is selected as which is when is , where is the variance of the both real and the imaginary part of the complex channel gain. The histogram of the product distribution , and the pdf curves in Fig.9 show exactly same behavior, hence we can make sure that the Eq.(38) totally represents the product distribution. Also, the product distribution for higher value displays more dispersive distribution as we have expected, hence it has bigger variance.
5.3 Coherent MIMO JRC Capacity Analysis
After determining the unintentional modulation distribution caused by the radiated signal towards communication direction, capacity formulas are evaluated for three different cases.
5.3.1 Case 1: AWGN Channel
Calculations are done for using the Eq.(22). Capacity of the communication direction fixed formula is evaluated when the radiated signal is constant at the level of . Results are displayed in Fig.10. As we have expected, applying no constraint, and , presents the worst capacity, which is also the ergodic channel capacity of the Rayleigh fading channel with known CSI, . Between these capacity curves in Fig.10, the techniques with a defined constraint with values, , are located. The capacity curves with higher values are getting closer to the no constraint case. More to that, increasing value with changing the probability always presents higher capacity due to its higher average power. On the other hand, smaller means higher search iterations, so higher search times during waveform generation process.
The capacity expressions in (25), (26) and (27) are calculated for three different values, as , and when is . Then, the differences between approximated expressions and general capacity formula are computed in Table.(2). If is greater than , the difference converges to significantly smaller values. Therefore, selecting as is enough to show the channel capacity behaviors.
5.3.2 Case 2: Slow Fading Rayleigh Channel
Before doing the corresponding calculations, the major design parameter needs to be determined to reach outage capacity results. must be selected as a number which will maximize the average rate in Eq.(55). For each communication method, we calculated the curves with respect to values. Capacity curves for Rayleigh fading is displayed in Fig.11 and for Rician fading in Fig.11. Outage probabilities for maximum rate values are displayed for each curves. These labeled values are used for the capacity with outage analysis are calculated for the slow fading case.
Capacity calculations are done for , and with using the Eq.(56). Results are displayed in Fig.11. Capacity of the communication direction fixed formula is evaluated when the radiated signal is constant at the level of . Since, we have considered the system in outage when symbols are received below the minimum SNR limit over a long period of time. More to that, the transmitter does not know instantaneous CSI. Then, the capacity with outage curves must display lower rate values. Capacity curves in Fig.11 shows the same behavior with AWGN case, however Rayleigh fading case presents lower capacity. On the other hand, Rician fading case displays more capacity due to its average power is derived by the shape parameter .
5.3.3 Case 3: Fast Fading Rayleigh Channel
In this case, we have assumed that channel has encountered all possible fades over a symbol period. Channel is characterized by the product distributions which are introduced and analyzed in the previous sections. This distribution has no truncated shape and it almost spans over all values zero to infinity. This results relatively small degradation on the capacity for all methods under fast fading conditions as seen on Fig.13. Capacity calculations are done for , and using the Eq.(57) and Eq.(58). Since all the capacity expressions are not in closed form, we have evaluated the integral with numerical methods using trapezoidal methods. These calculations are evaluated with MATLAB. Specifically, capacity curves for Rician fast fading case displays highest capacity by comparing all the other capacity curves.
6 Conclusion
The increasing interest on connected devices causes crowded spectrum. To solve this problem, various efforts on developing crafty and compact technologies are proposed as Joint Radar-Communication (JRC) systems. JRC ability can be established for coherent MIMO radar without disturbing the orthogonality and transmit beamforming requirements as in [4]. In this paper, the distribution characteristic of the unintentional modulation on signal amplitude is investigated. Then, the capacity expressions for different waveform generation methods under AWGN, slow and fast fading conditions are analyzed and evaluated. The results indicate that, using the null direction fixed waveform generation method with single constraint as where is equal to displays the best capacity curves for under all channel conditions. While reaching the capacity expressions, the pdf’s of the product of a double truncated Rayleigh r.v. and a Rayleigh/Rician distributed r.v. are expressed for the first time in the literature. Lastly, the resulting expressions given in this paper are also proven with the numerical evaluation.
Proposed MIMO radar waveform generation process in [4] under mobile scenarios may introduce computational complexity problems. To be more precise, the system must calculate the MIMO radar waveform in a very short period of time, since the direction of the receiver may change continuously due to the highly dynamic scenarios. The worst case may be that the MIMO waveform set changes from pulse to pulse. In order to reduce the computational cost, required number of iteration to generate MIMO waveform set can be decreased by loosing constraints. However, this will introduce some degradation on communication performance. This degradation is investigated in detail in [4]. MIMO waveform generation algorithm in [4] contains random process for each iteration. If the random permutation tables are generated off-line and stored to the system before the operation, the remaining process, i.e. summation and multiplication, will be fast and easy to implement on a system with a parallel processing capability.
Appendix A Proof of Amplitude Distribution of the Unintentional Modulation
The Rayleigh distribution is the distribution of the sum of a large number of coplanar vectors with constant amplitudes and uniformly distributed phases [10], [11]. The sum in Eq.(10) can be represented as a complex sum of random vectors as in Fig.14. This problem is very similar to the ”random-walk problem” in mathematical statistics. Hence, if presents a uniform distribution, will show Rayleigh distribution. Any phase distribution, , can be defined as a uniform distribution [11], if it satisfies the condition below,
[TABLE]
The first step is analyzing the phase component of , which is . can be written as a sum of two function as,
[TABLE]
where is the indexes of a random permutation sequence of for each sub-pulse and is a phase component of the -th antenna element and,
[TABLE]
where , , and are constants, equal to , , respectively. Let is one time selected randomly from and is a large number, . takes only different values, for each iteration these values are shuffled and summed with which is a linear function. Hence, takes values uniformly from the range . is a linear discrete function and takes M different values from the range .
Let , and are random variables which are selected from functions, , and respectively. Then, the distribution of can be written as,
[TABLE]
where () is convolution operation and probability density function of , , given as,
[TABLE]
The distribution of , can be expressed as,
[TABLE]
Let each value of is distinct, then . Then, is said to be each value of it is distinct due to is selected randomly, then becomes . The probability density function can be rewritten as,
[TABLE]
and rewriting Eq.(66),
[TABLE]
and since the is a phase component, it has periodic behavior as, . Therefore, is equal to . Therefore, we can expand dirac delta function with Fourier series expansion, (67) can be rewritten as,
[TABLE]
and after reordering Eq.(68),
[TABLE]
Since and are large numbers, we can always find a phase component couple for and , with a relation of,
[TABLE]
when . In order to prove this, can be written in terms of as,
[TABLE]
Hence, (69) can be rewritten as using the relation (71),
[TABLE]
Therefore, when , the probability density function approximates uniform distribution, , as,
[TABLE]
Since the phase distribution approximates uniform distribution, the distribution of sum of exponentials will follows Rayleigh distribution.
Appendix B Proof of Product Distribution of Truncated Rayleigh and Rayleigh Random Variables
There are two methods to reach product distribution. First is Mellin Convolution [15] which is defined as,
[TABLE]
where . Then, we can reach the pdf of the r.v. replacing the pdf functions with double truncated Rayleigh and Rayleigh distributions,
[TABLE]
where , , and by substituting using the transformation theorem we have,
[TABLE]
where is the generalized incomplete gamma function [12] as .
Second method uses Mellin transform to reach product distribution. Mellin transform of the product of two independent r.v. from different distributions is equal to the product of their Mellin transforms, . Mellin transform of a distribution can be given as [18],
[TABLE]
and the inverse Mellin transform is,
[TABLE]
where and the integration is along any path , such that exists. Then, Mellin transform of the r.v. can be expressed as,
[TABLE]
[TABLE]
where denotes the gamma function [14] and is the incomplete gamma integral [12], . Hence, Mellin transform of the product becomes,
[TABLE]
Then, by taking the inverse Mellin transform,
[TABLE]
product distribution can be defined with a contour integral above. If we replace by in the above equation, this does not affect the path of integration and the integral becomes,
[TABLE]
Generalized incomplete gamma functions have inverse Mellin transform representations from [13, Eq. 7.1] which can be given by,
[TABLE]
After applying the Eq.(84), can be rewritten as,
[TABLE]
Appendix C Proof of Product Distribution of Truncated Rayleigh and Rician Random Variables
Mellin Convolution in Eq.74 is also used for the derivation of the joint pdf of the product distribution of truncated Rayleigh and Rician r.v’s. We can reach the pdf of the r.v. replacing the pdf functions with double truncated Rayleigh and Rician distributions,
[TABLE]
where , . After applying the series representation of the which is , becomes,
[TABLE]
and by substituting using the transformation theorem we have,
[TABLE]
where is the generalized incomplete gamma function [12] as .
Appendix D A Numerical Method for the Calculation of Generalized Incomplete Gamma Functions
In [19], authors present a numerical method to compute generalized incomplete gamma functions based on Exponential Integrals with a relative precision of and an absolute precision of . This approach is also given in [20]. This method can be written as,
If , in , is an integer, :
[TABLE]
and,
[TABLE]
where
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] B. Paul, A. R. Chiriyath and D. W. Bliss, “Survey of RF Communications and Sensing Convergence Research,” in IEEE Access , vol. 5, pp. 252-270, 2017.
- 2[2] A. R. Chiriyath, B. Paul, G. M. Jacyna and D. W. Bliss, “Inner Bounds on Performance of Radar and Communications Co-Existence,” in IEEE Transactions on Signal Processing , vol. 64, no. 2, pp. 464-474, Jan.15, 2016.
- 3[3] I. Weiner, “High-SNR Channel Capacity for Communication Over Radar Waveforms,” in IEEE Trans. on Aerospace and Electronic Systems , vol. 55, no. 3, pp. 1261-1268, June 2019.
- 4[4] M. Arik, O. B. Akan “Realizing joint radar-communications in coherent MIMO radars,” in Elsevier Physical Communication , vol. 32, pp. 145-159, 2019.
- 5[5] H. Deng, Z. Geng and B. Himed, “MIMO Radar Waveform Design for Transmit Beamforming and Orthogonality,” in IEEE Trans. on Aerospace and Electronic Systems , vol. 52, no. 3, pp. 1421-1433, June 2016.
- 6[6] Z. Geng, H. Deng and B. Himed, “Waveform design for coherent MIMO radar radiation management and transmit-receive beam refocusing,” in Proc. of 2016 IEEE Radar Conference (Radar Conf) , Philadelphia, PA, 2016, pp. 1-5.
- 7[7] J. Li, A. Bose and Y. Q. Zhao, “Rayleigh flat fading channels’ capacity,” IEEE 3rd Annual Communication Networks and Services Research Conference (CNSR’05) , Halifax, NS, Canada, pp. 214-217, 2005.
- 8[8] A. Goldsmith and P. Varaiya, “Capacity of fading channels with channel side information,” in IEEE Trans. Information Theory , pp. 1986–1992, 1997.
