A Performance Evaluation of Signal Filtering with different Adaptive Algorithms on Doppler-Angular Adaptive Generalized Sidelobe Canceler (GSC)

In modern radar signal processing the Space-Time Adaptive (STAP) filtering establishes a significant part of statistical signal processing. Whenever there is requirement to process signals that result from operation in an environment of unknown statistics, the use of an adaptive filter offers an attractive solution to the problem as it usually provides a significant improvement in performance over the use of a fixed filter designed by conventional methods. In this paper we study the use of three different adapting algorithms for adaptive generalized sidelobe canceler of interference signal suppression which intercepted by a uniform linear phased array antenna. Results are plotted, discussed and compared to each other. The paper is divided into four main sections. The problem is defined in the introduction section. In the next section, the adaptive spatial-temporal processor is depicted as well the autocorrelation matrix R has been driven. Next, the filtering using the normalized least mean-square (NLMS) algorithm was discussed. Also, we study the use of normalized fractionally low-order moment algorithm (NFLOM) for adapting of the same uniform linear array antenna and it was explained in section 6. It was compared to the NLMS algorithm for rate of convergenc, minimum mean-square error and angular –doppler beam pattern. In the last section, The recursive least-square (RLS) algorithm was implemented to get the adaptive filtering weight vector and compared to previouse algorithms. Finally, we’ve concluded all the results discussed in each algorithm.


I. INTRODUCTION
Appling a temporal filtering and spatial beamforming combination on the signals intercepted by phased-array radar systems are denoted as Space-time adaptive processing (STAP). It utilizes an antenna composed of array of elements afterwards to obtain multiple pulses processed coherently it is passed through tapped-delay-lines, therefore it provides a greater capability to deduce the influence of jammers and clutters on the system while preserving desired signal target [1]. A thoroughly researching on STAP was employed and it is revealed to deliver gains of remarkable performance in target detection and interference suppression [2]. Various STAP processors manage frequent situations clutters, noises and jammers are complex Gaussian, that drives to solutions of tractable mathematics [4,24,25].one of significant solutions is the adaptive beamforming.
At different scenarios, ranging from radar and wireless communications to sonar is found to be applications for the adaptive beamforming technique. It is founded on a system where, an array of elements weights at attached filters are adjusted to achieve a selective spatial and spectral.  Attention here to linearly constrained minimum variance (LCMV) beamformer is directed, because of its efficiency to be employed as a generalized sidelobe canceller (GSC) which has an unconstrained optimization problem as an alternative of a constrained one. Instead of employing beamforming, techniques of adaptive filtering are applied at the GSC adaptive process structure to perform a deep interference nulling and high resolution by using multiple sensors placed in arrays and filter weights.
In next section the optimum processor for wideband interference cancellation is revised in terms of the standard direct processor system with a covariance matrix is known. The corresponding full-dimension generalized sidelobe canceller (GSLC) structure [3,4] of the processor is derived in Section 3. Two different methods referring to commonly International Journal of Engineering Research & Technology (IJERT) ISSN: 2278-0181 http://www.ijert.org used algorithmic families have been implemented for the adaptively estimated MSE optimized parameters based on the available set of data. On the basis of steepest descent stochastic approximation method which is known as the Least Mean Squared (LMS) family the first approach is utilized in Section 6. The latter is the Recursive Least Squares (RLS) approach which is referred as the Gauss-Newton stochastic approximation.
The equivalence between the the Constrained RLS (CRLS) and GSC-RLS was verified by Werner [5] leads to the GSC form by employing the Recursive Least Squares (RLS) algorithm, which is called GSC-RLS, was derived in Section 8. This paper is organized as follows. Section 2 describes Space-Time Adaptive Processing (STAP). Section 3 presents GSLC Processor Form. Section 4 presents Design of Blocking Matrix. Section 5 presents GSLC adaptation using the NLMS. Section 6 presents Simulation results of NLMS-GSC and NFLOM-GSC. Section 7 presents GSLC adaptation using proposed recursive least square (RLS) algorithm for STAP. Section 8 presents Simulation Results of RLS-GSC algorithm performance.

