【线性代数与矩阵论】矩阵的谱半径与条件数

矩阵的谱半径与条件数

2023年11月18日


文章目录

  • 矩阵的谱半径与条件数
    • 1. 矩阵的谱半径
    • 2. 谱半径与范数的关系
    • 3. 矩阵的条件数
    • 下链


1. 矩阵的谱半径

定义 A ∈ C n × n {A\in \mathbb C^{n \times n} } ACn×n λ 1 , λ 2 , ⋯ , λ n { \lambda_1,\lambda_2,\cdots,\lambda_n} λ1,λ2,,λn 是A的特征值,则称
ρ ( A ) = max ⁡ 1 ≤ i ≤ n ∣ λ i ∣ \rho(A)=\max_{1\le i\le n}| \lambda_i| ρ(A)=1inmaxλi
为矩阵 A {A} A 的谱半径。矩阵的谱指的是一个矩阵的特征值的集合。

定理 A ∈ C n × n {A\in \mathbb C^{n \times n} } ACn×n ,则

  1. ρ ( A H ) = ρ ( A T ) = ρ ( A ) {\rho(A^ \mathrm H)=\rho(A^ \mathrm T)=\rho(A)} ρ(AH)=ρ(AT)=ρ(A)
  2. ρ ( A k ) = [ ρ ( A ) ] k {\rho(A^k)=[\rho(A)]^k} ρ(Ak)=[ρ(A)]k
  3. A {A} A正规矩阵时, ρ ( A ) = ∣ ∣ A ∣ ∣ 2 {\rho(A)= || A ||_{2 }} ρ(A)=∣∣A2

2. 谱半径与范数的关系

A ∈ C n × n {A\in \mathbb C^{n \times n}} ACn×n λ {\lambda} λ A {A} A 的一个特征值, x {x} x A {A} A 属于 λ { \lambda} λ 的特征向量,则
A x = λ x Ax= \lambda x Ax=λx
C n × n { \mathbb C^{n \times n} } Cn×n 中任一矩阵范数 ∣ ∣ ⋅ ∣ ∣ m { || \cdot ||_{m }} ∣∣m 以及与它相容的向量范数 ∣ ∣ ⋅ ∣ ∣ v { || \cdot ||_{v }} ∣∣v ,有
∣ λ ∣ ⋅ ∣ ∣ x ∣ ∣ v = ∣ ∣ λ x ∣ ∣ v = ∣ ∣ A x ∣ ∣ v ≤ ∣ ∣ A ∣ ∣ ⋅ ∣ ∣ x ∣ ∣ v | \lambda| \cdot || x ||_{v }= || \lambda x ||_{ v}= || Ax ||_{ v} \le || A ||_{ } \cdot || x ||_{v } λ∣∣xv=∣∣λxv=∣∣Axv∣∣A∣∣xv
从而
∣ λ ∣ ≤ ∣ ∣ A ∣ ∣ | \lambda| \le || A ||_{ } λ∣∣A
于是有如下定理。

定理 ρ ( A ) ≤ ∣ ∣ A ∣ ∣ { \rho(A)\le || A ||_{ }} ρ(A)∣∣A ,其中 ∣ ∣ A ∣ ∣ { || A ||_{ }} ∣∣A A {A} A 的任一矩阵范数。
定理 A ∈ C n × n {A\in \mathbb C^{n \times n} } ACn×n ,则任意 ϵ > 0 { \epsilon>0} ϵ>0 ,必存在 C n × n { \mathbb C^{n \times n} } Cn×n 上矩阵范数 ∣ ∣ ⋅ ∣ ∣ m { || \cdot ||_{m }} ∣∣m ,使得
∣ ∣ A ∣ ∣ m ≤ ρ ( A ) + ϵ || A ||_{ m} \le \rho(A)+ \epsilon ∣∣Amρ(A)+ϵ
也就是虽然矩阵范数可能大于谱近谱半径的矩阵半径,却又总是存在无限接范数。数值分析中,谱半径可以认为是算子范数。


3. 矩阵的条件数

