International Journal of Performability Engineering, 2019, 15(1): 337-345 doi: 10.23940/ijpe.19.01.p34.337345

Mixing Matrix Estimation Algorithm for Underdetermined Instantaneous Mixing Model

Fan Shi,a and Chuang Liub

a Faculty of Maritime and Transportation, Ningbo University, Ningbo, 315211, China

b College of Information and Communication Engineering, Harbin Engineering University, Harbin, 150001, China

Corresponding authors: * E-mail address:shifan@nbu.edu.cn

Accepted: 2019-01-3   Online: 2019-01-1

Abstract

Estimating the mixing matrix is a research that focuses on underdetermined blind source separation. In order to get a more accurate estimated mixing matrix, we investigate a novel algorithm for mixing matrix estimation. Firstly, a new method for detecting single source points was introduced. Then, we reckoned the signal quantity and initial clustering centers by adopting an improved clustering method based on the potential function. Finally, the point density theory and the initial clustering centers were utilized to get more accurate clustering centers and estimated the mixing matrix. The simulation results illustrate that we can obtain the more accurate and stable estimation of the mixing matrix by using the proposed algorithm.

Keywords: underdetermined blind source separation; mixing matrix estimation; potential function; single source points

PDF (497KB) Metadata Related articles Export EndNote| Ris| Bibtex  Favorite

Cite this article

Fan Shi and Chuang Liu. Mixing Matrix Estimation Algorithm for Underdetermined Instantaneous Mixing Model. International Journal of Performability Engineering, 2019, 15(1): 337-345 doi:10.23940/ijpe.19.01.p34.337345

1. Introduction

As the mixing channel and the source signals are not established, how to get each source signal from mixed signals has been a worthy problem in science and engineering. Therefore, blind source separation (BSS) is put forward, and more and more researchers have been involved in research corresponding BSS algorithms. The BSS algorithm is now widely used in areas such as image processing [1-2], ship noise processing [3-4], bio-medical engineering [5-6] and wireless system [7-8]. In recent years, new methods are proposed unceasingly, and the theory of BSS is modified and improved. So far, BSS has been classed into standard blind source separation (SBSS), overdetermined blind source separation (OBSS) and underdetermined blind source separation (UBSS). For its extensive applicability, UBSS is now a hot topic in this field. Due to the limit of practical conditions, the number of sensors may be less than the sources. In this situation, how to solve the problem is relevant to the UBSS problem, which is even trickier than the other two problems. One of the most important methods to solve the SBSS problem and the OBSS problem is independent component analysis (ICA) [9-10]. Because the system of the UBSS model is irreversible, ICA is not suitable for solving this problem. Generally, the sparse component analysis (SCA)- based two-step method [11-12] is a fairly important method for UBSS problems. The first step of the method is to estimate a mixing matrix, and the next step is source signal recovery. Since an accurate estimation of the mixing matrix is the basis of signals recovery, the accuracy of the first step is very important. Therefore, the purpose of this study was to determine the estimation of the mixing matrix, the sparsity of signals and how to use the sparsity to determine the performance of the SCA methods. As a matter of fact, some signals are sparse in the time domain, but most signals are not. To obtain more sparse features of signals, the method of time-frequency (t-f) transformation is usually adopted because signals can present better sparsity in the t-f domain. Short-time Fourier transform (STFT) is a main t-f transform method for its implementation simplicity. Better sparsity leads to better linear clustering property. Combining this property and clustering algorithms, the mixing matrix is then generated with received signals. However, there will be some outliers if the signal in the t-f plane is not sufficiently sparse or noises exists. To solve the problem of outliers, the methods [13-17] based on single source point detection and single source region detection get more and more attention. Single source regions denote the t-f regions that only a source signal exists or only a source signal is dominant. Similarly, a single source point is denoted as the t-f point that only one source signal exists or only one source signal is dominant. These single source points make the observed signals present better linear clustering property. The better the performance of single source points detects, the higher the estimation accuracy is. Many scholars have researched the corresponding algorithms and have done a lot of work on single source region detection and single source point detection. From Literature [18], some t-f regions where only a source signal exists are first chosen by utilizing the t-f ratios, then are used to estimate the mixing matrix. Literature [19] extended the algorithm in Literature [18] to solve the problem of the time-delay mixing model. Lu [20] introduced an algorithm based on single source detection to estimate mixing matrix. First, the t-f region of observed signals is divided into many disjoint regions. Then, detect the regions where only a source signal occurs by using the method of eigenvalue decomposition. Finally, obtain the mixing matrix estimation with the hierarchical clustering. All of the above algorithms utilize single source regions. Based on the single source region algorithms, for the underdetermined instantaneous mixing model, scholars propose many single source point algorithms to get the estimation of the mixing matrix. Reju [21] investigated a method to detect the single source points by judging whether the direction of the real and imaginary parts of STFT coefficients of observed signals are the same, and gets the estimated mixing matrix utilizing a hierarchical clustering algorithm. Dong [22] studied an algorithm for estimating the mixing matrix. First of all, the real vector and the imaginary vector of t-f coefficients of observed signals are normalized. Then, the single source points are detected by measuring the distance of the normalized real vector and the normalized imaginary vector. Finally, the robustness clustering algorithm based on similarity is studied to get an estimated mixing matrix. Besides, considering single source points, many algorithms are used to estimate the mixing matrix in the underdetermined time-delay mixing model. The major contributions of this is a mixing matrix method with better performance for UBSS problems is investigated. Through proposing a new single source points detection algorithm and an improved clustering algorithm, the more accurate mixing matrices can be obtained.

