# Blind Source Separation for Complex-Valued Signals Using Generalized Autocorrelation

Authors: Xiaogang Tang; Sun’an Wang; Jiong Li
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

We introduce a new complex-valued blind source separation approach, based on generalized autocorrelations of sources, to improve the spectrum efficiency for the next-generation wireless communications system. The proposed algorithm considers the temporal structures of communication signals and the natural gradient-based method is used to optimize the demixing matrix. In addition, the local stability condition is proved. Simulation results are presented showing the superior performance of the proposed algorithm in the intersymbol interference of the estimated signals.


---

## Body

## 1. Introduction

For the next-generation mobile communications system (5G), spectrum efficiency is a critical performance metric, which aims at achieving more than 10 times spectrum efficiency compared with the current communications system (4G) [1]. Since radios must either transmit or receive on the same channel, but not simultaneously, in previous wireless communications systems, the frequency resource is not utilized sufficiently. Fortunately, an emerging technique, co-frequency and co-time full-duplex (CCFD) [2–4], is able to address this issue. CCFD enables radios to transmit and receive signals on the same channel simultaneously and, thus, theoretically, can double the spectrum efficiency. CCFD technique claims to be the most potential duplex scheme for the 5G network. But, there is a challenge which lies in the application of CCFD technique, i.e., mitigating the local self-interference (SI) [3, 4]. Since the transmit and receive antenna work on the same frequency band, traditional interference cancellation techniques are invalid. As our previous work shown [5], blind source separation (BSS) [6] has an excellent capability to address such issue.Since BSS is able to estimate original signals from their observed sensors’ signals without knowing both the mixing process and the sources, it has been one of the most attractive technologies in many fields, such as communications [7–9], speech signal processing [10, 11], and medical signal processing [12–14]. Depending on the application, the sources may be real or complex. This paper focuses on the separation of baseband communication signals which are usually complex-valued.Independent component analysis (ICA) is a common method to realize BSS, which can estimate demixing matrix and separate signals up to an arbitrary scaling and permutation ambiguities if the source signals are mutually independent and no more than one Gaussian source exists. And most of complex ICA algorithms exploit the non-Gaussianity or sample dependency of sources. Non-Gaussianity is often measured by high order statistics. Fourth-order moments is used in fourth-order blind identification (FOBI) [15], which is very efficient but requires that the sources have different kurtosis. However, sources have identical distribution in many applications, and then FOBI is invalid. Joint approximate diagonalization of eigenmatrices (JADE) [16] is based on fourth-order cumulants and performance suffers when the number of sources increases. The strong uncorrelating transform (SUT) [17] and the generalized uncorrelating transform (GUT) [18] estimate mixing matrix by jointly diagonalizing the covariance matrix and pseudo-covariance matrix of mixed signals. Also, it is important to note that, for noncircular sources, SUT can achieve separation of themselves as long as they have distinct noncircularity indices.The above methods make use of the non-Gaussianity property and neglect sample dependency. In this paper, we take the temporal characteristics of communication signals into account. We develop a new complex BSS algorithm motivated by generalized autocorrelations of the source signals. A natural gradient-based algorithm is derived for the generalized autocorrelations contrast function. Stability conditions are then derived in the complex domain for the generalized autocorrelations contrast function.This paper is organized as follows. First, the BSS problem is formulated in Section2. We then introduce the new contrast function based on the generalized autocorrelations of source signals and derive a natural gradient-based algorithm in Section 3. The performance of the algorithm is demonstrated with simulations in Section 4. Conclusions are drawn in Section 5.


## 2. Problem Formulation

Considering there areN sensors and N independent sources, the instantaneous linear mixtures of these sources are observed at the sensors:(1)xt=Astwhere st=s1t,s2t,…,sNtT (superscript T denotes transpose) is a vector of unknown zero mean and unit-variance source signals, xt=x1t,x2t,…,xNtT is the observation vector, and A is the unknown N×N mixing matrix. We impose another assumption on the source signals that the source signals have specific temporal structures, linear or nonlinear autocorrelations.In general, preprocessing operations of observed signals are needed before performing BSS algorithm. Two common preprocessing operations are removing mean and whitening. The whitening matrixQ can be obtained using the eigenvalue decomposition of covariance matrix Rx=ExtxtH (superscript H denotes conjugate transpose; E• denotes expectation).(2)Q=D-1/2VHwhere D and V are the eigenvalue matrix and the eigenvector matrix of the covariance matrix Rx. The observations are whitened by Q:(3)zt=QxtThen, Rz=EztztH=QARsQAH, due to Rs=I, Rz=I, and QA is an orthogonal matrix which imposes orthogonal constraint on the demixing matrix W. The sources can be estimated by(4)yt=Wztwhere yt=y1t,y2t,…,yNtT is the estimation of st.


