A Robust Complex α-Sigmoid Affine Projection Algorithm Under Non-Gaussian Noise

TL;DR

This paper introduces a new complex-valued adaptive filtering algorithm that performs better in noisy environments with correlated signals.

Contribution

The novel α-CSAP algorithm uses an α-Sigmoid cost function to improve performance in non-Gaussian noise.

Findings

The α-CSAP algorithm suppresses impulsive noise and reduces computational complexity.

Theoretical analysis provides a steady-state mean square deviation expression.

Simulations show superior performance in system identification and beamforming.

Abstract

To address the performance degradation of traditional adaptive filtering algorithms in environments with correlated input signals and non-Gaussian noise, this paper proposes a complex-valued affine projection algorithm based on the α-Sigmoid cost function (α-CSAP). The algorithm leverages the nonlinear characteristics of the α-Sigmoid function and implicitly achieves variable step-size updates by introducing a normalization factor, which effectively suppresses impulsive noise interference and avoids matrix inversion, thereby reducing computational complexity. Theoretical analysis derives the steady-state mean square deviation (MSD) expression for the algorithm. Simulation results demonstrate that the proposed α-CSAP algorithm exhibits superior performance compared to traditional complex adaptive filtering algorithms in both system identification and beamforming application scenarios.

Figures

Full text

1. Introduction

As a core technology in the field of digital signal processing, adaptive filtering, which does not require prior knowledge of the statistical characteristics of signals and noise and can dynamically adjust parameters to adapt to time-varying environments, has become a key means for radar systems to cope with complex signal and interference conditions [1,2,3,4]. Conventional fixed filters exhibit limited performance when dealing with non-stationary, nonlinear signals or unknown interference, whereas adaptive filtering, through iterative optimization algorithms such as Least Mean Square (LMS), Recursive Least Squares (RLS), and particle filtering, enables efficient denoising, target detection, and parameter estimation, and is widely applicable in engineering scenarios such as radar and communications [5,6,7].

In radar signal processing, adaptive filtering plays an indispensable and central role. Radar echoes are often affected by non-Gaussian noise such as spherically invariant clutter and pulse interference, making it difficult for traditional detection methods to achieve a satisfactory balance between detection probability and false alarm rate [8,9]. Adaptive filtering can effectively address such complex signal environments by dynamically adjusting filter weights or structures. For example, the adaptive matched filter can maintain a constant false alarm rate under unknown clutter covariance matrices, significantly improving radar target detection performance [10]; optimized particle filtering can achieve high-precision prediction of highly dynamic target trajectories, effectively handling tracking problems in nonlinear and non-Gaussian scenarios [11]. At the signal and data processing level, the RLS algorithm, with its fast convergence characteristics, is particularly suitable for processing non-stationary radar signals [12]. Furthermore, the integration of adaptive filtering with techniques such as digital beamforming and data fusion further enhances the anti-interference capability and target resolution accuracy of radar systems [13].

As a central direction in modern signal processing and machine learning, complex-valued adaptive filtering is of broad application value compared to real-valued filtering techniques, owing to its inherent suitability for complex signals prevalent in digital communications, power systems, radar, and other fields [14]. Traditional algorithms, often based on the Mean Square Error (MSE) criterion, are simple to implement but tend to become inaccurate or diverge in non-Gaussian impulsive noise environments. Moreover, they typically neglect the statistical properties of non-circular signals, leading to information loss [15]. The integration of information-theoretic criteria has led to algorithms based on complex correntropy and generalized complex correntropy, significantly enhancing robustness against non-Gaussian noise [16]. However, challenges remain, including the insufficient flexibility of single-kernel functions, incomplete utilization of error probability distributions, performance degradation under noisy input–output models, and the inherent trade-off between convergence speed and steady-state error [17]. In response, researchers have pursued optimization through various means, such as employing widely linear models to capture non-circular characteristics, introducing mixed-kernel functions to enhance adaptability, and designing fixed-point iterative or recursive forms to reduce complexity [18]. These efforts aim to better adapt these algorithms to complex engineering scenarios. Complementary to complex adaptive filtering, beamforming is a key technology in array signal processing. It enhances desired signals and suppresses interference by controlling the phase and amplitude of signals received by multiple sensors/antennas, demonstrating significant advantages in wireless communications, acoustic processing, and radar detection. Its core value lies in utilizing spatial diversity to improve signal quality and transmission efficiency, giving rise to specialized applications such as distributed collaborative beamforming [19], microphone array speech enhancement, and broadband beamforming with fractal antenna arrays [20]. Innovations in hybrid architectures and unsupervised methods further extend its applicability [21]. As important branches of adaptive signal processing, both complex adaptive filtering and beamforming face common challenges regarding adaptability in complex environments and balancing performance with complexity. Their technological integration and collaborative optimization have become important research directions driven by practical engineering needs.

