MATLAB矩阵求值
目录
方阵的行列式
矩阵的秩
矩阵的迹
矩阵的特征值和特征向量
向量和矩阵的范数
1、向量的3种常用范数
2、矩阵的范数
矩阵的条件数
方阵的行列式
在MATLAB中,求方阵A所对应的行列式的值的函数是det(A)
>> A=rand(4)A =0.8147 0.6324 0.9575 0.95720.9058 0.0975 0.9649 0.48540.1270 0.2785 0.1576 0.80030.9134 0.5469 0.9706 0.1419>> det(A)ans =-0.0261>> A=[1,2;3,4]A =1 23 4>> det(A)ans =-2
矩阵的秩
在MATLAB中,求矩阵秩的函数是rank(A)
>> A=[2,2,3;4,5,-6;7,8,9]A =2 2 34 5 -67 8 9>> rank(A)ans =3>> A=[1,2;1,2]A =1 21 2>> rank(A)ans =1
矩阵的迹
在MATLAB中,求矩阵迹的函数是trace(A)
>> A=[2,2,3;4,5,-6;7,8,9]A =2 2 34 5 -67 8 9>> trace(A)ans =16
矩阵的特征值和特征向量
在MATLAB中,计算矩阵A的特征值和特征向量的函数时eig(A)
>> A=[1,1,0.5;1,1,0.25;0.5,0.25,2];
>> eig(A)ans =-0.01661.48012.5365>> [V,v]=eig(A)V =0.7212 0.4443 0.5315-0.6863 0.5621 0.4615-0.0937 -0.6976 0.7103v =-0.0166 0 00 1.4801 00 0 2.5365>> [V,v]=eig(A,'nobalance')V =0.7212 0.4443 0.5315-0.6863 0.5621 0.4615-0.0937 -0.6976 0.7103v =-0.0166 0 00 1.4801 00 0 2.5365
向量和矩阵的范数
矩阵或向量的范数是用来度量矩阵或向量在某种意义下的长度。范数有多种方法定义,其定义不同,范数值也就不同。
1、向量的3种常用范数
设向量,他的3种范数定义如下
(1)向量1-范数:向量元素的绝对值之和。
(2)向量2-范数:向量元素平方和的平方根。
(2)向量-范数:所以向量元素绝对值中的最大值。
在MATLAB中,求这3种向量范数的函数如下:
norm(V,1)
norm(V,2)
norm(V,inf)
>> V=[-1,1/3,1];>> v1=norm(V,1)v1 =2.3333>> v2=norm(V,2)v2 =1.4530>> vinf=norm(V,inf)vinf =1
2、矩阵的范数
设A是一个mn矩阵,V是一个含有n个元素的列向量,矩阵A的范数可以定义为
,
上式只给出了矩阵范数的基本定义,未给出具体计算方法,完全按照上式是难以计算一个矩阵的某种具体范数的,从属于2种向量范数的矩阵范数计算公式如下。
(1)矩阵A的1-范数:所有矩阵列元素绝对值之和的最大值。
(2)矩阵A的2-范数:A'A矩阵的最大特征根的平方根。
(2)矩阵A的-范数:所有矩阵行元素绝对值之和的最大值。
在MATLAB中,求这3种矩阵范数的函数如下:
norm(A,1)
norm(A,2)
norm(A,inf)
>> A=[1,0,1,0,5;3,5,7,4,6]A =1 0 1 0 53 5 7 4 6>> a1=norm(A,1)a1 =11>> a2=norm(A,2)a2 =12.1738>> ainf=norm(A,inf)ainf =25
矩阵的条件数
在求解线性方程组 AX=b 时,一般认为,系数矩阵 A 中个别元素的微小扰动不会引起起解向量的很大变化。这样的假设在工程应用中非常重要,因为一般系数矩阵的数据是由实验数据获得的,并非精确值,但与精确值误差不大。
由上面的假设可以得出结论:当参与运算的系数与实际精确值误差很小时,所获得的解与问题的准确解误差也很小。
遗憾的是,上述假设并非总是正确的。
对于有的系数矩阵,个别元素的微小扰动会引起解的很大变化,在计算数学中,称这种矩阵为病态矩阵,而称解不因其系数矩阵的微小扰动而发生大的变化的矩阵为良性矩阵。当然,良性与病态是相对的,需要一个参数来描述,条件数就是用来描述矩阵的这种性能的一个参数。
矩阵 A 的条件数等于 A 的范数与A的逆矩阵的范数的乘积,即.这样定义的条件数总是大于1的。条件数越接近于1,矩阵的性能越好,反之,矩阵的性能越差。A 有3种范数,相应地可以定义3种条件数。
在 MATLAB 中,计算 A 的3种条件数的函数如下:
(1) cond(A,1):计算 A 的1-范数下的条件数.(2)cond(A) 或 cond(A,2) :计算A的2-范数下的条件数。
(3)cond(A,inf):计算A的-范数下的条件数。
>> A=[2,2,3;4,5,-6;7,8,9];
>> c1=cond(A)c1 =87.9754>> B=[2,-5,4;1,5,-2;-1,2,4];
>> c2=cond(B)c2 =3.7515矩阵B条件数比矩阵A的条件数更接近与1,因此,矩阵B的性能要好于矩阵A
MATLAB矩阵求值相关推荐
- 2.3 matlab矩阵求值(矩阵的行列式值、矩阵的秩、矩阵的积、矩阵的范数和矩阵的条件数)
1.方阵的行列式 把一个方阵看作一个行列式,并对其按行列式的规则求值,这个值就称为方阵所对应的行列式的值.det(A):求方阵A所对应的行列式的值. >> format rat; > ...
- 中南大学 科学计算与MATLAB语言 11矩阵求值
中南大学 科学计算与MATLAB语言 11矩阵求值 矩阵求值主要包括 矩阵的行列式值 矩阵的秩 矩阵的迹 矩阵的范数 矩阵的条件数 把一个方阵看作一个行列式,并对其按行列式的规则求值,这个值就称方阵所 ...
- matlab - 特殊矩阵、矩阵求值、稀疏矩阵
学习视频 1.特殊矩阵 1.1 通用特殊矩阵 format % 零矩阵(全0) 幺矩阵(全1) 单位矩阵 % zeros ones eye rand(生成0~1的随机元素) randn(生成均值为1, ...
- Matlab多项式求值(2)(roots和poly的使用)
Matlab多项式求值(2)(roots和poly的使用) #Matlab学习笔记 1.多项式求根 n次多项式具有n个根(包括实根和共轭复根) roots函数可以用于求解多项式的全部根,其调用格式为: ...
- java求矩阵_java实现的n*n矩阵求值及求逆矩阵算法示例
本文实例讲述了java实现的n*n矩阵求值及求逆矩阵算法.分享给大家供大家参考,具体如下: 先来看看运行结果: java版的写出来了,用的跟c语言相同的算法,然后看看能不能以后加个框做成程序: imp ...
- java n*n矩阵求值及求逆矩阵
java版的先写出来了,用的跟c语言相同的算法,然后看看能不能以后加个框做成程序: import java.math.*; import java.util.*; import java.text.* ...
- C语言 n*n矩阵求值及求逆矩阵
求出矩阵的值以及输出逆矩阵,英语不好,略拗口. 上代码: #include<stdio.h> #include<string.h> #include<math.h> ...
- MATLAB矩阵的值,迹,秩,范数,上三角矩阵,下三角矩阵,主对角线元素
设A为矩阵 det(A)求矩阵的值 trace(A)求矩阵的迹 rank(A)求矩阵的秩 norm(A)求矩阵的范数 norm(A,1)求矩阵的1范数 norm(A,inf)求矩阵的无穷范数 diag ...
- c语言逆值,C语言 n*n矩阵求值及求逆矩阵
求出矩阵的值以及输出逆矩阵,英语不好,略拗口. 上代码: #include #include #include int map1[110][110],i,j,k,l,m,n,p,vis[110][11 ...
最新文章
- Linux的Nginx四:功能|模块
- Cisco ASA 5500系列防火墙 Series Adaptive Security Appliances
- kd tree学习笔记 (最近邻域查询)
- 常用的色彩名、色彩值的对照表
- SAP License:SAP用户权限
- mysql卡住如何定位_MySQL 5.7中如何定位DDL被阻塞的问题
- 2012.02.13(rtsp)
- 计算机十个小技巧,Win10的10个实用小技巧,电脑小白轻松掌握,大大提升工作效率...
- 基于RFID定位技术的文物仓库管理--新导智能
- 两种操作botton的方法
- php惠新宸框架,首位国人惠新宸加入PHP官方开发组
- [画板]画PCB和玩植物大战僵尸
- Kaggle账号注册时验证码无法显示问题解决方法
- JSP页面查询显示常用模式 (附源代码)
- 经典网络架构学习-Transformer
- 解决英伟达CUDA和cuDNN下载过慢的问题
- oracle练习题(二)
- 牛津英语字典pdf下载_除了long time no see,你知道还有这些中式英语也进入了牛津字典吗...
- 中国大学MOOC C语言程序设计(大连理工大学) 课后编程题 第三周题解(个人向仅供参考)
- Liferay中的workflow
热门文章
- 一个计算器程序,要求用户输入两个整数和一个运算符,程序能够计算出两个数的相应加减乘除结果。(涉及到精确度的,一律保留两位小数,使用输出格式)
- Boosted Trees 介绍
- 剪映,最简单的PC端视频处理工具,大厂出品必属精品~
- 无限循环小数化分数、
- SQL习题“部门工资最高的员工”
- Linux云计算之shell 基本语法晚上干肺写出来的
- 凶手(算法 假设法)
- 微信小程序 TypeError: _typeof2 is not a function
- 手把手教你实现buffer(三)——接口及自动扩容
- ROS学习|SMACH状态机三(使用实例)