## 3. Proposed Algorithm

Now we generalize the framework in [19] for complex-valued signals. We start from an arbitrary nonlinear contrast function and prove that its extrema coincide with the sources. The generalized autocorrelation of estimated signal is calculated with the selected nonlinear function and the nonlinear function can be chosen freely. Our contrast function based on complex generalized autocorrelation functions is(5)Jwn=EGwnHzt2GwnHzt-τ2where G is a differentiable function which measures the generalized autocorrelation degree of the source signal, wn is the nth volume vector of demixing matrix W, and wn2=1, τ is a delay in time. Finding the extrema of a contrast is a well-defined problem only if the function is real. So we make our contrast functions operate on absolute values rather than complex values. Examples of choices of the differentiable function G are G1u=u, G2u=u2, G3u=log⁡cosh⁡u. For simplicity, the time index t is omitted in the following, i.e., zt=z, zt-τ=zτ.Theorem1 gives the local stability of the estimators, the proof of which is given in Appendix.Theorem 1.
Assume that the input data follows model (1). The prewhitening operation is performed on the observation signals as (3). Furthermore, assume that sn,snτ and sl,slτ are mutually independent and Esn2=1, Esn2=0. Then, the local maximum (respectively, minimum) of Jwn under the constraint wn=1 includes those rows of the inverse of the global mixing matrix QA, so the corresponding sources satisfy(6)Eα3-α1+2Resnsnτs1∗s1τ∗+snsnτ∗s1∗s1τgs12gs1τ2<0resp.>0where we define(7)α1=Es12gs12Gs1τ2+s1τ2Gs12gs1τ2(8)α2=2Es14g′s12Gs1τ2+2s12s1τ2gs12gs1τ2+s1τ4Gs12g′s1τ2(9)α3=Egs12Gs1τ2+Gs12gs1τ2+s12g′s12Gs1τ2+s1τ2Gs12g′s1τ2where function g′ is the derivative of g.
We now give the conjugate gradient algorithm for complex signals under model (1). The problem given in (5) can be written as the following Lagrangian function:(10)J0wn=Jwn+λwn2-1where the Lagrangian multiplier λ is a real number. The conjugate gradient of the contrast function J with respect to wn can be obtained as(11)∂J0wn∂wn∗=∂Jwn∂wn∗+λwn(12)∂Jwn∂wn∗=EgwnHz2GwnHzτ2wnTz∗z+GwnHz2gwnHzτ2wnTzτ∗zτFunction g is the derivative of G; z∗ designates the complex conjugate of z. Thus, the conjugate gradient-based update rule of wn can be written as(13)wn=wn-μ∂Jwn∂wn∗-RewnH∂Jwn∂wn∗wnwn=wnwnwhere μ>0 is the real-valued step-size.
It is known that, in general, relative gradient-based method converges slower than natural gradient-based method. To speed up the convergence, the natural gradient-based update rule of (13) can be written as(14)W=W-μ∂J0W∂WWHWW=WWHW-1/2


## 4. Simulation and Results

