## Target Detection Performance of Spectrum Sharing MIMO Radars

Awais Khawar, Ahmed Abdelhadi, and T. Charles Clancy

Abstract -Future wireless communication systems are envisioned to share radio frequency (RF) spectrum, with other services such as radars, in order to meet the growing spectrum demands. In this paper, we consider co-channel spectrum sharing between cellular systems and radars. We address the problem of target detection by radars that are subject to shape its waveform in a way that it doesnt cause interference to cellular systems. We consider a multiple-input multiple-output (MIMO) radar and a MIMO cellular communication system with K base stations (BS). We propose a spectrum sharing algorithm which steers radar nulls, by projecting radar waveform onto the null space of interference channel, towards a 'selected' BS, thus, protecting it from radar interference. This BS is selected, among K BSs, on the basis of guaranteeing minimum waveform degradation. We study target detection capabilities of this null-space projected (NSP) waveform and compare it with the orthogonal waveform. We derive the generalized likelihood ratio test (GLRT) for target detection and derive detector statistic for NSP and orthogonal waveform. The target detection performance for NSP and orthogonal waveform is studied theoretically and via Monte Carlo simulations.

few cases, due to the fear of harmful interference from these services to radar systems. Recently, in the United States (U.S.), the Federal Communications Commission (FCC), has proposed to use the 3550-3650 MHz band for commercial broadband use [6]. The incumbents in this band are radar and satellite systems. The Commission has proposed that incumbents share this band with commercial communication systems. The Commission's spectrum sharing initiative is motivated by many factors including the President's National Broadband plan, which called to free up to 500 MHz of federal-held spectrum by 2020 [7]; surge in consumers' demand for access to mobile broadband, which operators can't meet with current spectrum allocation; the report on efficient spectrum utilization by President's Council of Advisers on Science and Technology (PCAST), which emphasized to share 1000 MHz of government-held spectrum [8], and the low utilization of 3550-3650 MHz band by federal incumbents [9].

Index Terms -MIMO Radar, Null Space Projection, Cellular System Coexistence, Target Detection, GLRT, ML Estimation.

## I. INTRODUCTION

Spectrum sharing between wireless communication systems and radars is an emerging area of research. In the past, spectrum has been shared primarily between wireless communication systems using opportunistic approaches by users equipped with cognitive radios [1]. This type of spectrum sharing has been made possible with the use of spectrum sensing [2], or geolocation databases [3], or a combination of both in the form of radio environment maps (REM) [4]. Some recent efforts have explored co-channel sharing approaches among secondary network entities, please see [5] and reference therein. However, in contrast, co-channel spectrum sharing between wireless systems and radars has received little attention thus far because of regulatory concerns.

Spectrum policy regulators, in the past, have not allowed commercial wireless services in the radar bands, except in

Awais Khawar (awais@vt.edu) is with Virginia Polytechnic Institute and State University, Arlington, VA, 22203.

Ahmed Abdelhadi (aabdelhadi@vt.edu) is with Virginia Polytechnic Institute and State University, Arlington, VA, 22203.

T. Charles Clancy (tcc@vt.edu) is with Virginia Polytechnic Institute and State University, Arlington, VA, 22203.

This work was supported by DARPA under the SSPARC program. Contract Award Number: HR0011-14-C-0027. The views expressed are those of the author and do not reflect the official policy or position of the Department of Defense or the U.S. Government.

Distribution Statement A: Approved for public release; distribution is unlimited.

In the future, when radio frequency (RF) spectrum will be shared among many different systems, e.g., radars and cellular systems, it is important to access the interference scenario. Of course, radars will cause interference to communication systems and vice versa if proper interference mitigation methods and novel sharing algorithms are not employed. In a study conducted by National Telecommunications and Information Administration (NTIA), it was observed that in order to protect commercial cellular communication systems, from high power radar signal, large exclusion zones are required [9]. These exclusion zones cover a large portion of the U.S. where majority of the population lives and, thus, does not make a business case for commercial deployment in radar bands. In order to share radar bands for commercial operation we have to address the interference mitigation techniques at both the systems. In this work, we focus on interference caused by radar systems to communication system and propose methods to mitigate this interference.

The federal-commercial spectrum sharing is not a new practice. In fact, in the past, commercial wireless systems have shared government bands on a low transmit power basis, in order to protect incumbents from interference. An example of such a scenario is wireless local area network (WLAN) in the 5250-5350 MHz and 5470-5725 MHz radar bands [10]. So the Commission's latest initiative to share 3.5 GHz radar band with small cells, i.e. wireless base stations operating on a low power, is in harmony with previous practices [6].

## A. Related Work

On going research efforts have shown that there are numerous ways to share spectrum between radars and commu-

nication systems. Cooperative sensing based spectrum sharing approaches can be utilized where a radar's allocated bandwidth is shared with communication systems [11]-[13]. A joint communication-radar platform can be envisioned in which a spectrally-agile radar performs an additional task of spectrum sensing and upon finding of unused frequencies it can change its operating frequency. In addition of spectrum sharing such a setting can enable co-located radar and communication system platforms for integrated communications and radar applications [14]-[17]. Radar waveforms can be shaped in a way that they don't cause interference to communication systems [18][23]. Moreover, database-aided sensing at communication systems [24] and beamforming approaches at MIMO radars can be realized for spectrum sharing [25].

## B. Our Contributions

The problem of target estimation, detection, and tracking lies at the heart of radar signal processing. This problem becomes critically important when we talk about sharing radar spectrum with other systems, say cellular systems. The regulatory work going on in the 3.5 GHz band to share radar spectrum with commercial systems is the motivation of this work. The focus of this work is to study target detection performance of a radar that is subject to share its spectrum with a cellular system. We consider spectrum sharing between a MIMO radar and a cellular system with many base stations. In our previous work, we have addressed the problem of radar waveform projection onto the null space of interference channel, in order to mitigate radar interference to communication system, and the problem of selection of interference channel for projection when we have a cellular system with K base stations [19]. In this work, we consider the same sharing scenario but study the target detection performance of radar for the null-space projected (NSP) waveform and compare it with that of the orthogonal radar waveform. We use the generalized likelihood ratio test (GLRT) for target detection and derive detector statistic for NSP and the orthogonal waveform.

## C. Notations

Matrices are denoted by bold upper case letters, e.g. A , and vectors are denoted by bold lower case letters, e.g. a . Transpose, conjugate, and Hermitian operators are denoted by ( ) · T , ( ) · ∗ , and ( ) · H , respectively. Moreover, notations used throughout the paper are provided in Table I along with descriptions for quick reference.

## D. Organization

This paper is organized as follows. Section II discuss MIMO radar, target channel, orthogonal waveforms, interference channel, and our cellular system model. Moreover, it also discusses modeling and statistical assumptions. Section III discusses spectrum sharing between MIMO radar and cellular system and introduces sharing architecture and projection algorithms. Section IV presents the generalized likelihood ratio test (GLRT) for target detection and derives detector statistic for NSP and orthogonal waveform. Section V discusses numerical results and compares performance of NSP and orthogonal waveform. Section VI concludes the paper.

TABLE I TABLE OF NOTATIONS

| Notation     | Description                                            |
|--------------|--------------------------------------------------------|
| x ( t )      | Transmitted radar (orthogonal) waveform                |
| a ( θ )      | Steering vector to steer signal to target angle θ      |
| y ( t )      | Received radar waveform                                |
| R x          | Correlation matrix of orthogonal waveforms             |
| s UE j ( t ) | Signal transmitted by the j th UE in the i th cell     |
| L UE i       | Total number of user equipments (UEs) in the i th cell |
| K            | Total number of BSs                                    |
| M            | Radar transmit/receive antennas                        |
| N BS         | BS transmit/receive antennas                           |
| H i          | i th interference channel                              |
| r i ( t )    | Received signal at the i th BS                         |
| P i          | Projection matrix for the i th channel                 |

