MAPLE的使用_线性代数.ppt

线性代数和矩阵运算 MAPLE 的使用(三) 一向量及矩阵运算需调用线性代数软件包 linalg ? with(linalg ):?> v1:=vector([1,2,3]); v2:=vector([4,5,6]); > evalm(2 * v1-v2/5); > dotprod(v1,v2); > crossprod(v1,v2); 矩阵的输入有三种方式: array, matrix 和矩阵模块? A:=array(1..3,1..4,[[1,2,3,4], [5,6,7,8],[9,10,11,12]]); ? A:=array([[1,2,3,4],[5,6,7,8], [9,10,11,12]]); >A :=array([[1,2,3,4],[5,6,7],[9, 10]]); ? A:=matrix(3,4,[[1,2,3,4], [5,6,7,8],[9,10,11,12]]); ? A:=matrix([[1,2,3,4], [5,6,7],[9,10]]); 特殊矩阵> a:=array(1..3,1..3,sparse); > print(a); > b:=array(1..3,1..3,diagonal); > print(b); > c:=array(1..3,1..3,identity); > print(c); > with(linalg ): ? a:=array([[1,2,3],[2,3,4], [4,5,6]]); > swaprow(a,1,2); > swapcol(a,2,3); > delrows(a,1..1); > delcols(a,3..3); 向量及矩阵运算> mulrow(a,1,5); mulcol(a,3,x); > addrow(a,1,2,1); addcol(a,1,3,x); 倍法变换消法变换矩阵相关计算 a:=array([[1,2,3],[2,3,4], [4,5,6]]); ? b:=matrix([[2,-1,1],[1,1,1], [3,-2,1]]); > c:=matrix([[2,1],[1,2],[3,5]]); > evalm(a+b ); > evalm(a-b ); > evalm(4 * c); 矩阵相关计算> evalm(a &* c); > multiply(a,a,a); > transpose(c); > det(a ); > det(b ); > inverse(b); 矩阵相关计算> eigenvals(a ); > eigenvects(a ); > cond(b ); > jordan(a,'p '); > print(p);