In the remainder of this paper. Section 2 introduces the instantaneous mixing model about UBSS problems, and the theory of single source points detection is analyzed. Section 3 presents an improved clustering algorithm. In Section 4, we carry out several simulations, and the results and correlation analysis are shown. Section 5 gives out the conclusion.

2. Problem Formulation

The underdetermined instantaneous mixing model without noise interference canbe expressed as

$x(t)=As(t)$

where $x(t)={{[{{x}_{1}}(t),{{x}_{2}}(t),\cdots ,{{x}_{N}}(t)]}^{T}}$ denotes the mixed signal vector at the time point t. The mixing matrix is represented by $A=[{{a}_{1}},{{a}_{2}},\cdots ,{{a}_{M}}]\in {{R}^{N\times M}},\text{ }(M>N)$, in which,${{a}_{i}}$ denotes the mixing matrix’s column vector. $s(t)={{[{{s}_{1}}(t),{{s}_{2}}(t),\cdots ,{{s}_{M}}(t)]}^{T}}$ denotes the source signal vector at the time pointt. The Equation (1) can also be denoted as

$x(t)=\sum\limits_{i=1}^{M}{{{a}_{i}}{{s}_{i}}(t)}$

In order to get better sparsity, STFT is applied in this paper. The STFT of thenth mixed signal is expressed as

${{X}_{n}}(t, f)=\int_{-\infty }^{\infty }{{{x}_{n}}(\tau )h(t-\tau )}{{e}^{-j2\pi f\tau }}\text{d}\tau $

where $h(t)$ is the window function. The STFT of themth source signal can be written as

${{S}_{m}}(t, f)=\int_{-\infty }^{\infty }{{{s}_{m}}(\tau )}h(t-\tau ){{e}^{-j2\pi f\tau }}\text{d}\tau $

The t-f representation of the mixed signals is

$X(t, f)=AS(t, f)=\sum\limits_{i=1}^{M}{{{a}_{i}}{{S}_{i}}(t, f)}$

$X(t, f)={{[{{X}_{1}}(t, f),{{X}_{2}}(t, f),\cdots ,{{X}_{N}}(t, f)]}^{T}}$ and $S(t, f)={{[{{S}_{1}}(t, f),{{S}_{2}}(t, f),\cdots ,{{S}_{M}}(t, f)]}^{T}}$ are the mixed signal after STFT and the source signal after STFT, respectively. In this paper, two mixed signals are taken as an example to analyze the algorithms of this paper. The Equation (5) can also be written as

