A Phase-Only Optimization Null Control Method for FDA-MIMO Based on ADMM

: This paper investigates null control within the transmit–receive beampattern of Frequency Diverse Array-Multiple-Input and Multiple-Output (FDA-MIMO) systems, presenting a novel phase-only optimization approach for achieving null control in FDA-MIMO. We employ an alternating multiplier framework, which transforms the intricate and inherent constant modulus constraint and numerous amplitude constraints in optimization into more manageable projection problems. By employing a phase-only optimization strategy, the intricate hardware and computational burdens associated with null control in FDA-MIMO are effectively alleviated. The simulation results indicate that the algorithm proposed in this paper exhibits excellent null control ability while precisely maintaining constant modulus constraints, and it possesses an extremely high computational efficiency.


Introduction
Null control is a crucial technique in beamforming [1], widely employed in the field of radar systems.By controlling nulls in array antenna systems, it is possible to effectively suppress interference signals, thereby enhancing the performance of signal reception or transmission.Currently, null control has garnered extensive attention and research in practical applications, achieving significant advancements.For instance, in radar systems, null control can be utilized to suppress spurious targets or interference sources in echoes, thereby improving target detection and tracking performance [2].
Indeed, the problem of synthesizing beampatterns with widened nulls has been extensively investigated [3,4].A novel beamforming approach based on the Minimum Dispersion (MD) criterion has been designed in [5], which can limit the average power of interference in dynamic angular sectors to zero.A minimum variance beamformer with a predefined level of suppression on specified zero-angle sectors was proposed in [6], which utilizes a set of linear constraints, replacing the original quadratic constraints in the optimization problem, thereby reducing the computational complexity of the design.Additionally, an innovative robust null broadening adaptive beamformer was introduced in [7], which effectively widens interference nulls and enhances null depth through the use of projection transformation and diagonal loading techniques.The above methods are all aimed at null control in the spatial direction dimension, but it seems difficult to suppress the interference of the main lobe.
Frequency diverse arrays (FDAs) possess a unique distance-dependent characteristic [8]; therefore, FDAs can achieve precise null control in the distance-angle twodimensional plane, overcoming the limitation of traditional phased arrays that can only generate nulls in the angle dimension, giving FDAs a significant advantage in main lobe interference suppression [9].However, the time-varying nature of FDAs prevents them from transmitting beams due to achieving high-gain beams or deep nulls at constant distance-angle positions [10].An FDA-MIMO scheme was proposed, which allows the beamforming pattern to focus on addressing time-varying issues [11][12][13][14].Scholars have designed specialized receiving structures for FDA-MIMO [15][16][17][18][19] and conducted research on FDA-MIMO directional null control based on this [20][21][22][23][24].A null control method in the main lobe distance dimension using transmit weight vector design was proposed based on the standard FDA in [20], which can analytically set multiple nulls at an arbitrary distance dimension and control their width.However, this method neglects the time-varying impact on the beampattern, resulting in nulls that only briefly reside in spatial positions.Additionally, it fails to address the inherent coupling issue between the distance and angle, limiting target localization and interference performance.A null control technique based on a covariance matrix reconstruction was introduced based on the FDA-MIMO structure in [21], which widens nulls by placing virtual interference around actual interference locations.A method for designing a Preset Null Broadening Beamformer (PBNBF) by strategically placing artificial interference with appropriate power around the nulls in the equivalent transmit beampattern was proposed in [22], which widens the nulls and ensures effective suppression of deceptive interference.Lan et al. proposed a precise control method for nulls in the transmission-reception beampattern in [23].This method can calculate the power of virtual interference based on preset null levels and then construct a virtual covariance matrix, ultimately obtaining a closed-form solution for weight vectors using the adaptive array theory, which exhibits more precise null control capabilities compared to [21,23].The skew projection technique was further employed in [24], avoiding the mutual interference issues that arise when controlling multiple nulls simultaneously.
It can be observed that in the current research on null control in FDA beampatterns, the majority of studies are focused on FDA-MIMO.However, in the research presented in [21,23,24], all of them employ complex weighted vector optimization, which necessitates numerous high-dynamic amplifiers and intricate feed networks.In this scenario, the expansion of virtual array apertures leads to a higher dimensionality of optimization variables, making null control techniques for FDAs often associated with expensive computational costs and complex hardware requirements.Over the years, some researchers have considered using phase-only optimization methods to achieve antenna array beampattern nulling [25][26][27].S. T. Smith solved the optimal phase-only adaptive weighting problem by optimizing the signal-to-interference-plus-noise ratio (SINR) constrained by phase-only vectors [25].P. J. Kajenski proposed a phase-only method for generating notches in the beampattern of phased array antennas [26].Similarly, L. Pallottta et al. synthesized optimal space-time adaptive processing (STAP) using a phase-only optimization method [27].
Traditional methods need to handle both a significant number of constant modulus constraints and a large number of amplitude constraints simultaneously, whose dual requirement adds complexity to the optimization process.To address this challenge, this paper focuses on the research related to null control in FDA-MIMO transmit-receive beampatterns and proposes a method, which applies for FDA-MIMO null control based on the alternating direction method of multipliers (ADMM) [27] for phase-only optimization.The core idea behind this approach is to introduce the ADMM framework to transform complex constant modulus and amplitude constraints into simple projection problems, which reduces both hardware requirements and the computational complexity.
The remainder of this paper is organized as follows.Section 2 presents the FDA-MIMO transmit-receive signal model and principle.To improve the null depth and computational efficiency of FDA-MIMO, the ADMM and the double-layer alternating direction method of multipliers (DADMM) methods are investigated in Section 3. Simulation results are given in Section 4. Conclusions are drawn in Section 5.

