最小二乘解

定义:任意$A=A_{m \times n}$,方程$AX=b$可产生新方程$A^HAX=A^Hb$,叫$AX=b$的正规方程。

引理:正规方程组$A^HAX=A^Hb$一定有解(相容),且有特解$X_0=A^+b$(使$A^HAX=A^Hb$)

证明:

\[{A^H}A{X_0} = {A^H}A{A^ + }{X_0} = {A^H}{(A{A^ + })^H}{X_0} = {(A{A^ + }A)^H}{X_0} = {A^H}{X_0}\]

注:$A^HAX=A^Hb$有下面通解公式,$rank(A)=r$:

\[X = {X_0} + ({k_1}{Y_1} + ... + {k_{n - r}}{Y_{n - r}})\]

其中$X_0=A^+b$,$Y=k_1Y_1+...+k_{n-r}Y_{n-r}$

定义:任意$A=A_{m \times n}$可产生2个子空间

(1)$N(A)={X | AX=0, X \in C^n}  \subset C^n$,称为A的核空间

(2)$N(A)={AX |  X \in C^n}  \subset C^m$,称为A的像空间又称为值域

定理:

(1)$(A^+A)^2=A^+A$

(2)$(A+A)^2=A^+A$

(3)$AA^+ \ge 0$和$A^+A \ge 0$

 引理:若$AX=b$有解,则必有特解:$X_0=A^+b$

证明:设$X_1$为$AX=b$的任意解,$AX_1=b$,把$X_0=A^+b$代入方程:

\[A{X_0} = A({A^ + }b) = A({A^ + }A{X_1}) = A{A^ + }A{X_1} = A{X_1} = b\]

正交引理:

(1)任意$A=A_{m \times n} \in C^{m \times n}, b \in C^{m}$,$X_0=A^+b  \subset  C^n$,则$X_0 \bot N(A))$

证明:

\[(X,{X_0}) = X_0^HX = {({A^ + }b)^H}X = {({A^ + }A{A^ + }b)^H}X = {({A^ + }b)^H}{({A^ + }A)^H}X = {({A^ + }b)^H}{A^ + }(AX) = 0\]

(2)任意$A=A_{m \times n} \in C^{m \times n}, b \in C^{m}$,$X_0=A^+b  \subset  C^n$,则$(A{X_0} - b) \bot R(A)$

证明:

\[(A{X_0} - b,AX) = {(AX)^H}(A{X_0} - b) = {X^H}{A^H}(A{X_0} - b) = {X^H}({A^H}A{X_0} - {A^H}b) = 0\]

定义:

(1)若$AX=b$无解,则称$AX=b$为不相容方程。

(2)若$AX=b$有解,则称$AX=b$为相容方程。

定理:$\min ({\left| {AX - b} \right|^2}) = A{X_0} - b$,$X_0=AX_0-b$为最小值点

证明:

\[{\left| {AX - b} \right|^2} = {\left| {AX - A{X_0} + A{X_0} - b} \right|^2} = {\left| {A(X - {X_0})} \right|^2} + {\left| {A{X_0} - b} \right|^2} \ge {\left| {A{X_0} - b} \right|^2}\]

定理:若$AX=b$不相容,$X_0=A^+b$是一个最小二乘解,其他最小二乘解X适合条件$AX=AX_9 \Leftrightarrow A(X-X_0)=0$,通解为:

\[X - {X_0} = {t_1}{Y_1} + {t_2}{Y_2} + ... + {t_{n - r}}{Y_{n - r}} \Rightarrow X = {t_1}{Y_1} + {t_2}{Y_2} + ... + {t_{n - r}}{Y_{n - r}} + {X_0}\]

定理:若AX=b不相容,则$X_0=A^+b$恰是全体最小二乘解中的最小长度(2范数)解,被称为极小范数解。

\[{\left| X \right|^2} = {\left| {{X_0} + Y} \right|^2} = {\left| {{X_0}} \right|^2} + {\left| Y \right|^2} \ge {\left| {{X_0}} \right|^2}\]

定理:若$A=A_{m \times n}$与$X=X_{m \times n}$适合条件:$(1)AXA=A$,称X为A的一个减号逆或1号逆,记为$X=A^-$或$X=A^{(1)}$

定理:减号逆矩不唯一

求解$A^-$

(1)

\[A = {\left[ {\begin{array}{*{20}{c}}
{{I_r}}&0\\
0&0
\end{array}} \right]_{m \times n}} \Rightarrow {A^ - } = {\left[ {\begin{array}{*{20}{c}}
{{I_r}}&C\\
D&F
\end{array}} \right]_{n \times m}}\]

$C,D,F$的取值任意

(2)若$A=A_{m \times n}$:

\[PAQ = B = {\left[ {\begin{array}{*{20}{c}}
{{I_r}}&0\\
0&0
\end{array}} \right]_{m \times n}} \Rightarrow {A^ - } = Q{\left[ {\begin{array}{*{20}{c}}
{{I_r}}&C\\
D&F
\end{array}} \right]_{n \times m}}P\]

扩展:矩阵方程$AXB=D$:

(1)若$AXB=D$相容,则必有特解$X_0=A^+DB^+$和$Y_0=A^-DB^-$,使得$AX_0B=D, AY_0B=D$

(2)若$AXB=D$相容,则$X_0=A^+DB^+$是最小范数解

(3)若$AXB=D$不相容,则$X_0=A^+DB^+$是最小二乘解