$\left[ \begin{matrix} {{X}_{1}}(t, f) \\ {{X}_{2}}(t, f) \\ \end{matrix} \right]=\left[ \begin{matrix} {{a}_{11}} & {{a}_{12}} & \cdots & {{a}_{1M}} \\ {{a}_{21}} & {{a}_{22}} & \cdots & {{a}_{2M}} \\ \end{matrix} \right]\left[ \begin{matrix} {{S}_{1}}(t, f) \\ {{S}_{2}}(t, f) \\ \vdots \\ {{S}_{M}}(t, f) \\ \end{matrix} \right]$

The mixing matrix estimation accuracy is mainly relatedto the sparsity of signals. For the single source points detection algorithm, any two columns in the mixing matrix should not be correlated and some single source points must exist. If only one source signal ${{s}_{i}}$ presets at the time t, Formula (6) is rewritten as

$\left[ \begin{matrix} {{X}_{1}}(t, f) \\ {{X}_{2}}(t, f) \\ \end{matrix} \right]=\left[ \begin{matrix} {{a}_{1i}} \\ {{a}_{2i}} \\ \end{matrix} \right]{{S}_{i}}(t, f)$

For the single source point $({{t}_{p}},{{f}_{p}})$ only contain the source signal ${{s}_{1}}$, we can obtain

${{X}_{1}}({{t}_{p}},{{f}_{p}})={{a}_{11}}{{S}_{1}}({{t}_{p}},{{f}_{p}})={{a}_{11}}\left\{ \operatorname{Re}[{{S}_{1}}({{t}_{p}},{{f}_{p}})]+j\operatorname{Im}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}$
${{X}_{2}}({{t}_{p}},{{f}_{p}})={{a}_{21}}{{S}_{1}}({{t}_{p}},{{f}_{p}})={{a}_{21}}\left\{ \operatorname{Re}[{{S}_{1}}({{t}_{p}},{{f}_{p}})]+j\operatorname{Im}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}$

Based on the Equation (8) and the Equation (9), the following formula can be obtained

${{\left[ \frac{{{X}_{1}}({{t}_{p}},{{f}_{p}})}{{{X}_{2}}({{t}_{p}},{{f}_{p}})} \right]}^{2}}=\frac{a_{11}^{2}{{\left\{ \operatorname{Re}[{{S}_{1}}({{t}_{p}},{{f}_{p}})]+j\operatorname{Im}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}}^{2}}}{a_{21}^{2}{{\left\{ \operatorname{Re}[{{S}_{1}}({{t}_{p}},{{f}_{p}})]+j\operatorname{Im}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}}^{2}}}=\frac{a_{11}^{2}}{a_{21}^{2}}$
${{\left[ \frac{\left| {{X}_{1}}({{t}_{p}},{{f}_{p}}) \right|}{\left| {{X}_{2}}({{t}_{p}},{{f}_{p}}) \right|} \right]}^{2}}={{\left( \frac{\sqrt{a_{11}^{2}\left( {{\left\{ \operatorname{Re}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}}^{2}}+{{\left\{ \operatorname{Im}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}}^{2}} \right)}}{\sqrt{a_{21}^{2}\left( {{\left\{ \operatorname{Re}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}}^{2}}+{{\left\{ \operatorname{Im}[{{S}_{1}}({{t}_{p}},{{f}_{p}})] \right\}}^{2}} \right)}} \right)}^{2}}=\frac{a_{11}^{2}}{a_{21}^{2}}$
${{\left[ \frac{{{X}_{1}}({{t}_{p}},{{f}_{p}})}{{{X}_{2}}({{t}_{p}},{{f}_{p}})} \right]}^{2}}={{\left[ \frac{\left| {{X}_{1}}({{t}_{p}},{{f}_{p}}) \right|}{\left| {{X}_{2}}({{t}_{p}},{{f}_{p}}) \right|} \right]}^{2}}$

If two source signals ${{s}_{1}}$ and ${{s}_{2}}$ exist at the t-f point $({{t}_{q}},{{f}_{q}})$, we can denote the equations as

