CN109586728B - Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes - Google Patents

Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes Download PDF

Info

Publication number
CN109586728B
CN109586728B CN201811509899.4A CN201811509899A CN109586728B CN 109586728 B CN109586728 B CN 109586728B CN 201811509899 A CN201811509899 A CN 201811509899A CN 109586728 B CN109586728 B CN 109586728B
Authority
CN
China
Prior art keywords
signal
sparse
matrix
sampling
reconstruction
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
CN201811509899.4A
Other languages
Chinese (zh)
Other versions
CN109586728A (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 CN201811509899.4A priority Critical patent/CN109586728B/en
Publication of CN109586728A publication Critical patent/CN109586728A/en
Application granted granted Critical
Publication of CN109586728B publication Critical patent/CN109586728B/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/3059Digital compression and data reduction techniques where the original information is represented by a subset or similar information, e.g. lossy compression
    • H03M7/3062Compressive sampling or sensing

Abstract

A signal blind reconstruction method under a modulation broadband converter frame based on sparse Bayes is used in the technical field of reconstruction of compressed sensing signals. The invention solves the problem of poor reconstruction performance of the reconstruction method under the framework of the existing modulation broadband converter when the signal contains noise. Firstly multiplying an input sparse signal by a pseudorandom sequence, then carrying out low-speed sampling and filtering operation on the multiplied signal, then constructing an observation matrix, expressing the signal into a compressed sensing expression form, estimating the signal by adopting a sparse Bayesian method during recovery, iteratively solving a variance gamma of the input sparse signal by an EM (effective magnetic resonance) algorithm, and finishing the reconstruction of the sparse signal; under the condition that the signal-to-noise ratio is equal to-15 dB, compared with the prior method, the reconstruction method can reduce the steady state mean square error value by over 75 percent, and effectively improves the reconstruction performance. The invention can be applied to the field of reconstruction of compressed sensing signals.

Description

Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes
Technical Field
The invention belongs to the technical field of reconstruction of compressed sensing signals, and particularly relates to a signal blind reconstruction method under a modulation broadband converter framework.
Background
The CS theorists have proposed Analog Information Conversion (AIC) based on compressed sensing, one of which is mainly the sub-nyquist sampler based on random demodulation, and this scheme mainly aims at multi-tone signals with multiple discrete components, and has the advantages of low power consumption, simple circuit structure, etc., but its disadvantages are also obvious, in practice, because the frequency spectrum of the bandwidth is wide, it can be regarded as infinite discrete components, if a random demodulation scheme is adopted, the frequency spectrum should be discretized first, if the selected frequency resolution is too low, serious spectrum distortion will be caused, and if the selected resolution is higher, the storage and operation will be too heavy. In addition, studies have shown that random demodulation schemes are sensitive to signal models, and large errors can occur when the true sparse basis and the estimated sparse basis deviate.
In 2008, israeli proposed a modulation broadband converter (MWC), and derived and realized related formulas and signal reconstruction, and realized the construction work of hardware. Compared with other analog information converters, MWCs have the following main advantages: the selected signals are all frequency domain continuous signals which are common in life, the requirement on the matching degree of a signal model is not very high, and the existing low-speed ADC can be completely realized. The method has the advantages that the method has disadvantages, and the reconstruction method under the existing modulation broadband converter framework has certain defects, for example, the reconstruction condition of the existing reconstruction method is harsh, so that the reconstruction rate is low, the signal to noise ratio requirement on the signal is high, and the signal reconstruction performance is poor when the signal contains noise.
Disclosure of Invention
The invention aims to solve the problem that the reconstruction performance of the reconstruction method under the framework of the existing modulation broadband converter is poor when a signal contains noise.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a signal blind reconstruction method under a modulation broadband converter framework based on sparse Bayes comprises the following steps:
the method comprises the following steps: inputting a sparse signal x (t), and respectively connecting the input sparse signal x (t) with the pseudo random sequences of m channels of the modulation broadband converterColumn multiplication, the pseudo-random sequences of each channel being mutually orthogonal; then the multiplied result corresponding to each channel is passed through F 0 The sampling frequency of the' is sampled to obtain the sampling result of each channel, and the sampling result of each channel is filtered by a low-pass filter with the cut-off frequency fs/2 to obtain the sampling value y output by each channel i (n) a frequency domain DTFT;
step two: for the obtained sampling value y of each channel output i (n) performing windowing to obtain a windowed signal;
step three: adding Gaussian white noise to the signal subjected to windowing processing obtained in the step two to obtain a signal added with the Gaussian white noise;
step four: calculating an expression of the observation matrix A, and expressing the signal added with the Gaussian white noise obtained in the step three by using a compressed sensing model;
step five: and calculating the edge probability and the posterior probability of the signal added with the Gaussian white noise, applying a sparse Bayesian algorithm to the compressed sensing model obtained in the step four, and iteratively solving the variance gamma of the input sparse signal through an EM (effective-energy-efficient) algorithm to complete the reconstruction of the sparse signal.
The beneficial effects of the invention are: the invention provides a signal blind reconstruction method under a modulation broadband converter frame based on sparse Bayes, and particularly discloses a method for combining sparse Bayes with a modulation broadband converter, wherein the method comprises the steps of multiplying an input sparse signal by a pseudorandom sequence, then carrying out low-speed sampling and filtering operation on the signal obtained by multiplication, then constructing an observation matrix, representing the signal into a compressed sensing representation form, estimating the signal by adopting a sparse Bayes method during recovery, and iteratively solving a variance gamma of the input sparse signal by an EM (effective electromagnetic) algorithm to complete the reconstruction of the sparse signal; the invention combines the sparse Bayes with the modulation broadband converter framework to realize the reconstruction under the condition that the signal sparsity is unknown, and compared with the prior art, the reconstruction method can reduce the steady state mean square error value by more than 75 percent under the condition that the signal to noise ratio of the signal is-15 dB, thereby effectively improving the reconstruction performance.
Drawings
FIG. 1 is a time domain waveform of an original signal in the presence of noise;
FIG. 2 is a time domain waveform of a reconstructed signal in the presence of noise;
the abscissa of fig. 1 and 2 is time and the ordinate is signal amplitude;
FIG. 3 is a frequency domain waveform of an original signal in the presence of noise;
FIG. 4 is a frequency domain waveform of a reconstructed signal in the presence of noise;
the abscissa of fig. 3 and 4 is frequency and the ordinate is signal amplitude;
fig. 5 is a graph comparing the reconstruction performance using the method of the present invention with that using the conventional method.
The abscissa of fig. 5 is the signal-to-noise ratio and the ordinate is the steady state mean square error.
Detailed Description
The technical solution of the present invention is further described below with reference to the accompanying drawings, but not limited thereto, and any modification or equivalent replacement of the technical solution of the present invention without departing from the spirit and scope of the technical solution of the present invention shall be covered by the protection scope of the present invention.
The first embodiment is as follows: the blind signal reconstruction method based on the sparse Bayesian modulation broadband converter framework in the embodiment comprises the following steps:
the method comprises the following steps: inputting sparse signals x (t), and multiplying the input sparse signals x (t) with pseudo-random sequences of m channels of a modulation broadband converter respectively, wherein the pseudo-random sequences of each channel are mutually orthogonal; then the multiplied result corresponding to each channel is passed through F 0 The sampling frequency of the' is sampled to obtain the sampling result of each channel, and the sampling result of each channel is filtered by a low-pass filter with the cut-off frequency fs/2 to obtain the sampling value y output by each channel i (n) a frequency domain DTFT;
step two: for the obtained sampling value y of each channel output i (n) performing windowing to obtain a windowed signal;
step three: adding Gaussian white noise to the signal subjected to windowing processing obtained in the step two to obtain a signal added with the Gaussian white noise;
step four: calculating an expression of the observation matrix A, and expressing the signal added with the Gaussian white noise obtained in the third step by using a compressed sensing model;
step five: and calculating the edge probability and the posterior probability of the signal added with the Gaussian white noise, applying a sparse Bayesian algorithm to the compressed sensing model obtained in the step four, and iteratively solving the variance gamma of the input sparse signal through an EM (effective-energy-efficient) algorithm to complete the reconstruction of the sparse signal.
The second embodiment is as follows: the embodiment further defines the signal blind reconstruction method based on the modulation broadband converter framework based on the sparse bayes in the first embodiment, and the specific process of the first step is as follows:
the pseudo-random sequence of the ith channel of the modulated wideband converter is p i (t) the pseudo-random sequence p is known from Fourier transform i The specific expression of (t) is as follows:
Figure BDA0001900447760000031
wherein: l is a Fourier series, c il Is a Fourier coefficient, j is a complex unit, T P Is the period of the pseudorandom sequence, t is time;
obtaining a Fourier coefficient c according to the inverse Fourier transform il The expression of (a) is:
Figure BDA0001900447760000032
inputting a sparse signal x (t), multiplying the input sparse signal x (t) with the pseudo-random sequences of m channels of the modulation broadband converter respectively to obtain a frequency domain Y of the ith channel of a multiplied result i (f) The expression of (a) is:
Figure BDA0001900447760000033
wherein: f is the frequency domain;
by sampling frequency F 0 Sampling the obtained frequency domain of each channel to obtain a sampling result;
the sampling result is passed through a low-pass filter whose cut-off frequency is fs/2, and the high-frequency component of the filtered signal is filtered out and only the frequency of 0 to fs/2 is retained to obtain the filtered signal, i.e. the sampling value y outputted by the ith channel is obtained by filtering i The expression of frequency domain DTFT of (n) is:
Figure BDA0001900447760000041
wherein:
Figure BDA0001900447760000042
sample value y representing output of ith channel i Frequency domain of (n), T s Is F 0 Reciprocal of `, f p Is the frequency of the pseudorandom sequence, X (-) is the frequency domain of X (t), L 0 Is the smallest positive integer, L, that causes all non-zero values of the sparse signal x (t) to be contained in the low-pass filter 0 The expression of (a) is:
Figure BDA0001900447760000043
wherein: f. of nyq At the nyquist sampling rate.
The third concrete implementation mode: the embodiment further defines the signal blind reconstruction method under the framework of the modulation broadband converter based on the sparse bayes in the first embodiment, and the specific process of the fourth step is as follows:
and (4) expressing the signal added with the Gaussian white noise obtained in the step three by using a compressed sensing model as follows:
y(f)=Az(f) (5)
wherein: intermediate variable z (f) = [ z = 1 (f),...,z L (f)] T And z is i′ (f)=X(f+(i′-L 0 -1)f p ),1≤i′≤L,L=2L 0 +1; y (f) is a vector of length m, and
Figure BDA0001900447760000044
a is an observation matrix;
after sampling and filtering, the Fourier coefficient is represented by c il Is changed into c il ′:
Figure BDA0001900447760000045
Wherein: a is ik Is the value of the pseudo-random sequence of the ith channel, k =0,1, \8230;, L 0 -1;
Defining the integral term of equation (6) as d l
Figure BDA0001900447760000051
Wherein: intermediate variables
Figure BDA0001900447760000052
Then
Figure BDA0001900447760000053
The observation matrix a is then expressed as:
A=SFD (8)
wherein: the intermediate variable matrix F is L x L dimensional matrix, and the ith column of the matrix F is represented as F i″ =[θ 01*i″ ,...,θ (L-1)*i″ ] T Is representative of multiplying, -L 0 <i″<L 0 (ii) a D is a diagonal array of L rows and L columns, and the diagonal array D is in the form of
Figure BDA0001900447760000054
S is a symbol matrix with m rows and L columns, and the form of the symbol matrix S is
Figure BDA0001900447760000055
Therefore, the formula y (f) = Az (f) is converted into the form of formula (9):
Figure BDA0001900447760000056
the fourth concrete implementation mode: the embodiment further defines the signal blind reconstruction method under the framework of the modulation broadband converter based on the sparse bayes in the first embodiment, and the specific process of the fifth step is as follows:
step five, the edge probability of the signal after the Gaussian white noise is added is expressed as follows:
Figure BDA0001900447760000061
wherein: p (Y) ·s ;γ,σ 2 ) The representative parameters are gamma and sigma 2 Edge probability of p (Y) ·s |X ·s ;σ 2 ) Representative parameter is σ 2 Conditional probability of (A), p (X) ·s (ii) a γ) represents the prior probability of the parameter being γ;
Y ·s sampled value y for output of ith channel i At the s-th column of the frequency domain DTFT of (n), γ is the variance of the input sparse signal, σ 2 For the output sampled value y i (n) a variance; x ·s Is the s-th column of the input sparse signal;
Figure BDA0001900447760000062
is a matrix Y ·s The transpose of (a) is performed,
Figure BDA0001900447760000063
is a matrix sigma Y The inverse matrix of (d); edge probability covariance matrix sigma Y The expression of (a) is:
Σ Y =σ 2 I+AΓA T (11)
wherein: Γ = diag (γ), Γ being the diagonal matrix of the matrix γ; a is an observation matrix, I is an identity matrix,
posterior probability p (X) of signal after adding Gaussian white noise ·s |Y ·s ;γ,σ 2 ) Expressed as: the posterior probability is also Gaussian distribution
Figure BDA0001900447760000064
Wherein: p (X) ·s |Y ·s ;γ,σ 2 ) The representative parameters are gamma and sigma 2 A posterior probability of (d), mu ·s Is an average value, and
Figure BDA0001900447760000065
posterior probability covariance matrix
Figure BDA0001900447760000066
The expression of (a) is:
Figure BDA0001900447760000067
fifthly, randomly initializing the variance gamma of the input sparse signal and the output sampling value y i (n) variance σ 2
For the first iteration, by formula
Figure BDA0001900447760000068
To calculate the variance gamma of the input sparse signal corresponding to the first iteration 1 (ii) a Wherein:
Figure BDA0001900447760000069
||·|| 2 represents a norm;
intermediate variables
Figure BDA00019004477600000610
Σ ii Is composed of
Figure BDA00019004477600000611
Any element on the diagonal of (1);
by passingMaximum a posteriori probability to update
Figure BDA00019004477600000612
And mu ·s
Figure BDA00019004477600000613
Figure BDA00019004477600000614
For the second iteration, update Σ by equation (14) ii Reuse of updated sigma ii And mu ·s To calculate the variance gamma of the input sparse signal corresponding to the second iteration 2
Figure BDA0001900447760000071
Then according to the updated gamma 2 To update equations (14) and (15);
and so on until the variance gamma of the input sparse signal converges to a certain fixed point gamma * Stopping iteration, obtaining maximum posterior probability by using gamma value when stopping iteration, and obtaining mean value mu by using maximum posterior probability ·s Mean value of μ ·s Namely the most sparse data to be recovered, and the reconstruction of the sparse signal is completed.
The effects of the invention are explained with the attached drawings: as shown in fig. 1 and 2, fig. 1 and 2 are a time domain waveform diagram of an original signal and a time domain waveform diagram of a reconstructed signal in the presence of noise, respectively; as shown in fig. 3 and 4, fig. 3 and 4 are a frequency domain waveform diagram of an original signal and a frequency domain waveform diagram of a reconstructed signal, respectively, in the presence of noise;
the method can recover the sparse signal under the condition of unknown sparsity of the signal, as shown in fig. 5, compared with the traditional reconstruction method, under the condition that the signal-to-noise ratio is-15 dB, the reconstruction method can reduce the steady state mean square error value by over 75 percent, and the method can effectively reduce the sparsity of the signal and improve the reconstruction performance.
The fifth concrete implementation mode is as follows: the present embodiment further defines the signal blind reconstruction method under the framework of the modulation wideband converter based on the sparse bayes in the first or second embodiment, where F is 0 'the value range is [ fs/4,fs/2']。
The sixth specific implementation mode: in this embodiment, the signal blind reconstruction method under the framework of the modulation wideband converter based on the sparse bayes in the first embodiment is further defined, and in the second step, the sampling value y output by each channel is obtained i And (n) performing windowing, wherein a windowing function adopted in the windowing is a Hamming window. The fence effect of the signal is avoided.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art should be considered to be within the technical scope of the present invention, and the technical solutions and the inventive concepts thereof according to the present invention should be equivalent or changed within the scope of the present invention.

Claims (3)

1. A signal blind reconstruction method under a modulation broadband converter framework based on sparse Bayes is characterized by comprising the following steps:
the method comprises the following steps: inputting sparse signals x (t), and multiplying the input sparse signals x (t) with pseudo-random sequences of m channels of a modulation broadband converter respectively, wherein the pseudo-random sequences of each channel are mutually orthogonal; then the multiplied result corresponding to each channel is passed through F 0 The sampling frequency of the sampling frequency is sampled to obtain the sampling result of each channel, and the sampling result of each channel is filtered by a low-pass filter with the cut-off frequency of fs/2 to obtain the sampling value y output by each channel i (n) a frequency domain DTFT;
the specific process of the step one is as follows:
the pseudo-random sequence of the ith channel of the modulated wideband converter is p i (t) obtaining a pseudorandom sequence p by Fourier transform i The specific expression of (t) is as follows:
Figure FDA0003817041810000011
wherein: l is a Fourier series, c il Is a Fourier coefficient, j is a complex unit, T P Is the period of the pseudorandom sequence, t is time;
obtaining a Fourier coefficient c according to the inverse Fourier transform il The expression of (a) is:
Figure FDA0003817041810000012
inputting a sparse signal x (t), multiplying the input sparse signal x (t) with the pseudo-random sequences of m channels of the modulation broadband converter respectively to obtain a frequency domain Y of the ith channel of a multiplied result i (f) The expression of (a) is:
Figure FDA0003817041810000013
wherein: f is the frequency domain;
sampling the obtained frequency domain of each channel through a sampling frequency F to obtain a sampling result;
and the sampling result is passed through low-pass filter whose cut-off frequency is fs/2 to obtain filtered signal, i.e. filtering to obtain sampling value y outputted by ith channel i The expression of frequency domain DTFT of (n) is:
Figure FDA0003817041810000014
wherein:
Figure FDA0003817041810000015
sample value y representing output of ith channel i Frequency domain of (n), T s Is the reciprocal of F, F p Is the frequency of the pseudorandom sequence, X (-) is the frequency domain of X (t), L 0 All non-zero values of the sparse signal x (t) are included in the low-pass filteringMinimum positive integer of the device, L 0 The expression of (c) is:
Figure FDA0003817041810000021
wherein: f. of nyq At the nyquist sampling rate;
step two: for the obtained sampling value y of each channel output i (n) performing windowing to obtain a windowed signal;
step three: adding Gaussian white noise to the signal subjected to windowing processing obtained in the step two to obtain a signal added with the Gaussian white noise;
step four: calculating an expression of the observation matrix A, and expressing the signal added with the Gaussian white noise obtained in the third step by using a compressed sensing model;
the specific process of the step four is as follows:
and (3) expressing the signal added with the Gaussian white noise obtained in the step three by using a compressed sensing model as follows:
y(f)=Az(f) (5)
wherein: intermediate variable z (f) = [ z = 1 (f),...,z L (f)] T And z is i′ (f)=X(f+(i′-L 0 -1)f p ),1≤i′≤L,L=2L 0 +1; y (f) is a vector of length m, and
Figure FDA0003817041810000022
a is an observation matrix;
after sampling and filtering, the Fourier coefficient is represented by c il Is changed into c il ′:
Figure FDA0003817041810000023
Wherein: a is ik Is the value of the pseudo-random sequence of the ith channel, k =0,1, \8230;, L 0 -1;
Defining the integral term of equation (6) as d l
Figure FDA0003817041810000031
Wherein: intermediate variables
Figure FDA0003817041810000032
Then
Figure FDA0003817041810000033
The observation matrix a is then expressed as:
A=SFD (8)
wherein: the intermediate variable matrix F is a L-dimension matrix, and the ith' column of the matrix F is represented as F i″ =[θ 01*i″ ,...,θ (L-1)*i″ ] T ,-L 0 <i″<L 0 (ii) a D is a diagonal array of L rows and L columns, and the diagonal array D is in the form of
Figure FDA0003817041810000034
S is a symbol matrix with m rows and L columns, and the form of the symbol matrix S is
Figure FDA0003817041810000035
Therefore, the formula y (f) = Az (f) is converted into the form of formula (9):
Figure FDA0003817041810000036
step five: calculating the edge probability and the posterior probability of the signal added with the Gaussian white noise, applying a sparse Bayesian algorithm to the compressed sensing model obtained in the step four, iteratively solving the variance gamma of the input sparse signal through an EM (effective electromagnetic) algorithm, and completing the reconstruction of the sparse signal;
the concrete process of the step five is as follows:
step five, the edge probability of the signal after the Gaussian white noise is added is expressed as follows:
Figure FDA0003817041810000041
wherein: p (Y) ·s ;γ,σ 2 ) Representative parameters are γ and σ 2 Edge probability of p (Y) ·s |X ·s ;σ 2 ) The representative parameter is σ 2 Conditional probability of (A), p (X) ·s (ii) a γ) represents the prior probability of the parameter being γ;
Y ·s sample value y for output of ith channel i At the s-th column of the frequency domain DTFT of (n), γ is the variance of the input sparse signal, σ 2 For the output sampled value y i (n) variance; x ·s Is the s-th column of the input sparse signal;
Figure FDA0003817041810000042
is a matrix Y ·s The method (2) is implemented by the following steps,
Figure FDA0003817041810000043
is a matrix sigma Y The inverse matrix of (d); edge probability covariance matrix sigma Y The expression of (a) is:
Σ Y =σ 2 I+AΓA T (11)
wherein: Γ = diag (γ), Γ being the diagonal matrix of the matrix γ; a is an observation matrix, I is an identity matrix,
posterior probability p (X) of signal after adding Gaussian white noise ·s |Y ·s ;γ,σ 2 ) Expressed as:
Figure FDA0003817041810000044
wherein: p (X) ·s |Y ·s ;γ,σ 2 ) Representative parameters are γ and σ 2 A posterior probability of (d), mu ·s Is an average value, and
Figure FDA0003817041810000045
posterior probability covariance matrix
Figure FDA0003817041810000046
The expression of (c) is:
Figure FDA0003817041810000047
fifthly, randomly initializing the variance gamma of the input sparse signal and the output sampling value y i (n) variance σ 2
For the first iteration, by formula
Figure FDA0003817041810000048
To calculate the variance gamma of the input sparse signal corresponding to the first iteration 1 (ii) a Wherein:
Figure FDA0003817041810000049
||·|| 2 represents a norm;
intermediate variables
Figure FDA00038170418100000410
Σ ii Is composed of
Figure FDA00038170418100000411
Any element on the diagonal of (1);
updating by maximum a posteriori probability
Figure FDA00038170418100000412
And mu ·s
Figure FDA00038170418100000413
Figure FDA00038170418100000414
For the second iteration, update Σ by equation (14) ii Reuse of updated sigma ii And mu ·s To calculate the variance gamma of the input sparse signal corresponding to the second iteration 2
Figure FDA0003817041810000051
Then according to the updated gamma 2 To update equations (14) and (15);
and so on until the variance gamma of the input sparse signal converges to a certain fixed point gamma * Stopping iteration, obtaining maximum posterior probability by using gamma value when iteration is stopped, and obtaining mean value mu by using maximum posterior probability ·s Mean value of μ ·s Namely the most sparse data to be recovered, and the reconstruction of the sparse signal is completed.
2. The sparse Bayesian-based blind reconstruction method for signals under the framework of modulated wideband converters as claimed in claim 1, wherein F is 0 The value range of (b) is [ fs/4,fs/2]。
3. The method for blind reconstruction of signals under the framework of the sparse Bayesian-based modulation wideband converter as claimed in claim 1, wherein the sampling value y of each channel output obtained in the second step i And (n) performing windowing, wherein a windowing function adopted in the windowing is a Hamming window.
CN201811509899.4A 2018-12-11 2018-12-11 Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes Active CN109586728B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811509899.4A CN109586728B (en) 2018-12-11 2018-12-11 Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811509899.4A CN109586728B (en) 2018-12-11 2018-12-11 Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes

Publications (2)

Publication Number Publication Date
CN109586728A CN109586728A (en) 2019-04-05
CN109586728B true CN109586728B (en) 2022-10-25

Family

ID=65928939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811509899.4A Active CN109586728B (en) 2018-12-11 2018-12-11 Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes

Country Status (1)

Country Link
CN (1) CN109586728B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110896308B (en) * 2019-10-31 2023-09-12 中国工程物理研究院电子工程研究所 Single-tone signal reconstruction method
CN110943954B (en) * 2019-11-28 2020-10-23 北京理工大学 Greedy algorithm based mixing sequence optimization method for modulation broadband converter
CN111046329B (en) * 2019-12-10 2022-04-01 西安电子科技大学 Frequency hopping signal block sparse Bayesian reconstruction method based on edge maximum likelihood

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101464513A (en) * 2007-12-21 2009-06-24 中国电子科技集团公司第五十研究所 Non-continuous spectrum high-frequency radar range sidelobe suppression apparatus and method thereof
CN101662442A (en) * 2003-07-24 2010-03-03 科达无线私人有限公司 Method and system for communication in a multiple access network
CN102421105A (en) * 2012-01-12 2012-04-18 北京邮电大学 Method for modeling ultra wide-band (UWB) compressive sensing wireless channel
CN103401556A (en) * 2013-08-13 2013-11-20 哈尔滨工业大学 Comparator-based acquisition system and method for multiple parallel signals and signal reconstruction method
CN104767535A (en) * 2015-03-31 2015-07-08 电子科技大学 Low-complexity massive sparse signal reconstruction method
CA2942585A1 (en) * 2014-03-20 2015-09-24 Huawei Technologies Co., Ltd. Compressive sensing-based signal processing method and apparatus
CN105281779A (en) * 2015-11-04 2016-01-27 哈尔滨工业大学 Multiple observed value vector sparsity self-adaptive compressed sampling matching pursuit method
CN105515695A (en) * 2015-12-04 2016-04-20 哈尔滨工程大学 Compressed sampling signal detection method based on modulated wideband converter
WO2016134549A1 (en) * 2015-02-28 2016-09-01 华为技术有限公司 Compressed sampling method and apparatus
CN105933008A (en) * 2016-04-15 2016-09-07 哈尔滨工业大学 Multiband signal reconstruction method based on clustering sparse regularization orthogonal matching tracking algorithm
CN106093908A (en) * 2016-08-09 2016-11-09 西安电子科技大学 A kind of radar target detection method based on piecemeal segmentation AIC model
CN106663316A (en) * 2016-08-30 2017-05-10 深圳大学 Block sparse compressive sensing-based infrared image reconstruction method and system thereof

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101159136A (en) * 2007-11-13 2008-04-09 中国传媒大学 Low bit rate music signal coding method
KR101503201B1 (en) * 2011-01-10 2015-03-16 알까뗄 루슨트 Method and apparatus for measuring and recovering sparse signals
US8885783B2 (en) * 2011-07-29 2014-11-11 Telefonaktiebolaget L M Ericsson (Publ) Method and device for iterative blind wideband sampling
CN102833020B (en) * 2012-09-10 2015-01-28 杭州电子科技大学 Bayes compression broadband frequency spectrum detection method in cognitive radio network based on self-adaptive measurement
CN106878997B (en) * 2017-01-18 2020-08-14 北京邮电大学 Node sampling rate self-adaptive adjusting method based on distributed network

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101662442A (en) * 2003-07-24 2010-03-03 科达无线私人有限公司 Method and system for communication in a multiple access network
CN101464513A (en) * 2007-12-21 2009-06-24 中国电子科技集团公司第五十研究所 Non-continuous spectrum high-frequency radar range sidelobe suppression apparatus and method thereof
CN102421105A (en) * 2012-01-12 2012-04-18 北京邮电大学 Method for modeling ultra wide-band (UWB) compressive sensing wireless channel
CN103401556A (en) * 2013-08-13 2013-11-20 哈尔滨工业大学 Comparator-based acquisition system and method for multiple parallel signals and signal reconstruction method
CA2942585A1 (en) * 2014-03-20 2015-09-24 Huawei Technologies Co., Ltd. Compressive sensing-based signal processing method and apparatus
WO2016134549A1 (en) * 2015-02-28 2016-09-01 华为技术有限公司 Compressed sampling method and apparatus
CN104767535A (en) * 2015-03-31 2015-07-08 电子科技大学 Low-complexity massive sparse signal reconstruction method
CN105281779A (en) * 2015-11-04 2016-01-27 哈尔滨工业大学 Multiple observed value vector sparsity self-adaptive compressed sampling matching pursuit method
CN105515695A (en) * 2015-12-04 2016-04-20 哈尔滨工程大学 Compressed sampling signal detection method based on modulated wideband converter
CN105933008A (en) * 2016-04-15 2016-09-07 哈尔滨工业大学 Multiband signal reconstruction method based on clustering sparse regularization orthogonal matching tracking algorithm
CN106093908A (en) * 2016-08-09 2016-11-09 西安电子科技大学 A kind of radar target detection method based on piecemeal segmentation AIC model
CN106663316A (en) * 2016-08-30 2017-05-10 深圳大学 Block sparse compressive sensing-based infrared image reconstruction method and system thereof
WO2018039904A1 (en) * 2016-08-30 2018-03-08 深圳大学 Block sparse compressive sensing based infrared image reconstruction method and system thereof

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A Novel Signal Detection Scheme Based on Wavelet Packet Denoising and Cyclic Spectrum Theory;Yulong Gao;《2012 IEEE 12th International Conference on Computer and Information Technology》;20121231;全文 *
Sub-Nyquist Spectrum Sensing Based on Modulated Wideband Converter in Cognitive Radio Sensor Networks;XUE WANG;《IEEE Access》;20180725;全文 *
基于MWC的亚奈奎斯特采样理论和应用技术研究;孙伟朝;《中国优秀博士学位论文全文数据库信息科技辑》;20171115;全文 *
基于压缩感知的多频带信号重构算法设计及实现;史瑶;《中国优秀硕士学位论文全文数据库信息科技辑》;20180215;全文 *
高分辨直升机载旋转合成孔径雷达成像算法;李东;《系统工程与电子技术》;20131231;全文 *

Also Published As

Publication number Publication date
CN109586728A (en) 2019-04-05

Similar Documents

Publication Publication Date Title
CN109586728B (en) Signal blind reconstruction method under modulation broadband converter framework based on sparse Bayes
CN110096956B (en) Signal denoising method and device based on EEMD and permutation entropy second-order difference
Carrillo et al. Iterative hard thresholding for compressed sensing with partially known support
CN110909480B (en) Denoising method and denoising device for water turbine vibration signal
CN112731306A (en) UWB-LFM signal parameter estimation method based on CS and simplified FrFT
Vargas et al. Electrocardiogram signal denoising by clustering and soft thresholding
CN113238193B (en) SAR echo broadband interference suppression method based on multi-component joint reconstruction
CN109165432B (en) Magnetic resonance spectrum reconstruction method based on partial singular value sum
Vimala et al. Noise reduction based on double density discrete wavelet transform
CN111982489B (en) Weak fault feature extraction method for selectively integrating improved local feature decomposition
CN111582205B (en) Noise reduction method based on multi-resolution singular value decomposition model
CN110944336B (en) Time spectrum sensing method based on limited new information rate
CN109724693B (en) Fusion spectrum denoising method based on stationary wavelet
CN109558857B (en) Chaotic signal noise reduction method
CN111695444A (en) Radiation source individual feature extraction method based on wave atomic transformation
CN116502042A (en) Power quality disturbance denoising method based on variational modal decomposition and improved wavelet threshold
Kulkarni et al. Periodic signal denoising: An analysis-synthesis framework based on Ramanujan filter banks and dictionaries
Smith et al. An optimized de-noising scheme based on the null hypothesis of intrinsic mode functions
CN107612553B (en) Signal index regeneration nuclear sparse sampling method with random pulse positions
Kuang et al. Reconstructing signal from quantized signal based on singular spectral analysis
Novikov-Borodin Elimination of Systematic Distortions in the Signals of LTI Systems using Test Measurements
Yang et al. Decimations of intrinsic mode functions via semi-infinite programming based optimal adaptive nonuniform filter bank design approach
Liu et al. Converged Deep Framework Assembling Principled Modules for CS-MRI
CN115065367B (en) Low-complexity sparse signal sampling reconstruction system and method
CN115249028B (en) Blind deconvolution wireless communication signal reconstruction method based on sparse regularization constraint

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