MATLAB矩阵运算

MATLAB矩阵运算矩阵定义MATLAB以矩阵作为数据操作的基本单位,这使得矩阵运算变得非常简捷、方便、高效。矩阵是由m×n个数av(i=1,2,…,m;j=1,2,…,n)排成的m行n列数表,记成:若m=n,则该矩阵为n阶矩阵(n阶方阵)。矩阵创建在键盘上直接按行方式输入矩阵是最方便、最常用的创建数值矩阵的方法,尤其适合较小的简单矩阵。在用此方法创建矩阵时,应当注意以下几点:输入矩阵时要以“[]”为其标识符号,矩阵的所有元素必须都在括号内。 矩阵同行元素之间由空格(个数不限)或…

大家好,又见面了,我是你们的朋友全栈君。

矩阵定义

MATLAB以矩阵作为数据操作的基本单位,这使得矩阵运算变得非常简捷、方便、高效。矩阵是由m×n个数av (i=1,2,…,m; j = 1,2,…,n)排成的m行n列数表,记成:

     watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_15,color_FFFFFF,t_70,g_se,x_16

若m=n,则该矩阵为n阶矩阵(n阶方阵)。

矩阵创建

在键盘上直接按行方式输入矩阵是最方便、最常用的创建数值矩阵的方法,尤其适合较小的简单矩阵。在用此方法创建矩阵时,应当注意以下几点:

  • 输入矩阵时要以“[ ]”为其标识符号,矩阵的所有元素必须都在括号内。
  • 矩阵同行元素之间由空格(个数不限)或逗号分隔,行与行之间用分号或回车键分隔。矩阵大小不需要预先定义。
  • 矩阵元素可以是运算表达式。
  • 若“[ ]”中无元素,表示空矩阵。
  • 如果不想显示中间结果,可以用“;”结束。

矩阵生成

矩阵的生成除了直接输入法,还可以利用M文件生成法和文本文件生成法等。

1.利用M文件创建

       当矩阵的规模比较大时,直接输入法就显得笨拙,出差错也不易修改。为了解决这些问题,可以将所要输入的矩阵按格式先写入一文本文件中,并将此文件以m为其扩展名,即M文件。
       M文件是一种可以在 MATLAB环境下运行的文本文件,它可以分为命令式文件和函数式文件两种。在此处主要用到的是命令式M文件,用它的简单形式来创建大型矩阵。在MATLAB命令行窗口中输入M文件名,所要输入的大型矩阵即可被输入到内存中。
       M文件中的变量名与文件名不能相同,否则会造成变量名和函数名的混乱。运行M文件时,需要先将M文件sample.m复制到当前目录文件夹下,否则运行时无法调用。

2.利用文本创建
        MATLAB中的矩阵还可以由文本文件创建,即在文件夹(通常为work文件夹)中建立txt 文件,在命令行窗口中直接调用此文件名即可。

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16
         用户可以直接用函数来生成某些特定的矩阵,常用的函数如下:

  • eye(n):创建n×n单位矩阵。
  • eye(m,n):创建m×n的单位矩阵。
  • eye(size(A)):创建与A维数相同的单位矩阵。
  • ones(n):创建n×n全1矩阵。
  • ones(m,n):创建m×n全1矩阵。
  • ones(size(A)):创建与A维数相同的全1阵。
  • zeros(m,n):创建m×n全0矩阵。
  • zeros(size(A)):创建与A维数相同的全0阵。
  • rand(n):在[0,1]区间内创建一个n×n均匀分布的随机矩阵。
  • rand(m,n):在[0,1]区间内创建一个m×n均匀分布的随机矩阵。
  • rand(size(A)):在[0,1]区间内创建一个与A维数相同的均匀分布的随机矩阵。
  • compan(P):创建系数向量是P的多项式的伴随矩阵。
  • diag(v):创建一向量v中的元素为对角的对角阵。
  • hilb(n):创建n×n的Hilbert矩阵。
  • magic(n):生成n阶魔方矩阵。
  • sparse(A):将矩阵A转化为稀疏矩阵形式,即由A的非零元素和下标构成稀疏矩阵S。若A本身为稀疏矩阵,则返回A本身。

 MATLAB矩阵运算

1.矩阵元素的修改

 矩阵建立起来之后,还需要对其元素进行修改。下表列出了常用的矩阵元素修改命令:

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16
2.矩阵的变维

矩阵的变维可以用符号“:”法和reshape函数法。reshape函数的调用形式如下:

reshape(X,m,n)

将已知的矩阵变为m行n列的矩阵。

3.矩阵的变向

Rot(90):将A逆时针旋转90度

Rot(90,k):将A逆时针旋转90度*k,k可为正整数或负整数

Fliplr(X):将X左右翻转

flipud(X):将X上下翻转

flipdim(X,dim)q:dim=1时对行翻转,dim=2对列翻转

4.矩阵的抽取

对矩阵元素的抽取主要是指对角元素和上(下)三角阵的抽取。

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16