In applications such as system identification, echo cancellation, and signal prediction, the performance of traditional adaptive filtering algorithms is significantly affected when the input signal is correlated. Their core operation lies in adjusting filter coefficients in real time to adapt to dynamic environments [22]. The classic LMS algorithm is widely adopted due to its simplicity, but its convergence speed degrades notably under correlated inputs. RLS algorithms converge faster but face high computational complexity, making them less suitable for real-time requirements [23]. To balance convergence speed and computational complexity, Ozeki and Umeda first proposed the Affine Projection (AP) algorithm in 1984 [24]. By reusing past input data, it achieves faster convergence than LMS while maintaining lower complexity than RLS, representing a major breakthrough in adaptive filtering [25]. However, practical challenges such as impulsive noise, highly correlated inputs, and sparse system characteristics impose limitations on traditional AP algorithms. For example, in acoustic echo cancellation and hands-free vehicular communications, impulsive interference can degrade algorithm performance [26]. In sparse system identification, conventional AP algorithms fail to exploit sparsity, leading to inefficient convergence [27]. In massive data processing scenarios, indiscriminate data updates result in wasted computational resources. Consequently, extensive research has been dedicated to improving the AP algorithm. Various strategies have been developed to enhance its performance. Data-selective updating mechanisms have been introduced to reduce computational load while maintaining estimation accuracy. Robust AP algorithms based on M-estimation have been designed to improve resilience against impulsive noise [28]. Variants incorporating sparsity-inducing penalties have been proposed for sparse system identification [29]. Furthermore, approaches that avoid direct or indirect matrix inversion have been adopted to lower the computational complexity of AP-type algorithms [30]. Additional techniques, such as subband processing and variable step-size strategies, have also been employed to further optimize convergence behavior under correlated inputs [31]. These improved algorithms demonstrate advantages in applications like system identification, echo cancellation, and beamforming. Nonetheless, achieving an optimal balance among robustness, real-time performance, and accuracy in complex environments remains an open challenge. Furthermore, while most existing AP algorithms are confined to the real domain, their application in the complex domain is increasingly important due to the prevalence of complex signals in radar signal processing.

Recently, researchers proposed an adaptive projection algorithm based on an $[eqn]$ -sigmoid cost function, demonstrating a superior performance when dealing with strongly correlated input signals and non-Gaussian noise environments [32]. This provides a new approach for radar echo cancellation and channel estimation. Building upon this, we propose a novel complex-valued adaptive filtering algorithm based on the $[eqn]$ -sigmoid cost function. The steady-state MSD of the proposed algorithm is derived analytically. The performance of the proposed algorithm is compared with classical complex adaptive filtering algorithms in simulations of system identification and beamforming. The results confirm the effectiveness and robustness of the proposed algorithm. Furthermore, the proposed algorithm offers a new potential solution for beamforming applications.

The remainder of this paper is organized as follows: Section 2 introduces the proposed algorithm. Section 3 provides a theoretical analysis of the algorithm, including the derivation of its steady-state mean MSD expression and an analysis of its computational complexity. Section 4 verifies the theoretical analysis and provides simulation results to demonstrate the superiority and robustness of the proposed algorithm. Section 5 concludes the paper.