引理 P ∈ C n × n {P\in \mathbb C^{n \times n} } PCn×n ,若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ P ∣ ∣ < 1 { || P ||_{ }<1} ∣∣P<1 ,则 I − P {I-P} IP 可逆。
定理 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, δ A ∈ C n × n { \delta A\in \mathbb C^{n \times n} } δACn×n 。若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ A − 1 δ A ∣ ∣ < 1 { || A^{-1} \delta A ||_{ }<1} ∣∣A1δA<1 ,则

  1. A + δ A {A+ \delta A} A+δA 可逆
  2. ∣ ∣ ( A + δ A ) − 1 ∣ ∣ ≤ ∣ ∣ A − 1 ∣ ∣ 1 − ∣ ∣ A − 1 δ A ∣ ∣ {|| (A+ \delta A)^{-1} ||_{ }\le \frac{|| A^{-1} ||_{ }}{1- || A^{-1} \delta A ||_{ }}} ∣∣(A+δA)11∣∣A1δA∣∣A1
  3. ∣ ∣ A − 1 − ( A + δ A ) − 1 ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ≤ ∣ ∣ A − 1 δ A ∣ ∣ 1 − ∣ ∣ A − 1 δ A ∣ ∣ {\frac{||A^{-1}- (A+ \delta A)^{-1} ||_{ }}{|| A^{-1} ||_{ }} \le \frac{|| A^{-1} \delta A ||_{ }}{1- || A^{-1} \delta A ||_{ }}} ∣∣A1∣∣A1(A+δA)11∣∣A1δA∣∣A1δA 矩阵扰动后逆矩阵的相对误差小于右端式子

定义 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数,称
cond ( A ) = ∣ ∣ A ∣ ∣ ⋅ ∣ ∣ A − 1 ∣ ∣ \text{cond}(A)= || A ||_{ } \cdot || A^{-1} ||_{ } cond(A)=∣∣A∣∣A1
为矩阵 A {A} A条件数

推论 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, δ A ∈ C n × n { \delta A\in \mathbb C^{n \times n} } δACn×n 。若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ A − 1 ∣ ∣ ⋅ ∣ ∣ δ A ∣ ∣ < 1 { || A^{-1}|| \cdot || \delta A ||_{ }<1} ∣∣A1∣∣∣∣δA<1 ,则
∣ ∣ A − 1 − ( A + δ A ) − 1 ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ≤ ∣ ∣ A ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ 1 − ∣ ∣ A ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ = cond ( A ) ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ 1 − cond ( A ) ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ \frac{||A^{-1}- (A+ \delta A)^{-1} ||_{ }}{|| A^{-1} ||_{ }} \le \frac{|| A ||_{ }|| A^{-1} ||_{ }\frac{|| \delta A ||_{ }}{|| A ||_{ }}}{1-|| A ||_{ }|| A^{-1} ||_{ }\frac{|| \delta A ||_{ }}{|| A ||_{ }}}= \frac{\text{cond}(A) \frac{|| \delta A ||_{ }}{|| A ||_{ }}}{1-\text{cond}(A) \frac{|| \delta A ||_{ }}{|| A ||_{ }}} ∣∣A1∣∣A1(A+δA)11∣∣A∣∣A1∣∣A∣∣δA∣∣A∣∣A1∣∣A∣∣δA=1cond(A)∣∣A∣∣δAcond(A)∣∣A∣∣δA

定理 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, δ A ∈ C n × n { \delta A\in \mathbb C^{n \times n} } δACn×n b {b} b δ b ∈ C n { \delta b\in \mathbb C^n} δbCn 。若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ A − 1 ∣ ∣ ⋅ ∣ ∣ δ A ∣ ∣ < 1 { || A^{-1}|| \cdot || \delta A ||_{ }<1} ∣∣A1∣∣∣∣δA<1 ,则非齐次线性方程组
A x = b 与 ( A + δ A ) ( x + δ x ) = b + δ b Ax=b \,\,\, 与 \,\,\, (A+ \delta A)(x+ \delta x)=b+ \delta b Ax=b(A+δA)(x+δx)=b+δb
的解满足
∣ ∣ δ x ∣ ∣ v ∣ ∣ x ∣ ∣ v ≤ cond ( A ) 1 − cond ( A ) ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ ( ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ + ∣ ∣ δ b ∣ ∣ v ∣ ∣ b ∣ ∣ v ) \frac{|| \delta x ||_{ v}}{|| x ||_{v }}\le \frac{\text{cond}(A) }{1-\text{cond}(A) \frac{|| \delta A ||_{ }}{|| A ||_{ }}} \bigg( \frac{|| \delta A ||_{ }}{|| A ||_{ }}+\frac{|| \delta b ||_{ v}}{||b ||_{ v}} \bigg) ∣∣xv∣∣δxv1cond(A)∣∣A∣∣δAcond(A)(∣∣A∣∣δA+∣∣bv∣∣δbv)
其中 ∣ ∣ ⋅ ∣ ∣ v { || \cdot ||_{v }} ∣∣v C n { \mathbb C^n} Cn 上与矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ 相容的向量范数。
根据函数
f ( x ) = x 1 − x f(x)= \frac{x}{1-x} f(x)=1xx
的图像,当 cond ( A ) { \text{cond}(A)} cond(A) 很大,则说求逆或求解线性方程组是病态的。