II.
Space-Time Adaptive Processing (STAP) Signals intercepted by the antenna array are represented as 3D data of collected coherent processing interval (CPI) waveform at array elements. The data at the only desired range gate is then handled, which is represented as a CPI data slice. This part of CPI data contains a snapshot of NxM space-time data that consist of different components related to the information from the ith (i=1,2, ... ,N) element and jth (j=1,2, ... ,M) tap. Therefore, this two dimensional data structure of space-time comprises of tap information and space element. The samples is then stacked in a column to create vector of NM×l-dimensional. If a range cell of interest contains the target, then the received signal is formed of elements due to the target, white Gaussian noise, and the field interference signals (jammers sources and clutter): If the target does not exist, then the sample received will only consist of white noise, and interference. Then vector is: The covariance matrix of the input will be The NM×l-dimensional space-time steering vector Where ( ) is the steering temporal vector of M×1 dimension at the frequency of target signal and ( ) represent the target spatial steering vector of N×l-dimension at the direction delivered by spatial frequency of the targt vt. The operator of Kronecker product is noted by (.) ⊗ (.). Conveniencly, the steering vector in the space-time direction of interest is normalized and set to be The SINR unconstrained optimization provides vector that contains the optimum weights will be [6,7] = −1 (6) Considering the weight vector in (6) is providing target detection by utilizing a constrained optimization under assumption of the joint Gaussian. To maximize the output SINR, the vector of weights under optimal constrains that maintain the response of the space-time snapshot at targetdirection normalized, was originally given in [8,3] as following: The input data is applied to this weight vector, the output y is then formed as = = −1 −1 (8) and the output SINR be Where, | | 2 is a target waveform output power, and is the processor output power of the noise and it is given by As stated above the constrained optimization leads to the weight vector in (7). This weight vector is resulting in a beam, which is located in a specific direction of looking, spatially and temporally that is provided by the vector of steering s. That spatial-temporal processor structure is showed in Figure 2. The adaptive tapped delay line architecture required to perform STAP is shown in Fi. The weights w (adaptive vector), are employed at every one tap of the taps conducted at each spatial channel prior to the Final beamforming summation. Thus, N×M is the total adaptive coefficients number to be considered where N spatial Channels with M taps per channel are utilized. Note; τ is the taps in-between delay of time.
The minimum variance distortionless response (MVDR) scheme is usually applied for W the adaptation vector of weights in the Gaussian clutter environment for the output signal of second-order-moment minimization which is subjected to steering constraints Where the expectation operator used is E {·}, E {|y(t )| 2 }=W H Ruu W, and Ruu is the input vector U covariance matrix. S is the matrix of linear constraints set and h is the vector of desired response. Assume that, a constraint of a simple point can be selected as S = b( ωt )⊗ a( Θt ) and h=1, where a response of a unity gain is applied at the frequency Doppler ft and at the position of the target Θt. Assuming a full rank for the R (matrix of covariance)then the constrained minimization problem (11) has a well-known optimal solution [10]: Implementing (12) directly requires the knowledge of the input vector covariance matrix of the elements array and the sample-matrix-inversion (SMI) scheme is usually applied [2]. Instead, the decomposition of the vector of weights Wo into two modules orthogonally is employed: one is a vector of unconstrained adaptive weights Wa and the other is a fixed beamformer Wq. This decomposing type is recognized as the generalized sidelobe canceller (GSLC) where Wa could be adapted in iterative way by any least-square adaptive methods. Next section will point out this GSLC processor form and it's relying on the output signal to the desired one mean-square error optimization.

III. GSLC PROCESSOR FORM
The processor structure of the STAP in a direct way is reformed to a GSLC processor type. A GSLC leads to an unconstrained weight vector, where processor form of the STAP is transformed to a Wiener filter form. The GSLC processor structure delivers a deeper insight into problems of reducing dimension otherwise these problems appears to be possible. Considering operator "T" of the unitary matrix of dimension NM x NM as transformation operator of the detection priorly data. This operator structure is decomposed as follows: Where "s" is the conventional space-time beamformer of a dimension NM × 1 that defined in (5) and Bm is signal C × NM dimensional blocking matrix with full-raw rank which map the data input "u" into the null space of space-time beamformer "s", where C = (NM-1) Hence, Bm s = 0 (14) The waveforms intercepted at the spatio-temporal direction of interest is effectively blocked by the matrix Bm. Considering matrix Bm of a full-row rank to be a valid matrix of signal-blocking which fulfills (14) and provide T to be an invertible. such reformulation results in an output SINR is the same as the direct-form of the original processor of detection.