## II. SYSTEM MODEL

In this section, we introduce preliminaries of MIMO radar, point target in far-field, orthogonal waveforms, interference channel, and cellular system model. Moreover, we also discuss modeling and statistical assumptions along with RF environment assumptions used throughout the paper.

## A. Radar Model

The radar we consider in this paper is a colocated MIMO radar with M transmit and receive antennas and is mounted on a ship. The colocated MIMO radar has antennas that have spacing on the order of half the wavelength. Another class of MIMO radar is widely-spaced MIMO radar where elements are widely-spaced which results in enhanced spatial diversity [26]. The colocated radar gives better spatial resolution and target parameter identification as compared to the widelyspaced radar [27].

## B. Target Model/Channel

In this paper, we consider a point target model which is defined for targets having a scatterer with infinitesimal spatial extent. This model is a good assumption and is widely used in radar theory for the case when radar elements are colocated and their exists a large distance between the radar array and the target as compared to inter-element distance [28]. The signal reflected from a point target with unit radar crosssection (RCS) is mathematically represented by the Dirac delta function.

## C. Signal Model

Let x ( ) t be the signal transmitted from the M -element MIMO radar array, defined as

$$\text{in} \, \infty \, \text{raw} \, \text{any}, \, \text{un} \, \infty \, \text{as} \\ \text{ratio} \\ \text{in} \, \infty \, \text{raw} \, \text{any} \, \infty \, \text{in} \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \inno} \\ \text{general} \quad & \text{where} \, x _ { k } ( t ) e ^ { j \omega _ { c } t } \text{ is the baseband signal from the $k^{th}$ transmit}$$

where x k ( ) t e jω t c is the baseband signal from the k th transmit element, ω c is the carrier angular frequency, t ∈ [0 , T o ] , with

T o being the observation time. We define the transmit steering vector as

Then, the transmit-receive steering matrix can be written as

$$\text{vector as} & \text{in} \\ \text{a} _ { T } ( \theta ) \stackrel { \triangle } { = } \left [ e ^ { - j \omega _ { c } \tau _ { T _ { 1 } } ( \theta ) } \ e ^ { - j \omega _ { c } \tau _ { T _ { 2 } } ( \theta ) } \ \dots \ e ^ { - j \omega _ { c } \tau _ { T _ { M } } ( \theta ) } \right ] ^ { T }. & \text{in} \, \text{this} \\ \text{antenna} \colon \\ \text{Then, the transmit-receive steering matrix can be written as} &.$$

$$\mathbf A ( \theta ) \stackrel { \triangle } { = } \mathbf a _ { R } ( \theta ) \mathbf a _ { T } ^ { T } ( \theta ). \quad \quad \ \ ( 3 ) \quad \text{the } j$$

Since, we are considering M transmit and receive elements, we define a ( θ ) /defines a T ( θ ) /defines a R ( θ ) . The signal received from a single point target, in far-field with constant radial velocity v r , at an angle θ can be written as

$$\mathbf y ( t ) = \alpha \, e ^ { - j \omega _ { D } t } \mathbf A ( \theta ) \, \mathbf x ( t - \tau ( t ) ) + \mathbf n ( t ) \quad ( 4 ) \quad \text{where}$$

where τ ( ) = t τ T k ( )+ t τ R l ( ) t , denoting the sum of propagation delays between the target and the k th transmit element and the l th receive element, respectively; ω D is the Doppler frequency shift, and α represents the complex path loss including the propagation loss and the coefficient of reflection.

## D. Modeling Assumptions

In order to keep the analysis tractable we have made the following assumptions about our signal model:

- · The path loss α is assumed to be identical for all transmit and receive elements, due to the far-field assumption [29].
- · The angle θ is the azimuth angle of the target.
- · After compensating the range-Doppler parameters, we can simplify equation (4) as

$$\mathbf y ( t ) = \alpha \, \mathbf A ( \theta ) \, \mathbf x ( t ) + \mathbf n ( t ). \quad \quad ( 5 ) \quad \text{at th}$$

## E. Statistical Assumptions

We make the following assumptions for our received signal model in equation (5):

- · θ and α are deterministic unknown parameters representing the target's direction of arrival and the complex amplitude of the target, respectively.
- · The noise vector n ( ) t is independent, zero-mean complex Gaussian with known covariance matrix R n = σ 2 n M I , i.e. n ( ) t ∼ N c ( 0 M , σ 2 n M I ) , where N c denotes the complex Gaussian distribution.
- · With the above assumptions, the received signal model in equation (5) has an independent complex Gaussian distribution, i.e., y ( ) t ∼ N c ( α A ( θ ) x ( ) t , σ 2 n I M ) .

## F. Orthogonal Waveforms

In this paper, we consider orthogonal waveforms transmitted by MIMO radars, i.e.,

$$4 O \text{ radars, 1.e.,} \\ R _ { x } = \int _ { T _ { o } } x ( t ) x ^ { H } ( t ) d t = \mathbf I _ { M }. \quad ( 6 ) \quad S p \cdot \\ n \text{smission of orthogonal solnals oives MIMO radar} \ \text{snec}.$$

The transmission of orthogonal signals gives MIMO radar advantages in terms of digital beamforming at the transmitter in addition to receiver, improved angular resolution, extended array aperture in the form of virtual arrays, increased number of resolvable targets, lower sidelobes, and lower probability of intercept as compared to coherent waveforms [29].

## G. Communication System

In this paper, we consider a MIMO cellular system, with K base stations, each equipped with N BS transmit and receive antennas, with i th BS supporting L UE i user equipments (UEs). The UEs are also multi-antenna systems with N UE transmit and receive antennas. If s UE j ( ) t is the signals transmitted by the j th UE in the i th cell, then the received signal at the i th BS receiver can be written as

$$\text{ceiver can be written as} \\ \text{r} _ { i } ( t ) = & \sum _ { j } \text{H} _ { j } ^ { N ^ { \text{BS} } \times N ^ { \text{UE} } } s _ { j } ^ { \text{UE} } ( t ) + \text{w} ( t ) \quad \quad ( 7 ) \\ & \text{ for } 1 \leq i \leq \mathcal { K } \text{ and } 1 \leq j \leq \mathcal { L } _ { i } ^ { \text{UE} } \\ \dots$$

where w ( ) t is the additive white Gaussian noise.

## H. Interference Channel

In this section, we characterize the interference channel that exists between a MIMO cellular base station and a MIMO radar. In our paper, we are considering K cellular BSs that's why our model has H i , i = 1 2 , , . . . , K , interference channels, where the entries of H i are denoted by

$$\mathring { \ast } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \mathring { t } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \mathbb { t } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \

\mathring { t } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \. \quad \quad \quad$$

where h ( l,k ) i denotes the channel coefficient from the k th antenna element at the MIMO radar to the l th antenna element at the i th BS. We assume that elements of H i are independent, identically distributed (i.i.d.) and circularly symmetric complex Gaussian random variables with zero-mean and unitvariance, thus, having a i.i.d. Rayleigh distribution.

## I. Cooperative RF Environment

In the wireless communications literature, it is usually assumed that the transmitter (mostly BSs) has channel state information either by feedback from the receiver (mostly UEs), in FDD systems [30], or transmitters can reciprocate the channel, in TDD systems [30]. The feedback and reciprocity are valid and practical as long as the feedback has a reasonable overhead and coherence time of the RF channel is larger than the two-way communication time, respectively.

In the case of radars sharing their spectrum with communications systems one way to get CSI from communication systems is through feedback. Since, a radar signal is treated as interference at a communication system, we can characterize the channel as interference channel and refer to information about it as interference-channel state information (ICSI).

Spectrum sharing between radars and communications systems can be envisioned in two domains: military radars sharing spectrum with military communication systems, we call it Mil2Mil sharing; another possibility is military radars sharing spectrum with commercial communication systems, we call it Mil2Com sharing. In Mil2Mil sharing, ICSI can be acquired by radars fairly easily as both systems belong to military. In Mil2Com sharing, ICSI can be acquired by giving incentives

