抗差估计,有偏估计

1、 但简单易行。很小时此法很不可靠，故观测数常未知，和中识别粗差时，实际问题用鲁棒估计估计抗差估计稳健估计nRobusti3精品资料 分布的观测数据。用于定义完全服从正态经典平差模型）高斯马尔柯夫模型（120,:PDLDEBXLEBXL 在观测数据中出现在观测数据中出现0.2%的粗差时，最小二乘估的粗差时，最小二乘估值便失去了其最优性，但值便失去了其最优性，但0.2%的粗差概率完全正的粗差概率完全正常，特别是在现代的大数据量自动测量中。所以常，特别是在现代的大数据量自动测量中。所以(suy)经典平差适用的范围狭窄。经典平差适用的范围狭窄。精品资料抗差估计指导思想：在抗差能力和效率（指估值最优性）

2、中求抗差估计指导思想：在抗差能力和效率（指估值最优性）中求得最佳平衡。一般要求其效率达到经典平差效率的得最佳平衡。一般要求其效率达到经典平差效率的90%90%以上。以上。是在抗差的前提下谈效率。是在抗差的前提下谈效率。抗差估计实质：牺牲最小二乘估计的最优性，达到抵抗粗差污抗差估计实质：牺牲最小二乘估计的最优性，达到抵抗粗差污染的目的。染的目的。抗差估计的特点：当观测数据的实际分布偏离假定模型时的不抗差估计的特点：当观测数据的实际分布偏离假定模型时的不敏感性。其对子样分布要求不十分敏感性。其对子样分布要求不十分(shfn)(shfn)严格，只要子样严格，只要子样近似服从某一模型。近似服从某一模型

3、。若母体确实为正态时，抗差估计值无最小二乘估计值优良。若母体确实为正态时，抗差估计值无最小二乘估计值优良。精品资料最小二乘估计的优点：能够抵御大量随机小误差对最小二乘估计的优点：能够抵御大量随机小误差对参数估值的影响；估值无偏，方差最小。参数估值的影响；估值无偏，方差最小。估值的效率问题：可削弱大量小误差对参数估值的估值的效率问题：可削弱大量小误差对参数估值的影响。影响。抗差能力的标志：估值能容忍的粗差个数。抗差能力的标志：估值能容忍的粗差个数。抗差估计的适用范围：在确定性模型中有大量正确抗差估计的适用范围：在确定性模型中有大量正确的观测值存在，仅有少数几个是不正确的；统计的观测值存在，仅有少

4、数几个是不正确的；统计模型就不一定了，如果轨迹模型是你自己定的，模型就不一定了，如果轨迹模型是你自己定的，出界的点被认为是粗差而剔除，这是不正确的。出界的点被认为是粗差而剔除，这是不正确的。所以抗差估计适合所以抗差估计适合(shh)(shh)确定性模型而不适合确定性模型而不适合(shh)(shh)拟合模型。拟合模型。精品资料数据偏离正态分布的原因：数据偏离正态分布的原因：（1 1）有粗差（观测、记录、数据输入）有粗差（观测、记录、数据输入(shr)(shr)等）等）（2 2）数据组合与舍入误差）数据组合与舍入误差（3 3）就算数据中无粗差存在，但其分布仍有微小）就算数据中无粗差存在，但其分布仍

5、有微小明显的偏离正态趋势明显的偏离正态趋势（4 4）观测值之间并非完全独立）观测值之间并非完全独立 F F、t t、u u、2 2 四种检验方法由于都取决于正态分四种检验方法由于都取决于正态分布的母体，故对于偏离正态分布的数据检验是不布的母体，故对于偏离正态分布的数据检验是不可靠的。可靠的。举例：举例：0625. 5997. 4500. 5008. 5001. 5993. 4998. 4002. 5001. 5x术平均值为：其最小二乘估值，即算，观测值为精品资料高斯断言：如果最小二乘估计量不是最优估值的话，那么，高斯断言：如果最小二乘估计量不是最优估值的话，那么，观测列中必存在一种来自外界的、

