($\left\lceil x\right\rceil$) denotes the largest (smallest) integer not greater (smaller) than $x$. ${\left\langle\cdot\right\rangle_{N}}$ denotes the modulo-$N$ operation. $\delta(\cdot)$ denotes the delta function. Upper (lower) case boldface letters denote matrices (column vectors). The notation $\triangleq$ is used for definitions. Notations $\sim$ and $\propto$ represent “distributed as” and “proportional to”, respectively. We adopt ${{\mathbf{I}}}_{N}$ to denote the $N\times N$ dimensional identity matrix, and ${\mathbf{I}}_{N\times G}$ to denote the matrix composed of the first $G\left(\leq N\right)$ columns of ${\mathbf{I}}_{N}$. We adopt ${\mathbf{0}}$ to denote the all-zero vector or matrix. The superscripts $(\cdot)^{H}$, $(\cdot)^{T}$, and $(\cdot)^{*}$ denote the conjugate-transpose, transpose, and conjugate operations, respectively. The operator ${\mathsf{diag}\left\{{\mathbf{x}}\right\}}$ denotes the diagonal matrix with ${\mathbf{x}}$ along its main diagonal. We employ ${\left[{\mathbf{a}}\right]_{i}}$, ${\left[{\mathbf{A}}\right]_{i,j}}$ and ${\left[{\mathbf{A}}\right]_{i,:}}$ to denote the [$i$th]{} element of the vector ${\mathbf{a}}$, the $(i,j)$th element of the matrix ${\mathbf{A}}$ and the [$i$th]{} row of the matrix ${\mathbf{A}}$, respectively, where the element indices start with 0. ${\mathbb{C}}^{M\times N}$ (${\mathbb{R}}^{M\times N}$) denotes the $M\times N$ dimensional complex (real) vector space. ${{{\mathsf{E}}}\left\{\cdot\right\}}$ denotes the expectation operation. ${\mathcal{CN}\left( {\mathbf{a}}, {\mathbf{B}}\right) }$ denotes the circular symmetric complex Gaussian distribution with mean ${\mathbf{a}}$ and covariance ${\mathbf{B}}$. $\otimes$ and $\odot$ denote the Kronecker product and Hadamard product, respectively. ${\mathsf{vec}\left\{\cdot\right\}}$ represents the 