高等工程数学 —— 第一章 (2)矩阵的谱半径与条件数
第一章 (2)矩阵的谱半径与条件数
文章目录
- 第一章 (2)矩阵的谱半径与条件数
- 谱半径
- 收敛矩阵
- 矩阵级数
- 矩阵幂级数
- 矩阵的条件数及应用
谱半径
- 谱半径其实就是最大特征值
- 注意这里谱半径是小于等于矩阵的任意范数的。在求特征值比较麻烦的时候我们就可以用这条性质来估计谱半径的最大值。
- 当矩阵A为正规矩阵时,AH=AA^H = AAH=A,所以ρ(A)=∣∣A∣∣2\rho(A) = ||A||_2ρ(A)=∣∣A∣∣2。但是要注意谱半径不是矩阵范数,因为他不满足矩阵范数的性质。如下:
例1:
- 估计谱半径其实就是求矩阵范数中最小的。
例2:
例3:
收敛矩阵
- 矩阵A为收敛矩阵的充要条件为ρ(A)<1\rho(A) < 1ρ(A)<1。这里我是这样想的:矩阵的特征值就是某向量在矩阵定义的空间变换后长度的变化倍数。(可以去B站看下线性代数的本质)所以当最大特征值小于1时才收敛。
例:
- 这里A1A_1A1可以求F范数为1825<1\sqrt\frac{18}{25} < 12518<1
- 对于A2A_2A2发现常用的那几个范数都不小于1,所以我们求他的特征值。对于二维矩阵求特征值不必用∣λE−A∣=0|\lambda E - A| = 0∣λE−A∣=0来构建方程。形如:(abcd)\begin{pmatrix} a & b \\ c & d \end{pmatrix}(acbd)求特征值只需让求(λ−a)(λ−d)−bc=0(\lambda - a)(\lambda - d) - bc = 0(λ−a)(λ−d)−bc=0的解就好了。
矩阵级数
矩阵幂级数
- 由ρ(A)<∣∣A∣∣\rho(A) < ||A||ρ(A)<∣∣A∣∣可知,若存在∣∣A∣∣<R||A|| < R∣∣A∣∣<R,那么矩阵幂级数同样绝对收敛。
对于收敛半径RRR的求法如下:
例:
- 这里是将ak=ka_k = kak=k,Ak=15k(2132)kA^k = \frac{1}{5^k} \begin{pmatrix} 2 & 1 \\ 3 & 2 \end{pmatrix}^kAk=5k1(2312)k。我们也可以将ak=k5ka_k = \frac{k}{5^k}ak=5kk,Ak=(2132)kA^k = \begin{pmatrix} 2 & 1 \\ 3 & 2 \end{pmatrix}^kAk=(2312)k,这样的话R=5R = 5R=5、∣∣A∣∣F=18||A||_F = \sqrt{18}∣∣A∣∣F=18,所以R>∣∣A∣∣F>ρ(A)R > ||A||_F > \rho(A)R>∣∣A∣∣F>ρ(A)即绝对收敛
矩阵的条件数及应用
- 这里的范数可以是任意同等类型的范数
例:
特殊的:
例:
对于矩阵条件数一个比较重要的推论:
例:
条件数事实上表示了矩阵计算对于误差的敏感性。对于线性方程组Ax=b,如果A的条件数大,b的微小改变就能引起解x较大的改变,数值稳定性差。如果A的条件数小,b有微小的改变,x的改变也很微小,数值稳定性好。
例:
关于记号o,当x →a时,两个无穷小量α(x)、β(x)之间有记号α(x)=o[β(x)],就是说当x →a时,无穷小量α(x)关于β(x)是高阶无穷小,即当x →a时,α(x)/β(x)→0。
关于记号O,就是当x →a时,f(x) / g(x) 等于一个不为0的常数,就记作f(x) = O[g(x)]。
这里感觉取无穷范数更简单一点:
高等工程数学 —— 第一章 (2)矩阵的谱半径与条件数相关推荐
- 高等工程数学 —— 第一章 (1)距离与范数
前言 研一生活开始了,看了大家对我之前博客的鼓励让我知道写博客是一件多么有意义的事情.写这些让我遇见许多陌生的有缘人,有老骥伏枥的大叔.也有可爱温暖的学妹-- 这里将高等工程数学的笔记留给不爱吃香菜的 ...
- 高等工程数学 —— 第二章 (1) Smith标准型与Jordan标准型
高等工程数学 -- 第二章 (1) Smith标准型与Jordan标准型 文章目录 高等工程数学 -- 第二章 (1) Smith标准型与Jordan标准型 不变因子与Smith标准型 行列式因子 行 ...
- 高等工程数学 —— 第四章 (1)线性方程组的直接解法与广义逆矩阵求解矛盾方程组
高等工程数学 -- 第四章 (1)线性方程组的直接解法与广义逆矩阵求解矛盾方程组 文章目录 高等工程数学 -- 第四章 (1)线性方程组的直接解法与广义逆矩阵求解矛盾方程组 线性方程组的直接解法 Ga ...
- 高等工程数学 —— 第三章(2)奇异值分解和A的加号逆
高等工程数学 -- 第三章(2)奇异值分解和A的加号逆 文章目录 高等工程数学 -- 第三章(2)奇异值分解和A的加号逆 奇异值分解 广义逆矩阵 A + A^{+} A+的直接计算方法 奇异值分解计算 ...
- 【高等工程数学】南理工研究生课程 突击笔记2 谱半径,条件数
高等工程数学 突击笔记2 文章目录 高等工程数学 突击笔记2 写在前面 一.谱半径 什么是正规矩阵 矩阵级数敛散性判断 二.条件数 举个栗子 条件数求误差 总结 写在前面 提示:这里可以添加本文要记录 ...
- 【高等工程数学】南理工研究生课程 突击笔记1 距离与范数1
高等工程数学 突击笔记1 距离与范数1 文章目录 高等工程数学 突击笔记1 距离与范数1 写在开头 一.什么是范数 二.诱导范数 总结 写在开头 高等工程数学这个课真的很恶心,感谢B站UP ASH丶零 ...
- 【高等工程数学】南理工研究生课程 突击笔记3 不变因子与Jordan标准型
高等工程数学 突击笔记3 文章目录 高等工程数学 突击笔记3 一.标准型 λ矩阵 行列式因子D 不变因子d 初等因子 Jordan标准型 二.盖尔圆 特征值隔离 总结 第二章内容大致分成三个部分 标准 ...
- 【高等工程数学】南理工研究生课程 突击笔记4 幂迭代
幂迭代法,和逆幂迭代法 文章目录 幂迭代法,和逆幂迭代法 写在前面 一.幂迭代法 二.逆幂迭代法 三.规范化迭代方式 四.A分解 例 总结 写在前面 承接笔记3,先补一个盖尔圆的题目 如果特征值是复数 ...
- matlab工程数学实践,matlab高等工程数学作业实践报告
matlab高等工程数学作业实践报告 目 录1.非线性方程求根的数值解法房贷年利率---------------12.线性方程组的数值解法配置指定成分合金------------33.估计与检验铝合金 ...
最新文章
- jsp中空格字符怎么写_驰骋职场多年,如何对付Excel中的空格字符?
- glusterfs初步实践一(类似NFS)
- java long 对应mybati类型_修改 mybatis-generator 中数据库类型和 Java 类型的映射关系...
- 撑起12306网站,全靠这个世界第一的缓存框架!
- [HDU 3625] Examining the Rooms(第一类斯特林数)
- 英语语法学习2--句子的成分
- java桥_JAVA 桥模式
- 组合自定义控件的步骤详解
- 文本分类的一种对抗训练方法
- LintCode MySQL 1968. 查询首两个字母在 ‘Db‘ 和 ‘Dy‘ 之间的课程名称(REGEXP正则)
- Django自定义一个简单的中间件,并使用此中间件
- 【Elasticsearch】elasticsearch shard 分片
- 在阿里云主机的Debian操作系统上安装Docker
- 前端调用接口报错看不到报错响应时 console.dir
- vue基础之样式绑定(class,style)
- 有加密狗的软件怎样实现全网络电脑用_加密狗应用领域有哪些?为您揭开7大行业软件保护的奥秘...
- (转载)100+个程序员开发必备参考手册(在线及下载)
- 脑虎科技CEO彭雷专访│脑机未来,生命科学的无尽前沿
- VM虚拟机安装及安装Windows系统
- OpenStack详解、排错、经验总结
热门文章
- 文件is not found in the curren directory or on the Matlab path
- java判断字符串是否为数字、字母、汉字
- 【不专一的开发】UML(二)---行为图(状态图、活动图、序列图、协同图)
- 听肖邦Chopin: Waltz No.19 in A minor, Op.Posth的美好
- 计算机原理学习笔记——第一章 计算机的抽象和技术——2 程序的表象之下
- Win7系统IE浏览器提示网站还原错误怎么解决?
- MACBOOK 刻录系统盘及win7安装在mac上的步骤
- 《跨语言文本相似性检测》第一周—前期调研
- 数字可视化大屏边框制作
- Linux学习笔记-RH135之Apache的管理及优化