CN109995376B - Signal reconstruction method based on joint block sparse model - Google Patents

Signal reconstruction method based on joint block sparse model Download PDF

Info

Publication number
CN109995376B
CN109995376B CN201910350329.3A CN201910350329A CN109995376B CN 109995376 B CN109995376 B CN 109995376B CN 201910350329 A CN201910350329 A CN 201910350329A CN 109995376 B CN109995376 B CN 109995376B
Authority
CN
China
Prior art keywords
signal
block
reconstruction
sparse
observation matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910350329.3A
Other languages
Chinese (zh)
Other versions
CN109995376A (en
Inventor
高玉龙
轩启运
陈艳平
吴少川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201910350329.3A priority Critical patent/CN109995376B/en
Publication of CN109995376A publication Critical patent/CN109995376A/en
Application granted granted Critical
Publication of CN109995376B publication Critical patent/CN109995376B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/3057Distributed Source coding, e.g. Wyner-Ziv, Slepian Wolf
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention provides a signal reconstruction method based on a joint block sparse model, and belongs to the technical field of distributed compressed sensing. The method comprises the steps of firstly establishing a combined block sparse model based on a mixed support set model, then reconstructing a common part of signals by using the structural characteristics of the combined block sparse model based on the mixed support set model, then reconstructing a unique part of each signal one by using a BOMP algorithm, and finally adding reconstruction results of the common part and the unique part of the original signals to complete reconstruction of the original signals. The invention solves the problem of how to accurately reconstruct the original signal by a receiving end with low measurement value and low signal-to-noise ratio under the condition of multi-antenna and signal sparse coefficient blocking distribution. The invention can be used for signal reconstruction of a receiving end in an actual communication scene.

Description

Signal reconstruction method based on joint block sparse model
Technical Field
The invention relates to a signal reconstruction method, and belongs to the technical field of distributed compressed sensing.
Background
Compressed sensing is a novel signal sampling theory for finding sparse solution of an underdetermined linear system, under the premise of signal sparseness or sparse representation, discrete samples of a signal are obtained by sampling at a sampling rate far lower than the Nyquist theorem, and then the original signal is accurately reproduced by a nonlinear recovery algorithm. The distributed compressed sensing is based on the theory of compressed sensing, the compressed sensing of a single signal is extended to a multi-signal model, the correlation among a plurality of signals is fully utilized, and the condition of joint recovery is created, so that the number of measured values required for successful recovery is further reduced.
The model of distributed compressed sensing is shown in formula (1):
Figure BDA0002043716290000011
Figure BDA0002043716290000012
wherein X ∈ R JN ,Y∈R JM ,Φ∈R MNJ (ii) a X is a received signal, Y is measurement data, and phi is an observation matrix; r JN 、R JM 、R MNJ Respectively representing JN, JM and MNJ dimensional real number vector sets; j is the number of signals in the received signal, M is the number of measured values, and N is the length of the original signal.
The mixed support set model is one of distributed compressed sensing models, and a signal under the model comprises two parts, namely a common part with the same support set and an independent part specific to each signal, and can be described as follows:
Figure BDA0002043716290000013
Figure BDA0002043716290000014
Figure BDA0002043716290000015
wherein J belongs to {1, 2.., J },
Figure BDA0002043716290000016
and
Figure BDA0002043716290000017
common and unique independent portions of each signal, respectively; phi is an observation matrix;
Figure BDA0002043716290000018
and
Figure BDA0002043716290000019
there is a separate non-zero portion.
Figure BDA00020437162900000110
Is provided with
Figure BDA00020437162900000111
The non-zero values are randomly distributed in independent parts of the support set
Figure BDA00020437162900000112
C, removing;
Figure BDA00020437162900000113
with K (c) The number of non-zero values and all signals share the supporting set, i.e.: common part of supporting set
Figure BDA00020437162900000114
j∈{1,2,...J }. The supporting set of each signal is therefore:
Figure BDA0002043716290000021
definition of
Figure BDA0002043716290000022
The total sparsity is greater than the specific partial sparsity.
The multiband signal is a signal with a blocky distribution of frequency domains, which is common in the field of signal processing, and the mathematical model of the multiband signal is shown in equation (7).
Figure BDA0002043716290000023
Wherein sinc (x) = sin (π x)/π x, s (t) comprises L pairs of frequency bands, t is time, E i And τ i Energy and time offsets of the individual bands, i =1, \8230;, L; b is the maximum bandwidth, f i Is the carrier frequency of each frequency band and has a distribution interval of [ -f [ - ] nyq /2,f nyq /2],f nyq Is the Nyquist sampling rate.
In an actual communication scene, signals received by us usually have a joint block sparsity characteristic, and under the condition that multiple antennas and signal sparsity coefficients are distributed in a block mode, how a receiving end accurately reconstructs an original signal with a low measurement value and a low signal-to-noise ratio is a problem to be solved by us.
The current reconstruction algorithm based on block sparse signals is mainly divided into three types:
1. a Mixed l 2/l 1optimization (Mixed l 2/l 1optimization program, L-OPT) algorithm;
2. a Block-sparse Matching Pursuit (BMP) algorithm;
3. block-sparse Orthogonal Matching Pursuit (BOMP) algorithm.
The mixed l 2/l 1optimization algorithm is adopted as the reconstruction algorithm of the block sparse signal, so that the problems of high optimization complexity and difficult practical application exist; a block sparse matching tracking algorithm or a block sparse orthogonal matching tracking algorithm is adopted as a reconstruction algorithm of the block sparse signal, matching atoms are not changed after being determined, an over-matching phenomenon is easy to cause, and the method is not suitable for an actual communication scene.
Disclosure of Invention
The invention provides a signal reconstruction method based on a joint block sparse model, which aims to solve the problem of how to accurately reconstruct an original signal by a receiving end with a low measurement value and a low signal-to-noise ratio under the condition that multiple antennas are arranged and signal sparse coefficients are distributed in blocks.
The invention discloses a signal reconstruction method based on a joint block sparse model, which is realized by the following technical scheme:
step one, establishing a combined block sparse model based on a mixed support set model;
step two, initialization:
setting the iteration times l =1, and setting the vector selection index set omega to be null; let initialization residual r j,0 Is equal to the measured value y j Initializing the resulting sparse coefficient vector
Figure BDA0002043716290000024
Step three, sub-block selection:
selecting an atomic block with the largest average inner product value of all J residual errors in the observation matrix, and recording an index corresponding to the atomic block;
step four, sub-block reconstruction:
measured value y j Projecting in a vector space formed by the updated observation matrix D to obtain a measured value y j Coefficient beta linearly represented by each column vector in D j
Step five, residual error r is updated j,l (ii) a Updating sparse coefficient vectors
Figure BDA0002043716290000031
Step six, judging whether | | | r is satisfied j,l || 2 >ε||y j || 2 And l is less than K, if so, making l = l +1, and returning to the step three(ii) a If not, the Fourier transform matrix is used
Figure BDA0002043716290000032
As a result of reconstruction of the common portion of each original signal; epsilon is a residual error determination coefficient; k is the bulk sparsity;
step seven, respectively using BOMP algorithm to the residual errors, and reconstructing the specific part of each signal one by one;
and step eight, adding the common part of the original signal and the reconstruction result of the special part to complete the reconstruction of the original signal.
The most prominent characteristics and remarkable beneficial effects of the invention are as follows:
the invention relates to a signal reconstruction method based on a combined block sparse model, which takes the actual multi-antenna communication scene into consideration, combines a mixed support set model in distributed compressed sensing with a block sparse structure, establishes the combined block sparse model based on the mixed support set model closer to the actual communication scene, and then divides the reconstruction process into a public part and a special part to be respectively carried out; the method makes full use of the structural characteristics of the combined block sparse model based on the mixed support set model, and jointly utilizes the block sparse characteristic and the combined sparse information of the signals in the sub-block atom selection process, so that the sub-blocks containing non-zero points can be more accurately screened out, and the reconstruction accuracy is improved; the method can reconstruct the original signal with higher precision under the conditions of multi-antenna and signal sparse coefficient blocking distribution by utilizing the structural information of the model under the conditions of low measurement value and low signal to noise ratio. The reconstruction performance under the same measurement value number is obviously superior to the traditional BOMP algorithm, and the anti-noise performance is also superior to the traditional BOMP algorithm. In simulation experiments, the mean square error of reconstruction under the same measured value quantity is about 1.5dB lower than that of the traditional BOMP algorithm.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a diagram of a spectrum of an original signal X1 in an embodiment;
FIG. 3 is a diagram of the spectrum of an original signal X2 in an example;
FIG. 4 is a diagram of the spectrum of an original signal X3 in an example;
FIG. 5 is a diagram illustrating the effect of reconstructing the original signal X1 by using the method of the present invention in the embodiment;
FIG. 6 is a diagram illustrating the effect of reconstructing the original signal X2 according to the method of the present invention in the embodiment;
fig. 7 is a diagram illustrating the effect of reconstructing the original signal X3 by using the method of the present invention in the embodiment;
FIG. 8 is a graph comparing the reconstructed Mean Square Error (MSE) of the BOMP and the number of measurements M in accordance with the present invention;
FIG. 9 is a graph comparing the noise immunity of the method of the present invention with BOMP.
Detailed Description
The first embodiment is as follows: the present embodiment is described with reference to fig. 1, and the method for reconstructing a signal based on a joint block sparse model provided in the present embodiment specifically includes the following steps:
step one, establishing a combined block sparse model based on a mixed support set model;
step two, initialization:
setting the iteration times l =1, and setting the vector selection index set omega to be null; let initialization residual r j,0 Is equal to the measured value y j Initializing the obtained sparse coefficient vector
Figure BDA0002043716290000041
Step three, sub-block selection:
for an observation matrix containing BNum atom blocks, selecting an atom block with the largest average inner product value with all J residual errors, and recording an index corresponding to the atom block;
step four, sub-block reconstruction:
measured value y j Projecting in a vector space formed by the updated observation matrix D to obtain a measured value y j Coefficient beta linearly represented by each column vector in D j
Step five, updating residual error r j,l (ii) a Updating sparse coefficient vectors
Figure BDA0002043716290000042
Step six, judging whether | | | r is satisfied j,l || 2 >ε||y j || 2 If l is less than K, making l = l +1, and returning to the step three; if not, the Fourier transform matrix is used
Figure BDA0002043716290000043
As a result of reconstruction of the common portion of each original signal; stopping iteration when the convergence condition is met; wherein epsilon is a residual error judgment coefficient; k is the bulk sparsity;
step seven, after the public part of each signal is obtained, the BOMP algorithm is respectively used for the residual errors, and the specific part of each signal is reconstructed one by one;
and step eight, adding the reconstruction results of the common part and the specific part of the original signal to finish the accurate reconstruction of the original signal.
In order to fit actual communication scenes better, a combined block sparse model based on a mixed support set model is constructed, and a reconstruction method under the model is provided. The method fully utilizes the structural characteristics of the model, and jointly utilizes the block sparse characteristic and the joint sparse information of the signals in the sub-block atom selection process, so that the sub-blocks containing non-zero points can be screened out more accurately, and the reconstruction accuracy is improved.
The second embodiment is as follows: this embodiment is different from the first embodiment in that,
the combined block sparse model based on the mixed support set model in the first step specifically comprises:
Figure BDA0002043716290000051
Figure BDA0002043716290000052
Figure BDA0002043716290000053
wherein, X = [ X = 1 x 2 ... x J ] T J belongs to {1, 2.., J } for the received signal, and J is the number of signals in the received signal;
Figure BDA0002043716290000054
and
Figure BDA0002043716290000055
common and unique independent portions of each signal, respectively; phi is an observation matrix;
Figure BDA0002043716290000056
is provided with
Figure BDA0002043716290000057
A plurality of non-zero blocks;
Figure BDA0002043716290000058
with K (c) A non-zero block.
Other steps and parameters are the same as those in the first embodiment.
The third concrete implementation mode: the difference between this embodiment and the first or second embodiment is that the third step is specifically:
the solution model for distributed compressive sensing can be expressed as:
Figure BDA0002043716290000059
wherein the content of the first and second substances,
Figure BDA00020437162900000510
denotes x j K is the block sparsity;
according to the traditional BOMP algorithm based on the block sparse model, when sub-block atoms are selected, the principle that the inner product of original dictionary sub-blocks and signal residual errors is the largest is utilized, and when signals have the combined block sparse characteristic, the idea of synchronous orthogonal matching tracking algorithm in distributed compressed sensing can be referred to, and then the original sub-block with the largest average inner product value of all J residual errors in an observation matrix is selected:
Figure BDA00020437162900000511
wherein n is l Representing the index corresponding to the most matched atomic block selected by the iteration;<·>expressing the inner product, | | | · | | expresses a norm; phi j,n To observe the nth column in matrix Φ, N =1, \ 8230;, N; n is the length of the original signal (also the total number of columns of the observation matrix Φ); it can be seen from equation (9) that, in the sub-block atom selection process, the principle that the sum of the residual products of the original sub-blocks of the dictionary and all the J signals is the maximum is utilized, and as the J signals have the same common support set, the method is easier to find the most relevant atomic block than the traditional BOMP algorithm, and even if noise exists, the noise influence can be reduced by taking the sum of the residual products and the residual products of the J signals to be the maximum due to the random distribution of the noise, so that the anti-noise performance is obviously enhanced, the joint sparse information of the signals is fully utilized, and the probability of finding the optimal sub-block is improved.
Then record the index corresponding to the block atom:
Ω=[Ω,n l ] (10)
other steps and parameters are the same as those in the first or second embodiment.
The fourth concrete implementation mode: the present embodiment is different from the third embodiment in that β in the fourth step j The least squares method yields:
Figure BDA0002043716290000061
wherein D is an updated observation matrix composed of columns taken out of the original observation matrix.
Other steps and parameters are the same as those in the third embodiment.
The fifth concrete implementation mode: the difference between this embodiment and the fourth embodiment is that the specific process of updating the residual error in the fifth step is as follows:
Figure BDA0002043716290000062
other steps and parameters are the same as those in the fourth embodiment.
The sixth specific implementation mode: the difference between this embodiment and the fifth embodiment is that the specific process of updating the sparse coefficient vector in the fifth step is as follows:
Figure BDA0002043716290000063
the other steps and parameters are the same as those in the fifth embodiment.
The seventh embodiment: in this embodiment, the difference from the first embodiment is that the residual determination coefficient ∈ =0.01 in step six.
Other steps and parameters are the same as those in the first to sixth embodiments.
Examples
The following examples were employed to demonstrate the beneficial effects of the present invention:
the signal reconstruction method based on the joint block sparse model in the embodiment is carried out according to the following steps:
step one, establishing a combined block sparse model based on a mixed support set model:
Figure BDA0002043716290000064
wherein the content of the first and second substances,
Figure BDA0002043716290000071
is provided with
Figure BDA0002043716290000072
A plurality of non-zero blocks;
Figure BDA0002043716290000073
with K (c) A plurality of non-zero blocks;
step two, inputting: measured data Y = [ Y = [ ] 1 y 2 ... y J ] T Observation matrix phi, signal number J, block sparsity K and block division number BNum;
setting the iteration times l =1, and setting the vector selection index set omega to be null; let initialization residual r j,0 Is equal to the measured value y j Initializing the resulting sparse coefficient vector
Figure BDA0002043716290000074
Step three, sub-block selection:
selecting an atomic block with the maximum average inner product value of all J residual errors in the observation matrix, and recording an index corresponding to the atomic block;
Figure BDA0002043716290000075
step four, sub-block reconstruction:
measured value y j Projecting in a vector space formed by the updated observation matrix D to obtain a measured value y j Coefficient beta linearly represented by each column vector in D j (ii) a Obtained by the least squares method:
Figure BDA0002043716290000076
step five, residual error is updated
Figure BDA0002043716290000077
Updating sparse coefficient vectors
Figure BDA0002043716290000078
Step six, judging whether | | | r is satisfied j,l || 2 >ε||y j || 2 If l is less than K, making l = l +1, and returning to the step three; if not, the Fourier transform matrix is used
Figure BDA0002043716290000079
As a result of reconstruction of the common portion of each original signal; since the sparsity of the common part of the given signal set is K1 and the sparsity of the specific part is K2, K1 iterations are needed to find all the common parts.
Step seven, respectively using a BOMP algorithm to the residual errors, and reconstructing the specific parts of each signal one by one, wherein the iteration times of the process is 2X K2 due to the existence of noise, which is also the case of considering the poor anti-noise performance of the traditional BOMP algorithm;
and step eight, adding the reconstruction results of the common part and the specific part of the original signal to finish the accurate reconstruction of the original signal.
As shown in fig. 5, 6 and 7, which are diagrams of the effect of reconstructing the original signal in fig. 2, 3 and 4 by using the method of the present invention, it can be seen that the method of the present invention has a better reconstruction effect; FIG. 8 is a relationship between the reconstructed Mean Square Error (MSE) of the method (B-SOMP) of the present invention and the number of measured values M, and it can be seen that the performance of the method of the present invention gradually becomes better as the number of measured values increases (the mean square error gradually decreases as the number of measured values increases), and the average reconstructed mean square error at the same number of measured values is about 1.5dB lower than that of the conventional BOMP algorithm; fig. 9 is a comparison of the noise immunity of the method of the present invention (B-SOMP) and the conventional block orthogonal matching pursuit algorithm (BOMP), and it can be seen that the noise immunity of the proposed method is significantly better than the BOMP and the advantage is more significant at low snr.
The present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof, and it is therefore intended that all such changes and modifications be considered as within the spirit and scope of the appended claims.

Claims (2)

1. The signal reconstruction method based on the joint block sparse model is characterized by comprising the following steps:
step one, establishing a combined block sparse model based on a mixed support set model;
the joint block sparse model based on the mixed support set model specifically comprises the following steps:
Figure FDA0003956707060000011
Figure FDA0003956707060000012
Figure FDA0003956707060000013
wherein, X = [ X ] 1 x 2 ...x J ] T J belongs to {1, 2.., J } for receiving signals, and J is the number of signals in the received signals;
Figure FDA0003956707060000014
and
Figure FDA0003956707060000015
common and unique independent portions of each signal, respectively; phi is an observation matrix;
Figure FDA0003956707060000016
is provided with
Figure FDA0003956707060000017
A plurality of non-zero blocks;
Figure FDA0003956707060000018
with K (c) A plurality of non-zero blocks;
step two, initialization:
setting the iteration times l =1, and enabling a vector selection index set omega to be empty; let initialization residual r j,0 Is equal to the measured value y j Initializing the resulting sparse coefficient vector
Figure FDA0003956707060000019
Step three, sub-block selection:
selecting an atomic block with the largest average inner product value of all J residual errors in the observation matrix, and recording an index corresponding to the atomic block, wherein the specific steps are as follows:
selecting the atomic block with the largest average inner product value with all J residual errors in the observation matrix:
Figure FDA00039567070600000110
wherein n is l Representing the index corresponding to the original block with the maximum average inner product value of all J residual errors in the observation matrix; phi (phi) of j,n Is the nth column in the observation matrix phi; n is the length of the original signal,<·>expressing the inner product, | | | · | | expresses a norm;
then recording the index corresponding to the atomic block:
Ω=[Ω,n l ] (10);
step four, sub-block reconstruction:
measured value y j Projecting in a vector space formed by the updated observation matrix D to obtain a measured value y j Coefficient beta linearly represented by each column vector in D j
Beta is said j And is obtained by the least square method:
Figure FDA0003956707060000021
wherein D is an updated observation matrix composed of columns taken out of the original observation matrix;
step five, residual error r is updated j,l (ii) a Updating sparse coefficient vectors
Figure FDA0003956707060000022
The specific process of updating the residual error is as follows:
Figure FDA0003956707060000023
the specific process for updating the sparse coefficient vector comprises the following steps:
Figure FDA0003956707060000024
step six, judging whether | | | r is satisfied j,l || 2 >ε||y j || 2 If l is less than K, making l = l +1, and returning to the step three; if not, the Fourier transform matrix is used
Figure FDA0003956707060000025
As a result of reconstruction of the common portion of each original signal; epsilon is a residual error determination coefficient; k is the bulk sparsity;
step seven, respectively using BOMP algorithm to the residual errors, and reconstructing the specific part of each signal one by one;
and step eight, adding the common part of the original signal and the reconstruction result of the special part to complete the reconstruction of the original signal.
2. The method for reconstructing a signal based on a joint block sparse model as claimed in claim 1, wherein the residual error decision coefficient e =0.01 in step six.
CN201910350329.3A 2019-04-28 2019-04-28 Signal reconstruction method based on joint block sparse model Active CN109995376B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910350329.3A CN109995376B (en) 2019-04-28 2019-04-28 Signal reconstruction method based on joint block sparse model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910350329.3A CN109995376B (en) 2019-04-28 2019-04-28 Signal reconstruction method based on joint block sparse model