The proposed algorithm is compared with three competitive complex ICA algorithms: complex FastICA [20], EBM [21], and EASI [22]. FastICA is a Newton-based ICA algorithm which converges fast. EBM is a conjugate gradient-based algorithm. EASI is a relative gradient-based algorithm. In general, Newton method converges faster than conjugate gradient and relative gradient methods, and relative gradient method converges slowest. The performance index is defined as (15) which means the average intersymbol interference of the estimation sources.(15)PI=12N∑k=1N∑l=1NUkl2maxnUkn2-1+∑l=1N∑k=1NUkl2maxnUnl2-1where U=WQA is the combined separation-whitening-mixing matrix, and Ukl is the k,lth entry of U.This article considers the situation where two traditional and classic modulation signals are mixed, i.e., two 8PSK signals. Figure1 shows the constellation figures of two original 8PSK signals, their mixtures, and the estimations using the proposed algorithm employing the nonlinear function G1u=u. In this simulation, the two 8PSK are with the following parameters: symbol rate Rs=198kBps, raised cosine filter which has 0.5 roll-off factor, and sampling frequency 16Rs. The mixing matrix is randomly generated as(16)A=0.5559-0.4087i-0.8623+0.0475i-0.8571-1.1729i1.0418+0.9968iFrom Figure 1, we can see that the source signals are properly estimated. As most complex BSS algorithms, there are some ambiguities, including amplitude, order, and phase. The combined separation-whitening-mixing matrix is estimated as(17)U=-0.0012-0.0050i0.7640-0.6459i0.5893+0.8082i-0.0092-0.0167iFigure 1
Simulation result of the proposed algorithm for mixtures of two 8PSK signals, time delay,τ=1. From (a) to (c), the source signals 1-2, the mixtures 1-2, and the estimation signals 1-2.
(a)
Constellations of the source signals (b)
Constellations of the mixing signals (c)
Constellations of the estimation signalsIn order to exploit the influence of time delayτ to the proposed algorithm, the autocorrelation of the source signal is calculated ξτ=Esntsn∗t-τ with respect to different τ and the separation performance of the proposed algorithm with different time delay τ is simulated. Since the source signals are two 8PSK signals with the same modulation parameters, they have the same autocorrelation function. The autocorrelation function is shown in Figure 2 and the separation performance of the proposed algorithm is given in Figure 3. Form Figures 2 and 3, we can find that the separation results remain at a certain level when time delay τ<10 and then ξτ>0.5. For proposed algorithm with the three nonlinear functions, the proposed algorithm with G1u is sensitive to the time delay τ if τ>10; the influence of the time delay to the proposed algorithm with G3u is less than that with G2u. In the rest simulations, we set the time delay τ=1.Figure 2
Autocorrelation function of the source signal.Figure 3
Separation performance of the proposed algorithm at different time delayτ.Figure4 shows the convergence speed of different algorithms. The sources are two 8PSK signals with the same parameters as shown previously and the mixing matrix is randomly generated. It is clear that the convergence speed of EASI is slower than the other algorithm since EASI is a relative gradient-based algorithm with linear convergence speed. FastICA and EBM have a similar convergence speed and the convergence value of FastICA is lower than EBM, which means the residual intersymbol interference of EBM is higher than FastICA. For the proposed algorithm, we can see that its convergence speed is slightly slower than FastICA and EBM. It is because that the proposed algorithm is a natural gradient-based algorithm while FastICA is based on Newton method. In addition, it also can be seen that the convergence speed and convergence values of the proposed algorithm using the three different nonlinear functions are similar, and the convergence values of the proposed algorithm are lower than the two comparison algorithms.Figure 4
Comparison of average PIs over 100 independent runs for the four algorithms.Figure5 shows the averaged PI values for the separation of two 8PSK sources with different sample size. The mixing matrix is randomly generated in every individual experiment. From Figure 5, we observe that, with the increase of sample size, the PI values decrease. The proposed algorithm using the nonlinear function G1u is the best algorithm and shows almost the same performance as the proposed algorithm using nonlinear functions G2u and G3u when the sample size is bigger than 1600. FastICA and EBM show a similar performance. The performance of EASI is similar to EBM when sample size is bigger than 1600 even better than FastICA and EBM when sample size is smaller than 1600. However, as shown in Figure 2, the convergence speed of EASI is so slow that it would not be the first choice to perform sources separation.Figure 5
Comparison of average PIs of the four algorithms over 100 independent runs at different sample size for 8PSK signals.


## 5. Conclusion

In this paper, we address the complex blind source separation problem using generalized autocorrelations of the sources with temporal structures. A natural gradient-based algorithm is proposed and the stability condition is proved. In the simulations, we perform the proposed algorithm on the mixtures of two 8PSK signals. The simulation results show that although the convergence speed of the proposed algorithm is slightly slower than complex FastICA and EBM algorithms, its convergence value is smaller than the comparison algorithms, which means the averaged signal-to-interference ratio of the estimated signals is higher than the comparison algorithms.



---

*Source: 8076468-2018-07-12.xml*