VIII. DESIGN OF BLOCKING MATRIX
The Blocking Matrix (BM) structure, Bm acts significantly in the structure of the GSC, the computational complexity depends on its choice and in many cases, it strengths the overall system numerical stabilities [11]. Determining of a blocking matrix Bm and the proper quiescent vector s which satisfies the constraints required for the GSC structure. The suitable BM design can be acquired through employing one of two algorithms for directly obtaining an orthonormal Bm from the desired normalized steering vector s. The matrix Singular Value Decomposition (SVD) is invoked by the first algorithm which is more accurate. The SVD approach may be employed for any constraint matrix C and is not limited to broadside limitations case. Hence, it is a general strategy for finding BM [12]. The other approach utilizes the QR-decomposition, and is more rapidly. Each of these approaches are explained as well as is derived in depth in numerous numerical and linear algebra texts [13,14]. A third algorithm has been proposed to formulate the BM is the cascaded columns of differencing (CCD) [15].

a) Cascaded-Columns Differencing (CCD)
The proposed algorithm of CCD is offered to get the derivative constraints blocking matrix. The durability on front of the direction of interest signal errors is increased by these constraints through enlarging the targeted constraints' angular range. As the order of the derivative constraints is getting higher, the wider the beam pointing in the desired direction [15].At CCD approach, a blocking matrix is designed by S cascaded columns differencing operations of as in Fig [12]. The BM is formulated in a matrix form as follows: Where we have By i = M,iM-1.i.i.iM-S-1. It's clear that, if a waveform returns from a wide side, it won't be capable of passing this BM. The BM in broadside forms a zero response and with increasing S a lobe of wider width is obtained [12].
Much of the time, it may be easy to construct a nonorthonormal signal blocking matrix and the orthonormalization can be provided by implementing the procedure of Gram-Schmidt. One case of this is a system, which the angle and time array manifold is set by utilizing a primary separate vector steering of space-time. Accordingly, in regard to the array, the s normalized steering vector is engaged to a broadside point. In such cases, an appropriate input signals difference is simply taken to provide the nonorthonormal matrix of signal blocking which is the same as the array configuration for arrays of narrowband offered in [16].
The radar incoming signal u is transformed using the operator T in (13) incomes a vector ̃, which takes the form Where v is termed as the scalar-valued beamformed output.
Here, also the vector b with dimension (NM-1) is denoted as data vector of the noise subspace. Then vector ũ that provide the data transformed has its related covariance matrix R ũ , and it is given as: Where is the covariance of the noise subspace which is of dimension C x C is given by The Cx1-dimensional vector r bv is the output of the beamformer cross-correlation to data vector of the interference subspace and is expressed by Where the operator of the complex conjugate is represented by "*". The beamformer output variance is denoted as the scalar σ d 2 in (18) that is set to be

