If you wish to contribute or participate in the discussions about articles you are invited to contact the Editor

Weighted Least Square Solution (WLS): Difference between revisions

From Navipedia
Jump to navigation Jump to search
No edit summary
No edit summary
 
(16 intermediate revisions by the same user not shown)
Line 17: Line 17:




The measurements quality can be incorporated in the fitting criteria by introducing a symmetric, positive-definite weighting matrix <math>{\mathbf W}</math> and re-defining the "best-fit" as the condition <ref group="footnotes"> Notice that, if <math>{\mathbf {W}}</math> is a diagonal matrix, thence, the norm <math>\displaystyle \|{\mathbf {Y}}-{\mathbf {\hat{Y}}\|^{^{ 2}}_{_{\mathbf {W}}} =\sum{w_i(y_i-\hat{y}_i)^2}</math>, where the terms in the sums are weighted by the diagonal elements <math>w_i</math>. Thence, it follows that the measurements associated to large <math>w_i</math> values will contribute more to the adjustment.</ref>:
The measurements quality can be incorporated in the fitting criteria by introducing a symmetric, positive-definite weighting matrix <math>{\mathbf W}</math> and re-defining the "best-fit" as the condition <ref group="footnotes"> Notice that, if <math>{\mathbf W}</math> is a diagonal matrix, thence, the norm <math> \displaystyle \|{\mathbf {Y}-{\mathbf {\hat{Y}}\|^{^{ 2}}_{_{\mathbf {W}}}}} =\sum{w_i(y_i-\hat{y}_i)^2}</math>, where the terms in the sums are weighted by the diagonal elements <math>w_i</math>. Thence, it follows that the measurements associated to large <math>w_i</math> values will contribute more to the adjustment.</ref>:


:<math>
:<math>
Line 23: Line 23:
</math>
</math>


where the norm of the residuals vector <math>{\mathbf r}={\mathbf Y}-{\mathbf \hat{Y}}</math> to minimise is, now, associated to the scalar product defined by the weighting matrix <math>{\mathbf W}</math>.
where the norm of the residuals vector <math>{\mathbf{r}}={\mathbf{Y}}-{\mathbf {\hat{Y}}}</math> to minimise is, now, associated to the scalar product defined by the weighting matrix <math>{\mathbf W}</math>.


With this weighting criteria, estimator <math>\hat{\mathbf X}_{_{\mathbf W}}</math> and its covariance matrix <math>{\mathbf P}_{_{\mathbf W}}</math> are:
With this weighting criteria, estimator <math>\hat{\mathbf X}_{_{\mathbf W}}</math> and its covariance matrix <math>{\mathbf P}_{_{\mathbf W}}</math> are:

Latest revision as of 16:38, 3 April 2020


FundamentalsFundamentals
Title Weighted Least Square Solution (WLS)
Author(s) J. Sanz Subirana, J.M. Juan Zornoza and M. Hernández-Pajares, Technical University of Catalonia, Spain.
Level Advanced
Year of Publication 2011

In the least squares fit

[math]\displaystyle{ min \ |\mathbf {Y}-\mathbf {\hat{Y}}|^2 = min \left[\sum_{i=1}^n{\left ( y_i-\hat{y}_i \right )^2 }\right] \qquad \mbox{(1)} }[/math]
where [math]\displaystyle{ {\mathbf{\hat{Y}}=\mathbf{G}\,\mathbf{\hat{X}}} }[/math]

all measurements have been treated in the same way (see Parameters adjustment). Nevertheless in a real scenario, not all the measurements have the same quality, i.e., the same error. For instance, measurements at low elevation have larger multipath error than the measurements at high elevation. In the same way, the modelling errors due to the tropospheric or ionospheric mapping functions are larger at lower elevations.


The measurements quality can be incorporated in the fitting criteria by introducing a symmetric, positive-definite weighting matrix [math]\displaystyle{ {\mathbf W} }[/math] and re-defining the "best-fit" as the condition [footnotes 1]:

[math]\displaystyle{ min\|{\mathbf {Y}}-{\mathbf{\hat{Y}}}\|^{^{ 2}}_{_{\mathbf {W}}} \qquad \mbox{(2)} }[/math]

where the norm of the residuals vector [math]\displaystyle{ {\mathbf{r}}={\mathbf{Y}}-{\mathbf {\hat{Y}}} }[/math] to minimise is, now, associated to the scalar product defined by the weighting matrix [math]\displaystyle{ {\mathbf W} }[/math].

With this weighting criteria, estimator [math]\displaystyle{ \hat{\mathbf X}_{_{\mathbf W}} }[/math] and its covariance matrix [math]\displaystyle{ {\mathbf P}_{_{\mathbf W}} }[/math] are:


[math]\displaystyle{ \hat{\mathbf X}_{_{\mathbf W}}=({\mathbf G}^T\,{\mathbf W}\,{\mathbf G})^{-1}{\mathbf G}^T\,{\mathbf W}\,{\mathbf Y} \qquad \mbox{(3)} }[/math]


[math]\displaystyle{ {\mathbf P}_{_{\mathbf \Delta X_W}}=({\mathbf G}^T\,{\mathbf W}\,{\mathbf G})^{-1}{\mathbf G}^T \,{\mathbf W}\,\,{\mathbf R}\,\,{\mathbf W}\,{\mathbf G}({\mathbf G}^T\,{\mathbf W}\,{\mathbf G})^{-1} \qquad \mbox{(4)} }[/math]


Notes

  1. ^ Notice that, if [math]\displaystyle{ {\mathbf W} }[/math] is a diagonal matrix, thence, the norm [math]\displaystyle{ \displaystyle \|{\mathbf {Y}-{\mathbf {\hat{Y}}\|^{^{ 2}}_{_{\mathbf {W}}}}} =\sum{w_i(y_i-\hat{y}_i)^2} }[/math], where the terms in the sums are weighted by the diagonal elements [math]\displaystyle{ w_i }[/math]. Thence, it follows that the measurements associated to large [math]\displaystyle{ w_i }[/math] values will contribute more to the adjustment.