下链


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/638451.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

第11章_常用类和基础API拓展练习(字符串相关练习,日期时间API练习,比较器练习,其它API练习)

文章目录 第11章_常用类和基础API拓展练习字符串相关练习1&#xff1a;阅读题1、length说明2、阅读代码&#xff0c;分析结果3、阅读代码&#xff0c;分析结果4、阅读代码&#xff0c;分析结果5、阅读代码&#xff0c;分析结果6、阅读代码&#xff0c;分析结果7、阅读代码&…

【Linux】Linux基本操作(二):rm rmdir man cp mv cat echo

承接上文&#xff1a; 【【Linux】Linux基本操作&#xff08;一&#xff09;&#xff1a;初识操作系统、ls、cd、touch、mkdir、pwd 】 目录 1.rmdir指令 && rm 指令&#xff1a; rmdir -p #当子目录被删除后如果父目录也变成空目录的话&#xff0c;就连带父目录一…

Elasticsearch 字段更新机制

目录 不支持原地更新 更新过程 段合并&#xff08;Segment Merge&#xff09; 结论 不支持原地更新 Elasticsearch 不支持原地更新&#xff08;in-place update&#xff09;索引中的单个字段。由于 Elasticsearch 使用了不可变的倒排索引&#xff0c;一旦文档被写入&#x…

崩溃了!我说用attach进行问题定位,面试官问我原理

Arthas&#xff08;阿尔萨斯&#xff09;是一款开源的Java诊断和监控工具&#xff0c;可以在生产环境中进行实时的应用程序分析和故障排查。Arthas的实现原理主要基于Java Instrumentation API和Java Agent技术。 Java Agent 是 Java 编程语言提供的一种特殊机制&#xff0c;允…

多模态是什么意思,在生活工业中有哪些应用?

问题描述&#xff1a;多模态是什么意思&#xff0c;在生活工业中有哪些应用&#xff1f; 问题解答&#xff1a; 多模态&#xff08;Multimodal&#xff09;指的是同时利用多种不同模式或传感器的情境、系统或技术。在不同领域&#xff0c;多模态可以涉及到多种感官或信息源的…

第4周:综合应用和实战项目 Day 25-27: 模型调优和优化学习高级技巧

第4周&#xff1a;综合应用和实战项目 Day 25-27: 模型调优和优化学习高级技巧 在这个阶段&#xff0c;我们将专注于提高模型的性能&#xff0c;通过使用高级技巧如正则化、dropout、批标准化等。这些技术对于防止过拟合和提高模型的泛化能力非常重要。 重点学习内容&#xff…

