首页 > 其他分享 >【Matlab学习2.3】矩阵求值

【Matlab学习2.3】矩阵求值

时间:2023-01-31 20:57:40浏览次数:42  
标签:end 矩阵 2.3 Vmatrix Matlab 求值 范数 向量

方阵的行列式值

把一个方阵看作一个行列式,并对其按行列式的规则求值,这个值就称为所对应的行列式的值。

det(A):求方阵A所对应的行列式的值。 

例2.3.1:

验证det(A-1)=1/det(A)。

>> format rat
>> A = [1 3 2; -3 2 1; 4 1 2]
A =
       1              3              2       
      -3              2              1       
       4              1              2       
>> det(inv(A))
ans =
       1/11    
>> 1/det(A)
ans =
       1/11    

矩阵的秩

矩阵线性无关的行数或列数称为矩阵的秩。

rank(A):求矩阵A的秩。

例2.3.2:

求3~20阶魔方阵的秩。

>> for n = 3:20
r(n) = rank(magic(n));
end
>> bar(r)
>> grid on
>> axis([2,21,0,20])
>> [3:20;r(3:20)]

例2.3.2

矩阵的迹

矩阵的迹等于矩阵的对角线元素之和,也等于矩阵的特征值之和。

trace(A):求矩阵A的迹。

例2.3.3:

>> A=[1,3,2; -3,2,1; 4,1,2]
A =
       1              3              2       
      -3              2              1       
       4              1              2       
>> b = trace(A)
b =
       5       
>> t = sum(diag(A))
t =
       5     

矩阵的范数

矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度。 

向量的3种常用范数

向量1—范数:向量元素的绝对值之和。

$$ \begin{Vmatrix} V \end{Vmatrix}_1= \sum_{i=1}^n |v_i| $$

向量2—范数:向量元素绝对值的平方和的平方根。

$$\begin{Vmatrix} V \end{Vmatrix}_2 = \sqrt{\sum_{i=1}^n |v_i|^2}$$

向量∞—范数:所有向量元素绝对值中的最大值。

$$\begin{Vmatrix} V \end{Vmatrix}_\infty = \max_{1\leq i\leq n}{\lbrace|v_i|\rbrace}$$

在MATLAB中,求向量范数的函数为:

norm(V)norm(V,2):计算向量V的2—范数。

norm(V,1):计算向量V的1—范数。

norm(V,inf):计算向量V的∞—范数。

矩阵的范数

矩阵A的1—范数:所有矩阵列元素绝对值之和的最大值。

$$\begin{Vmatrix} A \end{Vmatrix}_1 = \max_{1\leq j\leq n}{\lbrace\sum_{i=1}^m|a_{ij}|\rbrace}$$

矩阵A的2—范数:A'A矩阵的最大特征值的平方根。

$$\begin{Vmatrix} A \end{Vmatrix}_2 = \sqrt{\lambda_1}$$

其中λ1为A'A的最大特征值。

矩阵A的∞—范数:所有矩阵行元素绝对值之和的最大值。

$$\begin{Vmatrix} A \end{Vmatrix}_\infty = \max_{1\leq i\leq m}{\lbrace\sum_{i=1}^m|a_{ij}|\rbrace}$$

Matlab提供了求3种矩阵范数的函数,其函数调用格式与求向量的范数的函数完全相同。

例2.3.4:

>> x = [2 0 1; -1 1 0; -3 3 0]
x =
       2              0              1       
      -1              1              0       
      -3              3              0       
>> n = norm(x)
n =
    5909/1251  
>> n = norm(x,1)
n =
       6       

矩阵的条件数

矩阵A的条件数等于A的范数与A的逆矩阵的范数的乘积。

条件数越接近于1,矩阵的性能越好,反之,矩阵的性能越差。

在Matlab中,计算矩阵A的3种条件数的函数是:

cond(A,1):计算A的1—范数下的条件数。

cond(A)cond(A,2):计算A的2—范数数下的条件数。

cond(A,inf):计算A的∞—范数下的条件数。

例2.3.5:

求2~10阶希尔伯特矩阵的条件数。

>> for n = 2:10
c(n) = cond(hilb(n));
end
format long
c'

例2.3.5

标签:end,矩阵,2.3,Vmatrix,Matlab,求值,范数,向量
From: https://www.cnblogs.com/crepuscule/p/17080143.html

相关文章

  • 【Matlab学习2.2】矩阵变换
    对角阵对角阵:只有对角线上有非零元素的矩阵。 数量矩阵:对角线上的元素相等的对角矩阵。 单位矩阵:对角线上的元素都为1的对角矩阵。提取矩阵的对角线元素diag(A):提取矩......
  • Matlab:读取、写入.txt
    写入txta=[1,2,3;4,5,6];saveC:\Users\Administrator\Desktop\a.txt-asciia参考:https://blog.csdn.net/hengyaha/article/details/78995291读取txt%对......
  • matlab练习程序(三体模型)
    最近在看三体电视剧,正好看到了计算三体数值解那一部分,就想起了上学时看三体,也用matlab实现了三体的运动模拟。不过当时是通过时域外推的方式实现的,不是很严谨。下面通过......
  • 下标引用函数调用和表达式求值
    前言:下标的用法非常非常的重要关乎到数组函数指针结构体等等一切,函数的调用也有很多细节表达式求值涉及优先级问题要注意1.下标引用、函数调用和结构成员注:一切下标又是从0......
  • tcp/ip详解之04 以太网和802.3封装
    1.以太网和802.3封装在TCP/IP世界中,以太网IP数据报的封装是在RFC894[Hornig1984]中定义的,IEEE802网络的IP数据报封装是在RFC1042[PostelandReynolds1988]中定义的。......
  • 【Matlab学习2.1】特殊矩阵
    通用性的特殊矩阵zeros函数:产生全0矩阵,即零矩阵。ones函数:产生全1矩阵,即幺矩阵。eye函数:产生对角线为1的矩阵。当矩阵是方阵时,得到一个单位矩阵。rand函数:产生(0,1)区间均匀......
  • 通过MATLAB实现雷达数字化正交解调、数字脉冲处理、固定目标对消、动目标检测MTD、恒
    1.算法描述发射信号为13位巴克码和线性调频混合调制的信号,线性调频的中心频率为30MHz,调频带宽为4MHz,每一位码宽为10微秒,发射信号的帧周期为1毫秒.该雷达具有数字化正交解调......
  • m基于MSOPSO多策略粒子群算法的目标优化matlab仿真
    1.算法描述粒子群优化算法(particleswarmoptimization,PSO)是一种源于对鸟群捕食行为的研究而发明的进化计算技术,最先由Eberhart博士和Kennedy博士提出,其主要想法是:将每......
  • m基于遗传算法的城市生活垃圾回收网络优化matlab仿真
    1.算法描述目前,城市生活垃圾成为困扰我国大多数城市健康发展的痼疾.与国外相关国家相比,我国垃圾回收处理模式存在资源化程度低,处理方式单一等弊端,很多城市在积极探索新......
  • 通过MATLAB实现雷达数字化正交解调、数字脉冲处理、固定目标对消、动目标检测MTD、恒
    1.算法描述       发射信号为13位巴克码和线性调频混合调制的信号,线性调频的中心频率为30MHz,调频带宽为4MHz,每一位码宽为10微秒,发射信号的帧周期为1毫秒.该雷达具......