A Waveform Design for MIMO Sensing on Two-Dimensional Arrays with Sparse Estimation

: With the increasing demand for target sensing in future mobile communication networks, integrated sensing and communications (ISAC) not only brings new system architecture and multi-functional integration, but also introduces new challenges. In order to further reduce the data volume of the sensing system signal processing, this paper proposes a waveform optimization design method and target parameter estimation technique for a compressed sensing multiple-input multiple-output (MIMO) radar based on a 2D antenna array. This scheme enables accurate estimation of target angles in the sparse time-domain sampling scenario and L-shaped antenna array. In this paper, extensive simulation results demonstrate the good performance of the proposed scheme for various compression ratios and signal-to-noise ratios (SNRs)


Introduction
At present, the devices used for communication and sensing are still in a state of independent separation or coexistence in adjacent frequency bands.Communication and sensing processes use different waveforms, spectra, antennas, processors, and other hardware equipment resources, resulting in low efficiency of overall resource utilization.With the continuous development of 5G technology and the global research layout of 6G technology, in order to meet the demands for high-precision wireless network connections such as the industrial internet of things, smart cities, and remote healthcare in the future, the network's requirements for positioning and sensing capabilities are gradually increasing [1].The integration, forming integrated sensing and communication (ISAC) technology, provides a new system or network architecture, enabling integration to achieve the symbiosis and reciprocity between communication and sensing through reasonable resource reuse design.Therefore, compared to traditional communication networks that only support communication capabilities, the ISAC poses new requirements and challenges for the management technology of radio resources and the design of network architecture [2].
Compared to traditional communication networks, the introduction of radar sensing capability will add signal processing processes for sensing in communication systems, resulting in additional computational resource requirements.For high range resolution radars, sensing requires a greater bandwidth compared to traditional communication.By introducing compressed sensing (CS) technology, sampling can be carried out at a lower sampling rate than Nyquist under the same bandwidth requirements, and the amount of data to be processed can be significantly reduced on the digital signal processing side, which can greatly alleviate the pressure on the hardware and computing power of the ISAC system [3].In addition, multiple-input multiple-output (MIMO) technology can greatly improve the overall performance of systems in both wireless communication and radar sensing fields.The design of MIMO technology for ISAC is still a popular research issue [4].Therefore, in order to meet the continuous evolution of the current 5G communication system for the ISAC and further improve the sensing performance in limited computing resources, this paper will focus on the optimal design of the sensing waveform for the MIMO radar based on the compressed sensing (CS) technology.The new design waveform can fully adapt to the compressed sensing system and reduce the data size for computation in the sparse signal processing process.In order to achieve sparse recovery for limited data with high probability, the sensing matrix, which can used for reconstruct the original signal of CS technology, needs to satisfy the restricted isometric properties (RIP) [5], which means that the sensing matrix should be as orthogonal as possible.Compressed sensing in the context of MIMO radars has been studied in [6,7]; it takes advantage of the sparsity of targets in the target space to perform multi-target detection and estimation through sparse recovery methods.Regarding a CS-based MIMO radar, waveform optimization can be widely achieved from multiple dimensions such as the frequency domain, time domain, spatial domain, and computing power domain.
In the frequency domain, there are several techniques to further improve the detection performance of CS-based MIMO radars.For example, significant range resolution gains can be achieved by using a stepping frequency scheme during transmission [8,9], in which optimal or near-optimal measurement matrices have been designed to improve the detection performance by minimizing the coherence and/or signal inference ratio (SIR) of the measurement matrix.Considering the spectrum policy of mobile communication systems, various operators of the wireless network will be strictly limited in frequency usage for ISAC in future 5G and 6G communication systems, making it difficult to achieve frequency hopping compression from the frequency domain dimension.
In the time domain, optimizing the waveform to improve the performance of the sparse recovery in CS-based MIMO radars has been investigated in several previous works.These optimization methods can compress the radar sampling signal in the time domain while ensuring sensing accuracy of the radar as much as possible, thereby reducing the complexity of the sampling and subsequent processing.To effectively reconstruct the target scene, a waveform design method to reduce the correlations between target responses has been introduced in [6].The waveform design problem considers beam-pattern design and optimizes phase hopping signals with the objective of minimizing coherence.In [10], the power allocation and waveform design problems are proposed for the widely separated and colocated setting using an orthogonal Hadamard sequence as the initialization.In [11], a new iteration algorithm was proposed for a joint design transmission waveform and sensing matrix, which then initialized the transmission waveform as a random complex vector.Study [12] proposed a formal algorithm for the joint design of radar waveforms and receiver filters of CS-based MIMO radar systems for extended target detection, and radar waveforms are designed based on Alltop sequences.It is not difficult to find that in order to explore the principle of waveform optimization under compressed sensing conditions, these studies first selected a relatively simple waveform to intuitively demonstrate the effectiveness of the method, but such ideal waveforms can not be directly used in practical applications.
Then, in the spatial domain, a waveform design approach based on given antenna locations where the design criterion focuses on the coherence of the sensing matrix has been presented in [13], and it also optimized the allocation of the total power budget among transmitters.The study considered a one-dimensional array structure, specifically assuming that transmitters and receivers are randomly distributed within a small area, forming a colocated CS-based MIMO radar system.A similar antenna configuration was used in [10].In [14], an optimal design of frequency hopping waveforms based on a linear array arrangement has been introduced.The waveform design problem can be simplified into a frequency code selection scheme, and then an iterative process was used to minimize the block coherence criterion for selection.Furthermore, study [6] considered the uniform linear array for TX and RX antennas.In order to verify the impact of the designed optimized waveforms on the sparse recovery performance, the above studies are all based on the assumption of one-dimensional arrays such as the uniform line array (ULA) for azimuth estimation.There is indeed little research on extending these methods to actual two-dimensional angle space.In the actual radar sensing process, it is necessary to simultaneously estimate the target's azimuth and elevation angles.
At last, in the computing power domain, a CS-based MIMO radar can realize the estimation of range, velocity, and angle.Compared to the basic functional requirements of range and speed, the estimation of direction of arrival (DOA) in the future wireless network will become more important.Therefore, for the angle estimation method of CS-based MIMO radar, in [15,16], the radar waveforms and receiver filters were jointly optimized with respect to a mutual coherence criterion.The Regularized Orthogonal Matched Basis Pursuit (ROMP) algorithm was used as the CS recovery algorithm in the studies.Moreover, the Orthogonal Matching Pursuit (OMP) algorithm was used in [6,11], and the Dantzig selector was used in [10].After research, it was found that the commonly used sparse recovery methods in previous works were greedy algorithms such as OMP or convex relaxation algorithms, but these methods may not perform well when dealing with highly correlated or high-dimensional signal problems.
To address the above issues, the contributions in this article are as follows: 1.
In the time domain for compressed optimization, compared to ideal orthogonal MIMO waveforms such as the Hadamard sequence and Alltop sequence as the initial waveform, this work will use a Hadamard-coded waveform as the initial waveform which is more suitable for millimeter wave radar; 2.
In the spatial domain for compressed optimization, unlike most previous studies that have considered one-dimensional linear arrays that can only estimate the azimuth information of the target, this paper extends the research to CS-based MIMO radar systems based on L-shaped arrays and redesigned waveform optimizations; 3.
In the computational power domain, due to the use of a two-dimensional array based MIMO radar scenarios in this article, sparse recovery belongs to the problem of high correlation and high-dimensional signals, and algorithms such as OMP will no longer be applicable.Therefore, in order to achieve higher resolution and accuracy, this article introduces the Sparse Bayesian Learning (SBL) algorithm for the sparse estimation of both azimuth and elevation angles.
The rest of the paper is organized as follows.In Section 2, the signal model and the expression of waveform design questions are discussed.In Section 3, the sparse recovery algorithm and optimization techniques have been introduced.The numerical experiments and performance validation of the proposed methods are presented in Section 4. Finally, Sections 5 and 6 contain the discussions and conclusions for this article.