Publications (2)

Publication Number Publication Date
CN109995376A CN109995376A (en) 2019-07-09
CN109995376B true CN109995376B (en) 2023-02-03

Family

ID=67135469

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910350329.3A Active CN109995376B (en) 2019-04-28 2019-04-28 Signal reconstruction method based on joint block sparse model

Country Status (1)

Country Link
CN (1) CN109995376B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111812644B (en) * 2020-08-29 2023-05-30 西安电子科技大学 MIMO radar imaging method based on sparse estimation
CN112202452A (en) * 2020-09-02 2021-01-08 北京电子工程总体研究所 Compressed sensing signal reconstruction method and system based on block sparsity and binary tree search
CN112653472A (en) * 2020-12-15 2021-04-13 哈尔滨工程大学 Dolphin whistle call signal reconstruction method based on block sparse compressed sensing

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101895297A (en) * 2010-07-30 2010-11-24 哈尔滨工业大学 Compressed sensing-oriented block-sparse signal reconfiguring method
CN101908889A (en) * 2010-07-30 2010-12-08 哈尔滨工业大学 Compressed sensing reconstructing method of sparse signal with unknown block sparsity
CN102405012A (en) * 2009-01-14 2012-04-04 高通股份有限公司 Distributed sensing of signals linked by sparse filtering
KR101172641B1 (en) * 2011-07-04 2012-08-08 한국과학기술원 Method and apparatus for compressed sensing with joint sparsity
CN104270158A (en) * 2014-03-26 2015-01-07 湘潭大学 A cooperative reconstruction method with adaptive sparsity
CN104333389A (en) * 2014-10-23 2015-02-04 湘潭大学 Adaptive threshold value iterative reconstruction method for distributed compressed sensing
CN104463927A (en) * 2014-12-12 2015-03-25 西安电子科技大学 Compressed sensing image reconstructing method based on multivariable observing and relevance constraining
CN105281779A (en) * 2015-11-04 2016-01-27 哈尔滨工业大学 Multiple observed value vector sparsity self-adaptive compressed sampling matching pursuit method
CN105515585A (en) * 2015-12-08 2016-04-20 宁波大学 Compressed sensing reconstruction method for signals with unknown sparseness
CN105763201A (en) * 2016-02-22 2016-07-13 广东顺德中山大学卡内基梅隆大学国际联合研究院 L1 norm minimization sparse reconstruction method and L1 norm minimization sparse reconstruction system based on weight resetting under low sampling rate
CN106374935A (en) * 2016-08-29 2017-02-01 电子科技大学 Block sparse signal reconstruction method adopting self-adaptive coupling parameters
CN106663316A (en) * 2016-08-30 2017-05-10 深圳大学 Block sparse compressive sensing-based infrared image reconstruction method and system thereof
CN107680036A (en) * 2017-08-15 2018-02-09 湖北工业大学 The joint sparse Vector Parallel method for reconstructing of network is stacked based on convolution depth
CN107977680A (en) * 2017-12-15 2018-05-01 武汉艾孚达信息技术有限公司 A kind of method of the sparse reconstruct of signal
CN108418769A (en) * 2018-01-17 2018-08-17 南京邮电大学 A kind of distributed compression perception adaptive method for reconstructing of degree of rarefication
CN109560841A (en) * 2018-12-13 2019-04-02 东北大学 Extensive mimo system channel estimation methods based on improved distributed compression perception algorithm

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7271747B2 (en) * 2005-05-10 2007-09-18 Rice University Method and apparatus for distributed compressed sensing
US8483492B2 (en) * 2005-10-25 2013-07-09 William Marsh Rice University Method and apparatus for signal detection, classification and estimation from compressive measurements
WO2012044380A2 (en) * 2010-06-18 2012-04-05 William Marsh Rice University Method and apparatus for compressive acquisition and recovery of dynamic imagery
US20120259590A1 (en) * 2011-04-11 2012-10-11 Jong Chul Ye Method and apparatus for compressed sensing with joint sparsity

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102405012A (en) * 2009-01-14 2012-04-04 高通股份有限公司 Distributed sensing of signals linked by sparse filtering
CN101895297A (en) * 2010-07-30 2010-11-24 哈尔滨工业大学 Compressed sensing-oriented block-sparse signal reconfiguring method
CN101908889A (en) * 2010-07-30 2010-12-08 哈尔滨工业大学 Compressed sensing reconstructing method of sparse signal with unknown block sparsity
KR101172641B1 (en) * 2011-07-04 2012-08-08 한국과학기술원 Method and apparatus for compressed sensing with joint sparsity
CN104270158A (en) * 2014-03-26 2015-01-07 湘潭大学 A cooperative reconstruction method with adaptive sparsity
CN104333389A (en) * 2014-10-23 2015-02-04 湘潭大学 Adaptive threshold value iterative reconstruction method for distributed compressed sensing
CN104463927A (en) * 2014-12-12 2015-03-25 西安电子科技大学 Compressed sensing image reconstructing method based on multivariable observing and relevance constraining
CN105281779A (en) * 2015-11-04 2016-01-27 哈尔滨工业大学 Multiple observed value vector sparsity self-adaptive compressed sampling matching pursuit method
CN105515585A (en) * 2015-12-08 2016-04-20 宁波大学 Compressed sensing reconstruction method for signals with unknown sparseness
CN105763201A (en) * 2016-02-22 2016-07-13 广东顺德中山大学卡内基梅隆大学国际联合研究院 L1 norm minimization sparse reconstruction method and L1 norm minimization sparse reconstruction system based on weight resetting under low sampling rate
CN106374935A (en) * 2016-08-29 2017-02-01 电子科技大学 Block sparse signal reconstruction method adopting self-adaptive coupling parameters
CN106663316A (en) * 2016-08-30 2017-05-10 深圳大学 Block sparse compressive sensing-based infrared image reconstruction method and system thereof
CN107680036A (en) * 2017-08-15 2018-02-09 湖北工业大学 The joint sparse Vector Parallel method for reconstructing of network is stacked based on convolution depth
CN107977680A (en) * 2017-12-15 2018-05-01 武汉艾孚达信息技术有限公司 A kind of method of the sparse reconstruct of signal
CN108418769A (en) * 2018-01-17 2018-08-17 南京邮电大学 A kind of distributed compression perception adaptive method for reconstructing of degree of rarefication
CN109560841A (en) * 2018-12-13 2019-04-02 东北大学 Extensive mimo system channel estimation methods based on improved distributed compression perception algorithm

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
Joint Sparse Recovery Method for Compressed Sensing With Structured Dictionary Mismatches;Zhao Tan等;《IEEE Transactions on Signal Processing》;20140729(第19期);4997-5008 *
Learning Joint-Sparse Codes for Calibration-Free Parallel MR Imaging;Shanshan Wang等;《IEEE Transactions on Medical Imaging》;20170829(第1期);251-261 *
Matching pursuits with time-frequency dictionaries;S.G. Mallat等;《IEEE Transactions on Signal Processing》;19931231(第12期);3397-3415 *
分布式压缩感知联合重建算法的研究;候肖兰;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20170115(第01期);I136-55 *
压缩感知增强型自适应分段正交匹配追踪算法;何雪云等;《信号处理》;20180925(第09期);1045-1052 *
基于压缩感知的多用户检测技术研究;杜睿;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20180215(第02期);I140-562 *
基于子空间的块稀疏信号压缩感知重构算法;付宁等;《电子学报》;20111015(第10期);2338-2342 *
基于联合块稀疏模型的随机调频步进ISAR成像方法;吕明久等;《电子与信息学报》;20180815(第11期);2614-2620 *
改进的压缩感知重构算法研究;刘馨月等;《青岛大学学报(自然科学版)》;20180215(第01期);61-68 *

