Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 484327, 12 pages

http://dx.doi.org/10.1155/2014/484327

Research Article

Blind Source Separation Based on Covariance Ratio and Artificial Bee Colony Algorithm

Lei Chen,¹ Liyi Zhang,1,2 Yanju Guo,³ Yong Huang,¹ and Jingyi Liang¹

¹School of Information Engineering, Tianjin University of Commerce, Tianjin 300134, China

²School of Electronic Information Engineering, Tianjin University, Tianjin 300072, China

³School of Information Engineering, Hebei University of Technology, Tianjin 300401, China

Received 19 November 2013; Revised 8 May 2014; Accepted 9 May 2014; Published 2 June 2014

Academic Editor: Dan Simon

Copyright © 2014 Lei Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The computation amount in blind source separation based on bioinspired intelligence optimization is high. In order to solve this problem, we propose an effective blind source separation algorithm based on the artificial bee colony algorithm. In the proposed algorithm, the covariance ratio of the signals is utilized as the objective function and the artificial bee colony algorithm is used to solve it. The source signal component which is separated out, is then wiped off from mixtures using the deflation method. All the source signals can be recovered successfully by repeating the separation process. Simulation experiments demonstrate that significant improvement of the computation amount and the quality of signal separation is achieved by the proposed algorithm when compared to previous algorithms.

1. Introduction

Blind source separation (BSS) is an advanced signal processing technique which aims to recover source signals from a received mixture with little prior information. It has received considerable attention in many research fields, such as, speech sound, image, communication, and biomedicine [1–4]. In the literatures, a lot of BSS algorithms have been proposed from different viewpoints, including non-Gaussianity [5], mutual information minimization [6], maximum likelihood [7], tensors [8], principle component analysis (PCA) [9], and neural networks [10] (for a detailed survey, see, e.g., [11–13]). Concerning BSS, the separation criteria and optimization method play the two most important roles and are detailed next. Separation criteria are proposed as the measure of separability and optimization method is used to find the optimum solution for objective function which comes from separation criteria.

Optimization methods based on stochastic gradient techniques are often used in conventional BSS algorithms, such as bigradient, back propagation, and natural gradient [14–16]. Nevertheless, these methods are inclined to fall into local optima unless suitable initial parameters and step length are given. Moreover, the nonlinear functions should be selected in some of these conventional BSS algorithms and it will also affect the separation precision.

The bioinspired intelligence optimization (BIO) method is a research area inspired by natural evolution. Compared with conventional optimization methods, its advantages are distinguished, enabling a lot of research fields to adopt it into their studies. In recent decades, it has been widely used in many fields such as feature selection, power system, data mining, inventory management, biomedical engineering [17, 18], and the performance enhancement of these applications had been confirmed.

In recent years, BIO methods have been given great importance in solving the BSS problem. For instance, genetic algorithms (GA) [19–21], particle swarm optimization (PSO) algorithm [22–25], and artificial bee colony (ABC) algorithm [26, 27] have been applied to identify the separation matrix or separation vectors within the BSS problem. Compared with conventional gradient-based approaches, the BIO-based approaches for BSS are characterized by high accuracy, efficiency, and robustness.

However, in these BIO-based BSS algorithms that have been proposed, most of the objective functions are based on fourth-order cumulant (kurtosis). The kurtosis of the signal mixtures is required to be computed repeatedly in the evolution process of each generation. Therefore, the computation complexity of these algorithms is high. In order to overcome the computation difficulties in BIO-based BSS algorithms, we proposed a novel BIO-based BSS algorithms using covariance ratio and modified ABC algorithm.

Our work differs from the above researches in that we use different separation criteria. The covariance ratio is used in the separation process and the modified ABC algorithm is used for solving the objective function based on it. The proposed algorithm differs from the BIO-based BSS algorithms using kurtosis in a way that the fixed two covariance matrices based on mixtures are computed in advance and they will be used in the separation process and without kurtosis of all the sample points from mixtures being computed in the evolution process of each generation. As a result, comparing with BIO-based BSS algorithms using kurtosis, the proposed algorithm can achieve blind separation efficiently with a high separation precision and a low computation complexity.

2. Blind Source Separation (BSS) Model

Let denote independent source signal vector that comes from signal sources. We can get observed mixtures under the circumstances of instantaneous linear mixture. Consider the following:

where is a nonsingular constant matrix. The task of the BSS algorithm is to recover the source signals from mixtures . The separation model of BSS algorithms can be divided into two categories: batch approach and extraction approach.

(1) Batch Approach. The separation model can be represented as follows:

where is a separation matrix and is the estimate of source signals.

(2) Extraction Approach. The separation model can be represented as follows:

where is th separation row vector, , and is the estimate of one source signal. source signals can be recovered sequentially by conducting this separation, otherwise known as extraction, process repeatedly. The task of the two categories of BSS algorithms is to find the separation matrix or separation row vector , respectively.

2.1. Objective Function Based on Kurtosis

Foremost, to perform blind source separation, initially we must define the objective function in BSS context. Most of the objective functions being used in previous BIO-based BSS algorithms are constructed with the idea that the output sources must be independent from their linear mixtures. For this purpose, a measure of independence must be utilized. Kurtosis is the classical measure of independence and it has been used in some BIO-based BSS algorithms.

These BIO-based BSS algorithms can also be divided into two categories: batch approach and extraction approach. In [23], a BIO-based BSS algorithm based on adaptive particle swarm optimization is proposed. The sum of the kurtosis absolute values is used as the objective function in this algorithm. It can be represented as follows:

where is the kurtosis of . Subsequently, in [27], a modified ABC algorithm is used to optimize this objective function in (4). The two algorithms are all belonging to batch approach and all the source signals can be recovered synchronously.

For extraction approach, the objective function can be represented as follows:

The task of the BIO algorithm in BSS is to find separation row vector who can maximize and one source can be recovered. In [25, 26], PSO and ABC are used, respectively, to optimize the objective function in (5). All the source signals can be separated out efficiently via multiple separations.

These BIO-based BSS algorithms have performed separation with high precision and excellent robustness. Unfortunately, the common problem with these BIO-based BSS algorithms is that their computation amount is large due to the high order cumulants of all the sample points from mixtures being calculated in each evolution process of a BIO algorithm. Therefore, we will propose one type of BIO-based BSS algorithm with different separation criteria and frame, and then, the computation difficulties in BIO-based BSS can be conquered.

2.2. Objective Function Based on the Covariance Ratio

The principle of most BSS algorithms is based on the independent component analysis (ICA) model. Alternatively, we can study the BSS problem from various viewpoints. As an example, many signals in nature such as speech signal, music signal, mechanical vibration signal, and seismic signal have temporal continuity. Therefore, we can separate outsource signals from these mixed signals using temporal structure information.

In [28], a seminal conjecture for blind separation is proposed as follows: the temporal predictability of a signal mixture is usually less than that of any of the source signals that contribute to it. Therefore, the measurement is defined to evaluate the temporal predictability of . The definition of is as follows:

where is the number of the sampling points. and are the short-term and long-term moving average value, respectively, in th separation. Consider the following:

where, , , , , (typically 100 times longer). and are defined as

where ; (6) can be rewritten as follows:

where is a matrix of long-term covariance between signal mixtures and is a matrix of short-term covariance. The task of the optimization algorithm is to find the separation row vector , which maximizing the objective function , and then, one source signal can be recovered.

Subsequently, the theoretical basis of Stone’s conjecture is discussed in [29] that Stone’s conjecture should be modified as follows: the temporal predictability of any linear mixture of the source signals is between the greatest temporal predictability and the smallest temporal predictability of the source signals.

Recently, Liu and Cheung [30] presented a generalized metric of those algorithms [28, 29], which are constructed on the idea that the metric value of any linear combination of source signals is less than the largest one of sources under a loose condition. Accordingly, the generalized metric is a logarithm of a ratio of two covariances, in which the numerator is the covariance of ’s transform, whereas the denominator is the covariance of . In this separation frame, BSS has been formulated as global optimization of a covariance ratio.

The general metric can be presented as follows:

The operators and may vary in types. Different BSS objective function can be obtained from different operators being selected. In this paper, we will use modified artificial bee colony algorithm to optimize the objective function from (10) and one novel BIO-based BSS algorithm will be proposed.

3. Artificial Bee Colony Algorithm

3.1. Basic ABC Algorithm

ABC algorithm is one type of BIO algorithm proposed by Karaboga, which simulates the foraging behavior of a bee colony [31]. It has been shown to be competitive to other BIO algorithms [32], such as, genetic algorithm (GA), particle swarm optimization (PSO), and ant colony optimization (ACO).

