Abstract:In this paper, we consider the block-sparse signals recovery problem in the context of multiple measurement vectors (MMV) with common row sparsity patterns. We develop a new method for recovery of common row sparsity MMV signals, where a pattern-coupled hierarchical Gaussian prior model is introduced to characterize both the block-sparsity of the coefficients and the statistical dependency between neighboring coefficients of the common row sparsity MMV signals. Unlike many other methods, the proposed method is able to automatically capture the block sparse structure of the unknown signal. Our method is developed using an expectation-maximization (EM) framework. Simulation results show that our proposed method offers competitive performance in recovering block-sparse common row sparsity pattern MMV signals.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to recover block - sparse signals with a common row - sparse pattern in the case of multiple measurement vectors (MMV). Specifically, the author focuses on the situation where in multiple measurement vectors, all signals share the same non - zero element positions (i.e., row - sparse pattern), and these non - zero elements appear in blocks.
### Problem Background
In practical applications, many signals not only have sparsity but also have additional structural characteristics, such as block - sparsity. Block - sparsity means that non - zero elements are not randomly distributed but concentrated in certain specific regions or blocks. This structural characteristic can be used to improve the performance of signal recovery. However, many existing methods mainly focus on the recovery problem of single measurement vectors (SMV) and ignore the common row - sparse pattern in multiple measurement vectors (MMV).
### Research Objectives
The objective of this paper is to develop a new method to recover block - sparse signals with a common row - sparse pattern. To this end, the author introduces a pattern - coupled hierarchical Gaussian prior model, which can not only capture the block - sparsity of coefficients but also describe the statistical dependence between adjacent coefficients.
### Method Innovation Points
1. **Pattern - coupled hierarchical Gaussian prior model**: This model characterizes the correlation between adjacent coefficients by introducing hyper - parameters, thereby better capturing the block - sparse structure of the signal.
2. **Automatically capture block - sparse structure**: Unlike existing methods, this method does not need to know the block partition information in advance but can automatically identify and utilize the block - sparse structure of the signal.
3. **Expectation - Maximization (EM) algorithm**: The author develops an effective algorithm based on the EM framework to estimate hyper - parameters and finally recover the signal.
### Experimental Results
Through simulation experiments, the author verifies the superior performance of the proposed method in recovering block - sparse signals. The experimental results show that this method can significantly outperform traditional sparse Bayesian learning methods (such as MSBL) under different noise levels, especially when the specific positions and sizes of non - zero blocks do not need to be known in advance.
### Conclusions
This paper proposes a new Bayesian method for recovering block - sparse signals with a common row - sparse pattern. By introducing a pattern - coupled hierarchical Gaussian prior model, this method can effectively capture the block - sparse structure of the signal and shows excellent performance in experiments.
### Formula Summary
- Model representation:
\[
y_l=\Phi x_l + v_l,\quad\forall l = 1,\ldots,L
\]
where \(y_l\in\mathbb{R}^{M\times1}\), \(\Phi\in\mathbb{R}^{M\times N}\), \(x_l\in\mathbb{R}^{N\times1}\), \(v_l\in\mathbb{R}^{M\times1}\) respectively represent the \(l\)-th measurement vector, sensing matrix, signal vector and noise vector.
- Gaussian prior distribution:
\[
p(x_{n\cdot})=\mathcal{N}(0,(\alpha_n+\beta\alpha_{n - 1}+\beta\alpha_{n+1})^{-1}B_2^{-1})
\]
- Posterior distribution:
\[
p(X|Y;\alpha,B_2,\lambda)\propto p(Y|X;\lambda)p(X;\alpha,B_2)
\]
- Update formulas:
\[
\mu=\lambda\Sigma(I\otimes\Phi)^T y
\]
\[
\Sigma=(\lambda I\otimes(\Phi^T\Phi)+B_2\otimes B_1)^{-1}
\]
\[
\alpha_i^{(t + 1)}=\frac{3L}{2\phi_i}