= (21)
Next, let (the vector of steering) is transformed by the operator T. The unit transformed, steering vector e1, is given by this operation and it is given by Then, these transformed coordinates has the optimal weight vector which is given by (as w0) [ Where, ̃ is that covariance matrix expressed in (18). The GSLC processor structure showed in Fig [18 , 19] is normally derived as the separating of the matrix operator T is applied, this processor leads for the upper branch a unity fixed weight and a C=NM-1 dimensional vector w as an adaptive weight vector in the lower branch. Regarding to the configuration shown in Fig the vector w in (23) with optimal value is provided by the Wiener solution. Then the optimum weight w defined as Where, Rb and rbv are computed in (19) and (20), respectively.

Fig. 4. full dimention GSLC processor
The GSLC structure processor employs the N×M weight vector in (23) by applying the splitting given in (17). The GSLC performance at the steady state in the operation mode of adaptive data is identical to the direct-form processor performance, but a lower dimension vector of the adaptive weights w is used by the GSLC. Therefore, a reduction of the computations required for the weight vector updating is obtained and hence, the canonical form of processor can be the GSLC processor. Employing the weight vectors defined in (7) and (23) on the data arrays, leads an identical output, and therefor utilizing (17) and (23), it may be defined as follows: Substituting (24) into (25) and (10) is used to the evaluate the mean-square value of y; the minimum noise power of the output is obtained as follows:

. Principal coordinates form of fulldimension GSLC processor
These normal coordinates with its GSLC form is illustrated in Error! Reference source not found.. In the adaptive data processing mode, the characteristics of the steady state of both array processors are the same, however they may perform differently. The GSLC output array defined by: Annotate that the minimum noise power at the output.
and that is the SINR, is conserved by the operator U.

IX.GSLC ADAPTATION USING THE NLMS
The structure Displayed in Fig [20] involves of two channels that is identified as lower channel and upper channel. The upper channel usually denoted as the fixed beam former and an adaptive section along with blocking matrix form the lower channel. Performance loss is suffered by the conventional STAP algorithm under the severe, impulsive clutter environment as a result of 2 causes: first is the receiving of snapshots where the outliers' probability is high; second is the snapshot covariance matrix with large eigenvalue spread. The FLOM approach is An algorithm to overcome these difficulties. Rather than the STAP output's variance minimization this approach reduces the p th order-moment [21,22].
The resolving of optimum weights where the cost function is minimized has no closed-form, but a gradient descent algorithm is presented. The same as the NLMS method, the iterative adaptation is applied to the NFLOM algorithm as We can have the NLMS algorithm by reducing the NFLOM algorithm using p=2 and when p=1 the NFLOM algorithm is reduced to the algorithm of normalized sign (NSA).

XI.SIMULATION RESULTS OF GSLC ADAPTATION USING NFLOM ALGORITHM
The offered NLMS-GSC and NFLOM-GSC algorithms' evaluation is demonstrated by two typical adaptive spacetime processor examples. The first experiment simulation scenario is adopted as a stationary narrowband target signal, sT(n), had a 0 dB power relative to noise power at background . It has 0• as arrival angle(AoA) relative to array antenna elements. The target Doppler frequency that is normalized is fixed at 0.3. The signal of the target is diminished by two stationary narrowband source of jamming sj1(n) sj2(n), in different intercepted angle.
The parameter θ is referred as the plane waveform incoming angle (with respect to the broadside) and f denotes the normalized frequency of the received signal. The equally spaced linear array elements of N = 15 form the input channels of the adaptive space-time processor. The distance between every two elements are half of the wavelength related to fmax, which is navigated to the target waveform direction. For each antenna element 10-delay elements are assigned to it. M = 10 is the coherent pulse interval (CPI) and a fixed range gate were employed for the space-time adaptive processing. The mean square error e(n) of the actual target signal and the STAP output, sT(n), (i.e.,iminiE[|sT(n)−e(n)| 2 i]i) was computed at every case through estimated instantaneous squared errors' averaging over a window exponentially decayed with instants of time memory equal. The NLMS-GSC and NFLOM-GSC algorithms have the learning curves that are showed in Fig GSLC-RLS method is built on the structure of the GSLC and the variable weight vector is adjusted by applying the Recursive Least Squares (RLS) method. The GSC-RLS method is described in details later. Using the RLS method which is a recursive solution of the least squares, the adjusted weight vector of GSLC can be adapted as [5,10]: B m is the matrix block of the signal, that is consisted of any set of basis that are orthogonal for the directing array of the target u t . Therefore, we have B m u t = 0, In general, the blocking matrix can be straightly acquired by utilizing the QR decomposition or the singular value decomposition (SVD) algorithms [42] as discussed in earlier section.
Where (n ) is the wanted signal of the RLS method and ( ) is the input signal. e(n ) is the a priori error, ( ) is the gain vector, the same as ( ) ( ) , ( ) is the inverse of the ensemble-averaged autocorrelation matrix, λ is the forgetting factor" of the RLS method.

XIII. SIMULATION RESULTS USING ADAPTATION RLS ALGORITHM
Computer simulations are used to drive the numerical results that are obtained here. The RLS-GSC algorithm performance is demonstrated using linear phased array of N = 15 equally spaced elements. The distance between elements is half of the operation frequency related wavelength. Every antenna element has Seven delay components assigned to it M = 10 and a range cell of interest was applied to (STAP) processor. 0 dB is the power of the target waveform which is used, relative to the noise at background and it has an angle of arrival 0 • degree referenced to the axis of the array. The target has a fixed Doppler frequency of 100Hz. Independent noises through taps of the CPI and the array elements were applied which has spectrum of a white Gaussian. A wideband jamming from Two sources represented at direction of arrival of −20 • and +50 • and 30 dB is the power of these jammers totally had also they had a full Doppler spectrum. Along with the jammers and noise, the array is impinged by many clutters from different DoAs. The clutters are between 180 • and −180 • and have a uniform distribution. These clutters has Doppler frequencies related to their DoAs  The adaptive filtering using the LMS algorithm at its simple structure delivers double the unknown parameters number of the scheme as a cost. Though, the algorithm has a rate of convergence heavily dependent on the spread of eigenvalue of the input correlation matrix. However, such a drawback is not suffered when the RLS-GSC algorithm is utilized but its complication of calculation is related to the squared of the unknown factors' number proportionally in the system. As the real time processing has to be achieved when using STAP, the algorithms of minimum cost of computations must be engaged. The built NFLM-GSC structure provides a lower cost of adapted signal processing; therefore, it is closer for adaptive space-time processing task implementation at real time.