MATLAB编程和工程应用——第2章-矩阵及运算.ppt

MATLAB编程与工程应用 教 师 ：曾晓东 电 话 E_mail : zengxiaodong@263.net QQ : 241747139 第2章 矩阵及其运算 2.1 矩阵的实现 2.2 矩阵的运算 2.3 矩阵分析 2.4 字符串 2.5 向量与多项式 2.1 矩阵的实现 矩阵赋值的方法 一、直接输入 二、通过语句和函数产生 三、建立在文件中 四、从外部的数据文件中装入 2.1 矩阵的实现 一、直接输入 a=1; b=2; c=3; x=[5 b c; a*b a+c c/b] x= 5.000 2.000 3.000 2.000 4.000 1.500 2.1 矩阵的实现 获取矩阵的大小 Size函数 [n,m] = size(A) N为A的行数，m为A的列数 Length函数 N = length(x) N为向量x中元素的个数 若x为矩阵，则n为该矩阵的行数和列数中的较大者，即 Length(A)=max(size(A)) 2.1 矩阵的实现 二、矩阵的元素 矩阵的元素可以是实数，也可是复数 B=[-1/3 1.3 sqrt(3) 1+2*i] B=[-0.3333 1.3000 1.7321 1+2.0000i} C=[A;[10,11,12]] C=1 2 3 4 5 6 7 8 9 10 11 12 2.1 矩阵的实现 序列 s1:s2:s3 S1：起始值 S2：步长,无则默认为1 S3：终值 y=1:2:10 y=1 3 5 7 9 y=2:3:10 y=2 5 7 10 y=1:5 y=1 2 3 4 5 2.1 矩阵的实现 三、特殊矩阵的实现 1. 单位矩阵函数 eye(n),eye(size(B)),eye(n,m) 2.1 矩阵的实现 5. 对角矩阵 diag(V) 2.1 矩阵的实现 8. linspace和logspace a=linspace(n1,n2,n) 在线性空间上，行矢量的值从n1到n2，数据个数为n，缺省n为100。 a=linspace(1,10,10) a= 1 2 3 4 5 6 7 8 9 10 a=logspace(n1,n2,n) 在对数空间上，行矢量的值从10n1到10n2，数据个数为n，缺省n为50。这个指令为建立对数频域轴坐标提供了方便。 a=logspace(1,3,3) a= 10 100 1000 2.1 矩阵的实现 9. reshape函数 用于建立数值矩阵。 Av=1:12 %产生12个元素的行向量av Bm=reshape(av,3,4) %利用向量Av创建(3x4)矩阵Bm 结果: Av=1 2 3 4 5 6 7 8 9 10 11 12 Bm= 1 4 7 10 2 5 8 11 3 6 9 12 2.2 矩阵的运算 一、矩阵的转置 运算符：’及 .’ 对于实矩阵用（’）符号或（.’）求转置结果是一样的； 对于含复数的矩阵， （’）将同时对复数进行共轭处理， （.’）则只是将其排列形式进行转置。 2.2 矩阵的运算 二、矩阵的四则运算和乘方 + ；- ；* ；\和/ ；^；.* ；.\ ；./；.^ 只有维数相同的矩阵才能进行加减运算。 注意只有当两个矩阵中前一个矩阵的列数和后一个矩阵的行数相同时，才可以进行乘法运算。a\b运算等效于求a*x=b的解；而a/b等效于求x*b=a的解。 只有方阵才可以求幂。 点运算是两个维数相同矩阵对应元素之间的运算，在有的教材中也定义为数组运算。 例exp2_1.m 2.2 矩阵的运算 三、矩阵的翻转 1. 矩阵的左右翻转 fliplr(A) 2.2 矩阵的运算 四、矩阵的超越函数(例exp2_2.m) 在MATLAB中exp、sqrt等命令也可以作用到矩阵上，但这种运算是定义在矩阵的单个元素上的，即分别对矩阵的每一个元素进行计算。 超越数学函数可以在函数后加上m而成为矩阵的超越函数，例如：expm,sqrtm。矩阵的超越函数要求运算矩阵为方阵。 五、逆矩阵与行列式计算 求逆： inv(A)； 求行列式： det(A) 要求矩阵必须为方阵 2.2 矩阵的运算 六、关系和逻辑运算 MATLAB的运算符有三种类型：算术运算符、关系运算符、逻辑运算符。它们的处理顺序依次为算术运算符、关系运算符、逻辑运算符。 1. 关系运算符(例exp2_3.m) 小于 大于 = 小于等于 = 大于等于 == 等于 ~= 不等于 问：如何判断两个矩阵相等？ 2.2 矩阵的运算 2. 逻辑运算符 注意：在处理逻辑运算时，运算元只有两个值