
(element-wise) and Kronecker products are and ⊗, respectively.
The vectorized form of a matrix Bis written as vec (B). The
s-dimensional all-ones vector, all-zeros vector, and the identity ma-
trix of size s×sare 1s,0N, and Is, respectively. The minimum
eigenvalue of Bis denoted by λmin(B). The real, imaginary, and
angle/phase components of a complex number are Re (·),Im (·),
and arg (·), respectively. vec−1
K,L (c)reshapes the input vector
c∈CKL into a matrix C∈CK×Lsuch that vec (C) = c.
2. MULTI-IRS-AIDED RADAR SYSTEM MODEL
Consider a colocated MIMO radar with Nttransmit and Nrreceive
antennas, each arranged as uniform arrays (ULA) with inter-element
spacing d. The MIRS platforms indexed as IRS1, IRS2,...,IRSM,
are implemented at stationary and known locations, each equipped
with Nmreflecting elements arranged as ULA, with element spac-
ing of dmbetween the antennas/reflecting elements of IRSm. The
continuous-time signal transmitted from the n-th antenna at time
instant tis xn(t). Denote the Nt×1vector of all transmit sig-
nals as x(t) = [xi(t)]Nt
i=1 ∈ΩNt, where the set of unimodular
sequences is Ωn=s∈Cn|s= [ejωi]n
i=1, ωi∈[0,2π]. The
steering vectors of radar transmitter, receiver and the m-th IRS
are, respectively, at(θ) = [1, ej2π
λdsinθ,...,ej2π
λd(Nt−1)sinθ]>,
ar(θ) = [1, ej2π
λdsinθ,...,ej2π
λd(Nr−1)sinθ]>, and bm(θ) =
[1, ej2π
λdmsinθ,...,ej2π
λdm(Nm−1)sinθ]>, where λ, is the carrier
wavelength and dand dmare usually assumed to be half the carrier
wavelength. Each reflecting element of IRSmreflects the incident
signal with a phase shift and amplitude change that is configured via
a smart controller [28]. We denote the phase shift vector of IRSm
by vm= [ejφm,1,...,ejφm,Nm]>∈CNm, where φm,k ∈[0,2π]
is the phase shift associated with the k-th passive element of IRSm.
Denote the angle between the radar-target, radar–IRSm, and
target-IRSmby θtr,θri,m , and θti,m, respectively. Denote target-
radar channel by Htr =ar(θtr)∈CNr×1; and radar-target by
Hrt =at(θtr )>∈C1×Nt. The LoS or radar-target-radar channel
matrix is Hrtr =ar(θtr )at(θtr )>∈CNr×NT. Analogously, for
the multi-IRS aided radar the NLoS channel matrices associated with
IRSmare defined as Hri,m =bm(θri,m)a>
t(θri,m)∈CNm×NT
for radar-IRSm;Hit,m =b>
m(θti,m)∈C1×Nmfor IRSm-
target; Hti,m =bm(θti,m)∈CNm×1for target-IRSm; and
Hir,m =ar(θri,m)b>
m(θri,m)∈CNr×Nmfor IRSm-radar paths.
The received signal back-scattered from a single target is the
superimposition of echoes from both LoS and NLoS paths as
y(t) = αrtr Hrtrx(t−τrtr)
+
M
X
m=1
αritr,m HtrHit,m ΦmHri,mx(t−τritr,m)
+
M
X
m=1
αrtir,m Hir,mΦmHti,mHrtx(t−τrtir,m)
+
M
X
m=1
αritir,m Hir,mΦmHti,mHit,mΦmHri,m
x(t−τritir,m) + w(t),∈CNr,(1)
where Φm= Diag (vm),α(·),m is the complex reflectivity which
depends on the target back-scattering coefficient and the atmospheric
attenuation, and w(t)∼ CN (0, σ2INt)denotes a stationary (ho-
moscedastic) additive white Gaussian noise (AWGN). In general,
the received signal may also have an additional inter-IRS interfer-
ence that should be included while accounting for the SNR. When
there is some blockage or obstruction between the radar and target,
we have αrtr '0,αritr,m '0and αrtir,m '0. We replace
αritir,m by αmfor notation brevity. The received signal becomes
y(t) =
M
X
m=1
αmHir,mΦmHti,mHit,mΦmHri,m
x(t−τritir,m) + w(t).(2)
Our goal is to design a radar system for inspecting a range cell
located at distance dtr with respect to (w.r.t.) the radar transmit-
ter/receiver for a potential target. Assume that the relative time
gaps between any two multipath signals are very small in compar-
ison to the actual roundtrip delays, i.e., τritir,m ≈τ0=2dtr
cfor
m∈ {1,...,M}, where cis the speed of light. We collect N
slow-time samples at the rate 1/Tsfrom the signal, at t=nTs,
n= 0,...,N −1. Hence, corresponding to the range-cell of inter-
est, the received signal vector is
y[n] =
M
X
m=1
αmHmx[n] + w[n],y[n]∈CNr×1,(3)
where x[n] = x(τ0+nTs)∈CNt×1,y[n]=[yi[n]]Nr
i=1, and we
define Hm=Hir,mΦmHti,mHit,mΦmHri,m ∈CNr×Nt. The
delay τ0is aligned on-the-grid so that n0=τ0/Tsis an integer [29].
Collecting all discrete-time samples for Nrreceiver antennas,
the received signal is the Nr×NmatrixY= [y[0],...,y[N−
1]] = PM
m=1 αmHmX+W,where X= [x[0],...,x[N−1]] ∈
CNt×Nand W= [w[0],...,w[N−1]] ∈CNr×N. Vectorizing as
y= vec (Y)yields
y=
M
X
m=1
αmvec (HmX) + vec (W) = ˜
X˜
Hα+˜
w,(4)
where ˜
X=X>⊗INr,˜
H= [ ˜
H1,..., ˜
HM],˜
Hm= vec (Hm)for
m∈ {1,...,M},˜
w= vec (W)and α= [αm]M
m=1. Given that
w(t)is AWGN in (1), it is easily observed that y∼ CN (µ,R),
where µ=˜
X˜
Hαand R=σ2INrN. Note that, since w(n)is a
stationary process and i.i.d. with σ2variance, through vectorization
and stacking all ensembles as one vector, the resulting process is still
stationary and i.i.d with the same variance.
Our goal is to show the effectiveness of placing MIRS plat-
forms in estimating the DoA of the target in the LoS path, i.e. θtr .
For simplicity, we consider a two-dimensional (2-D) scenario, where
the radar, IRS platforms and the target are in the same plane. Our
analysis can be easily extended to 3-D scenarios. The following re-
mark states that the estimation of DoAs in the NLoS paths, θti,m ,
for m∈ {1,...,M}is equivalent to an estimation of θtr.
Remark 1. Estimation of the vector of target DoAs, ζ= [θti,1,
...,θti,M ]>is equivalent to estimating scalar DoA parameter, θtr .
This follows because, given the locations of radar, IRS platforms
and potential target range in the 2-D plane, we have ζ= [θtr +
θ1,...,θtr +θM]>, where θmfor m∈ {1,...,M}are known.
3. UQP-BASED CRLB OPTIMIZATION
For an unbiased estimator of a parameter θtr (θ, hereafter), the vari-
ance of ˆ
θis lower bounded as E{(ˆ
θ−θ)(ˆ
θ−θ)H} ≥ CRLB(θ)[30].