数学实验第五讲 用mathematica的相应功能进行向量矩阵运算.ppt

数 学 实 验 授课老师： 所属院系：数计学院 E-mail ： 第五讲 用 Mathematica 的 相应功能进行向量、矩阵运算 用 Mathematica 的相应功能进行向量、矩阵运算 ? 获得表的元素 ? 向量和矩阵的输入 ? 表的维数和矩阵的加、减法 ? 向量和矩阵的乘法 ? 关于矩阵的几个常用函数 在 Mathematica 中，有序数组被称为 “表”。“表”既可以表示成集合，也 可以 表示成向量和矩阵。 Mathematica 中的许 多函数都可以作用在表上。 用 Mathematica 的相应功能进行向量、矩阵运算 ? 使用键盘输入一个表时，用 { } 将元素括 起，元素之间用逗号分隔。 例 1 ： 输入一组数据 0 ， 16 ， 64 ， 144 ， 256 ，并把这 个数组定义为变量 data 命令： data={0 ， 16 ， 64 ， 144 ， 256} 向量和矩阵的输入 例 2 ： 输入矩阵 M= 命令 : 2,-2}} 注意：矩阵的每一行用 { } 括起来 , 行与行之间用逗号分开。 2 5 -1 0 – 1 3 1 2 -2 例 3 ： 已知数列通项 ，请给出数列的前 10 项。 命令： Table[n^2,{n,1,10}] 2 n n x ? 例 4 ： 给出 30 以内的奇数。 命令： Table[n,{n,1,30,2}] 例 5 ： 生成四阶单位阵。 命令： IdentityMatrix[4] 例 6 ： 生成一个以 1,2,3,4,5 为对角元的对角矩阵， 并用 矩阵形式表示。 命令： DiagonalMatrix[{1,2,3,4,5}] MatrixForm[%] 向量和矩阵的输入 在 Mathematica 中获得表的元素的规则如下： ? 若 A 是一个向量，则 A[i] 表示向量的第 i 个元素。 ? 若 M 是一个 m 行 n 列矩阵，则用 M[[i]] 表示矩阵的 第 i 行。 ? 用 M[[i,j]] 表示第 i 行、第 j 列交叉点处的元素。 ? 用 Transpose[m][[j]] 表示 M 的第 j 列。 ? 用 ,j2}]] 表示取 M 的第 i1 、 i2 行， j1 、 j2 列构成的子矩阵。 获得表的元素 例 7 ： 构造一个 3*3 的矩阵，再取出它的元素。 命令： M=Array[a,{3,3}]; MatrixForm[%] M[[2]] M[[3,2]] Transpose[M][[3]] 3}]] 获得表的元素 ? 表的维数：用 Dimensions[list] 给出向量或矩阵的维数 例 8 ： 求向量 a=(1,2,3,4) 和矩阵 M= 命令： T={1,2,3,4} 5,6}} Dimensions[T] Dimensions[m] 1 2 3 4 5 6 表的维数和矩阵的加、减法 ? 矩阵的加、减法 在 Mathematica 中，矩阵可以表述成表，而相同维数 的表可以相加，它的和是两表对应元素相加所得的 同维的表。 例 9 ： {a1,a2,a3}+{b1,b2,b3} 例 10 ： m1=Array[a,{3,2}]; m2=Array[b,{3,2}]; MatrixForm[m1+m2] 表的维数和矩阵的加、减法 ? 向量的内积 命令格式： {a1,a2,a3}.{b1,b2,b3} ? 矩阵的乘积 例 11 ： 计算下列矩阵的乘积 命令： ,b2,b3}} ,e2}} m1 m2 注意： “ ”是 Mathematica 特有的，这种乘法不满足 交换律，当向量与矩阵相乘用“ ”时， Mathematica 能自动把向量看做行向量或列向量 a1 a2 a3 b1 b2 b3 c1 c2 d1 d2 e1 e2 ? ? ? 向量和矩阵的乘法 关于矩阵的几个常用函数 ? Inverse[M] : 求 M 的逆矩阵 ? Transpose[M] ：求 M 的转置矩阵 ? Det[M] ：方阵 M 的行列式 ? Eigenvalues[M] ：求矩阵 M 的特征值 例 12: (1). 求矩阵 的逆矩阵 (2). 求矩阵 的转置矩阵 (3). 求（ 2 ）中矩阵的行列式 (4). 求（ 2 ）中矩阵的逆矩阵 a b c d 1 2 3 4 5 6 7 8 9 关于矩阵的几个常用函数 例 13 ： 求方程组 的解 命令： 4,-7,6}} B={8,9,-5,0} Inverse[A].B//N 0 x 6 x 7 x 4 x 5 x 2 x x 2 9 x 6 x 3 x 8 x x