Fig. 1. Spectrum Sharing Scenario: A seaborne MIMO radar detecting a point target while simultaneously sharing spectrum with a MIMO cellular system without causing interference to the cellular system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0540v2-figure-2628dbd2-262e-41a4-a246-68c2821dd45e.jpg)

to commercial communication system. The biggest incentive in this scenario is null-steering and protection from radar interference. Thus, regardless of the sharing scenario, Mil2Mil or Mil2Com , we have ICSI for the sake of mitigating radar interference at communication systems.

## III. RADAR-CELLULAR SYSTEM SPECTRUM SHARING

After introducing our radar and cellular system models we can now discuss the spectrum sharing scenario between radar and cellular system. In our sharing architecture MIMO radar and cellular systems are the co-primary users of the 35503650 MHz band under consideration. In the following sections, we will discuss the architecture of spectrum sharing problem which is followed by our spectrum sharing algorithm.

## A. Architecture

We illustrate our coexistence scenario in Figure 1 where the maritime MIMO radar is sharing K interference channels with the cellular system. Considering this scenario, the received signal at the i th BS receiver can be written as

$$& \text{sigma at one } \imath \vdash \text{b.s. receiver can one written as} \\ & \text{r} _ { i } ( t ) = \mathbf H _ { i } ^ { N \text{BS} \times M } \mathbf x ( t ) + \sum _ { j } \mathbf H _ { j } ^ { N \text{BS} \times N ^ { \text{UE} } } s _ { j } ^ { \text{UE} } ( t ) + \mathbf w ( t ). \ \ ( 9 ) \quad \text{Proper} \quad \text{onto the} \\ & \text{The } \text{soal of the MIMO radar is to man } \mathbf v ( t ) \text{ onto the null}.$$

The goal of the MIMO radar is to map x ( ) t onto the nullspace of H i in order to avoid interference to the i th BS, i.e., H x i ( ) = t 0 , so that r i ( ) t has equation (7) instead of equation (9).

## B. Projection Matrix

In this section, we define the projection algorithm which projects radar signal onto the null space of interference channel H i . Assuming, the MIMO radar has channel state information of all H i interference channels, through feedback, in Mil2Mil or Mil2Com scenario, we can perform singular value decomposition (SVD) to find the null space and then construct a projector matrix. We proceed by first finding SVD of H i , i.e.,

$$\mathbf H _ { i } = \mathbf U _ { i } \Sigma _ { i } \mathbf V _ { i } ^ { H }. \quad \quad \ \ ( 1 0 ) \ \text{Onc}$$

Now, let us define

$$\widetilde { \Sigma } _ { i } \stackrel { \vartriangle } { = } \text{diag} ( \widetilde { \sigma } _ { i, 1 }, \widetilde { \sigma } _ { i, 2 }, \dots, \widetilde { \sigma } _ { i, p } ) \quad \text{cnanr} ( 1 1 )$$

where p /defines min( N BS , M ) and ˜ σ i, 1 &gt; σ ˜ i, 2 &gt; · · · &gt; σ ˜ i,q &gt; ˜ σ i,q +1 = ˜ σ i,q +2 = · · · = ˜ σ i,p = 0 are the singular values of H i . Next, we define where

$$\dots, \dots \dots \dots \int _ { i } ^ { \prime } \overset { \triangle } { \Sigma } _ { i } ^ { \prime } & \overset { \triangle } { = } \text{diag} ( \widetilde { \sigma } _ { i, 1 } ^ { \prime }, \widetilde { \sigma } _ { i, 2 } ^ { \prime }, \dots, \widetilde { \sigma } _ { i, M } ^ { \prime } ) \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \sum _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _{ i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i j } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime }\int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int_ i ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ * \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \end{ \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \int _ { i } ^ { \prime } \begin{array} { c } \int _ { i } \end{array}$$

Using above definitions we can now define our projection matrix, i.e.,

$$\widetilde { \sigma } _ { i, u } ^ { \prime } \stackrel { \triangle } { = } \begin{cases} 0, & \text{for $u \leq q$}, \\ 1, & \text{for $u > q$}. \end{cases} \quad ( 1 3 ) \\ \text{ove definitions we can now define our projection} \\., & \quad.$$

P i /defines V Σ V i ˜ ′ i H i . (14) In order to show that P i is a valid projection matrix we prove two results on projection matrices below.

Property 1. P i ∈ C M M × is a projection matrix if and only if P i = P H i = P 2 i .

Proof: Let's start by showing the 'only if' part. First, we show P i = P H i . Taking Harmition of equation (14) we have

Now, squaring equation (14) we have

$$\cdot \quad \begin{array} { c } \cdot & \cdot \\ \mathbf P _ { i } ^ { H } = ( \mathbf V _ { i } \widetilde { \Sigma } _ { i } ^ { \prime } \mathbf V _ { i } ^ { H } ) ^ { H } = \mathbf P _ { i }. \end{array} \quad & ( 1 5 ) \\ \intertext { a r i n g e q u a t i o n } \text{aring equation (14)} \, \text{ we have} \quad & \mathbf D ^ { 2 } = \mathbf V _ { i } \widetilde { \mathbf B } _ { i } \mathbf V _ { i } ^ { H } = \mathbf D _ { i } \quad & \mathbf V _ { i } \cdot \end{array}$$

P 2 i = V Σ V i ˜ i H i × V Σ V i ˜ i H i = P i (16) where above equation follows from V V H i i = I (since they are orthonormal matrices) and ( ˜ Σ ′ i ) 2 = ˜ Σ ′ i (by construction). From equations (15) and (16) it follows that P i = P H i = P 2 i . Next, we show P i is a projector by showing that if v ∈ range ( P i ), then P v i = v , i.e., for some w v , = P w i , then

$$P _ { i } \mathbf v = P _ { i } ( \mathbf P _ { i } \mathbf w ) = P _ { i } ^ { 2 } \mathbf w = P _ { i } \mathbf w = \mathbf v. \quad ( 1 7 )$$

Moreover, P v v null( P ), i -∈ i i.e.,

$$\text{ith} _ { \cdot } \quad P _ { i } ( P _ { i } \nu - \nu ) = P _ { i } ^ { 2 } \mathbf v - P _ { i } \mathbf v = P _ { i } \mathbf v - P _ { i } \mathbf v = 0. \quad ( 1 8 )$$

This concludes our proof.

Property 2. P i ∈ C M M × is an orthogonal projection matrix onto the null space of H i ∈ C N BS × M .

Proof: Since P i = P H i , we can write

H P i H i = U Σ V i ˜ i H i × V Σ V i ˜ ′ i H i = 0 . (19) The above results follows from noting that ˜ Σ Σ i ˜ ′ i = 0 by construction.

$$\text{sense, i.e.,} \\ i _ { \min } \stackrel { \triangle } { = } \arg \min _ { 1 \leq i \leq \mathcal { X } } \left | \left | \mathbf P _ { i } x ( t ) - x ( t ) \right | \right | _ { 2 } \quad ( 2 0 ) \\ \check { \mathbf P } \stackrel { \triangle } { = } \mathbf P _ { i _ { \min } }. \quad \quad ( 2 1 ) \\ \text{$two$ calculated our deviation matrix it is straight for}$$

In this paper, we are dealing with K interference channels. Therefore, we need to select the interference channel which results in least degradation of radar waveform in a minimum norm sense, i.e.,

Once we have selected our projection matrix it is straight forward to project radar signal onto the null space of interference channel via

$$\breve { x } ( t ) = \breve { \mathbf P } \, \mathbf x ( t ). \quad \quad \quad ( 2 2 )$$

The correlation matrix of our NSP waveform is given as

$$\i o n \, \text{matrix} x \, \text{or} \, \text{in} \, \text{$\nu$ waveform is given as} \\ \mathbf R _ { \check { x } } = \int _ { T _ { o } } \check { x } ( t ) \check { x } ^ { H } ( t ) d t & ( 2 3 ) \\ \text{longer identity and its rank denends union the rank}$$

which is no longer identity and its rank depends upon the rank of the projection matrix.

## C. Spectrum Sharing and Projection Algorithms

The process of spectrum sharing by forming projection matrices and selecting interference channels is executed with the help of Algorithms (1) and (2). First, at each pulse repetition interval (PRI), the radar obtains ICSI of all K interference channels. This information is sent to Algorithm (2) for the calculation of null spaces and formation of projection matrices. Algorithm (1) process K projection matrices, received from Algorithm (2), to find the projection matrix which results in least degradation of radar waveform in a minimum norm sense. This step is followed by the projection of radar waveform onto the null space of the selected BS, i.e, the BS to the corresponding selected projection matrix, and waveform transmission.

## Algorithm 1 Spectrum Sharing Algorithm

loop for i = 1 : K do Get CSI of H i through feedback from the i th BS. Send H i to Algorithm (2) for the formation of projection matrix P i . th from Algorithm

Receive the i projection matrix P i (2).

$$& \mathbb { A }. \\ & \text{end for} \\ & \text{End} \, i _ { \min } = \arg \min _ { 1 \leq i \leq \mathcal { X } } \left \| \left | \mathsf P _ { i } x ( t ) - x ( t ) \right \| _ { 2 }. \\ & \text{Set } \, \mathsf P = \mathsf P _ { i _ { \min } } \, \text{as the desired projector.} \\ & \text{Perform null space projection, i.e., } \, \breve { x } ( t ) = \, \mathsf P x |$$

Perform null space projection, i.e., ˘( ) = x t ˘ Px ( ) t . end loop

## Algorithm 2 Projection Algorithm

if H i received from Algorithm (1) then

Perform SVD on H i (i.e. H i = U Σ V i i H i )

Construct ˜ Σ i = diag ( ˜ σ i, 1 , σ ˜ i, 2 , . . . , σ ˜ i,p ) Construct ˜ Σ ′ i = diag ( ˜ σ ′ i, 1 , σ ˜ ′ i, 2 , . . . , σ ˜ ′ i,M ) Setup projection matrix P i = V Σ V i ˜ ′ i H i . Send P i to Algorithm (1).

end if

## IV. STATISTICAL DECISION TEST FOR TARGET DETECTION

In this section, we develop a statistical decision test for target illuminated with the orthogonal radar waveforms and the NSP projected radar waveforms. The goal is to compare the performance of the two waveforms by looking at the test decision on whether the target is present or not in the rangeDoppler cell of interest. We present a system-level architecture

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0540v2-figure-64914897-a0a5-48b0-9c94-331acfb1f9d7.jpg)

