CN106649198B - A kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network - Google Patents
A kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network Download PDFInfo
- Publication number
- CN106649198B CN106649198B CN201611018905.7A CN201611018905A CN106649198B CN 106649198 B CN106649198 B CN 106649198B CN 201611018905 A CN201611018905 A CN 201611018905A CN 106649198 B CN106649198 B CN 106649198B
- Authority
- CN
- China
- Prior art keywords
- signal
- wireless sensor
- sensor network
- dimension
- estimation
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
The invention discloses a kind of methods of higher-dimension signal reconstruction quality in detection wireless sensor network, are a kind of computational methods using force zero estimation and the performance indicator for detecting reconstructed signal after Linear Minimum Mean-Square Error Estimation.During signal reconstruction in wireless sensor network, wireless sensor collected signal can be indicated with the form of Fourier expansion.In view of the influence of noise in acquisition signal, need to rebuild signals and associated noises using estimator.When network size is larger, when signal dimension and larger wireless sensor quantity i.e. to be treated, the present invention provides a kind of methods for capableing of the easy probability density function for accurately obtaining generalized circular matrix characteristic value again, to the accurate detection performance index that can easily calculate system.
Description
Technical field
The present invention relates to a kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network, more particularly to one kind exists
Refer to the mean square error for being reconstructed signal is detected after Linear Minimum Mean-Square Error Estimation (LMMSE) method using force zero estimation (ZF)
Target computational methods belong to the technical field of signal reconstruction in wireless sensor network.
Background technology
During signal reconstruction in wireless sensor network, random generalized circular matrix plays important role.Wirelessly
Sensor collected signal can be indicated with the form of Fourier expansion.In view of the geography residing for each sensor
Position is different, and collected signal can be expressed as between a generalized circular matrix and a vector for indicating the signal spectrum
Product.In view of the influence of noise in acquisition signal, need to rebuild signals and associated noises using estimator.
Classical algorithm has force zero estimation and Minimum Mean Squared Error estimation, and original signal is being rebuild using both estimators
Afterwards, it needs to calculate its mean square error to the signal of reconstruction, so that it is determined that the accuracy of signal reconstruction.The calculating of mean square error needs
Solve the probability density function of the characteristic value for the generalized circular matrix mentioned before.
When the bandwidth for acquiring the dimension of signal, the number of sensor and collected signal is relatively low, generalized circular matrix
The probability density function of characteristic value can be estimated by the method for Monte Carlo numerical computations.However, working as above-mentioned parameter
When being worth larger, the calculation amount needed for the numerical computations of Monte Carlo is exponentially increased, it is therefore desirable to one kind can it is easy again accurately
The method for obtaining the probability density function of this matrix exgenvalue.
Invention content
In view of the above-mentioned problems, technical problem to be solved by the invention is to provide higher-dimensions in a kind of detection wireless sensor network
The method of signal reconstruction quality, the present invention provide it is a kind of can the easy probability for accurately obtaining generalized circular matrix characteristic value again it is close
The method for spending function, to the accurate detection performance index that can easily calculate system.
The present invention uses following technical scheme to solve above-mentioned technical problem:
The present invention provides a kind of method detecting higher-dimension signal reconstruction quality in wireless sensor network, the wireless sense network
It is made of m wireless sensor in network, each wireless sensor acquires the signal of d dimensions respectively, and the bandwidth of the signal indicates
For n, collected signal is expressed as y=VHa+nw, wherein y indicates that the measured value of signal, a indicate to need reconstructed signal,
nwIndicate that white noise signal, V indicate that a dimension is ndThe generalized circular matrix of × m.
When using force zero estimation or Minimum Mean Squared Error estimation method reconstruction signal, the mean square error difference of reconstructed results
ForWithWherein,For the variance of white noise signal, β=
nd/ m, Eλ{ } indicates to solve the desired value about λ, and λ is matrix β VVHCharacteristic value.
The probability density function f of λ is solved by the following methodλ(d, β, λ), specially:
1) constitution optimization equation:And meet conditional equation:Wherein, μpFor the p rank squares of λ;
2) by Gaussian quadrature rule the fixed point position f of the integral operation in optimization method and conditional equationλ(d,β,
λ) weighted sum of value shows, i.e.,:
Wherein, λjFor j-th of sampling defined in Gaussian quadrature rule
The coordinate of point, N are sampled point number, wjFor corresponding λjWeights;
3) solving-optimizing equation and conditional equation obtain corresponding to λjFλ(d,β,λj) optimal value, to fit fλ
(d,β,λ)。
According to obtained fλ(d, β, λ) can calculate the mean square error of reconstructed results, to weigh according to mean square error
Amount carries out the quality of signal reconstruction using force zero estimation or Minimum Mean Squared Error estimation method.
As the present invention further technical solution, use force zero method of estimation rebuild signal for:
As the present invention further technical solution, use Minimum Mean Squared Error estimation method rebuild signal for:Wherein, I indicates unit matrix.
As the further technical solution of the present invention, the element of generalized circular matrix V is:Its
In, Vs,tIndicate the s row t column elements of V, xtIndicate the geographical location of t-th of wireless sensor.
As the further technical solution of the present invention, a is a ndDimensional vector.
As the further technical solution of the present invention, y is a m dimensional vector.
The present invention has the following technical effects using above technical scheme is compared with the prior art:When the dimension of acquisition signal
Number, the number of sensor and collected signal bandwidth value it is larger when, technical solution using the present invention can it is easy again
Accurately obtain the probability density function of generalized circular matrix characteristic value.
Description of the drawings
Fig. 1 is the numerical value figure for the system mean square error being calculated in the specific embodiment of the invention.
Specific implementation mode
Below in conjunction with the accompanying drawings and specific embodiment is described in further detail technical scheme of the present invention:
In the wireless sensor network that one is made of m wireless sensor, each sensor is responsible for one d dimension of acquisition
The bandwidth of signal, the signal is expressed as n.Collected signal can be expressed as with the reconstructed signal of needs:Y=VHa+nw,
In, y indicates the measured value of signal, is a m dimensional vector;A indicates to need reconstructed signal, is a ndDimensional vector;nwIt indicates
White noise signal;V indicates that a dimension is ndThe s row t column elements of the generalized circular matrix of × m, this matrix can be expressed as:Wherein, xtIndicate the geographical location of t-th of wireless sensor.
When application force zero estimates (ZF), the signal of reconstruction is expressed as:The mean square error of the reconstructed results
Therefore difference can be written as:
As application Minimum Mean Squared Error estimation (LMMSE), the signal of reconstruction is expressed as:It should
Therefore the mean square error of reconstructed results can be written as:Wherein, I indicates unit matrix.
When the value of m, n, d are larger, it is assumed that ndIt is expressed as β=n with the ratio of md/ m then applies force zero estimation or minimum equal
The mean square error of square estimation error reconstruction signal further can be accurately estimated as:
Wherein, λ is matrix β VVHCharacteristic value.Therefore, probability density function (PDF) f of λλ(d, β, λ) needs a kind of side
Method is provided, wherein d and β is two parameter.
Before the probability density function (PDF) for calculating λ, the p ranks square of λ is it is known that be denoted as μp.The p rank squares of λ exist
A.Nordio,C.-F.Chiasserini,E.Viterbo,“Performance of linear field
reconstruction techniques with noise and uncertain sensor locations,”IEEE
It is discussed in detail in Trans.on Signal Processing, Vol.56, No.8, pp.3535-3547, Aug.2008..
Constitution optimization equation:And meet conditional equation:
By Gaussian quadrature rule the fixed point position f of the integral operation in optimization method and conditional equationλ(d,β,λ)
The weighted sum of value shows, i.e.,:
Wherein, λjFor j-th of sampling defined in Gaussian quadrature rule
The coordinate of point, N are sampled point number, wjFor corresponding λjWeights, also by Gaussian quadrature rule obtain.Therefore, excellent by solving
Change equation and conditional equation, corresponds to each λjThe f of sampled pointλ(d,β,λj) optimal value can be found out, so that fitting obtains fλ
(d,β,λ)。
In Simulation Test Environment, we test when d=1 or 4,And m=ndThe parameter setting of/β, corresponds to
The preceding 12 rank square of eigenvalue λ is known.Believe assuming that system is respectively adopted force zero estimation and is rebuild with Minimum Mean Squared Error estimation method
Number, we are respectively compared the system mean square error obtained by monte carlo method and are calculated with by method proposed by the invention
System mean square error.
The numerical value figure of system mean square error as shown in Figure 1, by observation chart 1 as can be seen that Minimum Mean Squared Error estimation
The performance of method is better than force zero method of estimation, because the target of least mean-square error is minimum mean square error.It is prior
It is that, as β value must increase, the mean square error result obtained with monte carlo method by the calculated mean square error of the present invention is very
It is close.But the invention enables the calculation amounts of mean square error to be substantially reduced, in Fig. 1, the point of " " and "×" indicates Monte Carlo
What the lines of the simulation result that method obtains, dotted line and solid line indicated is the simulation result of the method for the present invention.
The above, the only specific implementation mode in the present invention, but scope of protection of the present invention is not limited thereto, appoints
What is familiar with the people of the technology within the technical scope disclosed by the invention, it will be appreciated that expects transforms or replaces, and should all cover
Within the scope of the present invention, therefore, the scope of protection of the invention shall be subject to the scope of protection specified in the patent claim.
Claims (6)
1. a kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network, by m nothing in the wireless sensor network
Line sensor is constituted, and each wireless sensor acquires the signal of d dimensions respectively, and the bandwidth of the signal is expressed as n, collected
Signal is expressed as y=VHa+nw, wherein y indicates that the measured value of signal, a indicate to need reconstructed signal, nwIndicate white noise
Signal, V indicate that a dimension is ndThe generalized circular matrix of × m, which is characterized in that when using force zero estimation or least mean-square error
When method of estimation reconstruction signal, the mean square error of reconstructed results is respectivelyWithWherein,For the variance of white noise signal, β=nd/ m, Eλ{ } indicates to solve and close
In the desired value of λ, λ is matrix β VVHCharacteristic value;
The probability density function f of λ is solved by the following methodλ(d, β, λ), specially:
1) constitution optimization equation:And meet conditional equation:Wherein, μpFor the p rank squares of λ;
2) by Gaussian quadrature rule the fixed point position f of the integral operation in optimization method and conditional equationλ(d, β, λ) takes
The weighted sum of value shows, i.e.,:
Wherein, λjFor j-th of sampling defined in Gaussian quadrature rule
The coordinate of point, N are sampled point number, wjFor corresponding λjWeights;
3) solving-optimizing equation and conditional equation obtain corresponding to λjFλ(d,β,λj) optimal value, to fit fλ(d,β,
λ);
According to obtained fλ(d, β, λ) can calculate the mean square error of reconstructed results, to weigh use according to mean square error
Force zero is estimated or Minimum Mean Squared Error estimation method carries out the quality of signal reconstruction.
2. the method for higher-dimension signal reconstruction quality, feature in a kind of detection wireless sensor network according to claim 1
Be, use force zero method of estimation rebuild signal for:
3. the method for higher-dimension signal reconstruction quality, feature in a kind of detection wireless sensor network according to claim 1
Be, use Minimum Mean Squared Error estimation method rebuild signal for:Wherein, I indicates unit square
Battle array.
4. the method for higher-dimension signal reconstruction quality, feature in a kind of detection wireless sensor network according to claim 1
It is, the element of generalized circular matrix V is:Wherein, Vs,tIndicate the s row t column elements of V, xt
Indicate the geographical location of t-th of wireless sensor.
5. the method for higher-dimension signal reconstruction quality, feature in a kind of detection wireless sensor network according to claim 1
It is, a is a ndDimensional vector.
6. the method for higher-dimension signal reconstruction quality, feature in a kind of detection wireless sensor network according to claim 1
It is, y is a m dimensional vector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611018905.7A CN106649198B (en) | 2016-11-21 | 2016-11-21 | A kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611018905.7A CN106649198B (en) | 2016-11-21 | 2016-11-21 | A kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106649198A CN106649198A (en) | 2017-05-10 |
CN106649198B true CN106649198B (en) | 2018-11-02 |
Family
ID=58808872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611018905.7A Active CN106649198B (en) | 2016-11-21 | 2016-11-21 | A kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106649198B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102006113A (en) * | 2010-11-25 | 2011-04-06 | 西安电子科技大学 | Parallel MIMO (multiple input multiple output) signal detection method based on zero forcing predetection |
CN103107815A (en) * | 2013-01-09 | 2013-05-15 | 华为技术有限公司 | Method and device for signal reconstruction |
CN103546910A (en) * | 2013-10-15 | 2014-01-29 | 贵州师范大学 | Method for calculating lower bound of transmission capacity of mine laneway wireless sensor network |
CN104684041A (en) * | 2015-02-06 | 2015-06-03 | 中国科学院上海微系统与信息技术研究所 | Real-time wireless sensor network routing method supporting large-scale node application |
-
2016
- 2016-11-21 CN CN201611018905.7A patent/CN106649198B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102006113A (en) * | 2010-11-25 | 2011-04-06 | 西安电子科技大学 | Parallel MIMO (multiple input multiple output) signal detection method based on zero forcing predetection |
CN103107815A (en) * | 2013-01-09 | 2013-05-15 | 华为技术有限公司 | Method and device for signal reconstruction |
CN103546910A (en) * | 2013-10-15 | 2014-01-29 | 贵州师范大学 | Method for calculating lower bound of transmission capacity of mine laneway wireless sensor network |
CN104684041A (en) * | 2015-02-06 | 2015-06-03 | 中国科学院上海微系统与信息技术研究所 | Real-time wireless sensor network routing method supporting large-scale node application |
Non-Patent Citations (4)
Title |
---|
Adaptive Super-Exponentional Algorithms for Blind Deconvolution of MIMO Systems;Kiyotaka Kohno et al;《IEEE. ISCAS》;20041231;全文 * |
Sorted QR decomposition based detection for MU-MIMO LTE uplink;Chen S et al.;《2010 IEEE 71th Vehicular Technology Conference》;20101231;全文 * |
无线传感器水下监测网络稀疏采样和近似重构;邹志强等;《仪 器 仪 表 学 报》;20121231;第 33 卷(第 12 期);全文 * |
无线传感器网络中一种能量优化的运动单点信号重建算法;赵 明等;《湖南大学学报》;20091031;第36卷(第10期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106649198A (en) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stewart et al. | Estimating interchannel observation‐error correlations for IASI radiance data in the Met Office system | |
US10360472B2 (en) | System, method and computer-accessible medium for determining brain microstructure parameters from diffusion magnetic resonance imaging signal's rotational invariants | |
Zeng et al. | A Regional Gap-Filling Method Based on Spatiotemporal Variogram Model of $\hbox {CO} _ {2} $ Columns | |
CN103020478A (en) | Method for checking reality of ocean color remote sensing product | |
CN107196661B (en) | The signal reconfiguring method of vibration signal processing for rotor-support-foundation system | |
Kotwicki et al. | Improving area swept estimates from bottom trawl surveys | |
CN109212631A (en) | Satellite observation data three-dimensional variation assimilation method considering channel correlation | |
CN104298650B (en) | Multi-method fusion based Kalman filtering quantization method | |
Wang et al. | Equation-based InSAR data quadtree downsampling for earthquake slip distribution inversion | |
Moiseev et al. | Multi-core beamformers: derivation, limitations and improvements | |
Hu et al. | Fundamental comparison of time-domain experimental modal analysis methods based on high-and first-order matrix models | |
CN106485074B (en) | A kind of ocean thermohaline quarry sampling method based on adaptive sample rate | |
Jiang et al. | Effect of removing the common mode errors on linear regression analysis of noise amplitudes in position time series of a regional GPS network & a case study of GPS stations in Southern California | |
CN106649198B (en) | A kind of method of higher-dimension signal reconstruction quality in detection wireless sensor network | |
CN103888100B (en) | A kind of non-gaussian linear stochaastic system filtering method based on negentropy | |
Liu et al. | Efficient spatial modeling using the SPDE approach with bivariate splines | |
Cui et al. | PolInSAR coherence region modeling and inversion: The best normal matrix approximation solution | |
CN104407366A (en) | Pseudo-range smooth processing method | |
Refice et al. | On the use of anisotropic covariance models in estimating atmospheric DInSAR contributions | |
KR20190141321A (en) | Robust shrinkage range estimation method and system based on hampel and skipped filters | |
Kandaswamy et al. | Analytic sensing: Noniterative retrieval of point sources from boundary measurements | |
Pribić | Stochastic resolution analysis via a GLR test in radar | |
Kim et al. | Comparison of regression models for spatial downscaling of coarse scale satellite-based precipitation products | |
Rogass et al. | Performance of correlation approaches for the evaluation of spatial distortion reductions | |
Pedder et al. | Applications of EOF analysis to the spatial estimation of circulation features in the ocean sampled by high-resolution CTD soundings |
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 |