MATLAB矩阵运算

本文涉及的产品
实时数仓Hologres,5000CU*H 100GB 3个月
智能开放搜索 OpenSearch行业算法版,1GB 20LCU 1个月
检索分析服务 Elasticsearch 版,2核4GB开发者规格 1个月
简介: MATLAB矩阵运算

 

1. 矩阵的加减乘除和(共轭)转置

(1) 矩阵的加法和减法

  如果矩阵AB有相同的维度(行数和列数都相等),则可以定义它们的和A+B以及它们的差A-B,得到一个与AB同维度的矩阵C,其中Cij=Aij+BijAij-Bij.

     另外Matlab还支持任意一个矩阵A与一个标量s相加,结果为矩阵的每一个元素加减标量,得到一个与A同维度的新的矩阵,即A+s的各个元素为Aij+s.

(2) 矩阵的乘法

  如果矩阵A的列数等于矩阵B的行数,则可以将AB相乘,命令为A*B,得到一个新的矩阵CC的行数等于A的行数,列数等于B的列数. 由于矩阵的乘法不满足交换律,所以一般A*B不等于B*A.MATLAB矩阵运算文档

(3) 矩阵的张量积(tensor product)

  矩阵AB的张量积AB可以方便地用kron函数计算,即使用命令kron(A,B), 例如

image.gif编辑

(4) 矩阵的除法

  在MatLab中,有两个矩阵除法符号,左除\和右除/. 如果A是一个非奇异方阵(nonsingular square, 即满秩方阵),B的行数与A的行数相等,那么A\B=A-1B. 如果C的列数与A的列数相等,那么C/A=CA-1.

  从另一个角度来看,X=A\B是矩阵方程AX=B的解,X=C/A是矩阵方程XA=C的解. 如果b是一个行数与A的行数相等的列向量,则向量x=A\b是线性方程组 Ax=b的解. 且在矩阵方程AX=B中,A可以是一个m×n的矩阵,如果m=n则有唯一解;如果m<n则有多个解,Matlab会返回一个基础解;如果m>n则会返回一个最小方阵解.

(5) 矩阵的转置和共轭转置

  在Matlab中,矩阵的共轭转置用撇号表示,如果不需要对元素进行共轭运算,仅仅只对矩阵进行转置,则在撇号之前输入一个点号,即 .’ . 对于实数矩阵AA’A.’是相同的.

2. 矩阵元素操作运算

  矩阵的运算既可以是如前所述的正常的整体运算,也可以是矩阵对应的元素依次进行标量运算,也叫数组运算,即把矩阵看做是二维数组. 对矩阵进行数组运算后得到的结果是一个与参与运算的矩阵维度相同的新矩阵,.这种元素间的算术运算的前提是参与运算的两个矩阵的维数要相同.

     对于加法和减法,元素操作运算和矩阵运算没有差别,而对于乘、除和幂运算符,相应的数组运算符是在一般的算术运算符前面加上一个点号,如

+       -       .*       ./       .\       .^

其中,A./B 是指A中的元素除以B中相应的元素,即A./B 的第i行第j列的元素(A./B)ij=Aij/Bij,而(A.\B)ij=Bij\Aij. 这些元素运算符的使用例子如下所示:


image.gif编辑


image.gif编辑

  在Matlab中预定义的数学标准函数,如sin(x), abs(x)等都是基于对矩阵元素的运算. 如果函数f(x)是这样的一个函数,A是一个m×n的矩阵,其元素是aij ,那么 f(A)也是一个m×n的矩阵,其第i行第j列的元素为f(aij),例如


image.gif编辑

其中pi是Matlab的预定义变量,值为π,i也是预定义变量,表示复数的单位.

3. 常用的矩阵函数

  矩阵函数是指参数为矩阵的函数,函数结果可能是一个标量值也可能是一个函数或者向量. Matlab中常用的矩阵函数包括:

  (1) rank(A): 求矩阵A的秩,即A中线性无关的行数或者线性无关的列数.


image.gif编辑

  (2) det(A): 求矩阵A的行列式值.


image.gif编辑

  (3) inv(A): 如果A是一个非奇异(nonsingular)矩阵,则inv(A)返回A的逆矩阵. 另外还可以用左除A\eye(n)或右除eye(n)/A来计算A的逆,且在Matlab中用左除或右除来计算逆所花的计算时间比用inv函数要少,也比inv具有更好的容错性(error-detection properties).

  (4) dot(x,y): 求同维度的向量x和y的内积/点积. 若AB是两个具有相同维度的矩阵,则dot(A,B)是计算AB对应列的内积,结果是一个行向量,这个行向量的列数等于AB的列数. 例如


image.gif编辑

  (5) cross(x,y): 计算同维度的向量xy的叉积,结果是一个向量,其方向由右手定则决定,长度等于|x|*|y|sin<x,y>. 若AB是两个具有相同维度的矩阵,则cross(A,B)是计算AB对应列的叉积,结果是一个维度与AB相等的矩阵.

  (6) kron(A,B): 得到矩阵AB的张量积.

  (7) isequal(A,B): 如果矩阵AB是相同的,即具有相同的维数和相同的内容,则返回1.

  (8) isreal(A): 判断A是否是一个实矩阵,如果是则返回1,否则返回0.

  (9) trace(A): 计算方阵A的迹,即对角线元素之和.

  (10) eig(A): 计算方阵A的特征值,结果是一个列向量,向量中元素的个数等于特征值的个数,即A的维度(A的行数或列数).

  (11) [U,D]=eig(A): 计算方阵A的特征值和特征向量,得到两个方阵UD,其中D的对角线元素为A的特征值,U的列向量为A的特征向量(可能是未normalize的结果),例如


image.gif编辑

  (12) length(V): 求向量V的长度,即V的元素数量.

  (14) size(A): 若A是m行n列的矩阵,则返回行向量[m,n].

相关文章
|
算法
MATlAB运用——数值积分
MATlAB运用——数值积分
170 0
|
计算机视觉
专题七MATLAB符号计算-2
专题七MATLAB符号计算
88 0
专题七MATLAB符号计算-1
专题七MATLAB符号计算
85 0
|
2月前
|
算法
MATLAB符号计算
【10月更文挑战第9天】MATLAB不仅擅长数值计算,还具备强大的符号计算功能,支持代数运算、方程求解、微积分等。本文介绍如何使用MATLAB的符号工具箱进行符号变量定义、方程求解、微分积分及矩阵运算,并通过多个实际应用案例展示了其在机械系统、电路分析、经济优化和物理运动学等领域的应用。此外,文章还提供了符号计算的最佳实践和未来展望。
51 2
|
算法
MATLAB运用(4)-常微分方程的性态和初值问题
MATLAB运用(4)-常微分方程的性态和初值问题
115 1
|
机器学习/深度学习
专题二 MATLAB矩阵处理-1
专题二 MATLAB矩阵处理
110 0
|
存储 人工智能 NoSQL
专题二 MATLAB矩阵处理-2
专题二 MATLAB矩阵处理
104 0
|
机器学习/深度学习 算法
MATLAB实现傅里叶变换
MATLAB实现傅里叶变换
|
Windows
MATLAB的符号计算
MATLAB的符号计算
151 0
|
算法 索引
MATLAB矩阵运算(2)
1.2.14  特殊运算 1.矩阵对角线元素的抽取 函数  diag 格式  X = diag(v,k)   %以向量v的元素作为矩阵X的第k条对角线元素,当k=0时,v为X的主对角线;当k>0时,v为上方第k条对角线;当k0:抽取上方第k条对角线元素;k> v=[1 2 3]; >> x=...
950 0