Receiver

Transmitter

Fig. 2. Block diagram of spectrum sharing radar. The transmitter is modified to perform the functions of ICSI collection, projection matrix formation, interference channel selection, and radar waveform projection on to the selected interference channel for spectrum sharing. On the other hand, the receiver is a traditional radar receiver performing functions of parameter detection and estimation on radar returns.

of the spectrum sharing radar in Figure 2. In our architecture, the transmitter performs the functions of waveform generation, channel selection, and projection; and the receiver performs the functions of signal detection and estimation.

For target detection and estimation, we proceed by constructing a hypothesis test where we seek to choose between two hypothesis: the null hypothesis H 0 which represents the case when the target is absent or the alternate hypothesis H 1 which represents the case when the target is present. The hypothesis for a single target model in equation (5) can be written as

$$& \quad \text{when as} \\ - \quad & \quad \text{y} ($$

Since, θ and α are unknown, but deterministic, we use the generalize likelihood ratio test (GLRT). The advantage of using GLRT is that we can replace the unknown parameters with their maximum likelihood (ML) estimates. The ML estimate of α and θ are found for various signal models, targets, and interference sources in [29], [31] when using orthogonal signals. In this paper, we consider a simpler model with one target and no interference sources in order to study the impact of NSP on target detection in a tractable manner. Therefore, we present a simpler derivation of ML estimation and GLRT.

The received signal model in equation (5) can be written as

$$\mathbf y ( t ) = \mathbf Q ( t, \theta ) \alpha +$$

where

$$Q ( t, \theta ) = \mathbf A ( \theta ) \mathbf x$$

We use Karhunen-Lo` eve expansion for derivation of the log-likelihood function for estimating θ and α . Let Ω denote the space of the elements of { y ( ) t } { , Q ( t, θ ) } , and { n ( ) t } .

Moreover, let ψ z , z = 1 2 , , . . . , be an orthonormal basis function of Ω satisfying

$$\begin{array} {$$

where δ zz ′ is the Kr¨necker delta function. Then, the following o series can be used to expand the processes, { y ( ) t } { , Q ( t, θ ) } , and { n ( ) t } , as

$$y ( t ) = \sum _ { z = 1 } ^ { \infty } y _ {$$

$$\mathbf n ( t ) = \sum _ { z = 1 } ^ { z = 1 } \mathbf n _ { z } \psi _ { z } ( t ) \quad \quad ( 3 0 ) \quad & \text{Simi} \\.. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

$$\mathbf Q ( t, \theta ) = \sum _ { z = 1 } ^ { z = 1 } \mathbf Q _ { z } ( \theta ) \psi _ { z } ( t ) \quad \quad ( 2 9 ) \quad \text{where}$$

where y z , Q z , and n z are the coefficients in the KarhunenLo` eve expansion of the considered processes obtained by taking the corresponding inner product with basis function φ z ( ) t . Thus, an equivalent discrete model of equation (25) can be obtained as

$$y _ { z } = \mathbf Q _ { z } ( \theta ) \alpha + \mathbf n _ { z }, \ \ z = 1, 2, \dots \quad \ \ ( 3 1 )$$

For white circular complex Gaussian processes, i.e, E [ n ( ) t n ( t -τ ( ))] t = σ 2 n M I δ τ ( ( )) t , the sequence { n z } is i.i.d. and n z ∼ N c ( 0 M , σ 2 n M I ) . Thus, we can express the log-likelihood function as

$$L _ { y } = \\ L _ { y } ( \theta, \alpha ) = \sum _ { z = 1 } ^ { \infty } \left ( - M \log ( \pi \sigma _ { n } ^ { 2 } ) - \frac { 1 } { \sigma _ { n } ^ { 2 } } \left | \left | y _ { z } - \mathbf Q _ { z } ( \theta ) \alpha \right | \right | ^ { 2 } \right ). \quad \text{where $f$} \\ \text{$city$fun} \\ \mathcal { H } _ { 0 }, \text{res} ]$$

Maximizing with respect to α yields

$$L _ { y } ( \theta, \hat { \alpha } ) = \Gamma - \frac { 1 } { \sigma _ { n } ^ { 2 } } \left ( E _ { y y } - e _ { Q y } ^ { H } \mathbf E _ { Q Q } ^ { - 1 } e _ { Q y } \right ) \quad ( 3 3 ) \quad L _ { s } \\ \intertext { a h a }$$

where

$$\Gamma \stackrel { \triangle } { = } - M \log ( \pi \sigma _ { n } ^ { 2 } ) \text{ \quad \ \ } ( 3 4 )$$

$$1 & = - M \log ( \pi \sigma _ { n } ^ { - } ) & & ( 3 4 ) \\ E _ { y y } & \stackrel { \triangle } { = } \sum _ { z = 1 } ^ { \infty } \left \| y _ { z } \right \| ^ { 2 } & & ( 3 5 ) & & \text{where} \\ e _ { Q y } & \stackrel { \triangle } { = } \sum _ { z } ^ { \infty } Q _ { z } ^ { H } y _ { z } & & ( 3 6 ) & & \text{.}$$