In ABC algorithm, the position of a food source represents a possible solution of the optimization problem. The colonies of artificial bees are divided into two groups: employed bees and onlooker bees. The number of the bees colony is , the number of employed bees is , and the number of onlooker bees is . Under normal circumstances, .

At the initialization stage of ABC algorithm, let represents the th food source in the colony. And then, each food source is represented as follows:

where and . The employed bees are randomly assigned onto these food sources to evaluate their nectar amounts. The fitness value of the solution is proportional to the nectar amount of the food source in the position . Employed bees will search candidate food position in the neighbourhood of the old one. The position search equation is represented as follows:

where , has to be different from . By using a greedy selection mechanism, if the fitness value of is equal to or better than that of , will replace ; otherwise, is retained.

After all employed bees complete their neighbourhood searches, each onlooker bee will choose a food source depending on the probability value associated with that food source. The probability value with which is chosen by an onlooker bee can be calculated according to the following formula:

where is the fitness value of solution . Obviously, roulette selection strategy is used in it. Once the onlooker bee has selected one food source depending on the probability value , it will produce a modification according to (12). will be replaced by the modified food source if the modified food source has a better or equal nectar amount.

In ABC algorithm, if a position corresponding to one employed bee cannot be improved further through a predetermined number of cycles, the employed bee will become a scout bee and will start to produce a food source randomly according to (11). Meanwhile, the old food source is replaced by the new one.

To summarise, ABC algorithm employs three different steps: sending the employed bees onto their food sources and evaluating their nectar amounts; after sharing the nectar information of food sources, searching new food sources in the neighbourhood by the onlookers; determining the scout bees and then sending them randomly onto possible new food sources.

3.2. Modified ABC Algorithm

In order to enhance the optimization property of basic ABC algorithm, some efficient modification strategies have been proposed in [33–37]. The MABC algorithm [37] inspired by differential evolution (DE) is one simple yet efficient modified algorithm. In the BSS algorithm we proposed, the MABC algorithm will be used to optimize the objective function and the separation vector can be found. The two modifications in MABC algorithm can be briefly represented as follows.

(1) A Modified Search Equation Based on DE.  Differential evolution (DE) [38] strategy is an efficient evolutionary idea used in many hybrid optimization algorithms. In MABC, DE strategy is led into basic ABC algorithm to improve the convergence performance. Hence, the position search equation in neighbourhood is devised as follows:

 where the and are mutual integers randomly chosen from , and different from the base index ; is the best individual position with the best fitness value in the current colony; and . is a random number in the range . ABC/best/1 can drive the new candidate position only around the best solution of the previous iteration. Therefore, the proposed position search equation described by (14) can increase the exploitation ability of ABC.

(2) Initial Population Process. Population initialization is an important process in BIO algorithms because it can affect the convergence’s quality and speed of the final solution. In [37], instead of random initialization, chaotic maps have been used to initialize the population so that the search space information can be extracted to increase the population’s diversity and convergence speed. So the initialization approach based on chaotic system is used to generate initial population. Therefore, sinusoidal iterator is selected and its equation can be expressed as follows:

where is the iteration counter and is the maximum number of chaotic iterations. The mapped process can improve the ergodicity, randomness, and irregularity of the initial solution. The two modification strategies in MABC can prevent bees from falling into the local minimum and accelerate the evolution convergence process significantly. In the experimental studies [37], 28 benchmark functions were tested to show the superiority of MABC compared with other BIO algorithms.

Meanwhile, we will use the spherical coordinates transform in the initial population process according to information about the solution coming from objective function, and then, the convergence speed and search precision can also be improved. This modification will be explained in Section 4.

4. New BSS Algorithm Based on Covariance Ratio and MABC Optimization

4.1. Proposed Algorithm

In the proposed BSS algorithm, we will choose the operators and in the first place.

According to [39], let , . Equation (10) can be rewritten as follows:

Subject to the constraint . Next, the modified ABC algorithm is used to optimize the metric in (16), and then, all the source signals can be separated from mixtures sequentially. The proposed BSS algorithm will be described as follows.

Firstly, parameter coding process should be done according to the number of the unknown variables in the BSS algorithm. The variables needed to be found are the components of the separation row vector . That is, in blind separation problem with source signals, separation row vector is and the parameter coding for MABC is . The optimization problem can be represented as follows:

Subject to the constraint . Secondly, we use spherical coordinates transform to reduce the variable dimensions and search range of MABC according to the constraint. The transform principle can be represented as follows:

where ; set . Then, the objective function can be rewritten as follows:

The optimization problem for has been reformulated as optimization problem for angle vector . The variable dimensions decrease to and the searching range of MABC shrinks to . The initial population process will be carried out in this range according to chaotic maps in (15). As a result, the difficulties for searching artificial bee colony have been reduced and the searching time can be saved.

In the separation process, only two fixed covariance matrices are used repeatedly in evolution process of each generation, without high order cumulant (kurtosis) of all the sample points from mixtures being computed. Therefore, the computation amount of the separation is low.

In case of the same source signal being separated out repeatedly, deflation method is used to remove the source signal component separated out from mixtures. Let be the first signal separated using the algorithm we proposed and is the new mixtures with being removed from . The deflation process can be represented as follows:

There will be no component of in after the deflation process. It is only composed of . will be the new mixtures taking the place of in the next separation process and another source signal can be derived. Then, all the source signals can be separated out one by one through ongoing separation and deflation process.

4.2. Concrete Steps of the BSS Algorithm Proposed

Step 1: set the dimension and coding of artificial bee colony according to the number of source signals.

Step 2: initialize the population of artificial bees according to (15) in constraint range and they are divided into two sorts, which are employed bees and onlooker bees.

Step 3: compute the fitness value of the nectar source position located by each employed bee according to (18) and (19); seek out the best nectar position .

Step 4: each employed bee will search new nectar source positions in the neighbourhood of . If the better nectar source position is found before reaching the maximum search limit, replace the old nectar source position using the new one. Otherwise, abandon the old nectar source position and generate a new position according to (11).

Step 5: each onlooker bee will select one nectar source position according to the probability selection principle in (13) and search better nectar source position in its neighbourhood. If the better nectar source position is found before reaching the maximum search limit, replace the old nectar source position using the new one. Otherwise, abandon the old nectar source position and generate a new position according to (11).

Step 6: record the current best nectar source position . If the maximum evolution generation has not been met, return to step 4. Otherwise, transform , obtained from , into separation vector according to (18) and export separation signal according to (3).

Step 7: perform deflation process according to (20) and new mixtures can be obtained

Step 8: If all the source signals have been recovered, cease the computation process. Otherwise, return to step 1 and carry out the next separation process.

5. Experimental Results

To illustrate the performance of the proposed BSS algorithm, we report some experimental results on the separation of speech and music signals.

We compared the proposed algorithm with four BIO-based BSS algorithms presented in Section 2.1. We named them PsoBatch algorithm [23], PsoExtraction algorithm [25], AbcBatch algorithm [27], and AbcExtraction algorithm [26], respectively, according to their different combinations of optimization algorithm and separation category.

5.1. Source Signals and Parameter Settings

In the separation experiments, three source speech signals were selected from ICALAB for signal processing-benchmarks, taken from Cichocki Laboratory for Advanced Brain Signal Processing ICALAB Toolboxes website http://www.bsp.brain.riken.jp/ICALAB/. Three music source signals, piano, violin, and guitar sounds were selected from audiophile source website http://www.xrcd.com/.

The sensors and sources were set as 3 in the initialization. The elements of the 3 × 3 mixing matrix are random numbers distributed uniformly in the range . Consider the following:

The initialization values of the parameter for MABC in the proposed BSS algorithm were set as population size , , , , maximal generation for evolution ; the algorithm will be carried out in 30 independent runs. The parameter settings in comparison to algorithms were set as follows.

PsoBatch algorithm: population size , , , inertia factor linearly decreases from 0.9 to 0.4 and it will be adjusted adaptively with the change of fitness value, .

PsoExtraction algorithm: population size , , , inertia factor linearly decreases from 0.9 to 0.4 and .

AbcBatch algorithm: population size , , , , , , , and .

AbcExtraction algorithm: population size , , , , and .

5.2. Performance Measure

In order to measure the accuracy of the proposed algorithm, we evaluate it using three performance indexes: signal-to-noise ratio (SNR), absolute value of correlation coefficient (AVCC), and mean square error (MSE). They are, respectively, defined as follows.

The reconstruction index is defined as a signal-to-noise ratio (SNR) of the error [40]; that is,

The recovered source should be rescaled to the same energy level as their corresponding original sources. Furthermore, an absolute value of correlation coefficient (AVCC) [41] is used to evaluate the similar degree between recovered sources and original sources. Consider the following:

The mean square error (MSE) [21] is defined as

Higher SNR and AVCC and lower MSE represent better separating results, in other words, the separated signals will be more similar to the source signals.

5.3. Separation Results and Performance Analysis

All simulations are programmed using MATLAB version 2013 and executed on a personal computer equipped with an Intel core i5, 2.4 GHz processor, and 4 G RAM.

Firstly, we will evaluate the performance of the algorithm proposed using wave form. The source signals, mixed signals, and recovered signals using the proposed algorithm are shown in Figures 1(a), 1(b), and 1(c) and Figures 2(a), 2(b), and 2(c). Each signal has 8000 samples. From these figures, we can observe that the proposed algorithm performs separation efficiently.

(a) Sources , , and

(b) Mixtures , , and

(c) Recovered signals , , and

(d) Scatter plot

Figure 1: Separation results for speech signals.

(a) Sources , , and

(b) Mixtures , , and

(c) Recovered signals , , and

(d) Scatter plot

Figure 2: Separation results for music signals.

Figures 1(d) and 2(d) show the scatter plots of each recovered source signal versus three original source signals. The scatter plot of recovered signal versus the corresponding original source signal is almost a straight line, meaning that each source is recovered with very little interference from others. Then, we can also draw a conclusion that the proposed algorithm performs separation efficiently.

Secondly, we will illustrate the performance of the proposed algorithm using performance indexes: SNR, AVCC, and MSE. A number of 30 independent experiments have been carried out. The average experimental results of separating three speech signals and three music signals are summarized in Tables 1, 2, 3, and 4.

Table 1: SNR for different algorithms.

Signals
PsoBatch algorithm
PsoExtraction algorithm
AbcBatch algorithm
AbcExtraction algorithm
Proposed algorithm

Speech

31.4563
33.3178
31.7612
33.3171
31.5761

34.2591
30.3386
33.5527
30.3393
47.2315

33.9996
32.8233
33.6296
32.8136
35.7167

Music

18.5666
18.5669
18.2993
18.5663
25.8384

14.0873
20.8049
14.0664
20.8035
23.1312

24.7625
23.9757
24.3672
23.9760
50.9539

Table 2: AVCC for different algorithms.

Signals
PsoBatch algorithm
PsoExtraction algorithm
AbcBatch algorithm
AbcExtraction algorithm
Proposed algorithm

Speech

0.99947
0.99961
0.99729
0.99961
0.99959

0.99993
0.99973
0.99695
0.99973
0.99999

0.99980
0.99974
0.99470
0.99974
0.99986

Music

0.99223
0.99252
0.99128
0.99252
0.99857

0.99240
0.99591
0.99197
0.99591
0.99823

0.99895
0.99750
0.99828
0.99750
0.99999

Table 3: MSE for different algorithms.

Signals
PsoBatch algorithm
PsoExtraction algorithm
AbcBatch algorithm
AbcExtraction algorithm
Proposed algorithm

Speech

Music

Table 4: Average value of performance indexes.

Signals
Performance indexes
PsoBatch algorithm
PsoExtraction algorithm
AbcBatch algorithm
AbcExtraction algorithm
Proposed algorithm

Speech
SNR
33.2383
32.1599
32.9811
32.1567
38.1747

AVCC
0.99973
0.99969
0.99631
0.99969
0.99981

MSE

Music
SNR
19.1388
21.1158
18.9110
21.1153
33.3078

AVCC
0.99453
0.99531
0.99384
0.99531
0.99893

MSE

From Tables 1, 2, and 3, we can observe that the capability of the algorithm proposed is good. The SNR values between recovered signals and sources exceed 20 dB, the corresponding AVCC values exceed 0.998, and the MSE values are less than .

Meanwhile, the average value of different performance indexes for three way signals is given in Table 4. The results presented in Tables 1, 2, 3, and 4 show that the proposed algorithm provides a better quality of separation than other BIO-based BSS algorithms.

To illustrate the contribution of the improved artificial bee colony algorithm (MABC), the performance comparisons of MABC algorithm and other BIO algorithms (ABC, PSO) with the identical objective function (covariance ratio) have been given in Table 5 and Figures 3 and 4. The results show that MABC algorithm provides a better quality of separation and faster convergence than other BIO algorithms.

Table 5: Statistical results of different BIO algorithms.

Speech

Music

MABC
ABC
PSO

MABC
ABC
PSO

SNR

31.5761
31.2449
30.9816

25.8384
25.7068
24.2166

47.2315
45.9496
45.4138

23.1312
22.8939
21.4857

35.7167
35.1299
34.5042

50.9539
46.7180
43.6723

AVCC

0.99959
0.99956
0.99947

0.99857
0.99845
0.99689

0.99999
0.99997
0.99993

0.99823
0.99812
0.99778

0.99986
0.99984
0.99971

0.99999
0.99999
0.99793

MSE

(a) Speech

(b) Music

Figure 3: SNR for the three BIO algorithms.

(a) Speech

(b) Music

Figure 4: MSE for the three BIO algorithms.

The convergence curves for three different objective functions being optimized by MABC, respectively, are shown in Figures 5 and 6. We can draw a conclusion that the objective function based on covariance ratio results in better separation performance compared with other previous objective functions.

(a) Speech

(b) Music

Figure 5: SNR for different objective functions.

(a) Speech

(b) Music

Figure 6: MSE for different objective functions.

In order to further confirm the validation of the proposed algorithm, the performance comparisons between the proposed algorithm and some widely used algorithms like FastICA [42], extended Infomax [43], and JADE [44] are given in Table 6. The SNR, AVCC, and MSE for four algorithms are presented and the proposed algorithm performs with better results.

Table 6: Performance comparison between proposed algorithm and three widely used algorithms.

Speech

Music

FastICA
Extended Infomax
JADE
Proposed algorithm

FastICA
Extended Infomax
JADE
Proposed algorithm

SNR

30.3397
24.2091
31.3271
31.5761

16.4701
20.8582
19.3042
25.8384

32.8192
22.9595
32.4365
47.2315

11.6008
15.7295
17.4733
23.1312

32.1893
25.8227
32.2241
35.7167

16.0358
23.2738
37.3531
50.9539

AVCC

0.99953
0.99659
0.99955
0.99959

0.98688
0.99536
0.99345
0.99857

0.99974
0.99904
0.99981
0.99999

0.97552
0.99504
0.99446
0.99823

0.99961
0.99754
0.99971
0.99986

0.98633
0.99864
0.99990
0.99999

MSE

To illustrate the most distinguished advantage, we will compare the computation amount of the proposed algorithm with other BIO-based BSS algorithms. Therefore, we compare the computation time between different BIO-based BSS algorithms, under the same computer hardware platform and the same separation tasks, to evaluate the advantage of the proposed algorithm in computation amount.

The average computation time of 30 runs for separating three source signals is presented in Table 7. The proposed algorithm has outperformed the previous four BIO-based BSS algorithms in terms of the average computation time. From the above, we can draw a conclusion that the proposed algorithm can efficiently overcome the shortcoming of BIO-based BSS algorithms in computation amount. It has reduced the computation time for dozens of times and performed better separation than other algorithms.

Table 7: Average computation time for different algorithms.

Signals
PsoBatch algorithm
PsoExtraction algorithm
AbcBatch algorithm
AbcExtraction algorithm
Proposed algorithm

Speech
Max
17.56
10.89
33.05
8.82
0.86

Min
14.12
8.16
28.12
7.96
0.69

Ave
14.97
8.65
29.33
8.42
0.73

Music
Max
13.46
7.91
34.25
9.24
0.99

Min
12.74
7.39
27.44
8.00
0.71

Ave
12.91
7.51
30.65
8.55
0.78

6. Conclusion

In this paper, we have presented a new fast BIO-based BSS algorithm. In the proposed algorithm, the objective function for BSS based on covariance ratio is utilized and the modified ABC algorithm is used to solve the objective function. In the separation process, only two fixed covariance matrices are used repeatedly in evolution process of each generation, without kurtosis of all the sample points from mixtures being computed. The difficulty of computation amount for previous BIO-based BSS algorithms can be conquered. In order to prevent the same source signal being separated out repeatedly, deflation scheme is used to remove the source signal component separated out from mixtures after each separation process, and we can recover the source signals one by one. The results of the experiments show that the proposed algorithm outperforms the previous BIO-based BSS algorithm in terms of the computation amount and separation precision.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This study was supported by the National Science Foundation of China (no. 11127202).

References

Y. Zhang and Y. Zhao, “Modulation domain blind speech separation in noisy environments,” Speech Communication, vol. 55, no. 10, pp. 1081–1099, 2013.