原文地址:https://www.cnblogs.com/codeDog123/p/10223184.html

时间: 2024-10-07 10:37:47

最小二乘解的相关文章

白话机器学习-最小二乘解和无偏估计

机器学习之最小二乘法解和无偏估计 更新中 ----------------------------------------------- 照例广告时间,休息一下,马上回来 炼数成金课程报名5折优惠码:C181 性价比最高VPN:http://61.164.116.28:8080/?ad_str=49261 广告结束,没偷没抢,自力更生,力求赚回学费,现在学费已经消费近4000了,但我还没有被扣学费过,至少现在是这样. ---------------------------------------

5 最小二乘解

https://blog.csdn.net/kokerf/article/details/72437294 对于线性方程组,解的判别条件如下: 1. Ax=0Ax=0 总有解,至少有零解 2. Am×nx=0Am×nx=0 当r(A)=nr(A)=n,只有零解 当r(A)<nr(A)<n,有无穷多解 3. Am×nx=bAm×nx=b 当r(A)≠r(A|b)r(A)≠r(A|b),无解 当r(A)=r(A|b)=nr(A)=r(A|b)=n,有唯一解 当r(A)=r(A|b)=r<nr

《Python数据分析常用手册》一、NumPy和Pandas篇

一.常用链接: 1.Python官网:https://www.python.org/ 2.各种库的whl离线安装包:http://www.lfd.uci.edu/~gohlke/pythonlibs/#scikit-learn 3.数据分析常用库的离线安装包(pip+wheels)(百度云):http://pan.baidu.com/s/1dEMXbfN 密码:bbs2 二.常用库 1.NumPy NumPy是高性能科学计算和数据分析的基础包.部分功能如下: ndarray, 具有矢量算术运算和

【数理统计基础】 05 - 回归分析

参数估计和假设检验是数理统计的两个基础问题,它们不光运用于常见的分布,还会出现在各种问题的讨论中.本篇开始研究另一大类问题,就是讨论多个随机变量之间的关系.现实生活中的数据杂乱无章,够挖掘出各种变量之间的关系非常有用,它可以预估变量的走势,能帮助分析状态的根源.关系分析的着手点可以有很多,我们从最简单直观的开始,逐步展开讨论. 1. 一元线性回归 1.1 回归分析 如果把每个量都当做随机变量,问题的讨论会比较困难,或者得到的结论会比较受限.一个明智做法就是只把待考察的量\(Y\)看做随机变量,而

[Matlab]求解线性方程组

转自:http://silencethinking.blog.163.com/blog/static/911490562008928105813169/ AX=B或XA=B在MATLAB中,求解线性方程组时,主要采用前面章节介绍的除法运算符“/”和“\”.如: X=A\B表示求矩阵方程AX=B的解: X=B/A表示矩阵方程XA=B的解. 对方程组X=A\B,要求A和B用相同的行数,X和B有相同的列数,它的行数等于矩阵A的列数,方程X=B/A同理. 如果矩阵A不是方阵,其维数是m×n,则有: m=

小议线性回归

统计学习好玩,但不好学.无论是课本.讲义还是论文,都充斥着连篇累牍的公式,当然这是科学文献严密性的需要,只是对初学者来说,学习门槛太高.最近在看<统计学习要素>(The Elements of Statistical Learning),随手心得,仅供参考. 第三章-线性回归 如果我们已知了X和Y之间满足线性关系,在得到X和Y的数据之后,怎样估计它们之间的线性系数w呢? 这其实是在解一个线性方程组,xw = y.这个等式的意义是,用w作为x列向量的加权,得到y.但理论上,如果y并不在x的列空间

Scipy学习笔记 矩阵计算

Scipy学习笔记 非本人原创  原链接 http://blog.sina.com.cn/s/blog_70586e000100moen.html 1.逆矩阵的求解 >>>import scipy >>>from scipy import linalg >>>a=scipy.mat('[1 2 3;2 2 1;3 4 3]') >>>b=linalg.inv(a) >>>print b 输出结果 [[ 1.   3.

Python闲谈(二)聊聊最小二乘法以及leastsq函数

1 最小二乘法概述 自从开始做毕设以来,发现自己无时无刻不在接触最小二乘法.从求解线性透视图中的消失点,m元n次函数的拟合,包括后来学到的神经网络,其思想归根结底全都是最小二乘法. 1-1 “多线→一点”视角与“多点→一线”视角 最小二乘法非常简单,我把它分成两种视角描述: (1)已知多条近似交汇于同一个点的直线,想求解出一个近似交点:寻找到一个距离所有直线距离平方和最小的点,该点即最小二乘解: (2)已知多个近似分布于同一直线上的点,想拟合出一个直线方程:设该直线方程为y=kx+b,调整参数k

MATLAB基本命令

eye(n)创建n阶单位矩阵 zeros(n)创建n阶0方阵 rand(m,n)创建m*n阶元素为从0到1的均匀分布的随机数矩阵 round(A)对矩阵A中所有元素进行四舍五入运算 A^-1用幂运算求解矩阵A的逆 \左除,/右除, .(点号)群运算 已知:PA=B,AQ=B P=B/A(右除), Q=A\B(左除) U = rref(A)对矩阵A进行初等行变换,矩阵U为矩阵A的最简行阶梯矩阵 u = rref([a, eye(size(a))])通过初等行变换的方式得到a的逆 n = input