$$E _ { Q Q } ^ { - 1 } & \stackrel { \Delta } { = } \sum _ { z = 1 } ^ { \infty } Q _ { z } ^ { H } Q _ { z }. & \quad ( 3 7 ) & \quad. \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon$$

$$e _ { Q y } & \stackrel { \triangle } { = } \sum _ { z = 1 } ^ { z = 1 } Q _ { z } ^ { H } y _ { z } & \quad ( 3 6 ) & \quad \bullet \\ \longrightarrow & - 1 & \wedge \sum _ { z = 1 } ^ { \infty } \wedge H \wedge & \quad \dots$$

Note that, in equation (33), apart from the constant Γ , the remaining summation goes to infinity. However, due to the non-contribution of higher order terms in the estimation of θ and α the summation can be finite. Using the identity

$$\int _ { T _ { o } } \mathbf v _ { 1 } ( t ) \mathbf v _ { 2 } ^ { H } ( t ) d t = \sum _ { z = 1 } ^ { \infty } \mathbf v _ { 1 z } \mathbf v _ { 2 z } ^ { H } \, \int \quad \ \ ( 3 8 )$$

for v i ( ) = t ∑ ∞ z =1 v 1 z ψ z ( ) t , i = 1 2 , , equations (35)-(37) can be written as

$$- s s & J _ { T _ { o } } | | ^ { \nu \nu ^ { \prime } } | | ^ { \dots } & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \ e _ { \mathbf Q y } & \stackrel { \triangle } { = } & \int _ { T _ { o } } \mathbf Q ^ { H } ( t, \theta ) \mathbf y ( t ) d t & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\int _ { \substack { - \tilde { \psi } \\ \text{en as} } } \Delta ^ { \dagger } & \Delta ^ { \dagger } \Delta ^ { \dagger } \\ E _ { \text{yy} } & \Delta ^ { \dagger } \int _ { T _ { o } } \left | \left | \text{y} ( t ) \right | \right | ^ { 2 } d t \quad & \quad ( 3 9 ) \\ \text{e} _ { \text{Qy} } & \Delta ^ { \dagger } \int _ { T _ { o } } \mathbf Q ^ { H } ( t, \theta ) \text{y} ( t ) d t \quad & \quad ( 4 0 )$$

$$\mu \times & \quad J _ { T _ { o } } \stackrel { \sim } { \sim } \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \partial \, \alpha \, \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \alpha \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \beta \, \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \. \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \Delta \, \epsilon \, \, \,$$

Using the definition of Q ( t, θ ) in equation (26), we can write the f th element of e Qy as

$$[ e _ { Q y } ] _ { f } = \mathbf a ^ { H } ( \theta _ { f } ) \mathbf E ^ { T } \mathbf a ( \theta _ { f } ) \, \text{ \quad \ \ } ( 4 2 )$$

where

$$\mathbf E = \int _ { T _ { o } } y ( t ) x ^ { H } ( t ) d t. \quad \quad ( 4 3 ) \\ \intertext { a can writ a that a } \text{$\mathbf e_{\mathbf n}$ alamant of $\mathbf E_{\mathbf n}$} \text{$\mathbf e_{\mathbf n}$}$$

Similarly, we can write the fg th element of E QQ as

$$[ \mathbf E _ { \mathbf Q Q } ] _ { f g } = \mathbf a ^ { H } ( \theta _ { f } ) \mathbf a ( \theta _ { g } ) \mathbf a ^ { H } ( \theta _ { f } ) \mathbf R _ { x } ^ { T } \mathbf a ( \theta _ { g } ). \quad ( 4 4 )$$

Since, e Qy and E QQ are independent of the received signal, the sufficient statistic to calculate θ and α is given by E . Using equation (42)-(44) we can write the ML estimate in matrixvector form as

Then, the GLRT for our hypothesis testing model in equation (24) is given as

$$\text{vector form as} \\ L _ { y } ( \hat { \theta } _ { M L } ) = & \arg \max _ { \theta } \frac { \left | \mathbf a ^ { H } ( \hat { \theta } _ { M L } ) \mathbf E \mathbf a ^ { * } ( \hat { \theta } _ { M L } ) \right | ^ { 2 } } { M \mathbf a ^ { H } ( \hat { \theta } _ { M L } ) \mathbf R _ { x } ^ { T } \mathbf a ( \hat { \theta } _ { M L } ) } \cdot \quad ( 4 5 ) \\ \text{Then, the GLRT for our hypothesis testing model in equal}$$

$$L _ { \mathbf y } = \max _ { \theta, \alpha } \{ \log f _ { \mathbf y } ( \mathbf y, \theta, \alpha ; \mathcal { H } _ { 1 } ) \} - \log f ( \mathbf y ; \mathcal { H } _ { 0 } ) \overset { \mathcal { H } _ { 1 } } { \gtrless } \delta _ { \mathcal { H } _ { 0 } } ( 4 6 )$$

where f y ( y , θ, α ; H 1 ) and f ( y ; H 0 ) are the probability density functions of the received signal under hypothesis H 1 and H 0 , respectively. Hence, the GLRT can be expressed as

The asymptotic statistics of L θ ( ˆ ML ) for both the hypothesis is given by [32]

$$\swarrow \nu _ { 0 }, & \text{respectively. Hence, the GLKI can be expressed as} \\ L _ { \mathbf y } ( \hat { \theta } _ { \text{ML} } ) = & \arg \max _ { \theta } \frac { \left | \mathbf a ^ { H } ( \hat { \theta } _ { \text{ML} } ) \mathbf E \mathbf a ^ { * } ( \hat { \theta } _ { \text{ML} } ) \right | ^ { 2 } } { M \mathbf a ^ { H } ( \hat { \theta } _ { \text{ML} } ) \mathbf R _ { \mathbf x } ^ { T } \mathbf a ( \hat { \theta } _ { \text{ML} } ) } } \stackrel { \mathfrak { G } _ { \mathbf i } } { \gtrless } \delta. \quad ( 4 7 ) \\ & \text{The asymptotic statistics of } L ( \hat { \theta } _ { \text{ML} } ) \text{ for both the hypothesis}$$

$$\mathbf y [ 3 2 ] \\ L ( \hat { \theta } _ { \text{ML} } ) & \sim \begin{cases} \mathcal { H } _ { 1 } \colon \chi _ { 2 } ^ { 2 } ( \rho ), & \\ \mathcal { H } _ { 0 } \colon \chi _ { 2 } ^ { 2 }, & \end{cases}$$

where

- · χ 2 2 ( ρ ) is the noncentral chi-squared distributions with two degrees of freedom,
- · χ 2 2 is the central chi-squared distributions with two degrees of freedom,
- · and ρ is the noncentrality parameter, which is given by

$$\rho = \frac { | \alpha | ^ { 2 } } { \sigma _ { n } ^ { 2 } } | \mathbf a ^ { H } ( \theta ) \mathbf R _ { \mathbf x } ^ { T } \mathbf a ( \theta ) | ^ { 2 }. \text{ \quad \ \ } ( 4 9 )$$

For the general signal model, we set δ according to a desired probability of false alarm P FA, i.e.,

$$P _ { \text{FA} } = P ( L ( \mathbf y ) > \delta | \mathcal { H } _ { 0 } ) \text{ \quad \ \ } ( 5 0 )$$

$$\delta = \mathcal { F } _ { \chi _ { 2 } ^ { 2 } } ^ { - 1 } ( 1 - P _ { \text{FA} } )$$

where F -1 χ 2 2 is the inverse central chi-squared distribution function with two degrees of freedom. The probability of detection is given by

$$P _ { D } = P ( L ( \mathbf y ) > \delta | \mathcal { H } _ { 1 } ) & & ( 5 2 ) \quad \text{using} \\ \quad \text{of } M \text{.}$$

where F χ 2 2 ( ρ ) is the noncentral chi-squared distribution function with two degrees of freedom and noncentrality parameter ρ .