2. The Proposed α-CSAP Algorithm

As illustrated in Figure 1, the system identification model considered in this work is given by

$$[eqn]$$

where $[eqn]$ denotes the unknown system weight vector with $[eqn]$ being the filter length, $[eqn]$ is the complex-valued additive noise, and the input vector $[eqn]$ is defined as $[eqn]$ . Here, $[eqn]$ denotes the Hermitian transpose. The output of the adaptive filter is expressed as $[eqn]$ ; here, $[eqn]$ is the adaptive weight vector at time index $[eqn]$ . Accordingly, the instantaneous error signal is given by

$$[eqn]$$

In the AP algorithm framework, the projection order Q specifies the number of past input–error pairs used for coefficient adaptation. The stacked error vector is defined as $[eqn]$ and the input data matrix is constructed as $[eqn]$ . In order to enhance the convergence rate of the algorithm and effectively mitigate the adverse effects of non-Gaussian noise on the adaptation process, the $[eqn]$ -sigmoid function is employed as the cost function, as follows:

$$[eqn]$$

where the parameter $[eqn]$ governs the steepness of the cost function, thereby determining the convergence speed of the algorithm: a larger $[eqn]$ leads to a steeper function and faster convergence. The parameter $[eqn]$ adjusts the algorithm’s sensitivity to outliers. From a mathematical perspective, as the error $[eqn]$ increases, the curvature of the cost function gradually flattens, and its derivative with respect to the weights approaches zero, thereby achieving robust suppression of outliers.

Taking the Wirtinger derivative of $[eqn]$ with respect to the complex conjugate of the weight vector $[eqn]$ yields

$$[eqn]$$

where $[eqn]$ is a diagonal weighting matrix defined as

$$[eqn]$$

with

$$[eqn]$$

Under the gradient descent framework, a preliminary update expression for the algorithm is formulated as follows:

$$[eqn]$$

where $[eqn]$ is the step size.

In the derivation of the $[eqn]$ SAP algorithm, the introduction of the normalization factor $[eqn]$ is a key step to realize the $[eqn]$ -norm constraint. The squared Euclidean norm of the weight update vector is expressed as follows:

$$[eqn]$$

As can be observed from Equation (8), when the noise contains large outliers, the error $[eqn]$ will be significantly affected, which, in turn, substantially impacts the weight estimation and degrades the filter performance. The restriction of the filter update magnitude can effectively enhance the robustness of the filter. This can be expressed as follows:    

$$[eqn]$$

where $[eqn]$ is a positive scalar. Taking the partial derivative of Equation (6), and according to Equation (8), the diagonal elements of the weighting matrix $[eqn]$ have an upper bound $[eqn]$ , which leads to the following:

$$[eqn]$$

At this point, by introducing a time-varying parameter $[eqn]$ into the denominator of Equation (8) that shares the same order as the numerator, the upper bound for the squared norm of the squared Euclidean norm of the weight update vector is obtained, as follows:

$$[eqn]$$

The update formula with constraints is expressed as follows:

$$[eqn]$$

The proposed $[eqn]$ -CSAP algorithm is summarized in Algorithm 1. The introduction of the normalization factor can be understood from an implicit variable step-size perspective. In scenarios with large errors, the effective step size $[eqn]$ automatically decreases, thereby suppressing impulse noise interference and preventing filter divergence. In scenarios with small errors, the effective step size relatively increases, accelerating the convergence process. Moreover, traditional affine projection algorithms require the inversion or pseudo-inversion of the input matrix, which is prone to numerical instability due to ill-conditioned matrices. The normalization factor avoids matrix inversion operations, significantly reducing computational complexity. Algorithm 1: The $[eqn]$ -CSAP Algorithm

• **Require: **  $[eqn]$ , $[eqn]$ , $[eqn]$ , Q, $[eqn]$ , $[eqn]$ , $[eqn]$ , N

