正则化最小二乘与条件数(cond)

来源于自学《凸优化》和《矩阵分析与应用》笔记

  1. 正则化最小二乘

给定 AϵRm×n A ϵ R m × n , bϵRm b ϵ R m ,为函数F1和F2两个目标的优化问题, AxAxblabel,y A x 中 A 为 已 知 系 数 矩 阵 , x 表 示 要 求 的 系 数 参 数 , b 表 示 真 实 值 , l a b e l , 或 者 y 值

F1(x)=||Axb||22=xTATAx2bTAx+bTb F 1 ( x ) = | | A x − b | | 2 2 = x T A T A x − 2 b T A x + b T b

F2(x)=||x||22=xTx F 2 ( x ) = | | x | | 2 2 = x T x

也就是无约束的两准则问题,

minimizef0(x)=(F1(x),F2(x)) m i n i m i z e f 0 ( x ) = ( F 1 ( x ) , F 2 ( x ) )

选择 λ1>0,λ2>0 λ 1 > 0 , λ 2 > 0
则问题变为,
λTf0(x)=λ1F1(x)+λ2F2(x) λ T f 0 ( x ) = λ 1 F 1 ( x ) + λ 2 F 2 ( x )

=λT(λ1ATA+λ2I)x2λ1bTAx+λ1bTb = λ T ( λ 1 A T A + λ 2 I ) x − 2 λ 1 b T A x + λ 1 b T b

因为 Ax=b A x = b ,所以上式化为
xT(λATA+λ2I)xλ1bTAx x T ( λ A T A + λ 2 I ) x − λ 1 b T A x

x=(λ1ATA+λ2I)1λ1ATb=(ATA+μI)1ATb x = ( λ 1 A T A + λ 2 I ) − 1 λ 1 A T b = ( A T A + μ I ) − 1 A T b

其中 μ=λ2/λ1 μ = λ 2 / λ 1 ,
此时 μ μ 理解为F2相对于F1的权值

2 条件数,condition number
实际上,我们平时见到的普通最小二乘函数是这样的, Ax=b A x = b 也就是
minmize(Axb) m i n m i z e ∑ ( A x − b ) ,采用平方/二次的形式,或者叫欧式距离,是\

minmize(Axb)2 m i n m i z e ∑ ( A x − b ) 2

用矩阵表示也就是
minmize(Axb)(ATxTbT) m i n m i z e ∑ ( A x − b ) ( A T x T − b T )

其中的假设是,自变量x不含误差,只有因变量b存在服从正态分布的误差,但是很多时候,我们违反了假设,自变量也有误差,所以需要把上式化为,
(A+σA)(x+σx)=b ( A + σ A ) ( x + σ x ) = b

上式推导出
σx=(A+σA)1bA1b σ x = ( A + σ A ) − 1 b − A − 1 b

=[(A+σA)1A1]b = [ ( A + σ A ) − 1 − A − 1 ] b

=A1σA(x+σx) = − A − 1 σ A ( x + σ x )

由此得, ||σx||2||A1||2||σA||2||x+σx||2 | | σ x | | 2 ≤ | | A − 1 | | 2 | | σ A | | 2 | | x + σ x | | 2

||σx||2||x+σx||2=(||A||2||A||12)||σA||2||A||2 | | σ x | | 2 | | x + σ x | | 2 = ( | | A | | 2 | | A | | 2 − 1 ) | | σ A | | 2 | | A | | 2
,
向量X的相对误差与数值, cond(A)=||A||2||A1||2 c o n d ( A ) = | | A | | 2 ⋅ | | A − 1 | | 2
cond(A)Ak(A) c o n d ( A ) 为 矩 阵 A 的 条 件 数 , 有 时 也 用 k ( A ) 表 示

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值