Nothing Special   »   [go: up one dir, main page]

CN108228535A - A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion - Google Patents

A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion Download PDF

Info

Publication number
CN108228535A
CN108228535A CN201810001201.1A CN201810001201A CN108228535A CN 108228535 A CN108228535 A CN 108228535A CN 201810001201 A CN201810001201 A CN 201810001201A CN 108228535 A CN108228535 A CN 108228535A
Authority
CN
China
Prior art keywords
data
fusion
parameter
observation
estimated
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.)
Pending
Application number
CN201810001201.1A
Other languages
Chinese (zh)
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.)
Foshan Yi Jia Technology Co Ltd
Foshan University
Original Assignee
Foshan Yi Jia Technology Co Ltd
Foshan University
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 Foshan Yi Jia Technology Co Ltd, Foshan University filed Critical Foshan Yi Jia Technology Co Ltd
Priority to CN201810001201.1A priority Critical patent/CN108228535A/en
Publication of CN108228535A publication Critical patent/CN108228535A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Complex Calculations (AREA)

Abstract

The invention discloses a kind of optimal weighting parameter evaluation methods of unequal precision measurement data fusion, including obtaining first kind observation data and the second kind observation data, determine parameter to be estimated;Obtain first kind observation function;Obtain the second kind observation function;Structure first kind observation data merge majorized function with the second kind observation data;The fusion majorized function is solved, calculates optimal weighting value and to be estimated parameter of the first kind observation data with the second kind observation data.The calculating that the present invention passes through estimated bias and mean square error, the computational methods of optimum fusion weights when giving the processing of multiclass measurement data fusion, corresponding parameter estimation algorithm is established simultaneously, the Parameter fusion estimated value that can be further accurately calculated when realizing Data Fusion under unequal precision measurement data optimal weighting.The invention is used to merge the measurement data of unequal accuracy.

Description