${{X}_{1}}({{t}_{q}},{{f}_{q}})=\left\{ {{a}_{11}}\operatorname{Re}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{12}}\operatorname{Re}[{{S}_{2}}({{t}_{q}},{{f}_{q}}] \right\}+j\left\{ {{a}_{11}}\operatorname{Im}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{12}}\operatorname{Im}[{{S}_{2}}({{t}_{q}},{{f}_{q}}] \right\}$
${{X}_{2}}({{t}_{q}},{{f}_{q}})=\left\{ {{a}_{21}}\operatorname{Re}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{22}}\operatorname{Re}[{{S}_{2}}({{t}_{q}},{{f}_{q}}] \right\}+j\left\{ {{a}_{21}}\operatorname{Im}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{22}}\operatorname{Im}[{{S}_{2}}({{t}_{q}},{{f}_{q}}] \right\}$

The following equations are set

${{R}_{1}}({{t}_{q}},{{f}_{q}})={{a}_{11}}\operatorname{Re}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{12}}\operatorname{Re}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]$(15)
${{I}_{1}}({{t}_{q}},{{f}_{q}})={{a}_{11}}\operatorname{Im}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{12}}\operatorname{Im}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]$
${{R}_{2}}({{t}_{q}},{{f}_{q}})={{a}_{21}}\operatorname{Re}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{22}}\operatorname{Re}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]$
${{I}_{2}}({{t}_{q}},{{f}_{q}})={{a}_{21}}\operatorname{Im}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{22}}\operatorname{Im}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]$

The Equation (13) and the Equation (14) can also be denoted as

${{X}_{1}}({{t}_{q}},{{f}_{q}})={{R}_{1}}({{t}_{q}},{{f}_{q}})+j{{I}_{1}}({{t}_{q}},{{f}_{q}})$
${{X}_{2}}({{t}_{q}},{{f}_{q}})={{R}_{2}}({{t}_{q}},{{f}_{q}})+j{{I}_{2}}({{t}_{q}},{{f}_{q}})$

Based on the Equation (19) and the Equation (20), the following equations can be obtained

${{\left[ \frac{{{X}_{1}}({{t}_{q}},{{f}_{q}})}{{{X}_{2}}({{t}_{q}},{{f}_{q}})} \right]}^{2}}=\frac{{{\left\{ [{{R}_{1}}({{t}_{q}},{{f}_{q}})]+j[{{I}_{1}}({{t}_{q}},{{f}_{q}})] \right\}}^{2}}}{{{\left\{ [{{R}_{2}}({{t}_{q}},{{f}_{q}})]+j[{{I}_{2}}({{t}_{q}},{{f}_{q}})] \right\}}^{2}}}$
${{\left[ \frac{\left| {{X}_{1}}({{t}_{q}},{{f}_{q}}) \right|}{\left| {{X}_{2}}({{t}_{q}},{{f}_{q}}) \right|} \right]}^{2}}=\frac{{{[{{R}_{1}}({{t}_{q}},{{f}_{q}})]}^{2}}+{{[{{I}_{1}}({{t}_{q}},{{f}_{q}})]}^{2}}}{{{[{{R}_{2}}({{t}_{q}},{{f}_{q}})]}^{2}}+{{[{{I}_{2}}({{t}_{q}},{{f}_{q}})]}^{2}}}$

In order to make the Equation (21) to be equal to the Equation (22), the following equation needs to be satisfied

$\frac{{{R}_{1}}({{t}_{q}},{{f}_{q}})}{{{R}_{2}}({{t}_{q}},{{f}_{q}})}=\frac{{{I}_{1}}({{t}_{q}},{{f}_{q}})}{{{I}_{2}}({{t}_{q}},{{f}_{q}})}$

The Equation (23) can also be denoted as

