利用矩阵的逆(伪逆)与除法求解:

对于线性方程组Ax=b,若其为恰定方程组且A是非奇异的,则求x最明显的方法便是利用矩阵的逆,即x=A\b;若不是恰定方程组,则可利用伪逆来求其一个特解。

运用pinv命令求伪逆:

format rat %表示以有理数形式输出结果
A=[1,2,2,0;0,1,-2,-2;1,3,0,-2]
b=[1,2,3]' %'表示矩阵转置。当矩阵为复数时,求矩阵的共轭转置
x0=pinv(A)*b
Z=null(A,'r')

实例——矩阵特征值

演示求逆法与除法求解线性方程组在时间与精度上的区别

(tic和toc用来记录matlab命令执行的时间。tic用来保存当前时间,而后使用toc来记录程序完成时间。 两者往往结合使用)

对矩阵求逆。如果A是非奇异方阵,则B/A = B*inv(A),A\B = inv(A)*B。/表示右除,\表示左除)。

注意:使用inv时,必须对象为方阵。

A=1000*rand(1000,1000);
x=ones(1000,1);
b=A*x;
disp('利用矩阵的逆求解所用时间及误差为:')
tic %tic和toc用来记录matlab命令执行的时间。tic用来保存当前时间,而后使用toc来记录程序完成时间。 两 %者往往结合使用
y=inv(A)*b;%对矩阵求逆。如果A是非奇异方阵,则B/A = B*inv(A),A\B = inv(A)*B。/表示右除,\表示左 %除
t1=toc
error1=norm(y-x)
disp('利用除法求解所用时间及误差为:')
tic
y=A\b;
t2=toc
error2=norm(y-x)

利用矩阵的逆(伪逆)与除法求解相关推荐

  1. java 矩阵包求伪逆_Coursera-ML-AndrewNg-Notes

    第2周 [TOC] 四.多变量线性回归(Linear Regression with Multiple Variables) 4.1 多维特征 参考视频: 4 - 1 - Multiple Featu ...

  2. 机器人学笔记(3)矩阵的伪逆

    机器人学笔记(3)矩阵的伪逆 左逆矩阵与右逆矩阵 对于机器人而言,求解雅可比的矩阵可以通过机器人各个关节的速度求得末端的速度与角速度,而求解逆运动学时,经常会出现已知末端机器人末端速度,求机器人各个关 ...

  3. 矩阵的逆、伪逆、左右逆,最小二乘,投影矩阵

    主要内容: 矩阵的逆.伪逆.左右逆 矩阵的左逆与最小二乘 左右逆与投影矩阵 一.矩阵的逆.伪逆.左右逆 1.矩阵的逆 定义: 设A是数域上的一个n阶方阵,若在相同数域上存在另一个n阶矩阵B,使得: A ...

  4. numpy求矩阵的逆和伪逆

    我们可以使用np.linalg中的inv和pinv函数来求解矩阵的逆/伪逆. np.linalg.inv 对于可逆方阵M,我们使用下面这行代码求逆: np.linalg.inv(J(theta)) 示 ...

  5. (数学概念)矩阵的逆、伪逆、左右逆,最小二乘,投影矩阵

    主要内容: 矩阵的逆.伪逆.左右逆 矩阵的左逆与最小二乘 左右逆与投影矩阵 一.矩阵的逆.伪逆.左右逆 1.矩阵的逆 定义: 设A是数域上的一个n阶方阵,若在相同数域上存在另一个n阶矩阵B,使得: A ...

  6. 奇异矩阵,非奇异矩阵,矩阵的逆和伪逆

    奇异矩阵/非奇异矩阵 首先看矩阵是不是方阵,只有是方阵了,才有这两个概念. 方阵A的行列式等于零,记为|A|=0,A是奇异矩阵 方阵A的行列式不等于零,记为|A|~=0,A是非奇异矩阵 一些性质: 如 ...

  7. 矩阵伪逆的opencv实现

    1.矩阵的逆 定义: 设A是数域上的一个n阶方阵,若在相同数域上存在另一个n阶矩阵B,使得: AB=BA=I. 则我们称B是A的逆矩阵,而A则被称为可逆矩阵. 可逆条件: A是可逆矩阵的充分必要条件是 ...

  8. matlab求矩阵违逆,两种求矩阵伪逆的方法

    伪逆矩阵是逆矩阵的广义形式.由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在matlab里可以用函数pinv(A)求其伪逆矩阵.基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差, ...

  9. 清华大学公开课线性代数2——第6讲:伪逆

    此博客停止更新,迁移至SnailDove's blog,查看本文请点击此处,清华大学线性代数2笔记汇总:线性代数总结 笔记源自:清华大学公开课:线性代数2--第6讲:伪逆 **提示:**如果文中图片看 ...

  10. opencv求矩阵的逆_OpenCV 矩阵操作 加减乘除 求逆 行列式值 | 学步园

    GEMM 通用矩阵乘法 void cvGEMM( const CvArr* src1, const CvArr* src2, double alpha, const CvArr* src3, doub ...

最新文章

  1. leetcode--最小路径和--python
  2. python 抓取“煎蛋妹子图”并分页存在本地(普通版和多线程版)
  3. ceph与hdfs的比较_分布式存储中HDFS与Ceph两者的区别是什么,各有什么优势?
  4. C#创建WinForm应用程序的入口点
  5. SAP 电商云 Spartacus UI 的响应式 UI 实现细节
  6. 金士顿 8G u盘 红色 量产记录
  7. 当爬虫工程师遇到 CTF丨2021 年 B 站 1024 安全攻防题解
  8. nginx-exporter安装使用
  9. 维纳过程(Wiener Process)与高斯过程(Gaussian Process)
  10. Kubernetes 1.14.1快速升级
  11. logisim数据选择器_基于Logisim的CRC编码电路实现.pdf
  12. ectouch 概述
  13. smtplib 抄送邮件_Smtp发送及抄送邮件
  14. Python-Scrapy-抓取链家二手房信息
  15. 2022中元节前后几天不出门?前三天后三天不能出门是真的吗?
  16. xshell免费版安装
  17. C++不知算法系列之解析回溯算法中的人文哲学
  18. PLC(Power Line Carrier)电力线载波介绍
  19. gerrit git 邮箱不匹配的问题
  20. 芯片失效分析手段系列之X-RAY探伤

热门文章

  1. android编程微博的发送,基于Android手机微博客户端的开发.doc
  2. python英文情绪识别_Python能识别文字情绪?
  3. 计算机平均值的快捷键,Excel用快捷键和选项求平均值,且能一次对多行多列批量快速求平均值...
  4. 华为手机的nfc是什么功能_NFC功能是什么意思,华为手机NFC功能怎么用,NFC常用功能场景...
  5. shell脚本造数据
  6. C++训练:自由落体问题
  7. Java对接苹果账号授权登录
  8. 二元二次方程例题_二元二次方程练习题.doc
  9. 微信对账单功能开发(V2)
  10. 干货!自监督学习在推荐系统中的应用