• **Ensure: **  $[eqn]$

• 1:Initialize $[eqn]$

• 2:for  $[eqn]$ to N do

• 3:    Construct: $[eqn]$

• 4:    Construct: $[eqn]$

• 5:    Compute: $[eqn]$

• 6:    Compute: $[eqn]$

• 7:    Form: $[eqn]$

• 8:    for  $[eqn]$ to Q do

• 9:         $[eqn]$

• 10:    end for

• 11:    Construct: $[eqn]$

• 12:     $[eqn]$

• 13:end for

• 14:return  $[eqn]$

3. Theoretical Analysis

3.1. Unbiasedness in the Mean

This subsection employs a classical method from mean value analysis to conduct an approximate analysis of the mean unbiasedness of the proposed $[eqn]$ -CSAP algorithm under reasonable assumptions.

Define the weight error vector as

$$[eqn]$$

where $[eqn]$ denotes the unknown system weight vector. Substituting (12) into the above equation yields the update formula with constraints, which is expressed as:

$$[eqn]$$

Under the affine projection framework with projection order Q, let the stacked desired signal vector be $[eqn]$ and the stacked noise vector be $[eqn]$ . Thus, the stacked error vector can be expressed as

$$[eqn]$$

Substituting the above expression into the normalized update rule of the proposed $[eqn]$ -CSAP algorithm yields the following weight error recursion:

$$[eqn]$$

To facilitate the mean analysis, the following assumptions are adopted.

(A1) The input matrix $[eqn]$ is statistically independent of the weight error vector $[eqn]$ .(A2) In steady state, the diagonal weighting matrix $[eqn]$ varies slowly and can be approximated by

$$[eqn]$$

where $[eqn]$ .(A3) $[eqn]$ is weakly correlated with $[eqn]$ and $[eqn]$ , while $[eqn]$ is zero-mean and independent of both.

Under this assumption, it follows that $[eqn]$ . Then, taking the expectation on both sides of (16) yields

$$[eqn]$$

When the input signal is persistently exciting, the matrix inside the expectation is positive semi-definite. Therefore, for a sufficiently small step size $[eqn]$ , the only equilibrium point of (18) is

$$[eqn]$$

This result indicates that the proposed $[eqn]$ -CSAP algorithm is unbiased in the mean.

3.2. Steady-State MSD Analysis

3.2.1. Theoretical Value Derivation

In this subsection, the steady-state MSD performance of the proposed $[eqn]$ -CSAP algorithm is analyzed. Due to the presence of nonlinear error weighting and normalization terms, the analysis is carried out using the energy conservation approach commonly adopted in the mean-square analysis of AP algorithms.

The MSD is defined as

$$[eqn]$$

From the normalized update rule of the proposed algorithm, the weight error vector evolves as

$$[eqn]$$

where

$$[eqn]$$

Taking the squared Euclidean norm on both sides of (21), the expectation yields

$$[eqn]$$

Under steady-state conditions, $[eqn]$ , the MSD satisfies the following energy balance relation:

$$[eqn]$$

Define the a priori error vector as

$$[eqn]$$

Then,

$$[eqn]$$

Since $[eqn]$ is a diagonal matrix with entries $[eqn]$ , the numerator in (26) can be expanded as

$$[eqn]$$

Following the standard mean-square analysis of APs, it is assumed that under steady-state and small-step-size conditions, the correlation between the numerator and the normalization factor in (26) is weak. Thus, the following separation approximation is adopted:

$$[eqn]$$

Applying this approximation to (26), and noting that all Q components are statistically equivalent in steady state, yields

$$[eqn]$$

Define

$$[eqn]$$

Under steady-state conditions, the diagonal weighting function varies slowly and can be approximated by its average value, i.e.,

$$[eqn]$$

Therefore, the above expression reduces to

$$[eqn]$$

For the complex-valued error $[eqn]$ , where $[eqn]$ is a proper complex Gaussian random variable and z denotes Bernoulli–Gaussian impulsive noise, the weighting function operates on the magnitude of e only. By invoking Price’s theorem under the complex-valued setting, the following approximation holds:

$$[eqn]$$

Moreover, under steady-state conditions,

$$[eqn]$$

Combining the above results yields

$$[eqn]$$

The second-order term in (24) can be written as

$$[eqn]$$

Under steady-state and small-step-size conditions, this term can be approximated as

$$[eqn]$$

which is upper bounded by $[eqn]$ due to the boundedness of the weighting function.

Using the relations

$$[eqn]$$

the steady-state MSD is finally obtained as

$$[eqn]$$

The above expression provides an accurate approximation of the steady-state MSD.

3.2.2. Theoretical Validation

To verify the accuracy of the theoretical steady-state MSD model of the proposed $[eqn]$ -CSAP algorithm, the steady-state MSD expression (39) derived from the theoretical analysis was validated through simulations under a projection order of $[eqn]$ , with the step-size parameter $[eqn]$ varying within the range $[eqn]$ . The experimental setup was as follows: the filter length $[eqn]$ was set to three values, 16, 24, and 32, each matching the order of the unknown system; the number of iterations was $[eqn]$ , which was sufficient for the algorithm to reach steady state; and the results were averaged over 100 Monte Carlo trials. The input signal was a normalized complex-valued white noise, and the noise environment was modeled as Bernoulli–Gaussian mixture noise with an impulse probability of $[eqn]$ . Figure 2, Figure 3 and Figure 4 show the influence of the step-size parameter $[eqn]$ on the steady-state MSD for different filter lengths, with the projection order fixed at $[eqn]$ .

When $[eqn]$ (Figure 2), the theoretical MSD curve is slightly higher than the simulation data points, with an average error of 1.61 dB. As the filter length increases to $[eqn]$ (Figure 3), the theoretical and simulation results are in excellent agreement, verifying the accuracy of the derived expression under this configuration. When the filter length is further increased to $[eqn]$ (Figure 4), the predicted MSD values are slightly lower than the simulation results, with an average error of 1.35 dB. The slight deviations between the theoretical model and simulation results mainly originate from the necessary mathematical approximations adopted in the analysis, including the use of the separation assumption to handle the correlation of nonlinear terms, the steady-state averaging of the weighting function, and the approximate assumptions when applying Price’s theorem in the complex domain. Although these simplifications introduce systematic errors, they ensure the analytical tractability of the theoretical expressions. Experimental results demonstrate that the proposed theoretical model still accurately captures the overall trend of the steady-state MSD with respect to the step size, confirming the effectiveness of the analytical framework.

3.3. Computational Complexity Analysis

In this section, we conduct a comprehensive computational complexity analysis of the proposed $[eqn]$ -CSAP algorithm. The operation counting baseline follows complex arithmetic rules: complex addition requires two real additions, and complex multiplication requires four real multiplications and two real additions. For a length-n complex inner product, $[eqn]$ real additions are required. The exponential function $[eqn]$ , approximated using a second-order Taylor expansion, is estimated to need five multiplications and three additions.

Table 1 details the computational complexity of the $[eqn]$ -CSAP algorithm per iteration, where $[eqn]$ is the filter length and Q is the projection order.

Table 2 compares the computational complexity of different complex adaptive filtering algorithms. The algorithms under comparison include the Complex Least Mean Square (CLMS), the Complex Norm Least Mean Squares (CNLMS), the Complex Maximum Correntropy Criterion (CMCC), and the Complex Affine Projection (CAP) algorithms.

In terms of optimization potential, by employing fast recursive filtering (FRF) and an approximation of $[eqn]$ requiring only three multiplications, the complexity of $[eqn]$ -CSAP can be reduced to $[eqn]$ multiplications and $[eqn]$ additions.