6、未知的干扰因素所致。观测列中必存在一种来自外界的、未知的干扰因素所致。现已知，这种未知的干扰因素就是现已知，这种未知的干扰因素就是(jish)(jish)粗差。粗差。稳健估计的目标：稳健估计的目标：1 1、在采用的假定模型下，所估计的参数应具有最优或接、在采用的假定模型下，所估计的参数应具有最优或接近最优；近最优；2 2、如果实际模型与假定模型存在较小的偏差，则对应的、如果实际模型与假定模型存在较小的偏差，则对应的估计参数所受影响也较小；估计参数所受影响也较小；3 3、即使实际模型与假定模型有较大偏差，其参数估值的、即使实际模型与假定模型有较大偏差，其参数估值的性能也不应太差，不至于对估值产生

7、灾难性后果。性能也不应太差，不至于对估值产生灾难性后果。精品资料模型误差的产生和分类模型误差的产生和分类模型误差：模型与客观模型误差：模型与客观(kgun)(kgun)实际的误差，也分为粗实际的误差，也分为粗差、系统误差和偶然误差。差、系统误差和偶然误差。有粗差时用经典平差模型或无粗差时用抗差模型，都会有粗差时用经典平差模型或无粗差时用抗差模型，都会产生模型误差。产生模型误差。比例。的的数据在整个数据组中受污染（有粗差）部分污染率，干扰分布，主体分布，污染分布：xxFFF 1精品资料一、影响函数一、影响函数(hnsh)(hnsh) 影响函数影响函数(hnsh)(hnsh)是用来判断估计量对异常

8、值是用来判断估计量对异常值敏感程度的指标，即一个附加的观测值对估值的影响的敏感程度的指标，即一个附加的观测值对估值的影响的大小。大小。影响函数影响函数(hnsh)(hnsh)定义式：定义式： 在其它处处在：分布，其概率密度函数的（粗差）质量为粗差分布。在污染率。无粗差的主体分布。受污染的分布。其中：0111,00 xxxfxFFFFlinFFlinFlIFxx精品资料 污染率。观测值个数。粗差个数，估值。个粗差后的数据算得的用剔除个观测值算得的估值用全部的实用中，影响函数为：nsnssFnFnsFFFlIF,影响函数的重要用途：影响函数的重要用途：（1）若影响函数无界，则一个粗差可彻底破坏估计

9、量，此时该种方法）若影响函数无界，则一个粗差可彻底破坏估计量，此时该种方法(fngf)就不具有抗差性；就不具有抗差性；（2）大于某个限差的粗差应对平差结果不产生影响，即应设一个影响）大于某个限差的粗差应对平差结果不产生影响，即应设一个影响函数函数IF=0的误差界；的误差界；（3）影响函数可用图形表示，直观，重要。）影响函数可用图形表示，直观，重要。精品资料广义极大似然估计广义极大似然估计M M估计估计M M估计经典估计经典(jngdin)(jngdin)极大似然估计的推广，极大似然估计的推广，最接近传统的最小二乘估计。最接近传统的最小二乘估计。概率密度最大，即的组独立小误差联合出现以极大似然法

10、思想，一参数估值，误差的密度函数为，设其对应的，设有一组独立观测值XLXfLLLin21精品资料 211121minlnminlnmaxlnmaxiiiiiniiniiniinvPvvvvfLXfLXfLXfLXfLXfG最小二乘法中，的函数为残差数或估值函数，其称为目标函数、极值函，使其定义式广义化：上式中的差函数）代替（增长较慢的极小化残选用函数或：或：精品资料 0minmin1111niiniiiniiniivvvvvM的定义也可引申到其导数而由一类估计。式定义的估计由所以，也可以定义 iv iv由于由于(yuy) 或或 选择的不同，选择的不同，会得到不同的会得到不同的M估计法，其稳健性

