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
Gema.Cueto (talk | contribs) No edit summary |
Gema.Cueto (talk | contribs) No edit summary |
||
| (21 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 | 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> | ||
min | min\|{\mathbf {Y}}-{\mathbf{\hat{Y}}}\|^{^{ 2}}_{_{\mathbf {W}}} \qquad \mbox{(2)} | ||
</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
| Fundamentals | |
|---|---|
| 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
- ^ 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.