$\frac{{{a}_{11}}\operatorname{Re}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{12}}\operatorname{Re}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]}{{{a}_{21}}\operatorname{Re}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{22}}\operatorname{Re}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]}=\frac{{{a}_{11}}\operatorname{Im}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{12}}\operatorname{Im}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]}{{{a}_{21}}\operatorname{Im}[{{S}_{1}}({{t}_{q}},{{f}_{q}})]+{{a}_{22}}\operatorname{Im}[{{S}_{2}}({{t}_{q}},{{f}_{q}})]}$

The above equation is a condition. Based on this condition, the simpler condition can be denoted as

${{a}_{11}}{{a}_{22}}={{a}_{12}}{{a}_{21}}$

This condition doesn’t satisfy the hypothesis that any two columns are not correlated in the mixing matrix, so the Formula (21) can’t be equal to the Formula (22). Similarly, if more source signals exist at some t-f points, the result is similar. For detecting the single source points, we will take advantage of the following formula.

${{\left[ \frac{{{X}_{1}}(t, f)}{{{X}_{2}}(t, f)} \right]}^{2}}-{{\left[ \frac{\left| {{X}_{1}}(t, f) \right|}{\left| {{X}_{2}}(t, f) \right|} \right]}^{2}}=0$

In practice, the above condition is so strict that it is difficult to achieve it. To detect the single source points better, we relax Equation (26) as

${{\left[ \frac{{{X}_{1}}(t, f)}{{{X}_{2}}(t, f)} \right]}^{2}}-{{\left[ \frac{\left| {{X}_{1}}(t, f) \right|}{\left| {{X}_{2}}(t, f) \right|} \right]}^{2}}<{{\varepsilon }_{1}}$

where ${{\varepsilon }_{1}}$ is a positive value that is close to 0, through the above formula, we can detect the single source points.

3. Mixing Matrix Estimation based on the Improved Potential Function

To get the estimated mixing matrix, clustering algorithms are usually employed. If observed mixed signals are sparse enough, observed mixed signals have the linear clustering property which is determined by the mixing matrix column vector. When the observed signals are incomplete sparsity, the linear clustering property is not obvious. T-f transform is used to get more obvious linear clustering property. Get the clustering centers by using the clustering of the single source points. Then, an estimated mixing matrix can be achieved. In our work, we take advantage of the potential function to obtain the quantity of source signals, then the clustering centers are estimated by the point density theory. First, get the quantity of source signals and achieve a rough estimation of the clustering centers by utilizing the clustering method based on the potential function. The clustering centers are some initial values. More accurate clustering centers need subsequent processing. The improved potential function that makes a rough estimation of the clustering centers can be denoted as

$J(z)={{\sum\limits_{j=1}^{Q}{\left[ \exp \left( \frac{\sqrt{{{x}_{j}}{{z}^{T}}}}{{{\sigma }^{2}}} \right) \right]}}^{\gamma }}$

in which, ${{x}_{j}}$ represents the M-dimensional vector and it belongs to a data set $D=\{{{x}_{1}},{{x}_{2}},\cdots ,{{x}_{Q}}\}$ that includes the real part of each single source point. Qdenotes the single source point number. $\gamma$ is a parameter that has an impact on estimating the number of peaks and can be estimated with a correlation method. The formula of ${{\sigma }^{2}}$ is as follows

${{\sigma }^{2}}=\sum\nolimits_{j=1}^{Q}{\left\| {{x}_{j}}-\bar{x} \right\|}{}^{2}/Q$

where $\bar{x}=\sum\nolimits_{j=1}^{Q}{{{x}_{j}}/Q}$. We can change the value of z to get the local maximum of $J(z)$, which can be obtained when ${{x}_{j}}$ is equal to z. For eliminating some false peaks produced by noises and interferences, the following equation is adopted

$\frac{J({{z}_{i}})}{\max J({{z}_{i}})}\ge {{\varepsilon }_{2}},\text{ }i=1, 2,\cdots , Q$

where ${{\varepsilon }_{2}}$ denotes a threshold value and ${{z}_{i}}$ represents the clustering center before optimization. This formula is used to remove those small maxima. The number of peaks can be estimated through the peak number produced by the potential function. Meanwhile, the rough estimation of the clustering centers can be obtained through the above-mentioned potential function method.

