第二讲矩阵的设置优秀课件.ppt

《第二讲矩阵的设置优秀课件.ppt》由会员分享，可在线阅读，更多相关《第二讲矩阵的设置优秀课件.ppt（33页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第二讲矩阵的设置第1页，本讲稿共33页数值运算的功能创建矩阵矩阵运算多项式运算线性方程组数值统计线性插值函数优化微分方程的数值解第2页，本讲稿共33页一、命令行的基本操作1.创建矩阵的方法创建矩阵的方法直接输入法规则：矩阵元素必须用 括住 矩阵元素必须用逗号逗号或空格空格分隔 在 内矩阵的行与行之间必须 用分号分隔第3页，本讲稿共33页 矩阵元素可以是任何matlab表达式，可以是实数，也可以是复数，复数可用特殊函数I，j 输入 例：a=1 2 3;4 5 6 x=2 pi/2;sqrt(3)3+5i 第4页，本讲稿共33页注意：只要是赋过值的变量，不管是否在屏幕上显示过，都存储在工作空间中，

2、以后可随时显示或调用。变量名尽可能不要重复，否则会覆盖。当一个指令或矩阵太长时，可用续行第5页，本讲稿共33页2.用matlab函数创建矩阵空阵 matlab允许输入空阵，当一项操作无结果时，返回空阵。rand 随机矩阵eye 单位矩阵zeros 全部元素都为0的矩阵ones 全部元素都为1的矩阵第6页，本讲稿共33页 还有伴随矩阵、稀疏矩阵、魔方矩阵、对角矩阵、范德蒙等矩阵的创建，就不一一介绍了。注意：matlab严格区分大小写字母，因此a与A是两个不同的变量。matlab函数名必须小写。第7页，本讲稿共33页3.矩阵的修改 直接修改 可用键找到所要修改的矩阵，用键移动到要修改的矩阵元素上即

3、可修改。指令修改 可以用A(,)=来修改。第8页，本讲稿共33页例如例如a=1 2 0;3 0 5;7 8 9a=1 2 0 3 0 5 7 8 9a(3,3)=0a=1 2 0 3 0 5 7 8 0还可以用函数subs修改，matlab6.0还可用find函数修改。第9页，本讲稿共33页1.矩阵加、减（,）运算规则：相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。二、矩阵运算第10页，本讲稿共33页2.矩阵乘（）运算规则：A矩阵的列数必须等于B矩阵的行数标量可与任何矩阵相乘。a=1 2 3;4 5 6;7 8 0

4、;b=1;2;3;c=a*bc=14 32 23 第11页，本讲稿共33页d=-1;0;2;f=pi*df=-3.1416 0 6.2832 矩阵除的运算在线性代数中没有，有矩阵逆的运算，在matlab中有两种矩阵除运算第12页，本讲稿共33页 a p a 自乘p次幂 方阵方阵1的整数的整数3.矩阵乘方 an,ap,pa对于p的其它值,计算将涉及特征值和特征向量，如果p是矩阵，a是标量ap使用特征值和特征向量自乘到p次幂；如a,p都是矩阵，ap则无意义。第13页，本讲稿共33页 a=1,2,3;4,5,6;7,8,9;a2 ans=30 36 42 66 81 96 102 126 150当一

5、个方阵有复数特征值或负实特征值时，非整数幂是复数阵。第14页，本讲稿共33页 a0.5 ans=0.4498+0.7623i 0.5526+0.2068i 0.6555-0.3487i 1.0185+0.0842i 1.2515+0.0228i 1.4844-0.0385i 1.5873-0.5940i 1.9503-0.1611i 2.3134+0.2717i第15页，本讲稿共33页inv 矩阵求逆det 行列式的值eig 矩阵的特征值diag 对角矩阵 A 矩阵A的转置sqrt 矩阵开方4.矩阵的其它运算 第16页，本讲稿共33页例：求矩阵例：求矩阵A的特征向量和特征值的特征向量和特征值A