A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion
Technical field
The present invention relates to Data fusion technique field, more specifically to a kind of unequal precision measurement data fusion most Excellent weighting parameters evaluation method.
Background technology
In measurement data fusion treatment, the most typically fusion of different type, unequal accuracy data.When number will be observed After parameter model is expressed as, measurement data fusion problem can be converted to the Parameter Estimation Problem of regression model.Measurement data Fusion treatment be conducive to improve parameter Estimation precision.Wherein, different types of data or heterogeneous data refer to that observing data closes It is different in the functional relation of parameter to be estimated, so as to which its all-order derivative is also different.If to a variety of essences such as not involved in fusion treatment The observation data of degree, the different weighting scheme of these data, have a great impact to parameter estimation result, therefore for the essences such as not The weighting of degrees of data fusion is treated as improving the key technology of Parameter Estimation Precision.
For the parameter Estimation of linear regression model (LRM), Gauss-Markov theorems (Gauss-Markov theorem) give Unique optimal weighting principle of unequal precision measurement data, i.e. blending weight are only related with the precision of measurement data.
And for there are during the processing of the multiclass measurement data fusion of nonlinear model, optimum fusion weights are no longer by measurement number According to precision uniquely determine that but related with the structure of observation model simultaneously, i.e. the Gauss-Markov of linear model classics determines Reason is no longer set up.
But at present about the measurement data fusion there are nonlinear model during, always assume that it is all observation data for etc. Precision, that is, the random error for observing data is independent same distribution, so as to not consider to weight or directly use linear model The conclusion of Gauss-Markov theorems is weighted the multiclass measurement data for including linear and nonlinear model, so as to reduce The precision of data fusion.During for multiclass measurement data fusion parameter Estimation including nonlinear model, due to being all false If all observation data for equally accurate or directly merged using the precision for observing data to be weighted, therefore, for depositing In the multiclass unequal accuracy measurement data fusion of non-linear observation model, there are certain mistakes for Parameter fusion estimated result value Difference.
Invention content
The technical problem to be solved by the present invention is to:A kind of optimal weighting parameter of unequal precision measurement data fusion is provided to estimate Calculation method.
The present invention solve its technical problem solution be:
A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion, including:
First kind observation data and the second kind observation data are obtained, determine parameter to be estimated;
It determines the nonlinear function between the first kind observation data and parameter to be estimated, obtains first kind observation letter Number;
It determines the linear functional relation between the second kind observation data and parameter to be estimated, obtains the second kind observation letter Number;
Based on the first kind observation function and the second kind observation function, structure first kind observation data are seen with the second class The fusion majorized function of measured data;
Based on root mean squared error for channel estimation minimum criteria, the fusion majorized function is solved, and according to solving result, calculate Optimal weighting value and to be estimated parameter of the first kind observation data with the second kind observation data.
As a further improvement of the above technical scheme, the first kind observation function is as shown in expression formula 1, y1(ti)=f (ti,β)+ε1(ti), i=1 ..., m, wherein y1(ti) it is tiThe first kind observation data at moment, f (ti, β) and it is corresponding non- Systems with Linear Observation function, β are unitary parameter to be estimated, and β ∈ R, R are real number field, ε1(ti) (i=1 ..., m) it is that the first kind observes data Measurement random noise, independent same distribution is zero in mean value, and variance isNormal distribution,Data are observed for the first kind Measurement accuracy, m are the observation sample number that the first kind observes data.
As a further improvement of the above technical scheme, the second kind observation function is as shown in expression formula 2, y2(ti)=x (ti)β+ε2(ti), i=1 ..., k, wherein, y2(ti) it is tiThe second kind observation data at moment, x (ti) observed to be corresponding Matrix, ε2(ti) (i=1 ..., k) the second kind observation data measurement random noise, independent same distribution is zero in mean value, and variance isNormal distribution,For the measurement accuracy of the second kind observation data, k is the sample number of the second kind observation data.
As a further improvement of the above technical scheme, it is described fusion majorized function as shown in expression formula 3,WhereinData and institute are observed for the first kind The blending weight of the second kind observation data is stated,
As a further improvement of the above technical scheme, the fusion majorized function is solved, and according to solving result, is calculated Optimal weighting value and to be estimated parameter of the first kind observation data with the second kind observation data, this process specifically includes following Step:
Step A. setting weighting initial valuesMinimum is solved by the expression formula 3, is solved
The estimation mean square error that step B. calculates parameter beta to be estimated by expression formula 4 existsThe value at place
Wherein
Step C. solves minimumIt is right by expression formula 4Derivation simultaneously makes expression formula 4 Equal to 0, obtain
Step D. sets convergence criterion, described in judgementWhether set convergence criterion is met, if satisfied, then iteration Terminate, determineFor optimum fusion weights,Optimal estimation for parameter.
As a further improvement of the above technical scheme, convergence criterion described in step D is as shown in expression formula 5,Wherein τ is convergence threshold, takes τ=0.01.
As a further improvement of the above technical scheme, in step D, if describedSet convergence criterion is unsatisfactory for, It willIt is assigned toAnd return to step A, until iteration convergence, wherein,It is seen for the first kind The blending weight of measured data and the second kind observation data.
The beneficial effects of the invention are as follows:The present invention gives multiclass observation by the calculating of estimated bias and mean square error The computational methods of optimum fusion weights during Data Fusion, while corresponding parameter estimation algorithm is established, realize data The Parameter fusion estimated value that can be further accurately calculated during fusion treatment under unequal precision measurement data optimal weighting. The invention is used to merge the measurement data of unequal accuracy.
Description of the drawings
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment Attached drawing is briefly described.Obviously, described attached drawing is the part of the embodiment of the present invention rather than all implements Example, those skilled in the art without creative efforts, can also be obtained according to these attached drawings other designs Scheme and attached drawing.
Fig. 1 is the evaluation method flow chart of the present invention.
Specific embodiment
The technique effect of the design of the present invention, concrete structure and generation is carried out below with reference to embodiment and attached drawing clear Chu, complete description, to be completely understood by the purpose of the present invention, feature and effect.Obviously, described embodiment is this hair Bright part of the embodiment rather than whole embodiments, based on the embodiment of the present invention, those skilled in the art is not paying The other embodiment obtained under the premise of creative work, belongs to the scope of protection of the invention.
With reference to Fig. 1, the invention discloses a kind of optimal weighting parameter estimation side of unequal precision measurement data fusion Method, including:
First kind observation data and the second kind observation data are obtained, determine parameter to be estimated;
It determines the nonlinear function between the first kind observation data and parameter to be estimated, obtains first kind observation letter Number;
It determines the linear functional relation between the second kind observation data and parameter to be estimated, obtains the second kind observation letter Number;
Based on the first kind observation function and the second kind observation function, structure first kind observation data are seen with the second class The fusion majorized function of measured data;
Based on root mean squared error for channel estimation minimum criteria, the fusion majorized function is solved, and according to solving result, calculate Optimal weighting value and to be estimated parameter of the first kind observation data with the second kind observation data.
Specifically, the present invention gives the processing of multiclass measurement data fusion by the calculating of estimated bias and mean square error When optimum fusion weights computational methods, while establish corresponding parameter estimation algorithm, energy when realizing Data Fusion Enough Parameter fusion estimated values being further accurately calculated under unequal precision measurement data optimal weighting.
It is further used as preferred embodiment, in the invention specific embodiment, the first kind observation function As shown in expression formula 1, y1(ti)=f (ti,β)+ε1(ti), i=1 ..., m, wherein y1(ti) it is tiThe first kind observation at moment Data, f (ti, β) and for corresponding non-linear observation function, β is unitary parameter to be estimated, and β ∈ R, R are real number field, ε1(ti) (i= 1 ..., m) it is the measurement random noise that the first kind observes data, independent same distribution is zero in mean value, and variance isNormal state point Cloth,The measurement accuracy of data is observed for the first kind, m is the observation sample number that the first kind observes data.
It is further used as preferred embodiment, in the invention specific embodiment, the second kind observation function As shown in expression formula 2, y2(ti)=x (ti)β+ε2(ti), i=1 ..., k, wherein, y2(ti) it is tiThe second kind observation at moment Data, x (ti) for corresponding observing matrix, ε2(ti) (i=1 ..., k) the second kind observation data measurement random noise, it is independent It is zero with mean value is distributed in, variance isNormal distribution,For the measurement accuracy of the second kind observation data, k is the second class Observe the sample number of data.
It is further used as preferred embodiment, in the invention specific embodiment, the fusion majorized function is such as Shown in expression formula 3,WhereinFor the first kind The blending weight of data and the second kind observation data is observed,
It is further used as preferred embodiment, in the invention specific embodiment, solves the fusion optimization letter Number, and according to solving result, calculate the optimal weighting value of the first kind observation data and the second kind observation data and wait to estimate ginseng Number, this process specifically include the following steps:
Step A. setting weighting initial valuesMinimum is solved by the expression formula 3, is solved
The estimation mean square error that step B. calculates parameter beta to be estimated by expression formula 4 existsThe value at place
Wherein
Step C. solves minimumIt is right by expression formula 4Derivation simultaneously makes expression formula 4 Equal to 0, obtain
Step D. sets convergence criterion, described in judgementWhether set convergence criterion is met, if satisfied, then iteration Terminate, determineFor optimum fusion weights,Optimal estimation for parameter.
It is further used as preferred embodiment, in the invention specific embodiment, convergence criterion described in step D As shown in expression formula 5,Wherein τ is convergence threshold, takes τ=0.01.
Preferred embodiment is further used as, in the invention specific embodiment, in step D, if describedNo Meet set convergence criterion, it willIt is assigned toAnd return to step A, until iteration convergence, wherein,The blending weight of data and the second kind observation data is observed for the first kind.
A specific embodiment is lifted herein the invention is described in detail.
The present embodiment considers the measurement and estimation of parameter beta to be estimated to unitary, is illustrated using two classes observation data, a kind of It is non-Systems with Linear Observation function, another kind of is Systems with Linear Observation function.
First kind observation function is y1(ti)=f1(t1,β)+ε1(ti), wherein f1(ti, β) and=1+ (5+tiβ)0.1For about The nonlinear function of parameter beta to be estimated, ε1(ti) be mean value it is zero, variance 0.052Measurement random noise, wherein observation sample Number is:ti=0.05 × (i-1), i=1 ..., 300, i.e. observation sample m=300.
The second kind observation function is y2(tk)=f2(tk,β)+ε2(tk), whereinFor about waiting to estimate ginseng The linear function of number β, ε2(tk) be mean value it is zero, variance 0.012Measurement random noise, wherein observation sample number is:tk= 2+0.1 × (k-1), k=1 ..., 100, i.e. observation sample k=100.
If the true value of β is 10, observation data are generatedUtilize 100 systems of Monte Carlo simulation The estimated result of parameter beta to be estimated is counted, the estimated result that table 1 is given under various different blending weights compares.
When the weighting scheme using linear model is weighted fusion to above two class observation data, that is, the blending weight is taken to beWhen, the estimated value of obtained parameter beta is 10.064, and the mean square error of estimation is 0.054; And use the optimum fusion weights of iteration convergenceWhen, the estimated value of obtained parameter beta is 9.987, is estimated at this time Mean square error is minimum, is 0.015.
The better embodiment of the present invention is illustrated, but the invention is not limited to the implementation above Example, those skilled in the art can also make various equivalent modifications under the premise of without prejudice to spirit of the invention or replace It changes, these equivalent modifications or replacement are all contained in the application claim limited range.

