向量与矩阵的运算实验.ppt

编辑课件 数学实验 向量与矩阵运算 向量与矩阵的生成 向量与矩阵运算 向量的生成 直接输入: a=[1,2,3,4] 冒号运算符 a=[1:4] == a=[1, 2, 3, 4] b=[0:pi/3:pi] == b=[0, 1.0472, 2.0944, 3.1416] c=[6:-2:0] == c = [6, 4, 2, 0] 例： 从矩阵中抽取行或列 向量与矩阵的生成（续） 向量与矩阵运算 矩阵的生成 直接输入: A=[1, 2, 3; 4, 5, 6; 7, 8, 9] 由向量生成 由函数生成 通过编写m文件生成 例： x=[1,2,3];y=[2,3,4]; A=[x,y], B=[x;y] 例： C=magic(3) 常见矩阵生成函数 矩阵操作 提取矩阵的部分元素： 冒号运算符 A(:) A的所有元素 A(:,:) 二维矩阵A的所有元素 A(:,k) A的第 k 列， A(k,:) A的第 k 行 A(k:m) A的第 k 到第 m 个元素 A(:,k:m) A的第 k 到第 m 列组成的子矩阵 A(:) 与 A(:,:) 的区别 ? 如何获得由 A 的第一、三行和第一、二列组成的子矩阵？ 自己动手 矩阵操作 矩阵的旋转 fliplr(A) 左右旋转 flipud(A) 上下旋转 rot90(A) 逆时针旋转 90 度； rot90(A,k) 逆时针旋转 k×90 度 例： A=[1 2 3;4 5 6] B=fliplr(A) C=flipud(A) D=rot90(A), E=rot90(A,-1) 矩阵操作 矩阵的转置与共轭转置 ’ 共轭转置 .’ 转置，矩阵元素不取共轭 例： A=[1 2;2i 3i] B=A’ C=A.’ 点与单引号之间不能有空格！ 矩阵操作 改变矩阵的形状：reshape reshape(A,m,n): 将矩阵元素按 列方向 进行重组 重组后得到的新矩阵的元素个数 必须与原矩阵元素个数相等！ 矩阵操作 查看矩阵的大小：size size(A) 列出矩阵 A 的行数和列数 size(A,1) 返回矩阵 A 的行数 size(A,2) 返回矩阵 A 的列数 例： A=[1 2 3; 4 5 6] size(A) size(A,1) size(A,2) length(x) 返回向量 X 的长度 length(A) 等价于 max(size(A)) 矩阵基本运算 矩阵的加减：对应分量进行运算 要求参与加减运算的矩阵具有 相同的维数 例： A=[1 2 3; 4 5 6]; B=[3 2 1; 6 5 4] C=A+B; D=A-B; 矩阵的普通乘法 要求参与运算的矩阵满足线性代数中矩阵相乘的原则 例： A=[1 2 3; 4 5 6]; B=[2 1; 3 4]; C=A*B 矩阵基本运算 矩阵的除法：/、\ 右除和左除 若 A 可逆方阵，则 A\B == A 的逆左乘 B == inv(A)*B B/A == A 的逆右乘 B == B*inv(A) X=A\B == A*X=B X=B/A == X*A=B 通常，矩阵除法可以理解为 当 A 和 B 行数相等时即可进行左除 当 A 和 B 列数相等时即可进行右除 矩阵的乘方 A 是方阵，p 是正整数 A^p 表示 A 的 p 次幂，即 p 个 A 相乘。 若 A 是方阵，p 不是正整数 A^p 的计算涉及到 A 的特征值分解，即若 A = V*D*V-1 则 A^p=V*(D.^p)/V 矩阵的乘方 若 a 是标量，A 是方阵，且 [V,D] = eig(A)，则 a^A ＝ V*(a^D)/V 若 A, P 均是矩阵，则 A^P 无定义 若 a 是标量， 则 矩阵的 Kronecker 乘积 矩阵 Kronecker 乘积的定义 设A是n×m矩阵，B是p×q矩阵，则A与B的kronecker乘积为： Kronecker 乘积的性质 是 np×mq 矩阵；通常 任何两个矩阵都有 Kronecker 乘积 Matlab 中实现两个矩阵 Kronecker 相乘的函数为 kron(A,B) Kronecker乘积有时也称张量积 矩阵的数组运算 数组运算：对应元素进行运算 点与算术运算符之间