The proposed $[eqn]$ -CSAP algorithm provides a good trade-off between computational complexity and performance. Although it requires more operations than LMS-type algorithms ( $[eqn]$ ), it avoids the $[eqn]$ matrix inversion of CAP and remains robust to impulsive noise and colored signals. In engineering applications, $[eqn]$ is often much larger than Q, so the complexity of the $[eqn]$ -CSAP algorithm is dominated by the $[eqn]$ term, making it suitable for real-time implementation.

4. Simulation

4.1. System Identification Application

This section evaluates the system identification performance of the proposed $[eqn]$ -CSAP algorithm under various complex-valued non-Gaussian noise conditions, in comparison with several representative complex-valued adaptive filtering algorithms, including CLMS, CNLMS, CAP, and CMCC. The unknown system $[eqn]$ to be identified is represented by a complex impulse response vector of length $[eqn]$ , whose real and imaginary parts are independently drawn from a Gaussian distribution with zero mean and a variance of 0.5. The input signal $[eqn]$ is a first-order complex autoregressive (AR(1)) process, generated by filtering complex-valued white Gaussian noise through a system function $[eqn]$ , resulting in strong temporal correlation. The projection order is set to $[eqn]$ . For a fair comparison, all algorithmic parameters are adjusted to achieve the same initial convergence speed. Simulation results are averaged over 100 Monte Carlo trials, and the MSD learning curve is adopted as the performance metric. The MSD at the n-th iteration is calculated as follows:

$$[eqn]$$

To thoroughly assess the algorithm’s robustness, the following three distinct and representative models for complex-domain non-Gaussian noise are employed:

• 1.Impulsive Mixture Noise with Bernoulli–Gaussian Structure: The composite noise $[eqn]$ is defined as $[eqn]$ . Here, $[eqn]$ denotes the background complex Gaussian noise of variance $[eqn]$ . The impulsive element is represented by $[eqn]$ , having a larger variance $[eqn]$ . A Bernoulli random variable $[eqn]$ , distributed as $[eqn]$ with $[eqn]$ , governs the occurrence of impulses: $[eqn]$ signifies no impulse, while its complement $[eqn]$ indicates an impulse event. All variables $[eqn]$ , $[eqn]$ , and $[eqn]$ are mutually independent.
• 2.Contaminated Complex Gaussian Noise (CG Noise): This noise model is synthesized by summing two independent zero-mean complex Gaussian components. The background component has a variance of 0.008 for the real part and 0.002 for the imaginary part. An intermittent impulsive component, with a variance of 8 for the real part and 2 for the imaginary part, is superimposed with an occurrence probability of 0.01. The combined signal exhibits a heavy-tailed distribution characteristic due to the occasional high-power impulses.
• 3.Symmetric Complex $[eqn]$ -Stable Noise ( $[eqn]$ -Stable Noise): This model adopts a complex-valued symmetric $[eqn]$ -stable distribution, formulated by $[eqn]$ . The real component $[eqn]$ and imaginary component $[eqn]$ are generated as independent and identically distributed real $[eqn]$ -stable random variables, described by the distribution $[eqn]$ . The configuration employs the following: characteristic exponent $[eqn]$ , symmetry parameter $[eqn]$ , scale parameter $[eqn]$ , and location parameter $[eqn]$ .

Figure 5, Figure 6 and Figure 7 illustrate the MSD convergence curves of various algorithms under three typical non-Gaussian noise environments. As can be seen, under the adverse conditions of strongly correlated input signals accompanied by non-Gaussian noise, the conventional CLMS algorithm exhibits a clear divergent trend. In contrast, the CNLMS algorithm effectively overcomes the divergence issue of LMS caused by step-size mismatch in non-stationary or impulsive signals through input signal power normalization of the step size, thereby maintaining system stability. However, the CAP algorithm suffers from a significantly degraded convergence performance under the influence of non-Gaussian noise. Furthermore, even when compared with robust algorithms specifically designed to suppress impulsive interference, the proposed $[eqn]$ -CSAP algorithm still achieves the lowest steady-state MSD across all three noise scenarios. This is primarily attributed to the synergistic design of the $[eqn]$ -sigmoid cost function and the $[eqn]$ -norm constraint adopted in the $[eqn]$ -CSAP algorithm, which endows it with stronger anti-interference capability. In summary, the comparative results demonstrate that the proposed $[eqn]$ -CSAP algorithm exhibits superior robustness and a more accurate steady-state estimation performance when handling strongly correlated input signals and complex non-Gaussian noise.