Then, the distance from one point to another is used to subdivide the single source points. This approach divides these points by introducing the idea of density. If a point is surrounded by many other points, this point is considered to own high density. In this paper, we keep the points with high density and remove the points with low density in order to correct the clustering centers. After getting rid of outliers, the initial clustering centers can be correctedand the new clustering centers will be obtained. Finally, we will get the estimated mixing matrix.

The improved algorithm main steps are as follows:

Step 1 Calculate the t-f coefficients of observed signals through STFT.

Step 2 The new detection strategy is used to detect single source points.

Step 3 Obtain the initial clustering centers by using the potential function-based method.

Step 4 Estimate more accurate clustering centers and finally obtain the mixing matrix.

4. Simulation Results and Analysis

To evaluate the performance of our introduced method, four source signals and two observation signals are used for simulation experiments. Four source signals are speech signals described in [21]. Table 1 shows the parameters of specific simulations.

Table 1.   Parameters of simulations

Sampling numberSTFT sizeOverlappingWeighting function${{\varepsilon }_{1}}$${{\varepsilon }_{2}}$
1600001024896Hanning window0.00010.01

New window| CSV


The mixing matrix Ais set as

$A=\left[ \begin{matrix} 0.9920 & 0.8499 & 0.5317 & 0.1673 \\ 0.1247 & 0.5177 & 0.8472 & 0.9850 \\ \end{matrix} \right]$

STFT is applied to process the observed signals in order to obtain better estimation accuracy. Then, the single source points could be approximated through the single source detection. When there is no noise, the scatter plot of two mixed signals is considered. Because of a large number of t-f points will lead to high computational complexity, it is necessary to reduce some points closing to the origin. Figure 1 displays the scatter plot of two observed signals in t-f plane before single source points detection.

Figure 1.

Figure 1.   The scatter plot of two observed signals in t-f plane before detecting single source points


Figure 2.

Figure 2.   The scatter plot of two observed signals in t-f plane after single source point detection and removing the points near the origin


As shown in Figure 2, some outliers have been removed by single source detection.

The scatter plot of two observed signals in t-f plane with the sign changes and normalization can be shown in Figure 3.

Figure 3.

Figure 3.   The scatter plot of two observed signals in t-f domain with the sign changes and normalization


In Figure 3, these points have the clustering property but there are still some interferences. The potential function is introduced to process these points. Figure 4 shows the result processed by the potential function.

Figure 4.

Figure 4.   Three-dimensional plot of $J(z)$


From Figure 4, the number of peaks corresponds to the quantity of the source signals. Meanwhile, the positions of peaks are the estimated values of the mixing matrix. The initial clustering centers can help us get the mixing matrix, but the result is not satisfactory. A more accurate clustering result is needed, and the point density theory is adopted. Based on the density distribution of the points, some points are removed in order to get better performance. The two observed signals scatter plot in t-f plane with the sign changes and normalization after removing some points with low density are shown in Figure 5.

Figure 5.

Figure 5.   The scatter plot of two observed signals in t-f plane with the sign changes and normalization after removing some points with low density


As shown in Figure 5, a better clustering property is present after removing the points with low density. The corrected clustering centers are estimated based on the point density theory. Andfinally we can obtain the estimated mixing matrix

$\tilde{A}=\left[ \begin{matrix} 0.9915 & 0.8509 & 0.5320 & 0.1700 \\ 0.1296 & 0.5206 & 0.8468 & 0.9850 \\ \end{matrix} \right]$

Comparing the estimated mixing matrix with the original one, estimated values are close to the real values. It is illustrated that our studied algorithm can estimate the mixing matrix well. For the sake of illustrating the universality of the proposed algorithm, we carried out another simulation to test the proposed algorithm. Other conditions remain changeless. After the experiment, the original mixing matrix and the estimated one are as shown

$A=\left[ \begin{matrix} 0.5162 & 0.2015 & 0.7618 & 0.9905 \\ 0.8569 & 0.9763 & 0.6477 & 0.1377 \\ \end{matrix} \right]$