M. T. Özgen, E. E. Kuruoǧlu, and D. Herranz, “Astrophysical image separation by blind time-frequency source separation methods,” Digital Signal Processing, vol. 19, no. 2, pp. 360–369, 2009.

A. Ikhlef, K. Abed-Meraim, and D. Le Guennec, “Blind signal separation and equalization with controlled delay for MIMO convolutive systems,” Signal Processing, vol. 90, no. 9, pp. 2655–2666, 2010.

R. Romo Vázquez, H. Vélez-Pérez, R. Ranta, V. Louis Dorr, D. Maquin, and L. Maillard, “Blind source separation, wavelet denoising and discriminant analysis for EEG artefacts and noise cancelling,” Biomedical Signal Processing and Control, vol. 7, no. 4, pp. 389–400, 2012.

M. Kuraya, A. Uchida, S. Yoshimori, and K. Umeno, “Blind source separation of chaotic laser signals by independent component analysis,” Optics Express, vol. 16, no. 2, pp. 725–730, 2008.

M. Babaie-Zadeh and C. Jutten, “A general approach for mutual information minimization and its application to blind source separation,” Signal Processing, vol. 85, no. 5, pp. 975–995, 2005.

K. Todros and J. Tabrikian, “Blind separation of independent sources using Gaussian mixture model,” IEEE Transactions on Signal Processing, vol. 55, no. 7, pp. 3645–3658, 2007.

P. Comon, “Tensors: a brief introduction,” IEEE Signal Processing Magazine, vol. 31, no. 2, pp. 44–53, 2014.