FDA-MIMO Transmit-Receive Signal Model
Considering a standard FDA as shown in Figure 1, the transmitted signal superimposed on a far-field point target can be represented as where f 0 and ∆ f denote the reference carrier frequency and frequency increment, respectively; w T,m is the complex weight coefficient corresponding to the m-th element; f m represents the frequency of the m-th transmit array element; and τ m = r m /c represents the propagation delay from the m-th element to the target.ϕ(t) is the transmission baseband signal per unit energy, T p ϕ(t)ϕ * (t)dt = 1,where T p is the pulse duration.According to the narrowband assumption, ϕ(t − τ m ) ≈ ϕ(t − τ), where τ represents the propagation delay between the target and the reference element, Equation (1) can be approximately expressed as For ease of analysis, we first derive the case of a single-antenna receiver.Assuming the number of receiving element is 1, the signal at the receive antenna can be represented as where ξ represents the amplitude and phase variation caused by the electromagnetic wave propagation path and target scattering, (•) T represents the transpose operation, ⊙ represents for the Hadamard product operator, and w T = [w T,1 , w T,2 , ..., w T,M ], a t (t), a θ (θ), and a r (r) are defined as From (3), it can be observed that the time variable t simultaneously exists in the waveform correlation term ϕ(t − 2τ) and the beampattern correlation term w H T [a t (t) ⊙ a θ (θ) ⊙ a r (r)], the latter of which leads to the time-varying phenomenon of the standard FDA transmit beampattern.This phenomenon not only results in the inability to concentrate the transmitted energy but also renders traditional heterodyne filtering receiver designs unsuitable for FDAs [10,28].For this purpose, Gui et al. proposed a receiver model based on multi-channel mixing and matched filtering [18], as shown in Figure 2.  In the case of the single-antenna receiver, the output of the matched filter for the n-th channel can be expressed as where ξ = ξe j4π f 0 r c and a mn (r, where t r = t − 2τ.From (7), the time variable t has been separated into the ambiguity function A(t r , [m − n]∆ f ) after matched filtering, and the steering vector that influences the spatial beam direction is now independent of the time variable t.However, it is important to note that A(t r , [m − n]∆ f ) couples with the spatial term a mn (θ, r) and weight w * T,m through a summation operation, leading to space-time coupling that can cause range measurement biases, which can result in the peak of y n (t) not necessarily occurring at A(0, 0), thereby affecting the accurate measurement of the target's distance.As A(t r , f d ) is a function related to the envelope of the transmit waveform, it is necessary to carefully design the transmit waveform to satisfy the following decoupling condition: where A R (t r ) and A D ([m − n]∆ f ) are the time delay and Doppler tangent plane of the ambiguity function, respectively.Substituting ( 9) into (7), the n-th channel filter output is expressed as By applying weighting to y n (t r ) through the receiver beamformer, we can obtain where w F = [w F,0 , w F,1 , ..., w F,M−1 ] T is the receive weight vector of the beamformer and A θr = a θ (θ)a T r (r) represents the transmit-receive steering vector matrix.In fact, A D = {A D ([m − n]∆ f } M×M and A D can be regarded as the cross-correlation matrix of FDA's transmit signals. Equation ( 11) is a general FDA receiving signal model, which is applicable to both FDAs with coherent transmitted waveforms and with orthogonal transmit waveforms.Consider an FDA with orthogonal transmit signal, where A D = I M and where I M represents the diagonal matrix of M × M. In this case, ( 11) can be re-expressed as where represents the equivalent transmit steering vector and diag[•] represents the operation of taking diagonal elements.
In the multi-antenna situation, assuming that N elements are used for reception, the extension of ( 12) can be re-expressed as where 1 N represents a vector of ones with N dimensions, w F is the reception weight vector with MN dimensions, ⊗ denotes the Kronecker product operator, and b(θ) = [1, e j2π f 0 dsinθ/c , ..., e j2π f 0 (N−1)dsinθ/c ] T represents the steering vector of the receive array.

Problem Statement
After the treatment in Section 2, the normalized transmit-receive array factor of FDA-MIMO can be represented as In the case of an orthogonal transmit signal, the transmit weight coefficient w T only affects transmit power, so we do not consider optimizing w T .Without the loss of generality, let w T = 1 M and define u(θ, r) = a(θ, r) ⊗ b(θ); then, ( 14) can be expressed as In the phase-only optimization of null control, two general categories of null control problems are considered [26].The first type involves maximizing the main lobe gain while keeping the nulls in a specified region below a certain threshold value.The second type aims to minimize the null levels in a specified region while ensuring that the given main lobe gain is not less than a certain value.The corresponding optimization problems can be formulated as Q 1 and Q 2 : where ε n represents the desired null constraint level, Ω n represents the null area, and A 0 is the desired main lobe gain.Due to the presence of constant modulus constraints, both problems Q 1 and Q 2 are complex non-convex optimization problems and are difficult to deal with directly.Currently, there is limited research on phase-only null control based on FDAs.Classical semidefinite relaxation (SDR) techniques were employed to address the above problems in [29]; however, SDR techniques have two main drawbacks, as follows: Firstly, SDR transforms the original vector optimization problem of dimension MN × 1 into a matrix optimization problem of dimension MN × MN, leading to increased computational complexity.
Secondly, SDR discards the non-convex rank-one constraint, resulting in an approximate solution that may not satisfy the constant modulus constraint, necessitating further correction.
Additionally, compared to phased arrays, null control on a distance-angle 2D plane requires imposing more constraints, further increasing the computational burden of the SDR algorithm.To overcome these shortcomings, this paper introduces an efficient phase-only FDA null control technique based on the ADMM [30] in the subsequent discussion.
ADMM is a framework for solving separable convex optimization problems.It takes the form of a decomposition-composition process, where solutions to smaller local subproblems are combined to find a solution to the larger global problem.ADMM combines the advantages of dual decomposition and augmented Lagrangian methods, providing a decomposable algorithm with good convergence properties and fast processing speed.As a result, it has found widespread application in array signal processing and machine learning.Therefore, in this section, we utilize the ADMM method to solve problems Q 1 and Q 2 .

Solution for Q 1 Based on ADMM
It can be observed that due to the presence of absolute values, the objective function of problem Q 1 is not smooth.Therefore, we first introduce an auxiliary variable t that satisfies |t| = 1, and then we can obtain the following relationships for the objective function: where d p is the upper bound for |w H F u(θ 0 , r 0 )|, and when the magnitude of w F is 1, Equation (19) indicates that by updating t during the optimization process to satisfy ∠t = ∠w H F u(θ 0 , r 0 ), the objective function in Q 1 can be transformed into a typical quadratic programming (QP) form.After the transformation given in (19), Q 1 can be equivalently represented as where E k is a MN × MN diagonal matrix that is satisfied, as follows: Under the ADMM algorithm framework, introducing two auxiliary variables, v = w F and h j = w H F u(θ j , r j ), (θ j , r j ) ∈ Ω n , Q 1 can be further equivalently represented as where J represents the number of grid points after the discretization of the null region Ω n .Applying the ADMM algorithm framework to (23), the corresponding augmented Lagrangian function L p (w F , v, h, γ 1 , γ 2 ) can be expressed as where γ 1 and γ 2 are the dual variables, and 2 are penalty factors that affect the algorithm's convergence speed.Here, A J = [u(θ 1 , r 1 ), u(θ 2 , r 2 ), ..., u(θ J , r J )] and h J = [h 1 , h 2 , ..., h J ] T .The ADMM algorithm optimizes the variable (w F , v, h, γ 1 , γ 2 ) through alternating iterative updates, and the update at the g-th iteration can be expressed as According to (24), the closed-form update for the variable (w F , v, h, γ 1 , γ 2 ) at each iteration can be derived as follows: (1) Updating w g+1 F .Solving for w g+1 F involves finding the extreme of an unconstrained complex-valued function, which achieves its optimal value when the gradient is zero, satisfying the following: Clearly, w g+1 can be computed as where Importantly, after obtaining the update for w g+1 F , it is necessary to update the auxiliary variable t to maintain the equivalence of the problem; thus, the g-th iteration of t can be expressed as (2) Updating v g+1 .Solving for v g+1 can be equivalently transformed into the following problem: where Clearly, Equation ( 29) is an Euclidean projection problem, and its solution for v g+1 can be easily obtained as where ∅ represents element-wise division.
(2) Updating h g+1 .Similar to (29), the update for h g+1 is also an Euclidean projection problem, and its solution is are as shown in (24).The stopping criteria for the ADMM algorithm are controlled by the primal constraint residuals ε r1 , ε r2 , and the dual residual ε d , as defined in (32).Iterations are stopped when ε r1 , ε r2 , and ε d are smaller than their specified thresholds.
The stopping criteria is provided for the ADMM algorithm in [30], which are where ξ r1 , ξ r2 , and ξ d are the thresholds for the primal constraint residuals and dual residuals, and they can be represented as where ϵ abs represents the absolute error tolerance and ϵ rel represents the relative error tolerance.

Solution for Q 2 Based on DADMM
Similar to the approach in Q 1 , two auxiliary variables v and h = [h 1 , h 2 , ..., h J ] T are introduced to transform Q 2 into a more manageable form, as follows: The corresponding augmented Lagrangian function L p (w F , v, h, η, γ 1 , γ 2 ) can be expressed as The parameters in (36) are consistent with Equation (23).Notice the interdependence of h and η.Thus, {h, η} is iterated as a set of variables.According to the ADMM framework, the update form for variable (w F , v, {h, η}, γ 1 , γ 2 ) in the g-th iteration is given by (37).
According to Equation (37), we can derive the closed-form solution for the variable (w F , v, {h, η}, γ 1 , γ 2 ) at each iteration, as follows: (1) Updating w g+1 F .The update for w g+1 F is equivalent to solving the following suboptimal problem.
Equation ( 38) is a convex optimization problem under inequality constraints and can be solved using the Lagrange multiplier method in conjunction with the Karush-Kuhn-Tucker (KKT) conditions [31].However, the required algebraic operations are rather complex.Therefore, this paper introduces a double-layer ADMM (DADMM) framework to further simplify this problem.The specific approach involves introducing an auxiliary variable b 0 = w H F u(θ 0 , r 0 ) so that (38) can be equivalently transformed as Similarly, the augmented Lagrangian function for (39) can be derived as where κ is the dual variable and ρ 3 is the penalty factor.According to the ADMM variable update form, the iterative forms of w F , b 0 and κ can be obtained as Observing (41), we can see that the original inequality-constrained optimization problem (38) has been transformed into two simple unconstrained optimization problems and a straightforward Euclidean projection problem.Each step of the solution process has clear closed-form solutions, reducing the complexity of problem solving and improving efficiency.After a straightforward derivation, the iterative closed-form solutions for w where Once the inner ADMM iterations stop, the output w g+1 F is used for the outer ADMM, and the inner ADMM iteration stopping criterion is given by where (2) Updating v g+1 .Solving for v g+1 is likewise a simple Euclidean spherical projection problem, and its solution is consistent with (45).
According to , the problem of solving for {h g+1 , η g+1 } can be formulated as Observing (46), it can be seen that when the variable η is fixed, h g+1 can be regarded as a Euclidean projection problem, and the corresponding closed-form solution is given as where By combining (47) and ( 46), we can obtain an unconstrained optimization problem that depends only on η, as shown below.
Since the expression of h g+1 depends on min{|h g+1 j |, η}, let us consider different cases for η. , and problem (48) can be represented as It can be observed that problem (50) is a quadratic function with respect to η, and its minimum point η m can be expressed as Considering the constraint on η ≤ min(|h g+1 j |), η g+1 can be calculated as The objective function in ( 53) is also a quadratic function, and its minimum point η m in the unconstrained case is Considering the constraint on |h g+1 j p | ≤ η ≤ |h g+1 j p +1 |, the optimal value η g+1 can be calculated as Thus, the update for {h g+1 , η g+1 } can be completed through Equations ( 47), ( 49), (52), and (55), while the updates for the dual variables γ g+1 1 and γ g+1 2 are consistent with (37).In addition, the outer ADMM algorithm for solving problem Q 2 follows the same iteration stopping criteria as solving problem Q 1 in Equations ( 32)∼ (34).

Algorithm Summary and Computational Complexity Analysis
The specific execution steps of the proposed ADMM and DADMM methods are summarized in Algorithm 1 and Algorithm 2, respectively.In the ADMM method, the main computational workload is in the updates of w g+1 and t g+1 in step 2, while the other steps involve only vector multiplications and divisions, resulting in a computational complexity of O(M 3 N 3 + (J + 1)M 2 N 2 + (J + 2)MN).Assuming that the number of iterations required for problem Q 1 to converge is I Q 1 , the overall computational complexity is O(I Q 1 (M 3 N 3 + (J + 1)M 2 N 2 + (J + 2)MN)).In most cases, the value of I Q 1 depends on the desired level of precision in the solution.
In the DADMM method, the main computational workload includes the inner ADMM updates in step 2, resulting in a computational complexity of O(I I N (M 3 N 3 + (J + 1)M 2 N 2 + (J + 2)MN)), where I I N is the number of iterations required for inner ADMM to converge.Assuming that the number of iterations required for problem Q 2 to converge is I Q 2 , the total computational complexity for solving problem

Simulation Results
This section conducts several simulation experiments to validate the effectiveness of Algorithms 1 and 2. The simulation parameters are set as follows: the number of transmit elements M = 15, the number of receive elements N = 15, the carrier frequency f 0 = 10 GHz, the element spacing d = 0.015 m, and the frequency offset ∆ f = 5 kHz.The entire detection region is set as Ω = {(θ, r)| − 90 • ≤ θ ≤ 90 • , 15 km≤ r ≤ 45 km}, and the null region is defined as 18 km≤ r ≤ 22 km}, with sampling intervals in the range dimension and angle dimension set as 0.2 km and 0.5 • , respectively.The desired main lobe is directed toward (θ 0 , r 0 ) = (10 • , 30 km).Typical phase-only optimization methods based on SDR technology are used for a performance comparison, including the Semidefinite Relaxation with Linear Constraints (SDR−LC) algorithm and the Semidefinite Relaxation with Quadratic Constraints (SDR−QC) algorithm [26].In both SDR−LC and SDR−QC methods, a rank-one approximation method [32] is used to calculate the receive weight matrix w F .