$\tilde{A}=\left[ \begin{matrix} 0.5133 & 0.2047 & 0.7628 & 0.9904 \\ 0.8584 & 0.9772 & 0.6466 & 0.1385 \\ \end{matrix} \right]$

The above experiment proves that the above algorithm is able to estimate other mixing matrices effectively.

For determining the effectiveness of our studied algorithm, we make use of the normalized mean square error (NMSE) to evaluate the estimated mixing matrix, whose definition is

$NMSE=10\log \left( \frac{\sum\nolimits_{i, j}{{{({{{\tilde{a}}}_{ij}}-{{a}_{ij}})}^{2}}}}{\sum\nolimits_{i, j}{a_{i, j}^{2}}} \right)$

where ${{\tilde{a}}_{ij}}$ and ${{a}_{ij}}$ represent the (i, j)th element of the estimated mixing matrix and the (i, j)th element of the real one, respectively. Our proposed algorithm is compared with the other algorithms in Literature [21] or [22]. In a noisy case, the relative robustness of the algorithm is evaluated by adding Gaussian white noises. The curves of average NMSE for each algorithm obtained by 100 Monte Carlo trials are displayed in Figure 6.

Figure 6.

Figure 6.   The comparison of the proposed algorithm and other algorithms


As shown in Figure 6, our proposed algorithm is more robust when noises are present. The proposed algorithm is able to estimate the mixing matrix more accurately compared with other comparison algorithms.

5. Conclusions

An improved mixing matrix estimation algorithm was built in this paper as the number of the source signals is unknown, this method can estimate the mixing matrix more accurately and stably. A new single source point detection rule is proposed in t-f plane to get better sparsity of the observed signals. The clustering algorithm based on potential function was used to obtain the estimated value of the source signal number and initial clustering centers. (rewrite). This algorithm solves the problems that the clustering results are unstable, and the accuracy is low.

Acknowledgements

The work is funded by the National key research and development program (Grant No. 2016YFF0102806).

Reference

J.J. Yang and H. L. Liu ,

“Blind Identification of the Underdetermined Mixing Matrix based on K-Weighted HyperlineClustering, ”

Neurocomputing, Vol. 149, No. PB, pp. 483-489, 2015

[Cited within: 1]

M. T. Huang, C. H. Lee, C. M. Lin ,

“Blind Source Separation with Adaptive Learning Rates for Image Encryption, ”

Journal of Intelligent and Fuzzy Systems, Vol. 30, No. 1, pp. 451-460, 2015

[Cited within: 1]

J. Liu, S. E. Yang, S. C. Piao, Y. W. Huang ,

“Blind Source Separation of Ship-Radiated Noise using Single Observing Channel, ”

ACTA ACUSTICA, Vol. 36, No. 3, pp. 265-270, 2011

[Cited within: 1]

Z. C. Shan, C. S. Lin, Q. Xiang ,

“The Separation of Ship Signal based on Second OrderNonstationary Statistic, ”

SIGNALPROCESSING, Vol. 25, No. 6, pp. 973-976, 2009

[Cited within: 1]

F. Andreotti, J. Behar, S. Zaunseder, J. Oster, G. D. Clifford ,

“An Open-Source Framework for Stress-Testing Non-Invasive Foetal ECG Extraction Algorithms, ”

Physiological Measurement, Vol. 37, No. 5, pp. 627-648, 2016

[Cited within: 1]

S. Farashi ,

“Spike Detection using a Multiresolution Entropy basedMethod, ”

Biomedical Engineering-Biomedizinische Technik, Vol. 63, No. 4, pp. 361-376, 2018

[Cited within: 1]

M. G.S .Sriyananda, J. Joutsensalo, and T. Hamalainen ,

“Blind Source Separation based Interference Suppression Schemes for OFDM and DS-CDMA, ”

Telecommunication Systems, Vol. 4, No. 2, pp. 1-10, 2016

[Cited within: 1]

Z. C. Sha, Z. T. Huang, Y. Y. Zhou, F. H. Wang ,