11、也不同。估计法，其稳健性也不同。精品资料基本思想：基本思想：（1 1）平差仍采用经典）平差仍采用经典(jngdin)(jngdin)的最小二乘平差形的最小二乘平差形式；式； （2 2）每次平差后根据残差和有关参数构成下一步的）每次平差后根据残差和有关参数构成下一步的权函数；权函数；（3 3）迭代中止时相应的残差将直接指出粗差所在的）迭代中止时相应的残差将直接指出粗差所在的位置。平差后有：位置。平差后有： 保权区正常观测值保权区正常观测值 降权区非正常但可用的观测值降权区非正常但可用的观测值 除权区含粗差的观测值除权区含粗差的观测值精品资料 计估计就成为最小二乘估是常数阵，若权阵估计准则：其中每

12、次平差的权阵为阵形式：选权迭代法法方程的矩个观测值的残差次平差后第第个观测值的权函数第式中下标表示：令权函数为：MVPVVPVvPvPvPdiagVPlVPBxBVPBjiijjvvvvvPTiniiTTiiiiijmin021精品资料 代入。的初值用单位阵当观测值等权时，值符合限差要求为止。迭代到前后两次解的差选权迭代法的迭代公式IVPlxBVlVPBBVPBxvvvPkkkTkTkkikiikj111 v随着估计随着估计(gj)函数函数 选取的不同，构成了选取的不同，构成了不同的权函数形式，形成了不同的选权迭代法。不同的权函数形式，形成了不同的选权迭代法。精品资料权函数权函数 选取的要求：

13、选取的要求：（1 1）平差后粗差观测值的权应趋近于）平差后粗差观测值的权应趋近于0 0，其余，其余多余观测值的权趋近于多余观测值的权趋近于1 1；（2 2）迭代中止时，不含粗差的观测值的权应等）迭代中止时，不含粗差的观测值的权应等于于(dngy)(dngy)验前给定或验后方差估求的，平差验前给定或验后方差估求的，平差应回到通常的最小二乘法平差；应回到通常的最小二乘法平差；（3 3）权函数的选择应保证迭代尽快收敛。）权函数的选择应保证迭代尽快收敛。 VP 由于粗差的分布不同，不能象偶然误差一样有一个统一(tngy)的正态分布，有统一(tngy)的处理方法。所以不同的粗差分布对应了不同的处理方法。

14、精品资料1、Huber法2、一次范数最小法3、p范最小法4、丹麦(dn mi)法5、Hampel法精品资料1、Huber法2、一次范数最小法（L1估计(gj)）（中位数法） 212124222442222vvvvpvvsignvvvvvvvvv， kvvpvv1精品资料3、P范最小法（LP法）4、IGG法（周江文法）5、经典最小二乘法(chngf)（不具有抗差性） kvvpvvpp21 5 . 205 . 25 . 115 . 115 . 25 . 25 . 15 . 122vvvvvpvdvvvvv cvppvvpvv22精品资料1546372ABH3H1H2 如图，为模拟水准网，7个观测(

15、gunc)值配赋了随机误差，在第六条路线的观测(gunc)高差中附加了10mm的粗差。用各种算法结果列于下表。从表中看各种选权迭代法均有抗差性，而最小二乘法不具有抗差性。精品资料选权迭代法的缺陷：选权迭代法的缺陷：1 1、由于粗差的大小及位置未知，只能以残、由于粗差的大小及位置未知，只能以残差来研究，且目标函数差来研究，且目标函数选择选择(xunz)(xunz)成为成为残差残差v v的函数，这并不一定符合实际。的函数，这并不一定符合实际。2 2、选权迭代法中，第一次按最小二乘平差、选权迭代法中，第一次按最小二乘平差求得的残差受粗差的影响很大，由此将影响求得的残差受粗差的影响很大，由此将影响迭代

