码农知识堂 - 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)来判断
    奇异矩阵不是满秩

    参考链接

    参考链接

    参考链接

  • 相关阅读:
    QLineEdit 类(行编辑器)
    创建MySQL只读权限用户
    你是否感受到AI就在身边?
    linux运维基础2
    windows RocketMQ与可视化监控平台安装
    分类预测 | Matlab实现基于LFDA-SVM局部费歇尔判别数据降维结合支持向量机的多输入分类预测
    Git一台电脑同时绑定两个账号(如gitlab,github)
    数组(C语言)(详细过程!!!)
    【模型训练】YOLOv7不同颜色安全帽佩戴检测
    谁能想到Java技术能跟看小说一样学习,不知不觉就看完了,膜拜大老!
  • 原文地址:https://blog.csdn.net/weixin_44710960/article/details/126893675
  • 最新文章
  • 沪漂五周年了:我越来越迷茫了
    Agentic Skill Routing 实战:别再把所有 Skill 塞进 AI Agent 上下文
    MySQL-Seconds_behind_master的精度误差
    [MAF预定义ChatClient中间件-03]CachingChatClient——利用缓存省钱省时间
    AI的至暗历史:从万众期待到被政府撤资,AI的两次死亡徘徊
    Agent OS :五种驯服不确定性的范式
    PortSwigger SQL注入LAB11
    数据库即时编译JIT
    [Begin]AI Learn Data Day 0
    深度学习进阶(二十七)现代 LLM 的核心架构设计其二:SwiGLU
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号