Waveform Design for CS-Based MIMO Radar 2.1. Signal Model
In this paper, we consider a colocated MIMO radar system based on an L-shaped array, as shown in Figure 1.It is assumed that the L-shaped array consists of two mutually vertical uniform line arrays, where the X-axis is the transmitting array with M antennas and the Y-axis is the receiving array with N antennas, the reference antenna element (transmitting and receiving shared antenna) is at the coordinate origin.The transmitting and receiving antenna elements are separated by a distance of d t and d r , respectively.M transmitting antennas transmit orthogonal waveforms with the same bandwidth and center frequency.Lowercase bold letters denote vectors and capital bold letters denote matrices.Superscripts (•) T and (•) H denote, respectively, the transpose and Hermitian transpose.⊗ denotes the Kronecker tensor product.
Let us assume that there are K targets to be estimated in the far field of the antennas, ignoring the clutter.A clutter suppression filter can be used at the receiver if the clutter interference in the environment is large [17].Therefore, the following studies in this paper do not consider clutter.In the general case, the targets are located in different range cells and are moving.Therefore, we would require searching for targets in a three-dimensional (3D) space, i.e., range-speed-DOA space.For simplicity, we assume that all the targets are not moving and located in the same range cell, thus the only parameters that need to be estimated are the target azimuth angles θ k and elevation angles ϕ k , k = 1, . . ., K. We can be positioned in two dimensions by (θ k , ϕ k ), and θ k , ϕ k ∈ (0 • , 180 • ).Although, in the following, we focus on DOA estimation only, the results can be easily extended to the case of moving targets.

𝑥 𝑧
The k-th target To mathematically describe the radar signals, let x m = [x m (1), x m (2), . . ., x m (L)] T denote the transmitted waveform of the m-th transmitting antenna, where L denotes the number of samples in a pulse period.Then, the matrix of the transmitted waveforms corresponding to the transmitting array is denoted by X L×M = [x 1 , x 2 , . . . ,x M ] T , where each column corresponds to a transmitter and rows encode the time instances.Assuming that the angles between the k-th target and the X-axis and Y-axis are ρ k , γ k , respectively, from the geometric relationship in Figure 1, we can know that then, the transmit steering vector a M×1 (ρ k ) corresponding to X-axis and receive steering vector b N×1 (γ k ) corresponding to Y-axis of the k-th target can be given by The received signal (without matched filtering) at the n-th receive antenna can be represented as where X is an L × M matrix that contains the transmit waveforms as its columns, η k is the reflection coefficient (RCS) of the k-th target, and w n stands for the additive noise vector at the n-th receiver, which is modeled by a circularly symmetric complex Gaussian random vector with variance σ 2 .As a result, the target echo signal matrix of the L-shaped array MIMO radar can be expressed as where is the receive steering matrix, diag(η k ) is the diagonal matrix, and W is the additive noise matrix.At the n-th receive antenna, the L × 1 dimensional received signal is linearly compressed via the measurement matrix Φ to obtain a T × 1 dimensional observed signal, which can be expressed as where the size of Φ is T × L. Similarly, we can obtain the matrix form R = ΦU.
In general, the azimuth and elevation angles of each target are real numbers.However, in practice, due to quantization and other resolution-limiting factors, the estimated angles by the radar have finite precision and belong to a finite set.Therefore, we discretize the angle space into a linearly spaced grid with P × Q discrete angles, as shown in Figure 2. Define the azimuth angle space as θ1 , θ2 , . . ., θP , the elevation angle space as φ1 , φ2 , . . ., φQ , and normally P = Q.As a result, the angle between the target and the X-axis and Yaxis is discretized as ρ1 , ρ2 , . . ., ρPQ and γ1 , γ2 , . . ., γPQ , respectively.In addition, the discretization step is small enough so that the error of the nearest neighbor grid from the exact azimuth and elevation angles of the targets is tolerable.Hence, in the rest of the paper, we only consider the case where each target falls on some angle grid point.This way, Equation ( 6) can be rewritten as where Ψ n = e −j 2π λ (n−1)d r cos( γ1 ) Xa( ρ1 ), . . . ,e −j 2π λ (n−1)d r cos( γPQ) Xa ρPQ (8) is regarded as the sparse basis matrix for the n-th receiver, and s = s 1 , . . ., s pq , . . ., s PQ T is the K-sparse vector of target's reflection coefficients (PQ K), where s pq = η k , k -th target at angle θp , φq 0, no target at angle θp , φq .
Since the radar is colocated, the vector s is the same for all receive antennas.By concatenating the observation data from each receive antenna into a column vector, we construct y = r T 1 , . . ., r T N T , which can express as according to the formulation of compressed sensing, Θ = (ΦΨ 1 ) T , . . . ,(ΦΨ N ) T T denotes the sensing matrix.We let z = (Φw 1 ) T , . . . ,(Φw N ) T T , which denotes the noise vector.If the number of targets is significantly smaller than P × Q, then s PQ×1 can be regarded as a sparse vector, where the locations of non-zero elements indicate the target angles.Up to this point, the radar parameter estimation problem is transformed into a sparse recovery problem in compressed sensing, and s could be found as the sparse solution of Equation (10).Provided that s is sparse enough, a variety of CS methods are available for the sparse recovery of s, such as Basis Pursuit [18], Orthogonal Matching Pursuit [19,20], Regularized Orthogonal Matching Pursuit [21], the NESTA method [22], and so on.

Reconstruction Performance Analysis
In Section 2.1, we have transformed the problem of target parameter estimation in CS-based MIMO radars into a sparse recovery problem.Therefore, in order to improve the performance of radar systems, it becomes essential to increase the probability of successful sparse recovery.To achieve sparse recovery with high probability, the sensing matrix Θ M×N needs to satisfy the restricted isometry property (RIP) [5], which can be understood as follows: the sub-matrices composed of any K columns selected from the sensing matrix are approximately orthogonal to each other.Although RIP is a strong sufficient condition for stable recovery, it is an NP-hard problem to verify it for a specific matrix.In deterministic scenarios where a Θ M×N is given, it is common to evaluate the coherence of the matrix: where Θ l stands for the l-th column of Θ. µ(Θ) is known as the coherence of the sensing matrix (CSM) and reflects the lower bound on the performance of the sensing matrix; in other words, it reflects the recovery performance determined by the two columns of the sensing matrix with the worst orthogonality.For design purposes, it is desirable to keep µ(Θ) as small as possible to improve the upper limit of sparsity, the accuracy, and stability of target parameter estimation [23].

Waveform Design
RIP signifies that the sparse vector can be accurately recovered with a high probability if the sensing matrix Θ is orthogonal.This is impossible for a fat measurement matrix; however, we can optimize the sensing matrix to be as orthogonal as possible by minimizing the difference between the Gram matrix Θ H Θ and an identity matrix.The Gram matrix is a symmetric matrix, and the diagonal elements actually correspond to the auto-correlation coefficients of the column vectors in the sensing matrix, while the off-diagonal elements correspond to the cross-correlation coefficients between the different column vectors in the sensing matrix.Meanwhile, according to the definition of µ(Θ), it can be seen that µ(Θ) happens to be the maximum value of the absolute value of all the off-diagonal elements of the normalized Gram matrix.Previous work has studied waveform design methods that aim to minimize µ(Θ) [6], essentially also making the sensing matrix as orthogonal as possible.
Our design approach focuses on minimizing the coherence of the overall sensing matrix by optimizing the transmit waveforms, in order to improve the performance of CS-based MIMO radars [13].As a design criterion, we seek to minimize the difference between the Gram matrix and an identity matrix while maintaining a fixed total transmit power.The objective function is given by min where κ = P t × N, P t denotes the total transmit power, and the size of the identity matrix is PQ.The l-th column of the overall sensing matrix Θ equals where Ã( ρl ) = I L ⊗ a T ( ρl ), ⊗ is the Kronecker product, and x LM×1 = vec X T , it holds that x H x = P t .As a result, we can rewrite Equation ( 12) and then further formulate the following optimization problem to find waveforms: where the first term, denoted by SCSM, represents the sum of the square magnitude of the cross-correlation between different columns u l , u l in the sensing matrix.The second term aims to maintain the norm of all columns close to √ P t N. u i s are the columns of Θ that implicitly depend on x.The objective function in Equation ( 14) is not necessarily convex but can be transformed into a convex function.The derivation will be shown as follows.
For the sake of simplicity, we can simplify the inner product of u i s by using some properties of the Kronecker product, as shown in the following where is a constant that does not affect the optimization procedure, so it can be ignored.It is readily apparent that V ll is not a positive semidefinite (PSD) matrix unless l = l , which leads to a nonconvex objective function.To address this issue, we split V ll into Vr ll and Vi ll , where , as a result [10] x Clearly, both Vr ll and Vi ll are the Hermitian matrices.As is known, if all eigenvalues (denoted by λ) are greater than or equal to zero, a Hermitian matrix is positive semidefinite.Accordingly, we define b and d as non-negative real scalars that let b P t + λ min (Vr ll ) ≥ 0 and d P t + λ min (−2bGr ll ) ≥ 0.Then, the first term of Equation ( 16) can be rewritten as where C 1 is a constant that does not affect the objective function.In this way, Equation ( 17) is convex since Gr ll and Hr ll are PSD matrices.Performing the same procedure on x H Vi ll x 2 , we can obtain In the same way, the second term in the objective function of Equation ( 14) can be rewritten as Finally, the objective function is transformed into a convex function.However, we can see that the constraint of Equation ( 14) is a quadratic equality constraint, which is not a convex set.Thus, we replace it with a local affine approximation as follows [10] x where x is the estimate of x at the last iteration.In summary, the convex problem of optimizing waveforms with a total power constraint can be formulated as

Optimization Algorithm Based on Simulated Annealing
The solution to Equation ( 21) can be obtained with the Simulated Annealing algorithm (SA) [24].SA is a global optimization algorithm capable of conducting randomized searches in the solution space to find a globally optimal solution, or an approximation of it.SA is usually a good choice for problems where there may be multiple locally optimal solutions to the objective function.Consequently, it finds widespread application in various optimization problems, especially in the field of traditional radar waveform design [6].The process of the Algorithm 1 is shown as follows.

Algorithm 1 Simulated Annealing Algorithm
Input: the initial solution, i.e., initial waveform vector x 0 , the objective function y 0 .Output: the best solution, i.e., optimized waveform x * .
Step (2) : Generate a new neighboring solution x 1 by applying a small perturbation or modifiedtion to the current solution, and calculate the objective function value y 1 .
Step (3) : Verify the constraint and new solution.
if meet the constraint if y 1 < y 0 accept the new solution else accept the new solution in terms of the Metropolis criterion else discard the new solution end Step (4) : Update the temperature according to the cooling schedule.The acceptance probability of worse solutions decreases gradually as the temperature decreases.
Step (5) : Check if the stopping criteria are met.If met, terminate the iteration process and return the best solution as the final output, otherwise, repeat steps 2 to 4.

CS Recovery Algorithm Based on Sparse Bayesian Learning
In this paper, we use the Sparse Bayesian Learning (SBL) method [25,26] to recover the sparse vector s in Equation (10), which contains the angle information of all targets.SBL uses the sparse prior information of the signal to build a hierarchical probability model and perform sparse reconstruction by the maximum a posteriori (MAP) probability function to obtain the parameter estimation of targets.In the field of DOA estimation, SBL is an effective method.When applied to CS-based MIMO radar, it can estimate the target direction efficiently and has excellent performance in the case of limited observation data.In the case of a 2D array MIMO radar, the sparse recovery accuracy of the SBL algorithm is much better than the OMP algorithm.The process of the Algorithm 2 is shown as follows.

Algorithm 2 Sparse Bayesian Learning Algorithm
Input: the observation vector y, the sensing matrix Θ. Output: the estimated angles of targets.
Step (1) : Initialization.Initialize the sparse coefficient estimation vector randomly, representing the weights of received signals at different DOA angles.
Step (2) : Calculate the likelihood function of the current sparse coefficient vector.
Step (3) : Update the sparse coefficient vector by the Expectation-Maximization (EM) algorithm.E-step (3.1) : Calculate the posterior distribution of DOA angles, given the observed data.M-step (3.2) : Maximize the MAP function to update the sparse coefficient vector.
Step (4) : Check if the stopping criteria are met.If met, terminate the iteration process, otherwise, repeat steps 2 to 3.
Step (5) : DOA Estimation.Use the updated sparse coefficient vector to calculate the estimated angles of targets and return the result.