Also Published As

Publication number Publication date
CN109995376A (en) 2019-07-09

Similar Documents

Publication Publication Date Title
CN109995376B (en) Signal reconstruction method based on joint block sparse model
CN101908889B (en) Compressed sensing reconstructing method of sparse signal with unknown block sparsity
CN108322409A (en) Sparse OFDM channel estimation method based on generalized orthogonal matching pursuit algorithm
US20190266761A1 (en) System and method for sparse image reconstruction
CN108957390B (en) Arrival angle estimation method based on sparse Bayesian theory in presence of mutual coupling
CN107945129B (en) MRI image reconstruction method
CN103077510B (en) Multivariate compressive sensing reconstruction method based on wavelet HMT (Hidden Markov Tree) model
CN105277985A (en) OVT-domain seismic data regularization method based on image processing
CN112363110A (en) Grid-free single-bit DOA estimation method based on nested cross dipole array
CN107037409B (en) MIMO radar waveform separation method based on compressed sensing
CN102937706A (en) Magnetic resonance parallel imaging method and magnetic resonance imager
CN112130150B (en) Variable pulse repetition interval SAR imaging method based on compressed sensing
Liu et al. Calibrationless parallel imaging reconstruction for multislice MR data using low‐rank tensor completion
CN110954860B (en) DOA and polarization parameter estimation method
CN108132383B (en) Time-interleaved multi-harmonic signal undersampling method
CN112395546B (en) Finite information rate signal down-sampling and reconstruction method based on linear regular domain
CN109188327B (en) Magnetic resonance image fast reconstruction method based on tensor product complex small compact framework
CN109917330B (en) Arrival angle estimation method based on sparse orthogonal matching pursuit theory in the presence of phase errors
CN111830560B (en) Seismic data reconstruction method based on rank reduction algorithm
CN107403628B (en) Voice signal reconstruction method based on compressed sensing
CN109856673A (en) A kind of high-resolution Radon transformation data separating technology based on the weighting of dominant frequency iteration
Purnamasari et al. Recovery of radar's beat signal based on compressive sampling
Zhao et al. A high quality image reconstruction method based on nonconvex decoding
CN113009560B (en) Seismic data reconstruction method, device, equipment and storage medium
CN112712573B (en) Electromagnetic tomography image reconstruction method based on frequency spectrum characteristics

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant