矩阵论复习笔记:盖尔圆的隔离技巧
1.Gerschgorin定理
对于复数矩阵C={cij}n×nC=\{c_{ij}\}_{n\times n}C={cij}n×n,其矩阵特征值满足:
∣λ−ckk∣≤∑j≠k∣ckj∣∣λ−ckk∣≤∑j≠k∣cjk∣|\lambda-c_{kk}|\leq\sum_{j \neq k}|c_{kj}|\\ |\lambda-c_{kk}|\leq\sum_{j\neq k}|c_{jk}| ∣λ−ckk∣≤j=k∑∣ckj∣∣λ−ckk∣≤j=k∑∣cjk∣
且kkk个圆盘组成的联通区域包含kkk个矩阵特征值。
2.例题
对A=(−i0−11−i7i1020−20−190−6−9)A=\begin{pmatrix} -i&0&-1&1\\-i&7i&1&0\\2&0&-20&-1\\9&0&-6&-9\end{pmatrix}A=⎝⎜⎜⎛−i−i2907i00−11−20−610−1−9⎠⎟⎟⎞矩阵进行隔离其特征值(要求画图显示)。
解:先画出zzz域下的AAA的盖尔圆:
看图发现要隔离联通区域的话首先得缩小G4G_4G4,那么可以设DDD如下:
D=(111a)D = \begin{pmatrix} 1& & & \\ &1& & \\ & &1&\\& & & a\end{pmatrix} D=⎝⎜⎜⎛111a⎠⎟⎟⎞
其中0<a<10<a<10<a<1,那么隔离后的矩阵AAA为:
B=DAD−1=(−i0−11/a−i7i1020−20−1/a9a0−6a−9)B=DAD^{-1}=\begin{pmatrix} -i&0&-1&1/a\\-i&7i&1&0\\2&0&-20&-1/a\\9a&0&-6a&-9\end{pmatrix} B=DAD−1=⎝⎜⎜⎛−i−i29a07i00−11−20−6a1/a0−1/a−9⎠⎟⎟⎞
得到新的盖尔圆如下:
G1∗:∣z+i∣≤1+1aG2∗:∣z−7i∣≤2G3∗:∣z+20∣≤2+1aG4∗:∣z+9∣≤15aG_1^*:|z+i|\leq1+\frac{1}{a}\\ G_2^*:|z-7i|\leq 2\\ G_3^*:|z+20|\leq 2+\frac{1}{a}\\ G_4^*:|z+9|\leq15a G1∗:∣z+i∣≤1+a1G2∗:∣z−7i∣≤2G3∗:∣z+20∣≤2+a1G4∗:∣z+9∣≤15a
重点在于要使得G4∗G_4^*G4∗不与其他三个圆相交,可以得到:
G4∗与G1∗相离:∣9−i∣=82>1+1a+15aG4∗与G2∗相离:∣7i+9∣=130>2+15aG4∗与G3∗相离:∣11∣>2+1a+15aG_4^*与G_1^*相离:|9-i|=\sqrt{82}>1+\frac{1}{a}+15a\\ G_4^*与G_2^*相离:|7i+9|=\sqrt{130}>2+15a\\ G_4^*与G_3^*相离:|11|>2+\frac{1}{a}+15a\\ G4∗与G1∗相离:∣9−i∣=82>1+a1+15aG4∗与G2∗相离:∣7i+9∣=130>2+15aG4∗与G3∗相离:∣11∣>2+a1+15a
首先根据第一个式子发现尝试a=13a=\frac{1}{3}a=31发现条件成立,同时能满足后面两个式子,因此取a=13a=\frac{1}{3}a=31,得到隔离之后的盖尔圆如下:
矩阵论复习笔记:盖尔圆的隔离技巧相关推荐
- 矩阵论复习笔记-华中科技大学
来华科之后又的第一门有挂科风险的考试 也是复习时间时间最长的一门课 昨天考完试 赶紧记下来 再过一天可能就要忘记了 回顾考试涉及到的 填空题 矩阵范数 最小多项式 线性变换的特征值 K积 计算题 奇异 ...
- 矩阵论复习笔记:矩阵直积及其应用
1.矩阵直积及其应用类题目 1.1 常用矩阵直积知识回顾 直积常用性质: 1.A⨂B≠B⨂AA\bigotimes B \neq B\bigotimes AA⨂B=B⨂A 2.(A1+A2)⨂B= ...
- 【矩阵论】特征值的估计(上下界和盖尔圆)
前言:为什么不直接求特征值而是去估计特征值? 当我们遇到的不是书本上的3阶或4阶矩阵,而是高阶矩阵时(如图像中的256×256),我们再使用特征方程det(λI−A)=0\det(\lambda I ...
- python绘制盖尔圆并做特征值的隔离
本程序并非智能到直接运行隔离出所有特征值,而是需要高抬贵手,手动调节变换矩阵D的参数,以实现特征值的隔离.若期待直接找到能特征值隔离的D矩阵参数变化范围,怕足下要失望了,鄙人暂没有做到那一步,一是因为 ...
- 【Matlab 控制】绘制盖尔圆
Matlab 绘制盖尔圆 关于盖尔圆的原理可参考文章: [控制]盖尔圆盘定理 [数理知识]第4章-赋范线性空间与矩阵范数-<矩阵论>方保镕 根据输入矩阵,绘制盖尔圆 clear; clc; ...
- Opencv复习笔记
Opencv复习笔记 整理自:https://blog.csdn.net/u013162035/article/details/79645331 0.索引: [1]图像混合 [2]图像通道的分离和合并 ...
- 大学物理:第10章 静电场 复习笔记
大学物理:第10章 静电场 复习笔记 大学物理 大学物理:第10章 静电场 复习笔记 一级目录 二级目录 三级目录 一.电荷与库仑定律 1.电荷 2.库仑定律 3.静电力的叠加原理 二.电场与电场强度 ...
- 【北航软院+保研复习】计算机网络复习笔记
基于王道计网 1.0版本 北航软院大三上计网课程复习笔记 2.0版本 保研复习笔记 文章目录 第一章 1.1 计算机网络概述 1.1.1 计算机网络的概念 1.1.2 计算机网络的组成 1.1.4 计 ...
- 机原、机设的主干知识_简要复习笔记
声明:本人主要依据北航2022级硕士全国统考机械学院专业科目考试大纲对机原机设进行了知识梳理,仅个人复习使用,不作任何其他用途.如果写的有错误,欢迎指出. 在此之前,本人已对机原机设进行了较为系统的复 ...
最新文章
- python怎么打开文档_python打开怎么运行
- linux中yum源的配置和使用
- 22_Android中的本地音乐播放器和网络音乐播放器的编写,本地视频播放器和网络视频播放器,照相机案例,偷拍案例实现
- python 字符串分割_如何使用python语言split方法对不同字符串分割
- conv2d 公式_Pytorch 从0开始学(6)——Conv2d 详解
- python程序执行时间毫秒_如何使用python解析包含毫秒的时间字符串?
- 在Reporting Services (RDL)中自动生成大量列
- python查找多个关键词,在字符串列表中搜索多个关键字
- 持续集成jenkins工具介绍(一)
- 高淇Struts2.0教程之视频笔记(4)
- XMPP 即时通讯协议介绍
- 基于开源的BPM流程引擎VS天翎自研BPM流程引擎
- 华滋先生:互联网创业,加入社群是有用的吗?
- RF基础(一) RF内建函数库BuiltIn
- box-shadow四个边框设置阴影样式
- 谷歌学术、中国知网生成参考文献
- C#Windows学生信息管理系统
- 与PHP对抗招聘者垃圾邮件-概念证明
- win7系统下共享电脑网络
- TD-LTE获奖:守望5G 引领全球