码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • MATLAB----矩阵求逆的123!


    方阵求逆

    奇异矩阵:行列式为0的矩阵

    非奇异矩阵:行列式不为0的矩阵,即可逆矩阵

    对于方阵A,如果为非奇异方阵,则存在逆矩阵inv(A)

    对于奇异矩阵或者非方阵,并不存在逆矩阵,但可以使用pinv(A)求其伪逆

    非方阵

    如果不是方阵的话,会有个广义逆矩阵(伪逆),matlab里面用pinv()实现。
    基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差,pinv为pseudo-inverse的缩写:max(size(A))*norm(A)*eps。函数返回一个与A的转置矩阵A’ 同型的矩阵X,并且满足:AXA=A,XAX=X.此时,称矩阵X为矩阵A的伪逆,也称为广义逆矩阵。

    >> A=rand(2,3)
    A =
        0.9572    0.8003    0.4218
        0.4854    0.1419    0.9157
    >> B=pinv(A)
    B =
        0.5985   -0.0838
        0.7633   -0.5430
       -0.4355    1.2206
    >> A*B*A
    ans =
        0.9572    0.8003    0.4218
        0.4854    0.1419    0.9157
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10
    • 11
    • 12
    • 13

    接近奇异的矩阵求逆

    Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND =
    1.925555e-20.

    矩阵中的元素有的很接近,可能这样导致了矩阵接近奇异矩阵,出现上述警告。但是从数学上说该矩阵并非奇异矩阵,想求出准确的逆矩阵

    总结一下,就是首先给矩阵乘以一个常数,目的是使矩阵的行列式不至于太小,而后再用12楼的方法,用sym转化后计算。

    12楼的方法:

    现在的matla的符号计算能力不错!我的2013版。 可以开始时,转化为符号。 如 a=inv(sym([1e-29 1e-17;1e-16
    1e-3])); a =

    [ 1000000000000000000000000000000/9, -10000000000000000/9] [
    -100000000000000000/9, 10000/9]

    奇异矩阵的判断

    奇异矩阵一般可以通过求秩(rank)来判断
    奇异矩阵不是满秩

    参考链接

    参考链接

    参考链接

  • 相关阅读:
    485自收发专用芯片
    Flink---5、聚合算子、用户自定义函数、物理分区算子、分流、合流
    PTA 7-92 简化的插入排序
    Midjourney AI绘画中文教程详解(完整版)模型、命令、参数与各种高级用法
    mac 上 安装配置mat
    TortoiseGit图标覆盖的意义
    HedgeDoc的反向代理设置
    通过配置文件(.htaccess)实现文件上传
    每个开发都应该懂的正则表达式
    1.1.OpenCV技能树--第一单元--OpenCV简介
  • 原文地址:https://blog.csdn.net/weixin_44710960/article/details/126893675
  • 最新文章
  • 【JVM】编译执行与解释执行的区别是什么?JVM 使用哪种方式?
    用 Hashids 优雅解决 C 端自增 ID 暴露问题
    V8引擎 精品漫游指南--Ignition篇(上) 指令 栈帧 槽位 调用约定 内存布局 基础内容
    LLVM Pass快速入门(四):代码插桩
    milkup:桌面端 markdown AI续写和即时渲染
    基于项目工程构建SBOM(软件物料清单)的研究
    鸿蒙应用开发UI基础第二节:鸿蒙应用程序框架核心解析与实操
    .NET 中如何快速实现 List 集合去重?
    扣子Coze实战:从0到1打造抖音+小红书热点监控智能体
    浅谈数据访问层
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号