springboot(ssm电子产品销售网站 在线购物商城系

springboot(ssm电子产品销售网站 在线购物商城系统Java系统 开发语言&#xff1a;Java 框架&#xff1a;ssm/springboot vue JDK版本&#xff1a;JDK1.8&#xff08;或11&#xff09; 服务器&#xff1a;tomcat 数据库&#xff1a;mysql 5.7&#xff08;或8.0&#xff09;…

3.php开发-个人博客项目输入输出类留言板访问IPUA头来源

目录 知识点 : 输入输出 配置环境时&#xff1a; 搜索框&#xff1a; 留言板&#xff1a; 留言板的显示&#xff08;html&#xff09;&#xff1a; php代码显示提交的留言&#xff1a; 写入数据库 对留言内容进行显示&#xff1a; php全局变量-$_SERVER 检测来源 墨…

AMEYA360:航顺HK32AUTO39A—适用车载娱乐系统优化方案

车载娱乐系统&#xff0c;顾名思义最开始其功能主要是为驾驶员和乘客提供娱乐体验&#xff0c;但现在智能座舱时代到来&#xff0c;车载信息娱乐系统已成为“人-车-环境”交互的重要载体&#xff0c;是除了驾驶以外的其他任务的中心。 从功能模块上看&#xff0c;车载娱乐系统从…

【算法练习Day50】下一个更大元素II接雨水

​&#x1f4dd;个人主页&#xff1a;Sherry的成长之路 &#x1f3e0;学习社区&#xff1a;Sherry的成长之路&#xff08;个人社区&#xff09; &#x1f4d6;专栏链接&#xff1a;练题 &#x1f3af;长路漫漫浩浩&#xff0c;万事皆有期待 文章目录 下一个更大元素II接雨水单调…

批量终结ascp linux系统中,批量停止kill进程

假如我想批量kill ascp进程 #!/bin/bash# Get PIDs of all running ascp processes pids$(ps -ef | grep [a]scp | awk {print $2})# Loop through each PID and try to terminate the process gracefully for pid in $pids; doecho "Attempting to gracefully terminat…

Apache Wicket 9.10.0发布

Apache Wicket 是一个开源的面向 Java 组件的 Web 应用框架&#xff0c;为政府、商店、大学、城市、银行、电子邮件提供商等成千上万的 Web 应用和网站提供支持。 Wicket 的开发中使用了语义版本&#xff0c;因此与 9.0.0 相比&#xff0c;9.10 版本没有出现 API 中断。 9.10…

9.4 Lambda表达式

9.4 Lambda表达式 1 Lambda语法2. 基于Lambda实现函数式编程3. Stream流式处理 1 Lambda语法 2. 基于Lambda实现函数式编程 3. Stream流式处理

常见的点云下载地址/点云集合/点云库30个

1、pcl库自带教程所需的点云数据下载地址 Point Cloud Library - Browse /PCD datasets at SourceForge.net 另外一个地址 https://github.com/PointCloudLibrary/data 31、新增&#xff1a;VTK库的点云数据下载地址&#xff0c;包含.vtk\ply\pdb\vtp\tif\stl等格式的点云数据 …

【C++】结构体

目录 1.结构体基本概念 2.结构体的定义和使用 3.结构体数组 4.结构体指针 5.结构体嵌套结构体 6.结构体做函数参数 7.结构体中const使用场景 1.结构体基本概念 结构体属于用户 ---- 自定义的数据类型&#xff0c;允许用户储存不同的数据类型 2.结构体的定义和使用 语法…

【网站项目】基于SSM的273校园二手交易网站

&#x1f64a;作者简介&#xff1a;多年一线开发工作经验&#xff0c;分享技术代码帮助学生学习&#xff0c;独立完成自己的项目或者毕业设计。 代码可以私聊博主获取。&#x1f339;赠送计算机毕业设计600个选题excel文件&#xff0c;帮助大学选题。赠送开题报告模板&#xff…

cesium实现动态围栏

项目中使用到了cesium,需要实现动态的围栏的效果&#xff0c; 在网上也找了好多案例&#xff0c;通过着色器来实现效果&#xff0c;为此也有好多博主也附上了自己的代码&#xff0c;也许是因为使用方法不同&#xff0c;复制代码并修改依旧还是没有通过他们的方式实现效果【着色…

STM32标准库开发——PWM驱动代码

PWM驱动初始化代码 使能定时器二时钟 RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2,ENABLE);设置定时器时钟源 TIM_InternalClockConfig(TIM2);配置定时器二的时基单元 TIM_TimeBaseInitTypeDef TIM_TimeBaseInitStruct; TIM_TimeBaseInitStruct.TIM_ClockDivisionTIM_CKD_D…

vmware 安装Rocky-9.3系统

安装系统截图 安装完成&#xff0c;启动 查看版本和内核 开启远程登陆授权 1、编辑配置文件 #提升权限&#xff0c;输入su,并输入密码 su #编辑ssh文件开启root远程登陆 vi /etc/ssh/sshd_config找到以下内容&#xff1a;#PermitRootLogin prohibit-password 添加&#xff1a…

C语言——小细节和小知识11

一、自幂数 1、介绍 自幂数&#xff0c;也被称为阿姆斯特朗数&#xff0c;是一种特殊类型的数&#xff0c;在数学上具有一个有趣的性质&#xff1a;一个 n 位的正整数&#xff0c;其各个位上的数字的 n 次幂之和等于它本身。 这里是自幂数的定义步骤&#xff1a; 确定位数 (…