E. Oja and M. Plumbley, “Blind separation of positive sources using nonnegative PCA,” in Proceedings of the 4th International Symposium on Independent Component Analysis and Blind Signal Separation (ICA '03), pp. 11–16, Nara, Japan, April 2003.

W. L. Woo and S. S. Dlay, “Neural network approach to blind signal separation of mono-nonlinearly mixed sources,” IEEE Transactions on Circuits and Systems I, vol. 52, no. 6, pp. 1236–1247, 2005.

A. Cichocki and S. Amari, Adaptive Blind Signal and Image Processing, John Wiley & Sons, 2002.

A. Hyvarinen, J. Karhunen, and E. Oja, Independent Component Analsysis, John Wiley & Sons, 2001.

P. Comon and C. Jutten, Handbook of Blind Source Separation: Independent Component Analysis and Applications, Academic Press, 2010.

A. Cichocki and R. Unbehauen, “Robust neural networks with on-line learning for blind identification and blind separation of sources,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 43, no. 11, pp. 894–906, 1996.

S.-I. Amari, T.-P. Chen, and A. Cichocki, “Stability analysis of learning algorithms for blind source separation,” Neural Networks, vol. 10, no. 8, pp. 1345–1351, 1997.

S.-I. Amari, “Natural gradient works efficiently in learning,” Neural Computation, vol. 10, no. 2, pp. 251–276, 1998.

J. Kennedy, R. Eberhart, and Y. H. Shi, Swarm Intelligence, Elsevier, 2001.

X. S. Yang, Z. H. Cui, and R. B. Xiao, Swarm Intelligence and Bio-Inspired Computation Theory and Applications, Elsevier, 2013.

M. Li and J. Mao, “A new algorithm of evolutional blind source separation based on genetic algorithm,” in Proceedings of the 5th World Congress on Intelligent Control and Automation, pp. 2240–2244, Hangzhou, Zhejiang, China, June 2004.

S. Mavaddaty and A. Ebrahimzadeh, “Evaluation of performance of genetic algorithm for speech signals separation,” in Proceedings of the International Conference on Advances in Computing, Control and Telecommunication Technologies (ACT '09), pp. 681–683, Trivandrum, Kerala, India, December 2009.

J. M. Górriz, C. G. Puntonet, F. Rojas, R. Martin, S. Hornillo, and E. W. Lang, “Optimizing blind source separation with guided genetic algorithms,” Neurocomputing, vol. 69, no. 13–15, pp. 1442–1457, 2006.

S.-T. Hsieh, T.-Y. Sun, C.-L. Lin, and C.-C. Liu, “Effective learning rate adjustment of blind source separation based on an improved particle swarm optimizer,” IEEE Transactions on Evolutionary Computation, vol. 12, no. 2, pp. 242–251, 2008.

C.-Z. Zhang, J.-P. Zhang, and X.-D. Sun, “Blind source separation based on adaptive particle swarm optimization,” Systems Engineering and Electronics, vol. 31, no. 6, pp. 1275–1278, 2009.

T.-Y. Sun, C.-C. Liu, S.-J. Tsai, S.-T. Hsieh, and K.-Y. Li, “Cluster guide particle swarm optimization (CGPSO) for underdetermined blind source separation with advanced conditions,” IEEE Transactions on Evolutionary Computation, vol. 15, no. 6, pp. 798–811, 2011.

L. Chen, L. Zhang, Y. Guo, T. Liu, and Q. Li, “Sequential blind image extraction based on particle swarm optimization,” in Proceedings of the 6th International Conference on Natural Computation (ICNC '10), pp. 2697–2700, Yantai, Shandong, China, August 2010.

S. Mavaddaty and A. Ebrahimzadeh, “A comparative study of bees colony algorithm for blind source separation,” in Proceedings of the 20th Iranian Conference on Electrical Engineering (ICEE '12), pp. 1172–1177, Tehran, Iran, May 2012.

Y.-X. Zhang, X.-M. Tian, and X.-G. Deng, “Blind source separation based on modified artificial bee colony algorithm,” Acta Electronica Sinica, vol. 40, no. 10, pp. 2026–2030, 2012.

J. V. Stone, “Blind source separation using temporal predictability,” Neural Computation, vol. 13, no. 7, pp. 1559–1574, 2001.

S. Xie, Z. He, and Y. Fu, “A note on Stone's conjecture of blind signal separation,” Neural Computation, vol. 17, no. 2, pp. 321–330, 2005.

H.-L. Liu and Y.-M. Cheung, “On blind source separation using generalized eigenvalues with a new metric,” Neurocomputing, vol. 71, no. 4–6, pp. 973–982, 2008.

D. Karaboga and B. Basturk, “A powerful and efficient algorithm for numerical function optimization: artificial bee colony (ABC) algorithm,” Journal of Global Optimization, vol. 39, no. 3, pp. 459–471, 2007.

D. Karaboga and B. Akay, “A comparative study of artificial bee colony algorithm,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 108–132, 2009.

F. Kang, J. Li, and H. Li, “Artificial bee colony algorithm and pattern search hybridized for global optimization,” Applied Soft Computing Journal, vol. 13, no. 4, pp. 1781–1791, 2013.

F. Kang, J. Li, and Z. Ma, “An artificial bee colony algorithm for locating the critical slip surface in slope stability analysis,” Engineering Optimization, vol. 45, no. 2, pp. 207–223, 2013.

B. Akay and D. Karaboga, “A modified artificial bee colony algorithm for real-parameter optimization,” Information Sciences, vol. 192, pp. 120–142, 2012.

D. Karaboga and B. Akay, “A modified Artificial Bee Colony (ABC) algorithm for constrained optimization problems,” Applied Soft Computing Journal, vol. 11, no. 3, pp. 3021–3031, 2011.

W.-F. Gao and S.-Y. Liu, “A modified artificial bee colony algorithm,” Computers and Operations Research, vol. 39, no. 3, pp. 687–697, 2012.

R. Storn and K. Price, “Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces,” Journal of Global Optimization, vol. 11, no. 4, pp. 341–359, 1997.

Y.-M. Cheung and H. Liu, “A new approach to blind source separation with global optimal property,” in Proceedings of the IASTED International Conference on Neural Networks and Computational Intelligence, pp. 137–141, Grindelwald, Switzerland, February 2004.

P. Bofill and M. Zibulevsky, “Underdetermined blind source separation using sparse representations,” Signal Processing, vol. 81, no. 11, pp. 2353–2362, 2001.

D. Freedman, R. Pisani, and R. Purves, Statistics, W. W. Norton & Company, 2007.

A. Hyvärinen, “Fast and robust fixed-point algorithms for independent component analysis,” IEEE Transactions on Neural Networks, vol. 10, no. 3, pp. 626–634, 1999.

T.-W. Lee, M. Girolami, and T. J. Sejnowski, “Independent component analysis using an extended infomax algorithm for mixed subgaussian and supergaussian sources,” Neural Computation, vol. 11, no. 2, pp. 417–441, 1999.

J. F. Cardoso and A. Souloumiac, “Blind beamforming for non-Gaussian signals,” IEE Proceedings F: Radar and Signal Processing, vol. 140, no. 6, pp. 362–370, 1993.