16、的权函数迭代的权函数P(v)P(v)的选择的选择(xunz)(xunz)，可能导，可能导致错误的收敛。致错误的收敛。3 3、为避免选权迭代法的初值受、为避免选权迭代法的初值受LSLS平差法的平差法的影响，可以采用线性规划中的单纯形法进行影响，可以采用线性规划中的单纯形法进行初值确定。初值确定。精品资料 离原点最远的点上。内并且使最优解是在可行解区域约束条件：目标函数：解：间获得最大利润？何安排生产可在计划期可供资源如表，问：如元，元，乙产品得利润个甲产品得利润工，用料均为钢材，每和二车间的机械上加产品，均需要在一车间例：某厂生产甲乙两种约束函数：目标函数：数学模型：21212212121403

17、00,2005 . 07505 . 1600max403040300,minxxzxxxxxxxxxzXbAXXCxfT精品资料图解法的重要结论：图解法的重要结论：可行解的区域为凸多边形，其最优解若存在，一定可行解的区域为凸多边形，其最优解若存在，一定在某个极点（顶点）上。在某个极点（顶点）上。最优解：能使线性规划目标最优解：能使线性规划目标(mbio)达到极值的达到极值的可行解。可行解。 TTixxsssxsssxxsxsxxsxxsssxxzis002007506000,2005 . 07505 . 1600max00040303 , 2 , 12132103212132221121321

18、21方程的初始解为：约束条件：目标函数：准化：，使线性规划的模型标在上例中加入松弛变量精品资料单纯形凸类中的一种，在其内部任意两点间的连线仍处于图形的内部。单纯形法极点(jdin)迭代法。沿着凸多面体的棱向另一个极点(jdin)迭代，使目标函数的值逐次下降。)3(, 2 , 1, 0)2() 1 (max)min(221122222121112121112211nixbxaxaxabxaxaxabxaxaxaxcxcxcimnmnmmnnnnnn未知量并且要求满足：约束条件或目标函数：标准形式：精品资料从数学意义(yy)上讲，若有一组解1）满足（2）、（3）两式，称为基本可行解；2）同时满足（

19、1）、（2）和（3）三式，称为最优解。 个线性无关的列组成。中基底矩阵。由标准形式还可以化为标准形式的矩阵形式：mABXXbNXBXZXCXCXCXfXbAXXCXfNBNBNNBBTT0, 0min0min精品资料 行解，可停止迭代。，则已经取得了最优可若去寻找最优解；中进行换元，再迭代下、过对组更优的可行解，可通，一定存在一组比前一若判断：设解得目标函数：00,111rXXrNBCCrXNBCCbBCxfNBBNNBNB精品资料lbBBAVVXXXCVVXXlVVBXBXVVZVVXXVVXXVVVXXXxlxBVVTi1111000,minmin即，得数学模型为在非零解。互不独立，不能同

20、时存与、与，即：可正可负，故设、数均为非负值，但由于线性规划中要求参约束条件：目标函数：精品资料当观测当观测(gunc)(gunc)值是等权时，选权迭代法用值是等权时，选权迭代法用权函数进行平差；权函数进行平差；当观测当观测(gunc)(gunc)值不等权时，选权迭代法用值不等权时，选权迭代法用等价权进行平差。等价权进行平差。 式相同。权函数平差时的迭代公平差的时迭代公式与用权函数等价权定义：，或估计准则为：观测值，其权为对于独立不等权ijijjiijjiiiinvPvPpvvpPvpvpMppp0min,21精品资料xxxxTlTTxxiTQDBPBBPQPBBPBQPntnVPV 20 1

21、1 200参数方差近似公式：参数权逆阵近似公式：下的观测值个数。（含粗差观测值）后剩去除似公式：单位权中误差估计的近之所以称以上三式为近似公式，是因为近似地视等价权为常数矩之所以称以上三式为近似公式，是因为近似地视等价权为常数矩阵（其实等价权是随机阵（其实等价权是随机(su j)量，是残差的函数）。量，是残差的函数）。精品资料粗差作为一种模型误差粗差作为一种模型误差(wch)(wch)，可以从两种角度去，可以从两种角度去描述它：描述它：1 1）将粗差归入函数模型）将粗差归入函数模型数据探测法（也称均值漂数据探测法（也称均值漂移模型）移模型）2 2）将粗差归入随机模型）将粗差归入随机模型稳健估计

22、法（也称方差膨稳健估计法（也称方差膨胀模型）胀模型）的作用程度。反映了观测误差对残差称为平差的几何条件。常数阵的影响：对残差观测误差RPBBNIRRPBBNIVVTT11精品资料 测，未参与平差。个观测值为完全多余观，则第若；量，说明其无抗差能力个观测值无多余观测分，则第若好的控制网，有的分布范围：，且多余观测数量个观测值的多余观测分个主元为第阵中的第）称能发现粗差；时，即无需平差，也不，可见无多余观测，则）若无多余观测，即迹”特性，即秩阵是幂等阵，所以有“）产生影响；，将对所有的残差）某一观测误差相关的；是的影响，所以残差之间，它受到所有观测误差）对任意一个残差的描述：对irirrrrrrP

23、QrriiRVlBxtnrrRtrRrkRVvPQRRViiiiiniiVViiiiVV105 . 02 . 0.10:5004;32111精品资料假设观测值中仅有一个有粗差，用该法检测并剔除后，再假设观测值中仅有一个有粗差，用该法检测并剔除后，再建立新的平差系统重新建立新的平差系统重新(chngxn)平差后，再找出下一个平差后，再找出下一个粗差剔除，直到不含粗差。粗差剔除，直到不含粗差。 1) 10:20200uuPrrvQvvuuvEHiilivviviiiiiii接受域：很小时会影响判断统计量：检验法：已知，用若母体方差数据探测法的原假设：精品资料粗差的观测值。初剔除的有时并非是含）因各

