数值计算可以分为两类:矩阵运算和矩阵元素运算

3.1 矩阵函数和特殊矩阵

       矩阵代数的处理数组大部分以一维数组(向量)和二维数组(矩阵)为主。

常见的矩阵处理函数如下:

特殊矩阵:

稀疏矩阵:

稀疏矩阵的创建通常可以通过使用sparse和spdiags函数来进行:

  • sparse(i,j,s,m,n,nz_max)使用[i,j,s]来创建m*n维稀疏矩阵s,其中s为所有非零元素构成的向量,i,j分别为非零元素的行下标和列下标。
  • spdiags(B,d,m,n)可以抽取、创建或替代对角稀疏矩阵,其中m和n分别为矩阵的行和列的维数,d表示长度为p的整数向量,B为满矩阵。
  • 稀疏矩阵和满矩阵之间可以通过sparse和full函数进行转换。
  • spy(A)画出A矩阵

3-2  矩阵分析

  • 范数分析:norm()函数可以用来求解向量和矩阵的范数,

           如 n2=norm(x),ninf=norm(x,inf),nneg_inf=norm(x,-inf)分别求解二范数,无穷范数和负无穷范数

          求解大型矩阵的范数和稀疏矩阵时,可以考虑采用normest来求解。

  • 条件数分析:矩阵的条件数总是大于1,正交矩阵的条件数等于1,奇异矩阵的条件数为无穷大,

         病态矩阵的条件数是比较大的数据。 矩阵的条件数可以通过函数cond()求解。

3.3  线性方程组

下面介绍三种方程组:恰定方程组、超定方程组、欠定方程组。

  •  恰定方程组:方程组的个数和未知数的个数相等,分为奇异和非奇异两种。

%求解非奇异恰定方程组的解
>> A=[1 2 3;4 5 6;7 8 0];
>> y=[366;804;351]

%方法一:判断A的条件数、秩以及增广矩阵[A y]的秩
>> rank(A)

>>rank([A y])

>> cond(A)

>> x_1=inv(A)*y

>> %方法2 通过矩阵左除来求解线性方程组的解,LU分解,速度快
>> x_2=A\y

%求解奇异恰定方程组的解

>> A=[1 3 7;-1 4 4;1 10 18];
>> y=[6;4;15];
>> %方法1:通过矩阵左除来求解线性方程组的解
>> x_1=A\y

>> %方法2 判断A的条件数、A的秩以及增广矩阵[A y]的秩
>> det_A=det(A)

>> rank_A=rank(A)

>> rank_Ay=rank([A y])

>> cond_A=cond(A)

>> x_2=pinv(A)*y;
>> disp(x_2)

  •  欠定方程组:方程组的个数小于未知数的个数,存在无数个解。最小范数解,略。

方法1:左除法:x_1=A\y

方法2:伪逆矩阵法(范数最小)

x_2=pinv(A)*y

  •  超定方程组:方程组的个数大于未知数的个数,存在无数个解。最小二乘解,略。

3.4  矩阵分解

   主要包括Cholesky 分解、LU分解和QR分解

3.5 特征值分析(略)

