当前位置: 首页 > news >正文

a = Sw,其中a和w是向量,S是矩阵,求w等于什么?w可以写成关于a和S的什么样子的公式

给定公式:

a = S w a = S w a=Sw

其中:

  • a a a 是已知向量,
  • S S S 是已知矩阵,
  • w w w 是未知向量。

我们的目标是求解 w w w,即将 w w w 表示为 a a a S S S 的函数。


情况 1:矩阵 S S S 可逆

如果矩阵 S S S可逆矩阵,即 S S S 是方阵(行数等于列数)且行列式不为零,那么 S S S 存在逆矩阵 S − 1 S^{-1} S1。在这种情况下,可以直接求解:
w = S − 1 a \boxed{w = S^{-1} a} w=S1a

步骤:

  1. 计算 S S S 的逆矩阵 S − 1 S^{-1} S1

    • 仅当 S S S 为方阵且满秩(行列式不为零)时,逆矩阵 S − 1 S^{-1} S1 存在。
  2. 计算 w w w
    w = S − 1 a w = S^{-1} a w=S1a

示例:

设:
S = ( 2 1 0 1 ) , a = ( 5 1 ) S = \begin{pmatrix} 2 & 1 \\ 0 & 1 \end{pmatrix}, \quad a = \begin{pmatrix} 5 \\ 1 \end{pmatrix} S=(2011),a=(51)

计算 S − 1 S^{-1} S1
S − 1 = ( 1 2 − 1 2 0 1 ) S^{-1} = \begin{pmatrix} \frac{1}{2} & -\frac{1}{2} \\ 0 & 1 \end{pmatrix} S1=(210211)

求解 w w w
w = S − 1 a = ( 1 2 − 1 2 0 1 ) ( 5 1 ) = ( 1 2 × 5 − 1 2 × 1 1 × 1 ) = ( 2 1 ) w = S^{-1} a = \begin{pmatrix} \frac{1}{2} & -\frac{1}{2} \\ 0 & 1 \end{pmatrix} \begin{pmatrix} 5 \\ 1 \end{pmatrix} = \begin{pmatrix} \frac{1}{2} \times 5 - \frac{1}{2} \times 1 \\ 1 \times 1 \end{pmatrix} = \begin{pmatrix} 2 \\ 1 \end{pmatrix} w=S1a=(210211)(51)=(21×521×11×1)=(21)


情况 2:矩阵 S S S 不可逆或非方阵

如果矩阵 S S S 不可逆(行列式为零)或是一个非方矩阵(行数不等于列数),则不能直接计算逆矩阵。在这种情况下,可以使用广义逆矩阵(Moore-Penrose 伪逆)来求解。

求解步骤:

  1. 计算 S S S 的伪逆矩阵 S + S^{+} S+

    • S + S^{+} S+ 满足 S + = ( S T S ) − 1 S T S^{+} = (S^T S)^{-1} S^T S+=(STS)1ST,当 S T S S^T S STS 可逆时。
    • 对于更一般的情况,可以使用奇异值分解(SVD)来计算 S + S^{+} S+
  2. 计算 w w w
    w = S + a \boxed{w = S^{+} a} w=S+a

注意:

  • 当方程组有唯一解时, w = S + a w = S^{+} a w=S+a 给出唯一解。
  • 当方程组有无穷多解时, w = S + a w = S^{+} a w=S+a 给出最小范数解。
  • 当方程组无解(即 a a a 不在 S S S 的列空间内)时, w = S + a w = S^{+} a w=S+a 给出使 ∥ S w − a ∥ 2 \| S w - a \|_2 Swa2 最小的解。

示例:

设:
S = ( 1 2 3 4 5 6 ) , a = ( 7 8 9 ) S = \begin{pmatrix} 1 & 2 \\ 3 & 4 \\ 5 & 6 \end{pmatrix}, \quad a = \begin{pmatrix} 7 \\ 8 \\ 9 \end{pmatrix} S= 135246 ,a= 789

S S S 是一个 3 × 2 3 \times 2 3×2 的矩阵,无法直接求逆。

