If you wish to contribute or participate in the discussions about articles you are invited to contact the Editor
Weighted Least Square Solution (WLS)
| 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{ \displaystyle \|{\mathbf {Y}-{\mathbf {\hat{Y}}\|^ }[/math]
- [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.