4.2. Beamforming Application

This section evaluates the beamforming performance of the proposed $[eqn]$ -CSAP algorithm in impulsive noise environments. A hexagonal sensor array composed of 91 elements is employed, with an array radius set to three times the operating wavelength and an inter-element spacing of approximately half a wavelength. The operating frequency is 8 GHz, corresponding to a free-space wavelength of $[eqn]$ mm. Five independent narrowband signals with Quadrature Phase Shift Keying (QPSK) modulation are simulated, including one desired signal and four strong interferences. The direction of arrival (DOA) of the desired signal is at an azimuth of $[eqn]$ and an elevation of $[eqn]$ . The DOAs of the four interfering signals are $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ , respectively. Herein, the azimuth angle $[eqn]$ is defined clockwise with respect to the array normal, and the elevation angle $[eqn]$ denotes the angle between the incident signal and the array plane.

To simulate a realistic strong interference environment, high-intensity impulsive noise is introduced, which is modeled by a Bernoulli–Gaussian mixture distribution. The impulse occurrence probability is set to $[eqn]$ , and the variances of the corresponding Gaussian components are $[eqn]$ and $[eqn]$ , respectively. The Signal-to-Noise Ratio (SNR) and Interference-to-Noise Ratio (INR) are set to 20 dB and 40 dB, respectively.

The core parameters of the $[eqn]$ -CSAP algorithm are systematically tuned, as follows: the adaptive filter length $[eqn]$ equals the number of array elements (i.e., 91), the projection order is $[eqn]$ , and the step size is set to $[eqn]$ . Furthermore, the characteristic exponent $[eqn]$ and the scale parameter $[eqn]$ of the algorithm are configured as $[eqn]$ and 16, respectively. This configuration significantly enhances the stability and robustness of the algorithm under impulsive noise. For a fair comparison, the step sizes of the conventional algorithms are also adjusted accordingly, as follows: $[eqn]$ for the CLMS algorithm, $[eqn]$ for the CAP algorithm, and $[eqn]$ for the CMCC algorithm. This ensures that all algorithms exhibit similar initial convergence speeds under identical simulation conditions.

As shown in the beam pattern results in Figure 8, while maintaining the same mainlobe gain, the $[eqn]$ -CSAP algorithm achieves deeper nulls at the four interference directions of $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ , with null depths exceeding those of the other compared algorithms.

We consider the signal-to-interference-plus-noise ratio (SINR) as a metric to evaluate the performance of the beamforming system, which represents the ratio of the desired signal power to the interference-plus-noise power. As shown in Table 3, under the harsh condition of an input SINR as low as $[eqn]$ dB, the $[eqn]$ -CSAP algorithm elevates the output SINR to $[eqn]$ dB, achieving a remarkable improvement of $[eqn]$ dB. This performance enhancement significantly surpasses that of the other compared algorithms. These results fully demonstrate the superior capability of the $[eqn]$ -CSAP algorithm in interference suppression and overall beamforming performance.

5. Conclusions

This paper proposes a robust complex $[eqn]$ -CSAP suitable for non-Gaussian noise environments. By extending the affine projection algorithm based on the $[eqn]$ -Sigmoid cost function to the complex domain, the algorithm enhances robustness against non-Gaussian noise. The computational complexity is reduced by employing a normalization factor to avoid matrix inversion. In two typical application scenarios, system identification and beamforming, the proposed $[eqn]$ -CSAP algorithm demonstrates superior convergence characteristics and lower steady-state error compared to classical complex adaptive filtering algorithms such as CLMS, CAP, and CMCC. The results validate its effectiveness and engineering potential in complex-domain signal processing.