多维粗差定位定值LEGE法粗差探测MATLAB实现

多维粗差定位定值（LEGE）法粗差探测（MATLAB实现）。附带算例，通过输入系数矩阵、权阵。观测值向量来检测观测值中是否包含粗差

clear
n=19             %观测值个数n
t=10              %独立参数个数t
r=9            %粗差搜索计算要进行r-1次

B=double（load（‘B.txt‘））   %系数矩阵   （n*t）
d=double（load（‘d.txt‘））   %间接平差常数项  （n*1）
L=load（‘L.txt‘）   %真值向量      （n*1）
err=load（‘err.txt‘） %模拟误差
L=double（（L-err/1000））
%P=load（‘P.txt‘）;   %权   （n）
S=load（‘S.txt‘）
P=sum（S）/n./S;   %反距离定权
P=diag（P）         %权阵   （n*n）
P=eye（n）

X0=double（pinv（B）*（L-d））;  %因为参数估值X0是L的函数，因此不能手动输入，现利用广义逆矩阵求得
X0=[-178.8288;-7.5124;48.7509;91.1638;116.5219;134.3469;38.6390;-78.7160;230.5690;149.3715];

l=double（L-（B*X0+d））      %V=Bx-l   （n*1）
l=l*1000;
Nbb=double（B.‘*P*B）;       %Nbb   （t*t）
W=double（B‘*P*l）;          %W      （t*1）
x=double（Nbb\W）;      %参数改正数x   （t*1）R（:12）
x0=x
V=double（（B*x-l））;          %初始计算得到的改正数 （n*1） 单位mm
V0=V; 
Qvv=double（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3489  2019-06-03 15:53  LEGE2.m

     文件         74  2019-03-25 20:51  S.txt

     文件       1725  2019-03-25 20:39  B.txt

     文件        171  2019-03-25 20:38  d.txt

     文件        105  2019-05-09 15:08  err.txt

     文件        183  2019-03-26 16:01  L.txt

----------- ---------  ---------- -----  ----

                 5747                    6