Claims (7)

1. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion, which is characterized in that including:
First kind observation data and the second kind observation data are obtained, determine parameter to be estimated;
It determines the nonlinear function between the first kind observation data and parameter to be estimated, obtains first kind observation function;
It determines the linear functional relation between the second kind observation data and parameter to be estimated, obtains the second kind observation function;
Based on the first kind observation function and the second kind observation function, structure first kind observation data and the second kind observation number According to fusion majorized function;
Based on root mean squared error for channel estimation minimum criteria, the fusion majorized function is solved, and according to solving result, described in calculating The first kind observes the optimal weighting value of data and the second kind observation data and parameter to be estimated.
2. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion according to claim 1, special Sign is that the first kind observation function is as shown in expression formula 1, y1(ti)=f (ti,β)+ε1(ti), i=1 ..., m, wherein y1 (ti) it is tiThe first kind observation data at moment, f (ti, β) and for corresponding non-linear observation function, β is unitary parameter to be estimated, β ∈ R, R are real number field, ε1(ti) the measurement random noise of (i=1 ..., m) for first kind observation data, independent same distribution is in mean value It is zero, variance isNormal distribution,The measurement accuracy of data is observed for the first kind, m is the sight that the first kind observes data Survey sample number.
3. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion according to claim 2, special Sign is that the second kind observation function is as shown in expression formula 2, y2(ti)=x (ti)β+ε2(ti), i=1 ..., k, wherein, y2 (ti) it is tiThe second kind observation data at moment, x (ti) for corresponding observing matrix, ε2(ti) (i=1 ..., k) sight of the second class The measurement random noise of measured data, independent same distribution are zero in mean value, and variance isNormal distribution,For the second kind observation The measurement accuracy of data, k are the sample number of the second kind observation data.
4. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion according to claim 3, special Sign is, the fusion majorized function as shown in expression formula 3, WhereinThe blending weight of data and the second kind observation data is observed for the first kind,
5. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion according to claim 4, special Sign is, solves the fusion majorized function, and according to solving result, calculate the first kind observation data and the second kind observation The optimal weighting value and parameter to be estimated of data, this process specifically include the following steps:
Step A. setting weighting initial valuesMinimum is solved by the expression formula 3, is solved
The estimation mean square error that step B. calculates parameter beta to be estimated by expression formula 4 existsThe value at place Wherein
Step C. solves minimumIt is right by expression formula 4Derivation is simultaneously equal to expression formula 4 0, it obtains
Step D. sets convergence criterion, described in judgementWhether set convergence criterion is met, if satisfied, then iteration terminates, It determinesFor optimum fusion weights,Optimal estimation for parameter.
6. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion according to claim 5, special Sign is, convergence criterion described in step D as shown in expression formula 5,Wherein τ be convergence threshold, take τ= 0.01。
7. a kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion according to claim 6, special Sign is, in step D, if describedSet convergence criterion is unsatisfactory for, it willIt is assigned toAnd return to step A, directly Until iteration convergence, wherein,Melting for data and the second kind observation data is observed for the first kind Close weights.
CN201810001201.1A 2018-01-02 2018-01-02 A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion Pending CN108228535A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810001201.1A CN108228535A (en) 2018-01-02 2018-01-02 A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810001201.1A CN108228535A (en) 2018-01-02 2018-01-02 A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion

Publications (1)

Publication Number Publication Date
CN108228535A true CN108228535A (en) 2018-06-29

Family

ID=62642468

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810001201.1A Pending CN108228535A (en) 2018-01-02 2018-01-02 A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion

Country Status (1)

Country Link
CN (1) CN108228535A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110987068A (en) * 2019-11-28 2020-04-10 中国人民解放军陆军炮兵防空兵学院郑州校区 Data fusion method for multi-sensor integrated control system
CN111736626A (en) * 2020-06-22 2020-10-02 中国人民解放军国防科技大学 Stable missile path data processing method
CN113779502A (en) * 2021-08-20 2021-12-10 绥化学院 Image processing evidence function estimation method based on correlation vector machine

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412220A (en) * 2013-08-13 2013-11-27 无锡科技职业学院 Ship power station fault diagnosis method based on data fusion
US9572533B2 (en) * 2012-06-22 2017-02-21 Fitbit, Inc. GPS power conservation using environmental data

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9572533B2 (en) * 2012-06-22 2017-02-21 Fitbit, Inc. GPS power conservation using environmental data
CN103412220A (en) * 2013-08-13 2013-11-27 无锡科技职业学院 Ship power station fault diagnosis method based on data fusion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WANG JIONGQI 等: "Optimal Weight and Parameter Estimation of Multi-structure and Unequal-Precision Data Fusion", 《CHINESE JOURNAL OF ELECTRONICS》 *
周海银等: "基于多源信息融合的飞行器高精度定位方法", 《湖南大学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110987068A (en) * 2019-11-28 2020-04-10 中国人民解放军陆军炮兵防空兵学院郑州校区 Data fusion method for multi-sensor integrated control system
CN111736626A (en) * 2020-06-22 2020-10-02 中国人民解放军国防科技大学 Stable missile path data processing method
CN113779502A (en) * 2021-08-20 2021-12-10 绥化学院 Image processing evidence function estimation method based on correlation vector machine
CN113779502B (en) * 2021-08-20 2023-08-29 绥化学院 Image processing evidence function estimation method based on correlation vector machine