6、=1:3;2:4;3:5;（初始化矩阵（初始化矩阵A）V,D=eig(A)V=0.8277 0.4082 0.3851 0.1424 -0.8165 0.5595 -0.5428 0.4082 0.7339D=-0.6235 0 0 0 -0.0000 0 0 0 9.6235第17页，本讲稿共33页5.矩阵的一些特殊操作矩阵的变维 a=1:12;b=reshape(a,3,4)c=zeros(3,4);c(:)=a(:)矩阵的变向 rot90:旋转;fliplr:上翻;flipud:下翻矩阵的抽取 diag:抽取主对角线;tril:抽取主下三角;triu:抽取主上三角；第18页，本讲稿共33

7、页关系运算 关系符号意义=小于小于或等于大于大于或等于等于不等于第19页，本讲稿共33页例：a=magic(3);b=eye(3);a=magic(3);b=eye(3);a=b a=b ans=ans=1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1第20页，本讲稿共33页 数组运算指元素对元素的算术运算，与通常意义上的由符号表示的线性代数矩阵运算不同1.数组加减(+,-)a+b a-b6.矩阵的数组运算 对应元素相加减（与矩阵加减对应元素相加减（与矩阵加减等效）等效）第21页，本讲稿共33页2.数组乘除(，./，.）ab a，b两数组必须有相同的行 和列两数组相应元素

8、相乘。a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;a.*b（矩阵的数组相乘）ans=2 8 18 4 15 30 49 72 90 第22页，本讲稿共33页a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;a*b（矩阵相乘）ans=25 37 46 55 85 109 85 133 172 第23页，本讲稿共33页a./b=b.aa.b=b./aa./b=b.a 都是a的元素被b的对应元 素除a.b=b./a 都是b的元素被a的对应元 素除例:a=1 2 3;b=4 5 6;c1=a.b;c2=b./a c1=4.0000

9、2.5000 2.0000 c2=4.0000 2.5000 2.0000 给出a,b对应元素间的商.第24页，本讲稿共33页3.数组乘方(.)元素对元素的幂例:a=1 2 3;b=4 5 6;z=a.2z=1.00 4.00 9.00z=a.bz=1.00 32.00 729.00第25页，本讲稿共33页四、代数方程组求解matlab中有两种除运算左除和右除。对于方程ax=b，a 为anm矩阵，有三种情况：当n=m时，此方程成为“恰定”方程 当nm时，此方程成为“超定”方程 当nm时，此方程成为“欠定”方程 matlab定义的除运算可以很方便地解上述三种方程第26页，本讲稿共33页1.恰定方

10、程组的解方程ax=b(a为非奇异)x=a-1 b 矩阵逆两种解:x=inv(a)b 采用求逆运算解方程 x=ab 采用左除运算解方程 第27页，本讲稿共33页方程ax=ba=1 2;2 3;b=8;13;x=inv(a)*b x=ab x=x=2.00 2.00 3.00 3.00 =例:x1+2x2=8 2x1+3x2=13第28页，本讲稿共33页2.超定方程组的解方程 ax=b,mn时此时不存在唯一解。方程解(a a)x=a b x=(a a)-1 a b 求逆法 x=ab matlab用最小二乘法找一 个准确地基本解。第29页，本讲稿共33页例:x1+2x2=1 2x1+3x2=2 3x

11、1+4x2=3a=1 2;2 3;3 4;b=1;2;3;解1 x=ab 解2 x=inv(aa)a b x=x=1.00 1.00 0 0.00 =a x =b第30页，本讲稿共33页3.欠定方程组的解 当方程数少于未知量个数时,即不定情况,有无穷多个解存在。matlab可求出两个解：用除法求的解x是具有最多零元素的解是具有最小长度或范数的解，这个解是基于伪逆pinv求得的。第31页，本讲稿共33页 x1+2x2+3x3=1 2x1+3x2+4x3=2a=1 2 3;2 3 4;b=1;2;x=ab x=pinv(a)b x=x=1.00 0.83 0 0.33 0 -0.17=a x =b第32页，本讲稿共33页小 结 本节介绍了matlab语言的数值运算功能，通过学习应该掌握：如何创建矩阵、修改矩阵符号的用法矩阵及数组运算多项式运算线性方程组与微分运算第33页，本讲稿共33页