Given Null Depth Constraints to Maximize the Main Lobe Gain (Q 1 Based on ADMM)
In this case, the performance of different algorithms in solving the problem of maximizing the main lobe gain under given null depth constraints is analyzed and compared.The preset null depth level is −70 dB, i.e., ε n = 10 (−70/20) MN.Unless otherwise specified, the penalty factor is set to ρ 1 = ρ 2 = 100, and the absolute error tolerance and relative error tolerance are set to ϵ abs = 10 −4 and ϵ rel = 10 −4 , respectively.Figure 3 compares the transmit-receive beamforming patterns using different methods.According to Figure 3a, it can be observed that even without any optimization, the transmit-receive beampattern of the standard FDA-MIMO is decoupled in the range and angle dimension.According to Figure 3b-d, it can be seen that SDR−QC, SDR−LC, and the proposed Algorithm 1 all form lower nulls in the specified null region.To further demonstrate whether the null levels in the null region satisfy the preset null depth constraints, Figure 4 compares the projection of the null region in the direction patterns formed by different algorithms in the range dimension.It can be observed that due to the quadratic constraints used in SDR−QC, it can only describe the average null level, and the nulls formed by the SDR−QC algorithm exceed the constraint requirements in some regions, while SDR−LC and Algorithm 1 can satisfy the null depth constraints.However, the SDR−LC algorithm uses linear constraints to control the null levels, and achieving precise null control in the range-angle domain requires a large number of constraints, leading to extremely high computational complexity.Unlike SDR technology, Algorithm 1 in this paper adopts the ADMM framework, which transforms the constrained optimization problem into simple unconstrained optimization problems or projection problems, avoiding the need for complex matrix operations and significantly reducing computational complexity.
Table 1 provides a detailed performance comparison between different algorithms, where the violation of the constant modulus constraint (VCMC) is defined as || |w| − 1||.
From Table 1, it can be seen that Algorithm 1 requires only 0.4 s for solving problem Q 1 , which is a significant performance improvement compared to 10.2 s for SDR−QC and 842 s for SDR−LC.In terms of null depth, Algorithm 1 and SDR−LC achieved −69.76 dB and −69.68 dB, respectively, which are very close to the preset level requirements.The slight error Ω n is due to the discrete approximation of the null region and the limited precision of the solution.Regarding the main lobe gain loss, SDR−LC and Algorithm 1 have similar performance, with losses of 1.02 dB and 1.01 dB, respectively, while SDR−QC has a relatively smaller main lobe gain loss of 0.92 dB, which is a result of the trade-off between antenna main lobe gain and sidelobe performance.Examining the performance parameters of various algorithms in Table 1, it appears that the algorithm proposed in this paper holds an advantage solely in computational efficiency.However, it is noteworthy that Algorithm 1 exhibits significantly lower violation of the constant modulus constraint compared to the SDR−QC and SDR−LC algorithms.In fact, due to the highly dimensional receiving weight matrix in the FDA-MIMO, solutions obtained by SDR−LC and SDR−QC using rank-one approximation often contain numerous violations of the constant modulus constraint.Figure 5 compares the weights obtained with different algorithms, clearly revealing that both SDR−LC and SDR−QC have instances where the modulus values are not equal to 1.If the weight decreases, it will result in signal power loss [33], thereby reducing the effective range of the system.However, the constant modulus constraint can effectively reduce signal loss.It is easy to imagine that if the weights obtained with SDR−QC and SDR−LC are reset to 1, the performance of the beampattern will decrease significantly.The null performance of the beampattern after setting the weights to 1 is shown in Table 2.It can be observed that, under precise constant modulus constraints, SDR−QC and SDR−LC methods can only maintain null levels of −42.41 dB and −46.87 dB, respectively, significantly exceeding the prescribed null depth requirement of −70 dB.In contrast, Algorithm 1 proposed in this paper, by maintaining precise constant modulus constraints, preserves the original performance.In contrast, Algorithm 1 consistently ensures the constant modulus constraint, representing another advantage of the proposed algorithm over SDR−QC and SDR−LC.To illustrate the convergence process of Algorithm 1 in this paper, Figure 6 depicts the curves of the objective function and constraint variations.From Figure 6a, it is evident that the objective function rapidly decreases and then gradually increases, reaching an extremely low level within 30 iterations, demonstrating the excellent convergence capability of Algorithm 1.In Figure 6b, it can be observed that the constraint conditions are constantly approached as the iteration progresses, which is also the reason for the subsequent gradual increase in the objective function.Additionally, to validate the impact of different penalty factors on the convergence of the proposed algorithm, Figure 7 provides the curves of the objective function variation for different values of the penalty factors ρ 1 and ρ 2 .It can be observed that a smaller penalty factor leads to a faster convergence.However, it is important to note that the penalty factor should not be set too small as this may result in the algorithm failing to converge [34].In this section, an analysis and comparison of the performance of different algorithms in solving the problem of maximizing null depth under a given main lobe gain constraint were conducted.Unless otherwise specified, the penalty factors were set as ρ 1 = ρ 2 = 100 and ρ 3 = 50, and the main lobe gain constraint was denoted as A 0 = 0.9MN, indicating a main lobe loss of 0.92 dB. Figure 8 provides a comprehensive comparison of the beampatterns generated using different algorithms.It can be observed that the null depths formed by Algorithm 2 in this paper and the SDR−LC algorithm are comparable and superior to the SDR−QC method.Both Algorithm 2 and SDR−LC achieved precise control of null responses, which is crucial in certain radar-communication integration applications.On the other hand, the null region levels in SDR−QC exhibited noticeable fluctuations.
Table 3 provides a detailed comparison of the performance of different algorithms.The proposed Algorithm 2 still exhibits the highest computational efficiency in solving such problems, with a runtime of only 1.2 s, while SDR−QC and SDR−LC require 9.28 s and 560 s, respectively.Due to the presence of the main lobe gain constraint, SDR−QC, SDR−LC, and Algorithm 2 all achieve consistent main lobe gain performance.However, in terms of null depth, Algorithm 2 in this paper is slightly lower than SDR−QC and SDR−LC.As analyzed in Q 1 , this phenomenon is attributed to the disruption of constant modulus constraints by SDR−QC and SDR−LC methods.Figure 9 illustrates the modulus values of weights obtained with different algorithms.It can be observed that SDR−QC and SDR−LC have numerous violations of the constant modulus constraint, while Algorithm 2 in this paper consistently maintains the constant modulus constraint.Table 4 presents a performance comparison after projecting different algorithms into the feasible solution space.Due to the enforced constant modulus constraint, the SDR−QC and SDR−LC methods can only maintain null levels of −22.31 dB and −22.51 dB, respectively, whereas Algorithm 2 in this paper still maintains its original performance.It is worth mentioning that both the FDA-MIMO beampatterns synthesized by both the ADMM and the DADMM methods, as shown in Figures 3 and 8, show an increase in sidelobe regions.This is caused by the addition of null region, which introduces additional phase distortions or oscillations.For practical antenna applications, designers can add threshold conditions for the sidelobe level as constraints in problems Q 1 and Q 2 to meet different application requirements.Additionally, it is worth noting that random and unnecessary fluctuations during the matched filtering process can interfere with the phase of the ideal reference signal [35].Our theoretical derivations have neglected this impact.We plan to introduce phase noise in our future research and further study FDA-MIMO null control under phase noise conditions.