Similar Documents

Publication Publication Date Title
CN108228535A (en) A kind of optimal weighting parameter evaluation method of unequal precision measurement data fusion
Croon Using predicted latent scores in general latent structure models
Cheng et al. A modified Cobb–Douglas production function model and its application
Hanna et al. Data-efficient policy evaluation through behavior policy search
CN104102836B (en) A kind of quick robust state estimation method of electric system
CN106296434B (en) Grain yield prediction method based on PSO-LSSVM algorithm
CN111308896A (en) Nonlinear system self-adaptive optimal control method based on variable error
CN108960334A (en) A kind of multi-sensor data Weighted Fusion method
Xiao et al. Incentive mechanism design for federated learning: A two-stage stackelberg game approach
Tehranchi Uniform bounds for Black--Scholes implied volatility
CN109446614A (en) A kind of discrimination method and device of static reactive model parameter
CN104965506B (en) One kind is based on adjustable parameter Distributed Flight Control System real-time fault diagnosis method
CN106289727A (en) A kind of element damage from laser measuring method and device
CN105654053B (en) Based on the dynamic oscillation signal parameter discrimination method for improving constraint EKF algorithm
CN106599541B (en) A kind of structure and parameter on-line identification method of dynamic power load model
CN116956744A (en) Multi-loop groove cable steady-state temperature rise prediction method based on improved particle swarm optimization
CN109582915B (en) Improved nonlinear observability self-adaptive filtering method applied to pure azimuth tracking
CN110942346A (en) Dynamic estimation method for construction engineering cost based on particle swarm algorithm
CN102866629A (en) Dyanmic-static mixed nerve network modeling-based anti-interference control method for random system
CN116437290A (en) Model fusion method based on CSI fingerprint positioning
Wei et al. A Pearson-Type VII Distribution With Adaptive Parameters Selection-Based Interacting Multiple Model Kalman Filter
CN115685128A (en) Radar target tracking algorithm and electronic equipment under maneuvering target scene
CN109493065A (en) A kind of fraudulent trading detection method of Behavior-based control incremental update
CN108319571A (en) A kind of multi-Dimensional parameters evaluation method of multiclass isomery model optimum fusion
CN108009565B (en) Method for calculating design flood level of river network area under variable environment

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
CB02 Change of applicant information

Country or region after: China

Address after: 528000 Foshan Institute of science and technology, Xianxi reservoir West Road, Shishan town, Nanhai District, Foshan City, Guangdong Province

Applicant after: Foshan University

Applicant after: FOSHAN YOUYIJIA TECHNOLOGY Co.,Ltd.

Address before: 528000 Foshan Institute of science and technology, Xianxi reservoir West Road, Shishan town, Nanhai District, Foshan City, Guangdong Province

Applicant before: FOSHAN University

Country or region before: China

Applicant before: FOSHAN YOUYIJIA TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information