$$P _ { D } & = P ( L ( y ) > \delta | \mathcal { H } _ { 1 } ) & ( 5 2 ) & \text{                                                                                                                                                                                                        
\text{$P_{D} = 1 - \mathcal { F }_{\chi_{2}(\rho)}\left(\mathcal{F}_{\chi_{2}^{-1}(1-P_{FAP})\right)$} } & ( 5 3 ) & \text{                                                                                                                                                                                                  
\text{e $\mathcal{F}_{\chi_{2}(\rho)}$ is the noncentral chi-squared distribution func-} \quad \text{$null$ s}$$

## A. P D for Orthogonal Waveforms

For orthogonal waveforms R T x = I M , therefore, the GLRT can be expressed as

$$& \text{can be expressed as} \\ & L _ { \text{Orthog} } ( \hat { \theta } _ { \text{ML} } ) = \frac { \left | \mathfrak { a } ^ { H } ( \hat { \theta } _ { \text{ML} } ) \mathbf E \mathbf a ^ { * } ( \hat { \theta } _ { \text{ML} } ) \right | ^ { 2 } } { M \mathbf a ^ { H } ( \hat { \theta } _ { \text{ML} } ) \mathbf a ( \hat { \theta } _ { \text{ML} } ) } \stackrel { \mathfrak { H } _ { 1 } } { \gtrless } \delta _ { \text{Orthog} } \quad ( 5 4 ) \quad \frac { \text{Parame} } { \text{Radar} / \mathfrak { C } } \\ & \text{and the statistics of } L ( \hat { \theta } _ { \text{ML} } ) \text{ for this case is}$$

$$\ t n e \, \text{statistics of } L ( \vartheta _ { M L } ) \, \text{for } t n i s \, \text{case is } & & \text{computer} \\ L _ { \text{Orthog} } ( \hat { \theta } _ { M L } ) \sim & \begin{cases} \mathcal { H } _ { 1 } \colon \chi _ { 2 } ^ { 2 } ( \rho _ { \text{Orthog} } ), & \text{Comr} \\ \mathcal { H } _ { 0 } \colon \chi _ { 2 } ^ { 2 }, & \text{Wave} \\ & & \text{Wave} \end{cases}$$

and the statistics of L θ ( ˆ ML ) for this case is where

$$\rho _ { \text{Orthog} } = \frac { M ^ { 2 } | \alpha | ^ { 2 } } { \sigma _ { n } ^ { 2 } } \cdot \quad \quad \quad ( 5 6 ) \quad \text{Rac} \\ \text{Spec}$$

We set δ Orthog according to a desired probability of false alarm P PF-Orthog , i.e.,

$$\delta _ { \text{Orthog} } = \mathcal { F } _ { \chi _ { 2 } ^ { 2 } } ^ { - 1 } ( 1 - P _ { \text{PF-Orthog} } ) \quad \quad ( 5 7 ) \quad \text{Dopl}$$

and then the probability of detection for orthogonal waveforms is given by

$$\text{is given by} & \quad \text{$\frac{1}{2}$} \\ P _ { \text{D-Orthog} } = 1 - \mathcal { I } _ { \chi _ { 2 } ^ { 2 } ( \rho _ { \text{Orthog} } ) } \left ( \mathcal { I } _ { \chi _ { 2 } ^ { 2 } } ^ { - 1 } ( 1 - P _ { \text{PF-Orthog} } ) \right ). \quad ( 5 8 ) & \quad \text{$A.\ P erf} \\ n = \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi                                                                                                                                                                                                        
\text{$n$} \quad \text{$\Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde {                                                                                                                                                                                                        \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Psi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Psi } \tilde {                                                                                                                                                                                                 
\text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde {                         \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Psi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Psi } \quad \text{$\Phi } \quad \text{$\Phi } \ \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde{ \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde{ \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde {                                                                                               
\text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad\text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \ \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \quad \text{$\Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Psi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde {\Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Phi } \tilde { \Ph$$

## B. P D for NSP Waveforms

For spectrum sharing waveforms R T x = R T ˘ x , therefore, the GLRT can be expressed as

$$\text{GLK1 can be expressed as} \\ L _ { N S P } ( \hat { \theta } _ { \text{ML} } ) = \frac { \left | \mathfrak { a } ^ { H } ( \hat { \theta } _ { \text{ML} } ) \text{Ea} ^ { * } ( \hat { \$$

$$\text{ne status or } L ( \vartheta _ { M L } ) \text{ for this case is } & \quad \text{In} \\ L _ { N S P } ( \hat { \theta } _ { M L } ) \sim & \begin{cases} \mathcal { H } _ { 1 } \colon \chi _ { 2 } ^$$

and the statistics of L θ ( ˆ ML ) for this case is where

$$\rho _ { N S P } = \frac { | \alpha | ^ { 2 } } { \sigma _ { n } ^ { 2 } } | \mathbf a ^ { H } ( \theta ) \mathbf R _ { \mathbf x } ^ { T } \mathbf a ( \theta ) | ^ { 2 }. \quad \ \ ( 6 1$$

We set δ NSP according to a desired probability of false alarm P PF-NSP, i.e.,

$$\cdots, \quad \delta _ { N S P } = \mathcal { F } _ { \chi _ { 2 } ^ { 2 } } ^ { - 1 } ( 1 - P _ { P F \cdot N S P } ) \quad \quad ( 6 2 ) \quad \text{add} \quad \text{and} \, \, \,$$

and then the probability of detection for orthogonal waveforms is given by

$$\text{is given by} \\ P _ { \text{D-NSP} } = 1 - \mathcal { F } _ { \chi _ { 2 } ^ { 2 } ( \rho _ { \text{NSP} } ) } \left ( \mathcal { F } _ { \chi _ { 2 } ^ { 2 } } ^ { - 1 } ( 1 -$$

## V. NUMERICAL RESULTS

In order to study the detection performance of spectrum sharing MIMO radars, we carry out Monte Carlo simulation using the radar parameters mentioned in Table II. At each run of Monte Carlo simulation we generate K Rayleigh interference channels each with dimensions N BS × M , calculate their null spaces and construct corresponding projection matrices using Algorithm (2), determine the best channel to perform projection of radar signal using Algorithm (1), transmit NSP signal, estimate parameters θ and α from the received signal, and calculate the probability of detection for orthogonal and NSP waveforms.

TABLE II MIMO RADAR SYSTEM PARAMETERS

| Parameters                         | Notations   | Values          |
|------------------------------------|-------------|-----------------|
| Radar/Communication System RF band | -           | 3550 - 3650 MHz |
| Radar antennas                     | M           | 8, 4            |
| Communication System Antennas      | N BS        | 2               |
| Carrier frequency                  | f c         | 3.55 GHz        |
| Wavelength                         | λ           | 8.5 cm          |
| Inter-element antenna spacing      | 3 λ/ 4      | 6.42 cm         |
| Radial velocity                    | v r         | 2000 m/s        |
| Speed of light                     | c           | 3 × 10 8 m/s    |
| Target distance from the radar     | r 0         | 500 Km          |
| Target angle                       | θ           | ˆ θ             |
| Doppler angular frequency          | ω D         | 2 ω c v r /c    |
| Two way propagation delay          | τ r         | 2 r 0 /c        |
| Path loss                          | α           | ˆ α             |

## A. Performance of Algorithms (1) and (2)

In Figure 3, we demonstrate the use of Algorithms (1) and (2) in improving target detection performance when multiple BSs are present in detection space of radar and the radar has to reliably detect target while not interfering with communication system of interest. As an example, we consider a scenario with five BSs and the radar has to select a projection channel which minimizes degradation in its waveform, thus, maximizing its probability of detection of the target.

In Figure 3(a), we consider the case when N BS &lt; M . We show detection results for five different NSP signals, i.e, radar waveform projected onto five different BSs. Note that, in order to achieve a detection probability of 90%, we need 6 dB to 13 dB more gain in SNR as compared to the orthogonal waveform, depending upon which channel we select. Using Algorithms (1) and (2) we can select interference channel that results in minimum degradation of radar waveform and results in enhanced target detection performance with the minimum additional gain in SNR required. For example, Algorithms (1) and (2) would select BS#5 because in this case NSP waveform requires least gain in SNR to achieve a detection probability of 90% as compared to other BSs.

In Figure 3(b), we consider the case when N BS /lessmuch M . Similar to Figure 3(a) we show detection results for five

different NSP signals but now MIMO radar has a larger array of antennas as compared to the previous case. In this case, in order to achieve a detection probability of 90%, we need 3 dB to 5 dB more gain in SNR as compared to the orthogonal waveform. As in the previous case, using Algorithms (1) and (2) we can select interference channel that results in minimum degradation of radar waveform and results in enhanced target detection performance with the minimum additional gain in SNR required. For example, Algorithms (1) and (2) would select BS#2 because in this case NSP waveform requires least gain in SNR to achieve a detection probability of 90% as compared to the other BSs.

The above two examples demonstrates the importance of Algorithms (1) and (2) in selecting interference channel for radar signal projection to maximize detection probability and minimize gain in SNR required as a result of NSP of radar waveforms for spectrum sharing.

## B. Case 1: dim N [ H i ] = 2

In Figure 4, we plot the variations of probability of detection P D as a function of signal-to-noise ratio (SNR) for various values of probability of false alarm P FA. Each subplot represents the P D for a fixed P FA. We choose to evaluate P D against P FA values of 10 -1 , 10 -3 , 10 -5 and 10 -7 when the interference channel H i has dimensions 2 × 4 , i.e., the radar has M = 4 antennas and the communication system has N BS = 2 antennas, thus, we have a null space dimension of 'dim N [ H i ] = 2 '. When we compare the detection performance of two waveforms we note that in order to get a desired P D for a fixed P FA we need more SNR for NSP than orthogonal waveforms. For example, say we desire P D = 0 9 . , then according to Figure 4 we need 6 dB more gain in SNR for NSP waveform to get the same result produced by the orthogonal waveform.

## C. Case 2: dim N [ H i ] = 6

In Figure 5, similar to Figure 4, we do an analysis of P D against the same values of P FA but for interference channel H i having dimensions 2 × 8 , i.e., now the radar has M = 8 antennas and the communication system has N = 2 antennas, thus, we have a null space dimension of 'dim N [ H i ] = 6 '. Similar to Case 1, when we compare the detection performance of two waveforms we note that in order to get a desired P D for a fixed P FA we need more SNR for NSP than the orthogonal waveforms. For example, say we desire P D = 0 9 . , then according to Figure 5 we need 3.5 to 4.5 dB more gain in SNR for the NSP waveform to get the same result produced by the orthogonal waveform.

## D. Comparison of Case 1 and Case 2

As expected, when SNR increases detection performance increases for both waveforms. However, when we compare the two waveforms at a fixed value of SNR, the orthogonal waveforms perform much better than the NSP waveform in detecting target. This is because our transmitted waveforms are no longer orthogonal and we lose the advantages promised

(a) Probability of detection when N BS &lt; M . As an example we use 2 × 4 configuration. Note that 6 dB to 13 dB of additional gain in SNR is required to detect target with 90% probability, depending upon the NSP waveform transmitted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0540v2-figure-04938dca-6b47-4708-b1bc-7eb7d7d9df61.jpg)

(b) Probability of detection when N BS /lessmuch M . As an example we use 2 × 8 configuration. Note that 3 dB to 5 dB of additional gain in SNR is required to detect target with 90% probability, depending upon the NSP waveform transmitted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0540v2-figure-a97d3986-81d8-4282-9752-8f3f00ef024b.jpg)

Fig. 3. Performance of Algorithms (1) and (2): Using our spectrum sharing and projection algorithms, we can select interference channel for radar signal projection to maximize detection probability and minimize gain in SNR required as a result of NSP of radar waveforms. For example, Algorithms (1) and (2) select BS#5 and BS#2 for N BS &lt; M and N BS &lt; M cases, respectively, as they require minimum additional gain in SNR.

by orthogonal waveforms when used in MIMO radars as discussed in Section II-F, but, we ensure zero interference to the BS of interest, thus, sharing radar spectrum at an increased cost of target detection in terms of SNR.

In Case 1, in order to achieve a desired P D for a fixed P FA we need more SNR for NSP as compared to Case 2. This is because we are using more radar antennas, while the antennas at the BS remain fixed, in Case 2 which increases the dimension of the null space of the interference channel. This yields better detection performance even for NSP waveform. So, in order to mitigate the effect of NSP on radar performance one way is to employ a larger array at the radar transmitter.

Fig. 4. 'Case 1: dim N [ H i ] = 2 ': P D as a function of SNR for various values of probability of false alarm P FA, i.e., P FA = 10 -1 , 10 -3 , 10 -5 and 10 -7 . The interference channel H i has dimensions 2 × 4 , i.e., the radar has M = 4 antennas and the communication system has N BS = 2 antennas, thus, we have a null space dimension of 'dim N [ H i ] = 2 '. Note that we need 9 to 10 dB more gain in SNR for the NSP waveform to get the same result produced by the orthogonal waveform.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0540v2-figure-6860ed7b-33c3-4c91-ac13-3fc3313ffd5b.jpg)

Fig. 5. 'Case 2: dim N [ H i ] = 6 ': P D as a function of SNR for various values of probability of false alarm P FA, i.e., P FA = 10 -1 , 10 -3 , 10 -5 and 10 -7 . The interference channel H i has dimensions 2 × 8 , i.e., the radar has M = 8 antennas and the communication system has N = 2 antennas, thus, we have a null space dimension of 'dim N [ H i ] = 6 '. Note that we need 3.5 to 4.5 dB more gain in SNR for the NSP waveform to get the same result produced by the orthogonal waveform.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0540v2-figure-c0815acd-20fd-4e4a-953e-3678417ba719.jpg)

## VI. CONCLUSION

In the future, radar RF spectrum will be shared with wireless communication systems to meet the growing bandwidth demands and mitigate the effects of spectrum congestion for commercial wireless services. In this paper, we analyzed a similar spectrum sharing scenario between radars and cellular systems. We proposed a spectrum sharing scenario in which a MIMO radar is sharing spectrum with a cellular system. We proposed algorithms for interference channel selection and projection of radar waveform onto the selected interference channel in order to mitigate interference to the selected BS. We evaluated the detection performance of spectrum sharing MIMO radars. We formulated the statistical detection problem for target detection and used generalized likelihood ratio test to decide about the presence of target when using orthogonal waveforms and null-space projected (NSP) waveforms. We showed that by using our spectrum sharing and projection algorithms the radar can maximize target detection probability and minimize additional gain in SNR required to detect the target. We showed that when using the NSP waveforms, the detection performance degrades as compared to the orthogonal waveforms and we need more SNR to detect reliably. Our results showed that, about 6 dB of gain in SNR is required when N BS &lt; M and 3.5 to 4.5 dB of gain in SNR is required when N BS /lessmuch M when NSP waveforms are used instead of orthogonal waveforms for spectrum sharing. Our analysis showed that this degradation in performance can be mitigated by using a larger array at the MIMO radar transmitter.

## REFERENCES

- [1] S. Haykin, 'Cognitive radio: Brain-empowered wireless communications,' IEEE Journal on Selected Areas in Communications , vol. 23, pp. 201-220, Feb 2005.
- [2] T. Yucek and H. Arslan, 'A survey of spectrum sensing algorithms for cognitive radio applications,' IEEE Communications Surveys Tutorials , vol. 11, pp. 116-130, First 2009.
- [3] R. Murty, R. Chandra, T. Moscibroda, and P. V. Bahl, 'Senseless: A database-driven white spaces network,' IEEE Transactions on Mobile Computing , vol. 11, pp. 189-203, Feb. 2012.
- [4] Y. Zhao, L. Morales, J. Gaeddert, K. Bae, J.-S. Um, and J. Reed, 'Applying radio environment maps to cognitive wireless regional area networks,' in 2nd IEEE International Symposium on New Frontiers in Dynamic Spectrum Access Networks (DySPAN) , pp. 115-118, April 2007.
- [5] B. Gao, J. Park, and Y. Yang, 'Uplink soft frequency reuse for selfcoexistence of cognitive radio networks,' IEEE Transactions on Mobile Computing , vol. 13, pp. 1366-1378, June 2014.
- [6] Federal Communications Commission (FCC), 'FCC proposes innovative small cell use in 3.5 GHz band.' Online: http://www.fcc.gov/document/fcc-proposes-innovative-small-celluse-35-ghz-band, December 12, 2012.
- [7] Federal Communications Commission (FCC), 'Connecting America: The national broadband plan.' Online, 2010.
- [8] The Presidents Council of Advisors on Science and Technology (PCAST), 'Realizing the full potential of government-held spectrum to spur economic growth,' July 2012.
- [9] National Telecommunications and Information Administration (NTIA), 'An assessment of the near-term viability of accommodating wireless broadband systems in the 1675-1710 MHz, 1755-1780 MHz, 3500-3650 MHz, 4200-4220 MHz, and 4380-4400 MHz bands (Fast Track Report).' Online, October 2010.
- [10] Federal Communications Commission (FCC), 'In the matter of revision of parts 2 and 15 of the commissions rules to permit unlicensed national information infrastructure (U-NII) devices in the 5 GHz band.' MO&amp;O, ET Docket No. 03-122, June 2006.
- [11] L. S. Wang, J. P. McGeehan, C. Williams, and A. Doufexi, 'Application of cooperative sensing in radar-communications coexistence,' IET Communications , vol. 2, pp. 856-868, 2008.
- [12] S. S. Bhat, R. M. Narayanan, and M. Rangaswamy, 'Bandwidth sharing and scheduling for multimodal radar with communications and tracking,' in IEEE Sensor Array and Multichannel Signal Processing Workshop , p. 233236, 2012.
- [13] R. Saruthirathanaworakun, J. Peha, and L. Correia, 'Performance of data services in cellular networks sharing spectrum with a single rotating radar,' in IEEE International Symposium on a World of Wireless, Mobile and Multimedia Networks (WoWMoM) , pp. 1-6, 2012.
- [14] C. Rossler, E. Ertin, and R. Moses, 'A software defined radar system for joint communication and sensing,' in IEEE Radar Conference (RADAR) , pp. 1050-1055, May 2011.
- [15] R. Y. X. Li, Z. Zhang, and W. Cheng, 'Research of constructing method of complete complementary sequence in integrated radar and communication,' in IEEE Conference on Signal Processing , vol. 3, p. 17291732, 2012.
- [16] C. Sturm and W. Wiesbeck, 'Waveform design and signal processing aspects for fusion of wireless communications and radar sensing,' Proceedings of the IEEE , vol. 99, pp. 1236-1259, July 2011.
- [17] M. P. Fitz, T. R. Halford, I. Hossain, and S. W. Enserink, 'Towards simultaneous radar and spectral sensing,' in IEEE International Symposium on Dynamic Spectrum Access Networks (DYSPAN) , pp. 15-19, April 2014.
- [18] X. Chen, X. Wang, S. Xu, and J. Zhang, 'A novel radar waveform compatible with communication,' in International Conference on Computational Problem-Solving (ICCP) , p. 177181, 2011.
- [19] A. Khawar, A. Abdel-Hadi, and T. C. Clancy, 'Spectrum sharing between S-band radar and LTE cellular system: A spatial approach,' in 2014 IEEE International Symposium on Dynamic Spectrum Access Networks: SSPARC Workshop (IEEE DySPAN 2014 - SSPARC Workshop) , (McLean, USA), Apr. 2014.
- [20] A. Khawar, A. Abdel-Hadi, and T. C. Clancy, 'MIMO radar waveform design for coexistence with cellular systems,' in IEEE International Symposium on Dynamic Spectrum Access Networks: SSPARC Workshop (IEEE DySPAN 2014 - SSPARC Workshop) , (McLean, USA), Apr. 2014.
- [21] M. Ghorbanzadeh, A. Abdelhadi, and C. Clancy, 'A utility proportional fairness resource allocation in spectrally radar-coexistent cellular networks,' in Military Communications Conference (MILCOM) , 2014.
- [22] A. Khawar, A. Abdelhadi, and T. C. Clancy, 'On The Impact of TimeVarying Interference-Channel on the Spatial Approach of Spectrum Sharing between S-band Radar and Communication System,' in Military Communications Conference (MILCOM) , 2014.
- [23] H. Shajaiah, A. Khawar, A. Abdel-Hadi, and T. C. Clancy, 'Using resource allocation with carrier aggregation for spectrum sharing between radar and 4G-LTE cellular system,' in IEEE DySPAN , 2014.
- [24] F. Paisana, J. P. Miranda, N. Marchetti, and L. A. DaSilva, 'Databaseaided sensing for radar bands,' in IEEE International Symposium on Dynamic Spectrum Access Networks (DYSPAN) , pp. 1-6, April 2014.
- [25] H. Deng and B. Himed, 'Interference mitigation processing for spectrum-sharing between radar and wireless communications systems,' IEEE Transactions on Aerospace and Electronic Systems , vol. 49, no. 3, pp. 1911-1919, 2013.
- [26] A. M. Haimovich, R. S. Blum, and L. J. Cimini, 'MIMO radar with widely separated antennas,' IEEE Signal Processing Magazine , vol. 25, no. 1, pp. 116-129, 2008.
- [27] J. Li and P. Stoica, 'MIMO radar with colocated antennas,' IEEE Signal Processing Magazine , vol. 24, no. 5, pp. 106-114, 2007.
- [28] M. Skolnik, Radar Handbook . McGraw-Hill Professional, 3 rd ed., 2008.
- [29] J. Li and P. Stoica, MIMO Radar Signal Processing . Wiley-IEEE Press, 2008.
- [30] D. Tse and P. Viswanath, Fundamentals of Wireless Communication . Cambridge University Press, 2005.
- [31] I. Bekkerman and J. Tabrikian, 'Target detection and localization using mimo radars and sonars,' IEEE Transactions on Signal Processing , vol. 54, pp. 3873-3883, Oct 2006.
- [32] S. Kay, Fundamentals of Statistical Signal Processing: Detection Theory . Prentice Hall, 1998.
Awais Khawar Ahmed Abdelhadi T. Charles Clancy
Target Detection Performance of Spectrum Sharing MIMO Radars
Future wireless communication systems are envisioned to share radio frequency
(RF) spectrum, with other services such as radars, in order to meet the growing
spectrum demands. In this paper, we consider co-channel spectrum sharing
between cellular systems and radars. We address the problem of target detection
by radars that are subject to shape its waveform in a way that it does not
cause interference to cellular systems. We consider a multiple-input
multiple-output (MIMO) radar and a MIMO cellular communication system with $\mc
K$ base stations (BS). We propose a spectrum sharing algorithm which steers
radar nulls, by projecting radar waveform onto the null space of interference
channel, towards a `selected' BS, thus, protecting it from radar interference.
This BS is selected, among $\mc K$ BSs, on the basis of guaranteeing minimum
waveform degradation. We study target detection capabilities of this null-space
projected (NSP) waveform and compare it with the orthogonal waveform. We derive
the generalized likelihood ratio test (GLRT) for target detection and derive
detector statistic for NSP and orthogonal waveform. The target detection
performance for NSP and orthogonal waveform is studied theoretically and via
Monte Carlo simulations.