第1-2章节MATLAB概述幻灯片.pptVIP

方法3：在M文件输入矩阵的行列具体数值，存盘后在命令窗口中输入M文件的文件名，运行该文件就能自动生成矩阵。 如下述例子就是在M文件（ls1）存放的内容 exm=[1 2 3 -1 -2 -3 3 4 5] 在命令窗口中输入ls1 方法4: 利用外部文件装入到指定的矩阵 当数据量较多时.可借助于其它编辑的软件来实现 如:通过windows附件下的记事本进行编辑存盘后,使用命令:load [盘符][路径] 文件名.扩展名 如:load c:\lsq\lsq.txt 产生一个矩阵lsq,它的行列数由lsq.txt决定 方法5: 导入Excel数据文件中的数据该矩阵变量名为Excel数据文件中的文件名 方法为：执行主窗口下的FILE菜单下的Import Data功能项，按照向导提示初步完成。 2. 复数和复数矩阵 复数单位为i或j 复数表示: z=3+4i和 z=3+4j作用一样 复数矩阵表示 a=[12;34]+[56;78]*i a=[1+5i 2+6i ; 3+7i 4+8i] （二）矩阵元素 矩阵元素名称（下标列表） 一维：x(i) 二维： x(i,j) 矩阵元素可为任意的表达式,利用下标可直接输入或修改矩阵元素值或修正其大小 x=[-1.3,sqrt(3),(1+2+3)*4/5] 产生1*3矩阵 x(5)=abs(x(1)) x(4)值自动为0 矩阵连接(小矩阵构成大矩阵) a=[a ; r] 将同列数的a和r合并,构成a 大矩阵中抽取一小矩阵 b=a(1:3,:) 抽取a中1:3行所有元素构成b c=a(:,1:2) 抽取a中1:2列所有元素构成b 矩阵行列删除 利用空矩阵可从矩阵中删除指定的行和列 例如要删除b第二行,则使用命令: b(2,:)=[] 例如要删除b第二列,则使用命令: b(:,2)=[] （三）常数与算术运算符 （1）常数采用十进制数 （2） +、－、 * / (右除)、 \ (左除)、 ＾幂指数 （四） 函数 （1）MATLAB的内部函数 （2）MATLAB系统附带的各种工具箱中的M文件所提供的实用函数 （3）由用户自己增加的函数 二、基本矩阵和矩阵运算 1.产生基本矩阵的函数 Eye(n) 单位矩阵 格式: y= Eye(n) n*n的单位阵y y= Eye(m, n) m*n的单位阵y y= Eye(size(a)) 产生与A同维的单位阵y (2) Ones 全1列阵 格式: y= ones(n) n*n的单位阵y y= ones(m, n) m*n的单位阵y y=ones(size(a)) 产生与A同维的单位阵y (3) Zero 全0矩阵 (4) Rand 随机函数和阵列 2. 矩阵转置 格式：x’ 3.矩阵的行列式值 函数为：det（） 4.矩阵加,减 要求：A,B两矩阵整数相同: 矩阵与矩阵运算： 对应位量上的元素相加.减；矩阵与标量加、减: 矩阵中的每个元素都与标量进行加.,减运算 5.矩阵乘法 要求：前一点矩阵的列数等于后一点矩阵的行数 ① 矩阵与矩阵相乘 ② 矩阵与标量相乘 每个元素均乘以某一个标量 6.矩阵的除法 (左除\)、 (右除/) a\b=inv(a)*b ax=b ? x=a-1b b/a=b*inv(a) x*a=b ? x=b*a-1 注意： (1) inv(a) 求矩阵的逆，要求a 为方阵 (2) Matlab采用IEEE算法,即使A是奇异(A的行列式值为0)照样进行,但会给出“Waring: Matrix is singular to working”,当矩阵A为病态时, Matlab算法会造成计算误差可能很大,“Waring:Matrix is badly scaled to working precision” 应用: （1）解静定方程的解:ax=b 方法1: x=inv(a)*b 方法2: x=a\b （2）求解的方程是超静定的(方程的个数多于未知数的个数) 用Householder直接求解方程的最小二乘解 举例:方程的