matlab对称正定矩阵,对称阵与正定矩阵 · Issue #18 · cyfile/Matlab-base-toolbox · GitHub...
正定性和主子式的关系(借助正定阵的主元特性进行说明)
只保留A的前m行和前m列,得到的新矩阵叫A的m阶主子矩阵,其行列式叫主子式.
以前我有一个错误的思路说明这个问题:
利用矩阵A的特征值分解,通过分块矩阵的方式表示矩阵的主子阵,
再利用A的特征值,表示主子式的值
因为有 A=P*D*inv(P)
去掉A的第k行和第k列,相当于去掉P的第k行的同时去掉D的第k行,第k列.
进一步对A同时去掉行号列号相同的多行多列,相当于对V,P做相应的销行销列变化,
设新的A,D,P,分别表示为B,C,W
则有 B=W*C*W'.
进而有 |B|=|W|*|C|*|W'|=|W|^2*|C|
上式说明A的剩余矩阵B的行列式的符号完全和D中剩余特征值的积(|C|)的符号相同.
错误原因:实际上无法通过分块矩阵表示主子阵
正确说明思路
利用矩阵A的主元分解,同样利用分块的方式表示矩阵的主子阵.
因为 x'Ax = x'C*BAB'*C'x =x'CFC'x,中的B,C都是三角阵,
在取B,C的主子阵后,进行运算可以正确得到A的主子阵.
即,通过可以通过对主元分解的各个矩阵进行分块从而正确表示A的主子阵
再利用主元值表示主子式的值.即可说明主子式和矩阵正定的关系.
但是直接这样说明比较麻烦.
可以直接应用主元公式(矩阵A的第n个主元为相邻的主子式之商 |An|/|An-1|)
因为正定阵主元都大于0,
所以 A0, |An|/|An-1| (||是行列式,不是绝对值)都大于0
即各阶主子式 A0 到 |An| 都大于0
反之 如果各阶主子式 A0 到 |An| 都大于0
也可以得到 A0, |An|/|An-1| (||是行列式,不是绝对值)都大于0
即 A的各个主元也都大于0 ,A是正定的
A的各阶主子式大于0 ⇔ A为正定矩阵
A的奇数阶主子式小于0,偶数阶主子式大于0 A的所有特征值小于0⇔ A为负定矩阵
其他思路
因为A正定,所以无论x如何取值都有x'Ax >0
令向量x的某些位置恒为0,上式依然满足
再进行x'Ax运算时,可以看到A中与x非零位置进行运算的部分,就是A的主子式部分
即A的任何主子阵都正定
这个思路只能说明必要性(A正定推主子矩阵正定),不易说明充分性
因为顺序主子阵和原矩阵A共用相同主元.所以主子阵继承A的(正)定性.
同理因为A的主元都能在顺序主子阵的主元中找到,所以顺序主子阵都是正定的话,A正定
这个思路只能说明顺序主子阵
matlab对称正定矩阵,对称阵与正定矩阵 · Issue #18 · cyfile/Matlab-base-toolbox · GitHub...相关推荐
- 频域波束形成matlab,SONAR_code 波束形成技术,包括相移 ,频域 matlab代码,以及相关实验的ppt 263万源代码下载- www.pudn.com...
文件名称: SONAR_code下载 收藏√ [ 5 4 3 2 1 ] 开发工具: matlab 文件大小: 25848 KB 上传时间: 2015-05-03 下载次数: 21 提 供 ...
- 通过有限差分和matlab矩阵运算直接求解一维薛定谔方程,通过有限差分和MATLAB矩阵运算直接求解一维薛定谔方程.doc...
通过有限差分和MATLAB矩阵运算直接求解一维薛定谔方程 通过有限差分和MATLAB矩阵运算直接求解一维薛定谔方程 文章编号:1672-8785(2010)03-0042-05 通过有限差分和MATL ...
- 【 MATLAB 】两个序列的卷积和运算的MATLAB实现(2)
已知下面两个序列: 求这两个序列的卷积. 求卷积的函数是conv,但是使用这个函数有个问题,就是下标问题,也就是求卷积之后的元素值的位置.因此,我们必须要定一个起始点和一个结束点. 方法: 是两个有限 ...
- matlab求二元函数极值算法_最优化计算与matlab实现(3)——进退法
参考资料 <精通MATLAB最优化计算(第二版)> 数值实现 Matlab 2019a 目录 石中居士:最优化计算与Matlab实现--目录zhuanlan.zhihu.com 进退法 ...
- matlab中的libsvm怎么录入数据啊,LibSVM在MATLAB中使用时的几个问题
在科研中需要用到支持向量机(Support Vector Machines, SVM)来进行分类,而目前比较成熟的用于实现SVM的软件包则首推LibSVM.LibSVM目前的版本已经能直接在MATLA ...
- matlab数学建模可应用到第几章,《MATLAB在数学建模中的应用(第2版)》
<MATLAB在数学建模中的应用(第2版)> 卓金武 (编者)) 基本信息 •出版社: 北京航空航天大学出版社; 第2版 (2014年9月1日) •丛书名: MATLAB开发实例系列图书 ...
- matlab中函数或变量无法识别怎么办_用MATLAB巧解微分方程实例分析
点"考研竞赛数学"↑可每天"涨姿势"哦! MATLAB巧解微分方程实例分析 王少华 西安电子科技大学 微分方程求解难, 字母一堆看着烦. 写错数字一时爽, 一直 ...
- github issue 搜索_回顾 2020 年 GitHub 的大事件,你知道多少?
作者:HelloGitHub-小鱼干 这里是 HelloGitHub 出品的年度盘点系列,本期我们将盘点 GitHub 在 2020 发生的大事件,回顾一下今年 GitHub 给我们带来了那些惊喜.故 ...
- matlab使用_重磅!哈工大、哈工程无法使用 MATLAB 软件
近日,有消息称:哈工大.哈工程的师生们已经无法使用 MATLAB 软件了. 随后,有哈工大的学生表示收到了正版MATLAB 软件取消激活的通知,学校在与 MATLAB 软件的开发公司 MathWork ...
最新文章
- Chapter 2 Open Book——8
- JAVA_OA(八):springMVC对JDBC的操作小项目b
- 错误计算机怎么打开,电脑开机出错怎么回事
- shell 进入hadoop_Hadoop关于HDFS的基本操作(Shell命令)
- 独立游戏:四人团队如何100天做到百万用户
- 通过javascript 函数的调用方式来理解this
- java过滤器命名_java 过滤指定后缀文件
- 《深入浅出数据分析》读后感
- ofd文件怎么复制文字
- python数字转中文字符_python数字金额转换为中文大写金额
- python蓝屏代码_死机、卡顿、蓝屏,Python部分的老江湖告诉我的一些超等反常代码...
- 鸟哥的linux私房菜学习笔记《三十七》deamon与服务
- 浅谈屏幕适配 dp dip sp dpi ppi px sp
- 机器学习基础——损失函数与风险函数
- 车联网真要来了?全球汽车网络安全市场将实现跨越式增长
- 关于经典蓝牙和低功耗蓝牙的区别
- 荣峰科技携手互动未来开启“行销宝”战略合作
- 保姆级教程 | Java 8 安装及环境变量配置
- Python简单实现ATM自动存取款机
- 通用算法 - [树结构] - 红黑树