“Frequency-Hopping Signals Sorting based on Underdetermined Blind Source Separation, ”

IET Communications, Vol. 7, No. 14, pp. 1456-1464, 2013

[Cited within: 1]

G.R. Naik and D. K. Kumar ,

“An Overview of Independent Component Analysis and its Applications, ”

International Journal of Computing and Information Sciences,Vol. 35, No. 1, pp. 63-81, 2011

[Cited within: 1]

X. He, F. He, T. Zhu ,

“Large-Scale Super-Gaussian Sources Separation using Fast-ICA with Rational Nonlinearities, ”

International Journal of Adaptive Control and Signal Processing, Vol. 31, No. 3, pp. 379-397, 2017

[Cited within: 1]

S. Qin, J. Guo, C. Zhu ,

“Sparse Component Analysis using Time-Frequency Representations for Operational Modal Analysis, ”

Sensors, Vol. 15, No. 3, pp. 6497-6519, 2015

[Cited within: 1]

Y. Zhong, X. Wang, L. Zhao, R. Y. Feng, L. P. Zhang, Y. Y. Xu ,

“Blind Spectral Unmixing based on Sparse Component Analysis for Hyperspectral Remote Sensing Imagery, ”

Isprs Journal of Photogrammetry and Remote Sensing, Vol. 119, pp. 49-63, 2016

[Cited within: 1]

Y. Li, S. I. Amari, A. Cichocki, D. W.C. Ho, and S. L. Xie ,

“Underdetermined Blind Source Separation based on Sparse Representation, ”

IEEE Transactions on Signal Processing,Vol. 54, No. 2, pp. 423-437, 2006

[Cited within: 1]

Yilmaz and S. Rickard ,

“Blind Separation of Speech Mixtures Via Time-Frequency Masking, ”

IEEE Transactions on Signal Processing, Vol. 52, No. 7, pp. 1830-1847, 2004

[Cited within: 1]

A. Aissa-El-Bey, N. Linh-Trung, K. Abed-Meraim, A. Belouchrani, and Y. Grenier ,

“Underdetermined Blind Separation of Nondisjoint Sources in the Time-Frequency Domain, ”

IEEE Transactions on Signal Processing, Vol. 55, No. 3, pp. 897-907, 2007

J. Thiagarajan, K. Natesan Ramamurthy, A. Spanias ,

“Mixing Matrix Estimation using Discriminative Clustering for Blind Source Separation, ”

Digital Signal Processing, Vol. 23, No. 1, pp. 9-18, 2013

F.C. Feng and M. Kowalski ,

“Revisiting Sparse ICA from a Synthesis Point of View: Blind Source Separation for over and Underdetermined Mixtures, ”

Signal Processing, Vol. 152, pp. 165-177, 2018

F. Abrard and Y. Deville ,

“A Time-Frequency Blind Signal Separation Method Applicable to Underdetermined Mixtures of Dependent Sources, ”

Signal Processing, Vol. 85, No. 7, pp. 1389-1403, 2005

Q. Guo and G. Q. Ruan ,

“A Complex-Valued Mixing Matrix Estimation Algorithm for Underdetermined Blind Source Separation, ”

Circuits Systems and Signal Processing,Vol.37,No. 8,pp. 3206-3226,2018

[Cited within: 1]

F. Lu, Z. Huang, W. Jiang ,

“Underdetermined Blind Separation of Non-Disjoint Signals in Time-FrequencyDomain based on Matrix Diagonalization, ”

Signal Processing, Vol. 91, No. 7, pp. 1568-1577, 2011

[Cited within: 2]

V. G. Reju, S. N. Koh, I. Y. Soon ,

“An Algorithm for Mixing Matrix Estimation in Instantaneous Blind Source Separation, ”

Signal Processing, Vol. 89, No. 9, pp. 1762-1773, 2009

[Cited within: 1]

T. Dong, Y. Lei, J. Yang ,

“An Algorithm for Underdetermined Mixing Matrix Estimation, ”

Neurocomputing, Vol. 104, pp. 26-34, 2013

[Cited within: 1]

/