24、残差相关，故最差。）一次只能发现一个粗数据探测法缺点：检验粗差。个观测值重新平差，再的剔除粗差后，再将剩余参数个数。其中统计量：检验：）母体方差未知时，用21111122200ntrvpPVVtntntQvttiiiTvviiii精品资料34 为什么要研究病态方程： 当误差(wch)方程为病态时，即使观测数据服从正态分布，其最小二乘估值也不理想，甚至很差，平差结果中方差虽然最小，但方差的值却很大，即平差精度很差，而且解也相当的不稳定。精品资料1 1、矩阵、矩阵(j zhn)(j zhn)的的条件数条件数TxxybyybxxxbAx110001. 00,0001. 220001. 1111022

25、20001. 11112121方程组的解变为：其中有微小变化时，有当常数项其精确解为，即例：有方程组一、病态(bngti)问题与条件数若系数阵若系数阵A A或常数项或常数项b b的微小变化，会引起方程组的解的微小变化，会引起方程组的解x x有巨大变化，有巨大变化，则这种方程组称为则这种方程组称为“病态方程组病态方程组”。A A称为病态矩阵。称为病态矩阵。精品资料AAbbAAAAAAxxbbxxAAxbAbbAAxxbAxxAbbAxbAxbbxxAxbbAbAx11111211得解的相对误差为：即有多大误差？均有误差时，解、）设所以：（范数特性）即有多大误差？，导致解有误差非奇异，设正常，）中

