Cond(A)称作矩阵A的条件数,为矩阵A的范数与A的逆矩阵的范数的乘积
定义 
在MATLAB中,计算矩阵A的3种条件数的函数是:
(1) cond(A,1) 计算A的1—范数下的条件数。
(2) cond(A)或cond(A,2) 计算A的2—范数数下的条件数。
(3) cond(A,inf) 计算A的 ∞—范数下的条件数。
(关于范数知识请看我博客中关于norm函数的相关介绍)
举例
输入矩阵:A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19];
命令:cond(A,1) 结果:ans = 92.1325
命令:cond(A,2) 结果:ans =48.3117
命令:cond(A,inf) 结果:ans =68.6582
在matlab中,判断矩阵的奇异性一般用函数cond,求矩阵的条件数,即矩阵最大奇异值与最小奇异值的比值。当矩阵的条件数很大时,说明矩阵趋近于奇异
矩阵的条件数
条件数事实上表示了矩阵计算对于误差的敏感性。对于线性方程组Ax=b,如果A的条件数大,b的微小改变就能引起解x较大的改变,数值稳定性差。如果A的条件数小,b有微小的改变,x的改变也很微小,数值稳定性好。它也可以表示b不变,而A有微小改变时,x的变化情况。
比如线性方程组
的解是(x,y)=(2,1),
的解是(x,y)=(-3.999,4.000)
可见b很小的扰动就引起了x很大的变化,这就是A矩阵条件数大的表现。
一个极端的例子,当A奇异时,条件数为无穷,这时即使不改变b,x也可以改变。奇异的本质原因在于矩阵有0特征值,x在对应特征向量的方向上运动不改变Ax的值。如果一个特征值比其它特征值在数量级上小很多,x在对应特征向量方向上很大的移动才能产生b微小的变化,这就解释了为什么这个矩阵为什么会有大的条件数,事实上,正规阵在二范数下的条件数就可以表示成 abs(最大特征值/最小特征值)。

matlabnbsp;cond函数nbsp;矩阵的条件数相关推荐

  1. 条件数cond函数(matlab)

    条件数 为什么计算矩阵的条件数 如何调用matlab内置函数 为什么计算矩阵的条件数 矩阵的条件数是判断矩阵"病态"程度的一个指标:用于衡量线性方程组的解对数据误差的敏感性,它反映 ...

  2. matlabnbsp;sumsqr函数nbsp;求矩阵的平方和

    matlab中的sumsqr函数 求矩阵的平方和 例子 a = 0     1     2     3     4     5 >> sumsqr(a) ans = 55 >> ...

  3. 2.3 matlab矩阵求值(矩阵的行列式值、矩阵的秩、矩阵的积、矩阵的范数和矩阵的条件数)

    1.方阵的行列式 把一个方阵看作一个行列式,并对其按行列式的规则求值,这个值就称为方阵所对应的行列式的值.det(A):求方阵A所对应的行列式的值. >> format rat; > ...

  4. 线性代数知识回顾:矩阵的秩,矩阵的范数,矩阵的条件数,矩阵的特征值和特征向量

    一.矩阵的秩 1.定义: 矩阵线性无关的行数或列数称为矩阵的秩 补充: 线性代数中的线性相关是指: 如果对于向量α1,α2,-,αn, 存在一组不全为0的实数k1.k2.-.kn, 使得:k1·α1+ ...

  5. 数值计算---求希尔伯特矩阵的条件数

    这几天数值计算老师交给我们一个课程设计,计算希尔伯特矩阵的条件数,观察其随维数的变化情况. 下面是程序,主要用到幂法和反幂法. #include <iostream> #include & ...

  6. matlab cond函数

    文章出处:点击打开链接 [功能简介]求矩阵的条件数.矩阵的条件数用于衡量线性方程组的解对数据误差的敏感性,它反映出矩阵求逆及线性方程组解的精确程度. [语法格式] 1.c=cond(X) 求矩阵X的2 ...

  7. 矩阵论(六):矩阵的条件数

    矩阵论专栏:专栏(文章按照顺序排序) 参考资料: 线性代数基础知识系列:1.2.3.4.5 矩阵分解-从Schur分解.特征值分解EVD到奇异值分解SVD(下) 矩阵的正定性 矩阵的条件数用于界定一个 ...

  8. cond--求矩阵的条件数

    [功能简介]求矩阵的条件数.矩阵的条件数用于衡量线性方程组的解对数据误差的敏感性,它反映出矩阵求逆及线性方程组解的精确程度. [语法格式] 1.c=cond(X) 求矩阵X的2-范数的条件数,即X的最 ...

  9. 数值计算·第二集:矩阵的条件数(Matlab版)

    条件数的倒数:rcond(A):A为矩阵,rcond(A)为A的1范数的条件数的倒数的 估计值.如果A的条件数越好,那么其值在1.0附近:反之,则在无穷小附近. %%矩阵的条件数 A = [11,2, ...

最新文章

  1. hibernate一级缓存和二级缓存的区别[转]
  2. Spring之IoC总结帖
  3. 使用Convirt 2.0.1管理虚拟机环境
  4. php搜索文件名,php实现按文件名搜索文件的远程文件查找器
  5. python自动发邮件运行正常就是收不到邮件是为什么_python stmp module 163邮箱发送邮件不成功...
  6. 计算机二级链表,计算机二级c语言上机考试——结构体与链表(3页)-原创力文档...
  7. python整形魔法_python 魔法方法
  8. 【BZOJ2460】元素,线性基+贪心
  9. 计算机云客户端技术指标,云服务器技术指标
  10. 根据 List 对象中的某个字段进行排序
  11. OpenGL--纹理贴图基础
  12. 模块“dm.dll”已加载,但对 DllRegisterServer 的调用失败,错误代码为 0x80020009
  13. 【原】豆瓣电台桌面版插件开发——搜索插件
  14. win10远程桌面连接凭据怎么设置_如何解决Win10远程桌面提示你的凭证不工作?...
  15. 数据可视化专属配色方案
  16. go 语言环境安装 WIMDOWS + LINUX 系统
  17. 服务器raid卡维修,服务器Raid卡故障|服务器RAID卡坏恢复数据
  18. tcl文件调用c语言,TCL与c/c++的互相调用
  19. Oracle sql 对多个sql count的值再求和
  20. sudo进入root权限时,sudo :unable to resolve host me(此为你的Linux登录用户名): Connection timed out.

热门文章

  1. 机器学习算法优缺点及其应用领域
  2. 向量的点乘与叉乘公式
  3. 耦合性、内聚度的介绍
  4. 求弹性模量和泊松比计算题_弹性模量及泊松比计算公式
  5. 1-8记录联合查询与注入逻辑
  6. 吃鸡又一重大更新,只有跑毒才能拯救你
  7. 单片机通过FTP私有协议升级步骤
  8. 仪征警方“猎狐行动”取得重大战果一外逃越南嫌疑人落网
  9. ArcGIS基础实验操作100例--实验60三维模型转体模型
  10. c语言strsep,C/C++ 字符串分割: strtok 与 strsep 函数说明(示例代码)