Simulation Configuration
First of all, let us discuss the initial waveform x 0 .In addition to the general requirements for radar waveforms such as good resolution, low sidelobes, etc, MIMO radar waveforms should also possess good orthogonality.However, strictly speaking, the ideal orthogonal MIMO waveform does not exist [27].To be closer to practical radar systems, we use a Hadamard-coded waveform as the initial waveform, unlike the ideal orthogonal waveforms used in other literature.We consider a colocated MIMO radar system based on an L-shaped array, as shown in Figure 1, where d t = d r = λ 2 and the number of transmit antennas is M = 8.The Hadamard coding waveform is a kind of slow-time code division multiple access (CDMA) waveform.For the scenario with 8 transmit antennas, each antenna transmits a chirp signal with the same bandwidth and carrier frequency.The initial phases of eight consecutive pulse signals transmitted by eight antennas are modulated by a 8 × 8 Hadamard matrix.
The parameter settings of the simulation are given in Table 1.The transmitted waveform is an orthogonal Hadamard-coded waveform with unit power.The received signal is corrupted by zero-mean complex Gaussian noise.The measurement matrix Φ is chosen as a Gaussian random matrix with T = round(τ × L), and τ denotes the compression ratio.There are three targets on the angle grid [20 It should be noted that, in CS-based scenarios, the angular grid can range from 0 • to 360 • .However, there might be some prior information about the DOA of targets, or the antennas are not omnidirectional, so only part of the angles are illuminated.In these cases, the angle grid can be limited to a small range.Therefore, to keep the computational complexity low, we consider here a small grid located in the middle of the angle space for the performance evaluation.The reflection coefficients η k s come from a complex Gaussian distribution with zero mean.Each simulation is run 2000 times based on the Monte Carlo method.In each run, three targets are randomly placed on the angle grids.[0.4, 0.5, 0.6, 0.7, 0.8] Signal-to-noise ratio [0, 3, 6, 9, 12, 15]  The number of runs per simulation 2000

Simulation Results
In this section, simulation experiments are performed to illustrate the effectiveness of the proposed waveform design scheme.Our main focus is to evaluate the incoherence of the sensing matrix and the CS-based multi-target detection and parameter estimation performance for the MIMO radar in the above configuration.
(1) Incoherence: First, analyze the incoherence of the sensing matrix.Using the Gram matrix mentioned in Equation (12), calculate the absolute value of normalized coherence coefficients between the different columns of the sensing matrix.All normalized coherence coefficients are counted in steps of 0.01, and then histograms are plotted.The statistical result is shown in Figure 3.It can be seen that normalized coherence coefficients are more concentrated around 0 after optimization.Since the coherence between the target returns from different search cells is large in the 2D array, the maximum value of normalized coherence coefficients, i.e., µ(Θ), is close to 1. Therefore, we use the average normalized coherence coefficient denoted by ACSM to measure the incoherence of the sensing matrix, which is more intuitive.The ACSMs before and after optimization at different compression ratios are given in Table 2.It can be seen that the ACSM is lower after optimization.These two simulation experiments show that the optimized sensing matrix has better incoherence, i.e., better orthogonality, and is therefore more likely to satisfy the RIP criterion, which provides improved sparse reconstruction performance.
(2) Target Parameter Estimation: To verify the effectiveness of the optimized waveform, we examine the root mean squared error (RMSE) of the target scene recovery (i.e., angle estimation) for different signal-to-noise ratio (SNR) values and different compression ratios τ.In the cases of M = 8, N = 20 with τ = 0.5, 0.6, Figures 4 and 5 demonstrate the RMSE curves of the angle estimation produced by the initial waveform (Hadamard coded waveform), the proposed waveform optimization scheme, and the initial waveform without compression for different SNRs, respectively.From the figure, we can see that the performance improvement due to the optimized waveform is evident.Under the same SNR condition, the optimized waveform has a smaller target parameter estimation error, indicating that the effectiveness and robustness of the CS-MIMO radar are improved.The RMSE of angle estimation decreases as the SNR increases.Moreover, the performance improvement due to the optimized waveform is more obvious when the SNR is low, indicating that when the noise interference is high, the waveform optimization scheme can well compensate for the performance loss caused by noise.Although the target parameter estimation performance in compressed perception scenarios is consistently worse than in general scenarios without compression, the smaller amount of data to be processed in compressed perception reduces the signal processing complexity of the system.Thus, a better balance between the accuracy and complexity of the system can be achieved by our waveform optimization scheme.Figure 6 shows recovery results for different compression ratios in the conditions of M = 8, N = 20 and SNR = 15.From the figure, we can clearly see that the optimized waveform has a better target parameter estimation performance than the initial waveform, regardless of the compression ratio.It is worth noting that the performance improvement is particularly significant at the low compression ratio.In conclusion, we have learned that applying compressive sensing to MIMO radar systems offers numerous benefits.However, the reduction in data volume may lead to a loss of accuracy during parameter estimation.In such cases, using the waveform optimization scheme can improve the performance of CS-based MIMO radar systems.The case of M = 8, τ = 0.6 with a different number of receiving antennas is shown in Figure 7.One can see that the waveform optimization scheme has improved RMSE performance as compared to the initial waveform.Again, increasing the number of receiving antennas can improve the target parameter estimation performance of the CS-based MIMO radar when all other conditions are the same.The performance gain with the increase in the number of receiving antennas is more prominent at low SNR.This is because using more receiving antennas effectively increases the array aperture.However, it is important to note that the number of antennas and apertures cannot be increased without limitations, otherwise the practical requirements such as hardware cost, amount of data for signal processing, and size of the radar cannot be met, especially in the field of automotive radar.In this context, our proposed waveform optimization scheme is particularly applicable as it achieves a significant improvement in the performance of the CS-based MIMO radar at a relatively low cost.

Discussion
For the sake of simplicity, the objective function in Equation ( 14) only pertains to the DOA.However, in the general scenario where targets are situated in different range cells and are moving, a range-dependent time delay will be introduced into the waveforms, and a Doppler shift will manifest in the phase term.In this case, when formulating the signal model, the sparse basis matrix has to be represented based on a three-dimensional space and the targets require to be searched in the 3D space.The method proposed in this paper can be easily extended to accommodate the aforementioned model involving range, speed, and angle, but the computational complexity will increase significantly.Therefore, it becomes crucial to propose a more efficient optimization algorithm in future work.

Conclusions
We propose a novel waveform design scheme for colocated CS-based MIMO radars.Under the constraint of fixed total transmit power, we optimize the transmit waveform to minimize the difference between the Gram matrix of the sensing matrix and an identity matrix, thus achieving maximum orthogonality of the sensing matrix.Unlike the onedimensional arrays considered in most of the previous studies, we consider a MIMO radar system based on an L-shaped array with the capability of estimating both the azimuth and elevation angles of the targets.It is more widely used in real-life applications.Furthermore, we use the SBL algorithm for sparse recovery, which is more suitable for 2D array MIMO radar systems.In simulation experiments, we optimize a Hadamard-coded waveform that can be used in practical application scenarios.The simulation results show that under the condition of a 0.5 compression rate of the received data, the MIMO sensing performance has an overall improvement of more than 10% in the range of 0 dB∼15 dB SNR, compared to the non-optimized waveform, and under the condition of 15dB SNR, the MIMO sensing performance has an overall improvement of more than 10% in the range of 0.4∼0.8compression ratio, compared to the non-optimized waveform.It proves that our waveform design scheme can significantly improve the performance of the radar system.

Figure 1 .
Figure 1.Colocated MIMO radar system based on an L-shaped array.

Figure 6 .
Figure 6.The recovery results for different compression ratios with M = 8, N = 20 and SNR = 15.

Figure 7 .
Figure 7.When M = 8 and τ = 0.6, the recovery results for different values of SNR with different receiving antennas.