26、，讨论：在方程组精品资料定义定义(dngy(dngy) ) 。中最大和最小的特征值分别是正定对称矩阵和其中，数为阵的条件阵），（如法方程的为正定实对称矩阵时，当感程度。组的解对原始数据的敏条件数。它刻画了方程的为矩阵，称乘积异阵矩阵的条件数：对非奇AAAAcondANAAAAAminmaxminmax121不稳定模型不稳定模型(mxng)(mxng)：观测数据很小的误差会引起待估参数很大的：观测数据很小的误差会引起待估参数很大的误差。所以病态方程也是不稳定模型误差。所以病态方程也是不稳定模型(mxng)(mxng)。精品资料2 2、病态、病态(bngti)(bngti)性程度的衡性程度的衡量方

27、法量方法。时，有严重的复共线性中等强度复共线性；时，有时，有弱复共线性；在复共线性；时，可以认为不存中，的特征值通常的判断标准，模糊的说法。很接近于零”是一个很个复共线性关系。但“中就有多少于零，设计矩阵有多少个特征值很接近法矩阵、特征分析法01. 005. 001. 01 . 005. 01 . 0iiiiiNBNa解释：复共线性解释：复共线性 复共线性，指的是平差参数之间具有近似相关关系，反映在误差方程的设计矩阵上，就是列向量间的某些数据(shj)列可以由其余的数据(shj)列近似（非精确）地线性表示。 在最小二乘平差中，“复共线性”就是指“病态性”。精品资料 方差分解比方法。条件指标法、

28、。中有几个复共线性关系判定设计矩阵条件数法的缺点是不能正取舍。准则应根据实际情况修左右。所以对上述致在快速定位中，条件数大处理实际应用中，如数据提下得到的。但在测量对数据中心标准化的前呈病态。这个指标是在，系统时存在严重的复共线性时没有复共线性；一般认为、条件数法CTVDPcBKKNNNcondKb13minmax110GPS10001003 3、病态方程、病态方程(fngchng)(fngchng)产生原因产生原因1、参数选取原因。（参数近似相关或过度参数化）、参数选取原因。（参数近似相关或过度参数化）2、观测原因。（样本为局部采样或接近重复采样）、观测原因。（样本为局部采样或接近重复采样）

29、3、模型选择原因。（模型建立、模型选择原因。（模型建立(jinl)的方法不同，其病态程度不同）的方法不同，其病态程度不同）4、计算方面原因。（计算方法要稳定，计算机字节长度应长一些）、计算方面原因。（计算方法要稳定，计算机字节长度应长一些）精品资料4 4、病态方程、病态方程(fngchng)(fngchng)最小二乘估值的最小二乘估值的性质性质 病态方程处理的观测值可以是正态分布，但其LS估值并不理想，甚至很差。 虽然LS估计的方差(fn ch)在线性无偏类中是最小，但数值却很大，并表现得相当不稳定。 常用均方误差MSE来评价病态情形下参数的估值质量。 的减小。部分，换取方差部分大偏差有偏估计

30、实质：适当增。估值，即参数估值将不再是无偏用解病态方程法得到的的一个良好估值了。不再是估值很大，此时值较小，会导致的最小特征。而若法矩阵的无偏估值，即是真值上式的条件为由均方误差公式：122120120212 2001xxExxLSxMSENxxBBtrxxEQtrxxxxExMSEitiiTxxT精品资料问题的适定性：问题的适定性：人们根据已获取的观测数据和物理规律，列出的数学模型，当这些模型具人们根据已获取的观测数据和物理规律，列出的数学模型，当这些模型具有下述性质：有下述性质： 1 1、解存在、解存在(cnzi)(cnzi)；、解唯一；、解唯一；、解稳定。、解稳定。则这个问题称为适定性问