Matlab--从入门到精通(Chapter3 矩阵运算)的更多相关文章

  1. Matlab 从入门到精通 Chapter11 文件读取I/O

    11.1 工作空间数据读取 将工作空间的变量保存为文件,可以使用save命令.  save('filename') 将文件保存在当前目录下,文件名为filename.mat save('filenam ...

  2. 《Matlab实用案例》系列Matlab从入门到精通实用100例案例教程目录(持续更新)

    目录 1. 专栏简介 2. 专栏地址 3. 专栏目录 1. 专栏简介 2. 专栏地址 「 刘一哥与GIS的故事 」之<Matlab使用案例> 3. 专栏目录 [MATLAB统计分析与应用1 ...

  3. Matlab从入门到精通 Chapter5 数据可视化--

    5-1 图形绘制示例 >> x2=-17:0.02:3; >> y2=1./((x2+3).^2+1)+1./((x2+9).^2+4)+5; >> subplot ...

  4. CUDA从入门到精通

    http://blog.csdn.net/augusdi/article/details/12833235 CUDA从入门到精通(零):写在前面 在老板的要求下.本博主从2012年上高性能计算课程開始 ...

  5. CUDA从入门到精通 - Augusdi的专栏 - 博客频道 - CSDN.NET

    http://blog.csdn.net/augusdi/article/details/12833235 CUDA从入门到精通 - Augusdi的专栏 - 博客频道 - CSDN.NET CUDA ...

  6. MyBatis从入门到精通(1):MyBatis入门

    作为一个自学Java的自动化专业211大学本科生,在学习和实践过程中"趟了不少雷",所以有志于建立一个适合同样有热情学习Java技术的参考"排雷手册". 最近在 ...

  7. 【halcon教程资料】全网汇总如何快速、高效率学习机器视觉从入门到精通

    我以八年的视觉工程师开发的工作经验告诉你,你不要再因为学习halcon发愁了,我接触过很多学习halcon的小白,并不是不愿意学,而是不知道怎么快速.高效率的学习精通,一天天的过去了,对学习halco ...

  8. 【CUDA开发】CUDA从入门到精通

    CUDA从入门到精通(零):写在前面 在老板的要求下,本博主从2012年上高性能计算课程开始接触CUDA编程,随后将该技术应用到了实际项目中,使处理程序加速超过1K,可见基于图形显示器的并行计算对于追 ...

  9. Simulink仿真入门到精通(十九) 总结回顾&自我练习

    从2019年12月27到2020年2月12日,学习了Simulink仿真及代码生成技术入门到精通,历时17天. 学习的比较粗糙,有一些地方还没理解透彻,全书梳理总结: Simulink的基础模块已基本 ...

  10. matlab快速入门

    matlab快速入门 1矩阵 生成矩阵 ​ % 直接法 a = [1,2,3;4,5,6;7,8,9]; % 冒号一维矩阵 a = 开始:步长:结束,步长为1可省略 b = 1:1:10; % 1,2 ...

随机推荐

  1. FrameLsyout

    <FrameLayout xmlns:android="http://schemas.android.com/apk/res/android" android:layout_ ...

  2. RePlace函数

    <html><body> <script type="text/javascript"> var str="Visit qq1!&qu ...

  3. windows下Keras框架搭建

    1. 安装Anaconda https://mirrors.tuna.tsinghua.edu.cn/anaconda/archive/ conda info来查询安装信息 conda list可以查 ...

  4. POJ 3624 Charm Bracelet【01背包】

    解题思路:直接套公式就能做的01背包, for(i=1;i<=n;i++) { for(v=w[i];v<=m;v++) f[i,v]=max(f[i,v],f[i-1,v-w[i]]+d ...

  5. 12、Camel: Content-Aware and Meta-path Augmented Metric Learning for Author Identification----作者识别

    摘自:https://blog.csdn.net/me_yundou/article/details/80459341 具体看上面链接 一.摘要: 这篇文章主要介绍的是作者识别(author iden ...

  6. BZOJ 2342 [SHOI2011]双倍回文 (回文自动机)

    题目大意:略 先建出$PAM$ 因为双倍回文串一定是4的倍数,所以找出$PAM$里所有$dep$能整除4的节点 看这个串是否存在一个回文后缀,长度恰好为它的一半,沿着$pre$链往上跳就行了 暴跳可能 ...

  7. webpack配置相关的页面异常

    原文:https://www.cnblogs.com/Hsong/p/9023341.html 前言 在团队协作开发中,为了统一代码风格,避免一些低级错误,应该设有团队成员统一遵守的编码规范.很多语言 ...

  8. [Ubuntu]安装中文输入法

    安装了英文版的Ubuntu12.04之后.使用ctrt+space无法切换到中文,才反应过来. 于是查询了下怎么在Ubuntu下安装和配置中文输入法. 第一步安装语言包: 在左側的菜单条中点击 sys ...

  9. volatile 和 mutable 关键字

    经常接触,但是过一段时间可能又忘了.做个记录. volatile是表示变量易变,不要放缓存,每次实际取,尤其是多线程. mutable表示一个const 类或者数据结构里面,某个字段是可以改变的.

  10. POJ 3696

    这里面的一个转换的小技巧很重要,把888...8转换成(10^x-1)/9*8.神来之笔,佩服. 这样有(10^x-1)/9*8=L*p得10^x-1=L*p*9/8,设m=9*L/gcd(L,8). ...