2.3 matlab矩阵求值(矩阵的行列式值、矩阵的秩、矩阵的积、矩阵的范数和矩阵的条件数)
1、方阵的行列式
把一个方阵看作一个行列式,并对其按行列式的规则求值,这个值就称为方阵所对应的行列式的值。det(A):求方阵A所对应的行列式的值。
>> format rat;
>> a = [1 2 3;-3 -5 4;-5 9 1]a =1 2 3 -3 -5 4 -5 9 1 >> det(inv(a))ans =-1/231 >> 1/det(a)ans =-1/231
2、矩阵的秩
for n = 3:20r(n) = rank(magic(n));
end
bar(r)
grid on
axis([2,21,0,20])
n=3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
r =3 3 5 5 7 3 9 7 11 3 13 9 15 3 17 11 19 3
①奇数阶魔方阵秩为η,即奇数阶魔方阵是满秩矩阵;②重偶数阶魔方阵秩为n2+2(n是2的倍数,但非4的倍数);③双重偶数阶魔方阵秩均为3(阶数是4的倍数)。
3、矩阵的迹
矩阵的迹等于矩阵的对角线元素之和,也等于矩阵的特征值之和。trace(A):求矩阵A的迹。
>> a = [1 2 3;-3 -5 4;-5 9 1]a =1 2 3-3 -5 4-5 9 1>> b = trace(a)b =-3>> b = sum(diag(a)) *可先提取矩阵的对角线元素,再进行求和b =-3
4、矩阵和向量的范数
(1)向量的3种范数
①向量1—范数:向量元素的绝对值之和。
②向量2—范数:向量元素绝对值的平方和的平方根。
③向量∞—范数:所有向量元素绝对值中的最大值。
(2)矩阵的3种范数
①矩阵A的1—范数:矩阵列元素绝对值之和的最大值。
②矩阵A的2—范数:A’A矩阵的最大特征值的平方根。
③矩阵A的oo—范数:所有矩阵行元素绝对值之和的最大值。
矩阵和向量的范数求法一致:
①norm(V)或norm(V,2):计算2—范数。
②norm(V,1): 计算1—范数。
③norm(V,inf):计算∞—范数。
例
>> a = [1 2 3;-3 -5 4;-5 9 1]a =1 2 3-3 -5 4-5 9 1>> b = norm(a,1) *矩阵1范数b =16>> b = norm(a) *矩阵2范数b =10.9225>> b = norm(a,inf) *矩阵∞范数b =15
5、矩阵的条件数
矩阵A的条件数等于A的范数与A的逆矩阵的范数的乘积。
条件数越接近于1,矩阵的性能越好,反之,矩阵的性能越差。
在MATLAB中,计算矩阵A的3种条件数的函数是:
①cond(A,1):计算A的1—范数下的条件数。
②cond(A)或cond(A,2):计算A的2—范数下的条件数。
③cond(A,inf):计算A的○o—范数下的条件数。
for n = 2:10r(n) = cond(hilb(n));
end
format long
r'
ans =1.0e+13 *00.0000000000019280.0000000000524060.0000000015513740.0000000476607250.0000014951058640.0000475367356910.0015257575566630.0493153404551011.6025028168113180.0000000000011000.0000000000003000.0000000000013000.0000000000009000.0000000000015000.0000000000003000.0000000000017000.0000000000011000.0000000000019000.000000000000300
希尔伯特矩阵的条件数不断增大,矩阵性能变差
2.3 matlab矩阵求值(矩阵的行列式值、矩阵的秩、矩阵的积、矩阵的范数和矩阵的条件数)相关推荐
- python矩阵求导代码_搞点枯燥的公式推导:深度学习中的矩阵求导
最近看到一篇极视角转的文章算法推导核心!一次性梳理清楚,是时候搞定矩阵求导了!,想到前不久面试爱奇艺的时候一面的小哥一直让我手推全连接的公式推导,以及用Python+Numpy将过程实现,感觉自己对深 ...
- MATLAB中求开环传递函数的幅值裕度、相位裕度、截止频率的margin()函数用法
一.幅值裕度.相位裕度 相位裕度γ就是Bode图中幅值为0dB时对应的相位加上180°,如图中0dB时的相位是-99.1°,所以该系统的相位裕度就是 -99.1°+180°=80.9°.所谓相位裕度的 ...
- 一个F范数对矩阵求导例子
==PART1 == 1. 问题 如何求下列表达式中的未知参数W?求偏导? minf(W)=minw∥(XW∘D˘)B∥F2(0)\min f(W)=\min_w \|(XW \circ \bre ...
- 矩阵求导与BP的证明的建议
前言 在有些博客推导神经网络的BP时,涉及到多次矩阵求导运算,尤其是反向传播时候,求的梯度结果被转置了,比如假设最后一层的输出为 y=σ(w⋅x+b) y=\sigma\left(w\cdot x+b ...
- 算法推导核心!一次性梳理清楚,是时候搞定矩阵求导了!附参考资料
点击上方"Datawhale",选择"星标"公众号 第一时间获取价值内容 1. 前言与感谢 以前就对矩阵求导有心理阴影,到底要不要转置,哪个矩阵放在前面,哪个矩 ...
- 线性代数之 矩阵求导(3)标量、向量求导的快速记忆
线性代数之 矩阵求导(3)标量.向量求导的快速记忆 前言 基本约定 引例 标量对标量求导 标量对向量求导 向量对向量求导 包含两个变量的求导 总结 扩展 前言 上一次记录了矩阵求导的基本法则和公式,并 ...
- 中南大学 科学计算与MATLAB语言 11矩阵求值
中南大学 科学计算与MATLAB语言 11矩阵求值 矩阵求值主要包括 矩阵的行列式值 矩阵的秩 矩阵的迹 矩阵的范数 矩阵的条件数 把一个方阵看作一个行列式,并对其按行列式的规则求值,这个值就称方阵所 ...
- MATLAB矩阵求值
目录 方阵的行列式 矩阵的秩 矩阵的迹 矩阵的特征值和特征向量 向量和矩阵的范数 1.向量的3种常用范数 2.矩阵的范数 矩阵的条件数 方阵的行列式 在MATLAB中,求方阵A所对应的行列式的值的函数 ...
- matlab - 特殊矩阵、矩阵求值、稀疏矩阵
学习视频 1.特殊矩阵 1.1 通用特殊矩阵 format % 零矩阵(全0) 幺矩阵(全1) 单位矩阵 % zeros ones eye rand(生成0~1的随机元素) randn(生成均值为1, ...
最新文章
- 百度搜索引擎优化指南3.0_深圳网站搜索引擎排名优化电话,百度优化排名费用_华阳网络...
- android的AIDL的调用
- win7系统如何添加环境变量
- 微信小程序WebSocket使用案例
- 企业级负载均衡LVS集群——DR模式下的(加权)轮询调度器、DR模式下的健康检测(ldirectord)
- SSM框架医院信息管理系统-患者医生考勤管理挂号管理-病房管理-科室管理- (idea开发javaweb-php-asp.netC#-j2ee-springboot)
- QQ客户端管理聊天记录漫游功能 | 解决经常弹窗消息漫游安全验证,QQ服务独立密码的问题
- 使用 Python 脚本执行国密 sm2 加解密
- 九连环解法( 基于递归 )
- 群论(Burnside引理和Polya定理)
- 范数(简单的理解)、范数的用途、什么是范数
- Win Server 2012发布!回顾Server发展史
- Unity声音 Sound (圣典篇)
- 什么是微信防火墙_【听课】第4节 什么是“合并单元”(MU)?智能站中的合并单元有什么作用?...
- android项目重构经验总结
- matlab运行后没反映,这个matlab程序运行没反应是怎么回事?在线急等
- 网安基础学习之“文件上传漏洞原理与实现”
- BIM Revit2016下载安装教程
- halcon sobel 边缘检测 sobel_dir
- 应用程序接口(API)安全的入门指南
热门文章
- 黑马程序员C++ 第三阶段 提高篇(1)
- 使用ffmpeg搭建HLS直播系统
- 今年Gartner超融合软件魔力象限有点“魔性”啊
- vcenter 7.0 安装VMware-vRealize-Log-Insight-8.1.1.0
- 141 分割字符串的方法
- WPF工程出现编译错误:error CS0426: The type name 'Views' does not exist in the type 'Calibration'
- 协同过滤推荐算法和基于内容推荐算法的区别?
- 如何清除FTP自动登录状态
- 电机震动噪声(NVH)入门笔记
- Xen虚拟化架构简介