31、题。则这个问题称为适定性问题。不适定性：不适定性： 不满足上面三个条件中的任意一个或多个。不满足上面三个条件中的任意一个或多个。不适定问题通常是病态的，但病态问题不一定就是不适定问题。不适定问不适定问题通常是病态的，但病态问题不一定就是不适定问题。不适定问题通常是求方程的稳定近似解。题通常是求方程的稳定近似解。5 5、什么样的方程可能、什么样的方程可能(knng)(knng)是是病态的？病态的？1）行列式的值很大或很小（如某些行、列近代相关）；2）元素间相差大数量级，且无规则；3）主元消去过程中出现小主元；4）特征值相差大数量级。精品资料1 1、病态方程、病态方程(fngchng)(fngch

32、ng)的截断奇异的截断奇异值解法值解法奇异奇异(qy)值分解技术值分解技术(Singular Value Decomposintion Technique，简记为，简记为SVD法法)均为正交矩阵。、为半正定的对角阵；式中阵可分解为时，对）当（进行奇异值分解：下面对的广义逆。是为的最小二乘最小范数解是误差向量。得是设计矩阵，已经单位化）：的权阵测值向量设有观测方程（式中观VUVUAAtnppArankAAALAxxeAeLxAPLttTtnnntnLS,min)(11n1n1ttn精品资料 minmax221212121,:. 0),min(,000AcondAVUvvvVuuuUVUAAAAt

33、nARpdiagDDtniiiTiipptn数与奇异值的关系为：为长方阵时，得其条件：奇异值与条件数的关系见后面算例。阵的计算和阵按列划分，为和将。的关系为：的特征值或与矩阵奇异值阵全部的非零奇异值。是且其中：阵的分块形式为精品资料TTTTTTiTTTTTTvvvdduuVUAAuuDAVUUAArankvVVVvvvVvvvAAddAAAAAA32121212111112132132121321000031313102121626161010003212121212121212110031062216121611101012)(,313131,02121,6261611,3,0,1,3211

34、110101110101的奇异值分解为从而可得），则的个数所占的，（令别为的正交单位特征向量分且。的奇异值为故的特征值为，得解：因为的奇异值分解。例：求矩阵精品资料 。近似秩亏的线性方程组法可解算满秩、秩亏和看出，由，而的通解为：线性方程组组的数值稳定的方法。、特别是线性病态方程分解法是求线性方程组，可见奇异值超过时，奇异值的变化不会有扰动此特性说明当矩阵，有，则对均属于与若、奇异值分解的扰动。SVDVUAtnpdiagDDUVVUALAxLAxEEAEAEAtptnREAATpntTTpptn11121111122,min,000, 2 , 13精品资料件数，提高了稳定性。很强的约束，降低了

35、条这相当于舍去了相关性及其相关的特征向量，的奇异值，舍去小于值截断原则：选择一个阈程的截断奇异值法解：步对其截断，得病态方在第得解为：的奇异值分解式可写为相应地，阵的奇异值为：则，且有的特征值为设inTinTitiitTnpinTitiitLSpiTiiiTpppiLuvxTLuvLAxuvUVAAApiA11111111111111221121,0, 2 , 1通过截断，适当去除（通过截断，适当去除（t-T）个大误差项，恢复了一些解的主要特性）个大误差项，恢复了一些解的主要特性(txng)，但也丧失了一些解的精确性。，但也丧失了一些解的精确性。精品资料47 均为单位阵）内容岭估计。（以下时，

36、正则化估计也称为正则化矩阵。当满足正则化参数光滑函数。RIRRTikhononvxxRxPVVxTT02min2 2、病态、病态(bngti)(bngti)方程的正则方程的正则化解法化解法）式：的正则化准则作用于（，此时可取的特征值单调地趋向于若上式病态，则法矩阵）（有误差方程：101TikhonovNlxBVPlBIPBBxTtT1得正则化参数解为： 可见 ，正则化方法的核心是通过附加“全部或部分参数（或其改正数）加权平方和极小”的条件，增加约束，补充（先验）信息，来克服(kf)不适定性，使解唯一且稳定。精品资料48 QIIINQNQQxxQNQQxMSEMxQxxExBiasxBiasxxBlElIPBBQlxBVPlBQxINQtttTnTTt其中考虑了：阵为评定精度的均方误差矩为的偏差值有参