Conclusions
This paper focuses on the research of null control in the transmit-receive beampattern of phased array antennas.It begins by deriving the signal model for the phased array receiver, establishing the model for null control in the transmit-receive beampattern of the FDA-MIMO system.Based on this, a phase-only optimization approach for null control in the FDA-MIMO system is proposed using the ADMM algorithm.This method addresses the challenges of the hardware complexity and computational load associated with the high dimensionality of FDA-MIMO null control by adopting a phase-only optimization strategy to reduce hardware complexity.Subsequently, the ADMM framework is employed to solve this constant modulus constraint optimization problem.Numerical experimental results demonstrate that the proposed algorithm exhibits a superior null control capability, precise constant modulus constraint retention ability, and extremely high computational efficiency.In our future work, we will test the effectiveness of the algorithm using real recorded interference signals to further validate its practical applicability.

Figure 1 .
Figure 1.Uniform FDA with linear frequency offset.
θ) = e j2π(m d f 0 sinθ c −n 2∆ f r c ) , with * representing the convolution operator.The definition of the ambiguity function A(t r , [m − n]∆ f ) is given as

6 :
Calculate the primal and dual residuals ε

6 :
Calculate the primal and dual residuals ε

Figure 4 .
Figure 4. Comparison of projections of null regions in beampatterns on the distance dimension.(a) SDR−QC FDA.(b) SDR−LC FDA.(c) Proposed method.

Figure 5 .
Figure 5.Comparison of modulus values for receiving weights obtained with different algorithms.

Figure 6 .
Figure 6.The curves depicting changes in the objective function and constraints.(a) The evolution of the objective function.(b) The evolution of constraint changes.

Figure 7 .
Figure 7.Comparison of objective function variation curves under different penalty factors.

Figure 8 .
Figure 8. Beampattern comparison of different null control algorithms.(a) SDR−QC beampattern projection in distance-angle domain.(b) SDR−QC beampattern projection in distance domain.(c) SDR−LC beampattern projection in distance-angle domain.(d) SDR−LC beampattern projection in distance domain.(e) Proposed method beampattern projection in distance−angle domain.(f) Proposed method beampattern projection in distance domain.

Figure 9 .
Figure 9.Comparison of modulus values for receiving weights obtained with different algorithms.

Table 1 .
Performance comparison of different null control algorithms.

Table 2 .
Performance comparison of different algorithms after projection onto feasible solution space.

Table 3 .
Performance comparison of different null control algorithms.

Table 4 .
Performance comparison of different algorithms after projection onto feasible solution space.