计算 S S S 的伪逆 S + S^{+} S+

  1. 计算 S T S S^T S STS
    S T S = ( 1 3 5 2 4 6 ) ( 1 2 3 4 5 6 ) = ( 35 44 44 56 ) S^T S = \begin{pmatrix} 1 & 3 & 5 \\ 2 & 4 & 6 \end{pmatrix} \begin{pmatrix} 1 & 2 \\ 3 & 4 \\ 5 & 6 \end{pmatrix} = \begin{pmatrix} 35 & 44 \\ 44 & 56 \end{pmatrix} STS=(123456) 135246 =(35444456)

  2. 计算 ( S T S ) − 1 (S^T S)^{-1} (STS)1
    ( S T S ) − 1 = 1 35 × 56 − 4 4 2 ( 56 − 44 − 44 35 ) ≈ ( 1.4 − 1.1 − 1.1 0.875 ) (S^T S)^{-1} = \frac{1}{35 \times 56 - 44^2} \begin{pmatrix} 56 & -44 \\ -44 & 35 \end{pmatrix} \approx \begin{pmatrix} 1.4 & -1.1 \\ -1.1 & 0.875 \end{pmatrix} (STS)1=35×564421(56444435)(1.41.11.10.875)

  3. 计算 S + = ( S T S ) − 1 S T S^{+} = (S^T S)^{-1} S^T S+=(STS)1ST
    S + = ( 1.4 − 1.1 − 1.1 0.875 ) ( 1 3 5 2 4 6 ) = ( − 0.2 0 0.2 0.15 0 − 0.15 ) S^{+} = \begin{pmatrix} 1.4 & -1.1 \\ -1.1 & 0.875 \end{pmatrix} \begin{pmatrix} 1 & 3 & 5 \\ 2 & 4 & 6 \end{pmatrix} = \begin{pmatrix} -0.2 & 0 & 0.2 \\ 0.15 & 0 & -0.15 \end{pmatrix} S+=(1.41.11.10.875)(123456)=(0.20.15000.20.15)

  4. 计算 w = S + a w = S^{+} a w=S+a
    w = ( − 0.2 0 0.2 0.15 0 − 0.15 ) ( 7 8 9 ) = ( ( − 0.2 ) × 7 + 0 × 8 + 0.2 × 9 0.15 × 7 + 0 × 8 + ( − 0.15 ) × 9 ) = ( 0.4 0 ) w = \begin{pmatrix} -0.2 & 0 & 0.2 \\ 0.15 & 0 & -0.15 \end{pmatrix} \begin{pmatrix} 7 \\ 8 \\ 9 \end{pmatrix} = \begin{pmatrix} (-0.2) \times 7 + 0 \times 8 + 0.2 \times 9 \\ 0.15 \times 7 + 0 \times 8 + (-0.15) \times 9 \end{pmatrix} = \begin{pmatrix} 0.4 \\ 0 \end{pmatrix} w=(0.20.15000.20.15) 789 =((0.2)×7+0×8+0.2×90.15×7+0×8+(0.15)×9)=(0.40)


总结

  • S S S 可逆时 w = S − 1 a w = S^{-1} a w=S1a
  • S S S 不可逆或非方阵时 w = S + a w = S^{+} a w=S+a,其中 S + S^{+} S+ S S S 的伪逆矩阵。

注意事项

  • 求逆矩阵的条件:只有当矩阵 S S S 满足可逆条件(方阵且满秩)时,才可以直接求逆。
  • 计算伪逆的方法:对于任意矩阵 S S S,都可以使用奇异值分解(SVD)来计算其伪逆 S + S^{+} S+
  • 数值计算:在实际计算中,建议使用数学软件(如 MATLAB、NumPy 等)来计算逆矩阵或伪逆矩阵,以避免手算过程中的误差。
http://www.lryc.cn/news/449554.html

相关文章:

  • 多线程事务管理:Spring Boot 实现全局事务回滚
  • Vue3 中集成海康 H5 监控视频播放功能
  • Linux: eBPF: libbpf-bootstrap-master 编译
  • 1.1.4 计算机网络的分类
  • 周家庄智慧旅游小程序
  • 【在Linux世界中追寻伟大的One Piece】命名管道
  • 如意控物联网项目-ML307R模组软件及硬件调试环境搭建
  • 大模型分布式训练并行技术(九)-总结
  • uniapp view设置当前view之外的点击事件
  • 【Mybatis篇】动态SQL的详细带练
  • 【MyBatis-Plus】 学习记录 常用功能及代码生成器使用
  • HalconDotNet实现OCR详解
  • 手搓一个Agent#Datawhale 组队学习Task3
  • 基于SpringBoot+Vue+MySQL的在线酷听音乐系统
  • 大数据实时数仓Hologres(一):Hologres 简单介绍
  • 【鸿蒙HarmonyOS NEXT】数据存储之分布式键值数据库
  • 基于springboot+小程序的儿童预防接种预约管理系统(疫苗1)(源码+sql脚本+视频导入教程+文档)
  • 计算物理精解【8】-计算原理精解【5】
  • 【Linux】 tcp | 解除服务器对tcp连接的限制 | 物联网项目配置
  • 如何隐藏Windows10「安全删除硬件」里的USB无线网卡
  • 【QT Quick】基础语法:导入外部JS文件及调试
  • 【质优价廉】GAP9 AI算力处理器赋能智能可听耳机,超低功耗畅享未来音频体验!
  • 用Flutter几年了,Flutter每个版本有什么区别?
  • 解决Qt每次修改代码后首次运行崩溃,后几次不崩溃问题
  • 语言的变量交换
  • 【muduo源码分析】「阻塞」「非阻塞」「同步」「异步」
  • 顶顶通呼叫中心中间件-机器人话术挂机后是否处理完成事件
  • Springboot Mybatis 动态SQL
  • ORM的了解
  • 关于大模型的10个思考