5.矩阵的加法运算

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16

6.矩阵的乘法运算 

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16

 7.矩阵的点乘运算

 点乘运算指将两矩阵中相同位置的元素进行相乘运算,将积保存在原位置组成新矩阵。

8.矩阵的左除运算

线性方程组D*X=B,如果D非奇异,即它的逆矩阵inv(D)存在,则其解用MATLAB表示为:

X=inv(D)*B=D\B

符号“\”称为左除,即分母放在左边。

左除的条件:B的行数等于D的阶数(D的行数和列数相同,简称阶数)。

9.矩阵的右除运算

线性方程组D*X=B,如果D非奇异,即它的逆矩阵inv(D)存在,则其解用MATLAB表示为:

X=B*inv(D)=B/D

符号“/”称为右除,即分母放在右边。

右除的条件:B的列数等于D的阶数(D的行数和列数相同,简称阶数)。

10.矩阵的幂运算

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16

11.矩阵的逆

MATLAB矩阵运算

 

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBAc2hybHhscWw=,size_20,color_FFFFFF,t_70,g_se,x_16    

12.范数

MATLAB矩阵运算

MATLAB矩阵运算

13.奇异值分解

奇异值分解(SVD)是现代数值分析(尤其是数值计算)的最基本和最重要的工具之一,因此在实际工程中有着广泛的应用。

所谓的SVD分解指的是将mxn矩阵A表示为3个矩阵乘积形式: USV^T,其中U为mxm酉矩阵,V为n×n酉矩阵,S为对角矩阵,其对角线元素为矩阵A奇异值且满足S1>=S2>=..>=Sr>Sr+1=..=Sn,r为矩阵A的秩。在 MATLAB中,这种分解是通过SVD命令来实现的。

SVD的命令调用格式:
MATLAB矩阵运算
 

 

 

                                                                                                                                                              

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/153005.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 如何获取手机gps数据_获取手机gps数据

    如何获取手机gps数据_获取手机gps数据不借助QXDM,直接通过手机抓取qxdm日志。主要用于测试GPS的时候使用,因为测试GPS大都需要在室外进行(室外卫星信号良好),并且GPS问题一般需要获取QXDM日志,在室外不能连接PC端使用QXDM客户端工具。adbshellmkdir/sdcard/diag_logsadbpushgnss.cfg/sdcard/diag_logs/Diag.cfg_GNSS_Generic_2020.cfgadbshelldiag_mdlog-f/sdcard/diag_logs/Di

    2022年10月2日
    3
  • 小米8手机相册中的图片怎么识别文字?[通俗易懂]

    小米8手机相册中的图片怎么识别文字?[通俗易懂]小米8大家都买了吗?最近新出的还是很好的,小米8手机相册中的图片怎么识别?对于用小米手机的人来说很简单,下面分享一个方法。1在手机上输入PDF阅读器,然后开始识别相册图片的文字2打开就是这样的界面,找到倒数第2个的小功能3选择拍照识别中的相册,接下来添加图片4这张是…

    2022年5月18日
    83
  • 向navicat中导入数据库时出现错误_sqlserver导入sql文件

    向navicat中导入数据库时出现错误_sqlserver导入sql文件在Navicat导出的 或者别的sql文件,在使用Navicat导入时候 出现异常失败报错问题。搜索了很多资料查看,发现是没有解决掉的。最后无意间想起使用 MySql 直接使用命令导入尝试,发现可行的简单粗暴,直接打开你的MySql 登录以后 选择 要导入的数据库use 数据库名称;source 文件的绝对路径;完事 ,坐等~…

    2022年8月19日
    27
  • 国家标准《电子计算机机房设计规范》(GB50174-93)_计算机机房建设标准规范

    国家标准《电子计算机机房设计规范》(GB50174-93)_计算机机房建设标准规范dgtool:GB50174-2008电子计算机机房设计规范目次1总则(1)2术语(2)3机房分级与性能要求(6)3.1机房分级(6)3.2性能要求(6)4机房位置及设备布置(7)4.1机房位置选择(7)4.2机房组成(7)4.3设备布置(8)5环境要求(9)5.1温度、相对湿度及空气含尘浓度(9)5.2噪声、电磁干扰、振动及静电(9)6建筑与结构(1o)6.1一般规定(1o)6.2人流、物…

    2022年10月2日
    4
  • 做淘宝的经验_阿里社招offer审批通过率

    做淘宝的经验_阿里社招offer审批通过率如何轻松拿到 淘宝前端 offer

    2022年4月20日
    38
  • css filter属性不生效问题处理

    css filter属性不生效问题处理明天是清明节假期 感谢疫情以来第一个假期就让我学到了一个棘手问题解决方案由于清明节公祭日 公司的网站都要改为灰度样式 明明就是一个很简单的问题 再是在样式表里增加这个一个样式即可 答案很简单 html filter grayscale 100 webkit filter grayscale 100 moz filter grayscale 100 ms fi

    2025年11月10日
    3

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号