线性代数之超定方程使用最小二乘求解

在超定方程中, 线性方程组一般不能被完全精确满足,因此需要使用最小二乘法求解参数 XXX。

假设线性方程组可以写成:

AX≈B AX \approx B AX≈B

其中,AAA 为系数矩阵,XXX 为待求参数向量,BBB 为观测值向量。

定义残差向量为:

r=AX−B r = AX - B r=AX−B

最小二乘法的目标是使残差平方和最小,即:

X=arg⁡min⁡X∥AX−B∥2 X = \arg\min_X \|AX - B\|^2 X=argXmin∥AX−B∥2

为了便于推导,定义目标函数:

J(X)=∥AX−B∥2 J(X) = \|AX - B\|^2 J(X)=∥AX−B∥2

向量二范数的平方可以写成内积形式:

J(X)=(AX−B)T(AX−B) J(X) = (AX - B)^T(AX - B) J(X)=(AX−B)T(AX−B)

展开可得:

J(X)=XTATAX−XTATB−BTAX+BTB J(X) = X^T A^T A X - X^T A^T B - B^T A X + B^T B J(X)=XTATAX−XTATB−BTAX+BTB

由于 XTATBX^T A^T BXTATB 和 BTAXB^T A XBTAX 都是标量,并且二者相等: (计算两者矩阵维度)

XTATB=BTAX X^T A^T B = B^T A X XTATB=BTAX

所以目标函数可以化简为:

J(X)=XTATAX−2XTATB+BTB J(X) = X^T A^T A X - 2X^T A^T B + B^T B J(X)=XTATAX−2XTATB+BTB

下面对 J(X)J(X)J(X) 关于 XXX 求导。

首先有:

∂∂X(XTATAX)=2ATAX \frac{\partial}{\partial X} (X^T A^T A X) = 2A^TAX ∂X∂(XTATAX)=2ATAX

这是因为 ATAA^TAATA 是对称矩阵,即:

(ATA)T=ATA (A^TA)^T = A^TA (ATA)T=ATA

其次有:

∂∂X(−2XTATB)=−2ATB \frac{\partial}{\partial X} (-2X^TA^TB) = -2A^TB ∂X∂(−2XTATB)=−2ATB

最后,BTBB^TBBTB 中不含有待求参数 XXX,因此:

∂∂X(BTB)=0 \frac{\partial}{\partial X} (B^TB) = 0 ∂X∂(BTB)=0

所以目标函数的导数为:

∂J(X)∂X=2ATAX−2ATB \frac{\partial J(X)}{\partial X} = 2A^TAX - 2A^TB ∂X∂J(X)=2ATAX−2ATB

最小二乘问题要求目标函数取得最小值,因此令导数为零:

2ATAX−2ATB=0 2A^TAX - 2A^TB = 0 2ATAX−2ATB=0

化简可得:

ATAX=ATB A^TAX = A^TB ATAX=ATB

这就是最小二乘法对应的正规方程。

当 ATAA^TAATA 可逆时,可以在等式两边左乘 (ATA)−1(A^TA)^{-1}(ATA)−1,得到:

X=(ATA)−1ATB X = (A^TA)^{-1}A^TB X=(ATA)−1ATB

因此,正规方程的本质就是:对最小二乘目标函数 J(X)=∥AX−B∥2J(X)=\|AX-B\|^2J(X)=∥AX−B∥2 关于待求参数 XXX 求导,并令导数为零后得到的方程。

相关推荐
树獭哥7 分钟前
量化金融入门:从数据规则到随机游走与凯利公式
人工智能·算法·金融·量化金融
联盟分享专家8 分钟前
2026独立站引流推广:低成本获取网站精准流量的7个渠道
人工智能
G智爱AI14 分钟前
2026工作汇报PPT AI工具实测|主流省时办公工具横向对比
人工智能·powerpoint
迪康coolmu15 分钟前
当大模型遇上终端安全—AI驱动的智能运维实践
运维·人工智能·安全
财迅通Ai16 分钟前
科源制药携手腾讯云推进AI办公迭代升级
人工智能·云计算·腾讯云·科源制药
DS随心转APP21 分钟前
实测 AI 导出鸭实操效果,依托 AI 直接生成 word 功能横向比对五类导出办法,梳理办公文档转换最优路径
人工智能·ai·word·deepseek·ai导出鸭
前沿在线30 分钟前
2026 WRC世界机器人大会| 黑芝麻智能展台亮点
人工智能·ai·大模型
QYRdata36 分钟前
28.7%年复合增速锚定AI模型安全赛道,2026-2032年行业步入高速扩容新阶段
人工智能·安全
皮卡丘不断更37 分钟前
Vibe Coding 做完原型后:给网页项目加上任务、验证和人工确认
人工智能·软件工程
Leo.yuan44 分钟前
FineDataLink 5.0 vs Great Expectations:数据质量治理工具深度横评,国产低代码方案能否替代开源标杆?
大数据·人工智能