最小二乘平差

1、基本观测

观测方程

L = A X +\Delta  (1)

最小二乘

V = AX + l    (2 )   l = -L +\Delta

最优估计

VTPV =MIN    (3)

V中含有X,

(AX-L)TP(AX-L) =MIN

(AX-L)TPV/DX  + VTP(AX-L)/DX = ATPV + (VTPA)T = 2ATPV

蓝色求导用到矩阵求导

Y = AXT(C)  ;DY/DX =AC ;

红色求导用到矩阵求导

Y = CAX     DY/DX =(CA)T ;

令其等于0

ATPV = ATP(AX-L) =0;

X=( ATPA)-1ATPL

2、条件平差

模型

A (L+V) +A0 =0;

A V +W =0;

W =AL+A0

最小二乘

VTPV =MIN

拉格朗日乘数法求极值 ,K为联系数向量。K = - (AP-1AT)-1W

\Phi = VTPV -2KT(AV+W)

对V求一阶导  用到矩阵求导公式

Y = XTA DY/DX = A  ;

Y =AX  DY/DX =AT

\Phi/dV = VT/DVPV + VTPV/DV -2KTAV/DV =PV + (VTP)T - 2(KTA) T = 2PV - 2ATK =0

PV = ATK

V = P-1ATK = P-1 AT (-AP-1AT)-1W = - P-1AT(AP-1AT)-1W 

1、求矩阵A  、W

2、法方程系数AP-1AT

3、求K

4、求V

5、平差V+L

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值