资源描述：
《maple数学软件3》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数学软件3.2矩阵1.用array定义和输入格式为：A:=array(1..m,1..n)定义m*n的空矩阵A:=array(1..m,1..n,list)A:=array(list)如:A:=array(1..3,1..4);B:=array(1..3,1..4,[[1,2,3,4],[5,6,7,8],[9,10,11,12]]);C:=array([[1,2,3,4],[5,6,7,8],[9,10,11,12]]);1、矩阵的定义和输入2.用matrix定义和输入格式为：A:=matrix(m,n)定义m*n的空矩阵A:=matrix(m,n,list)A:=matrix(lis
2、t)3.用矩阵模板定义和输入矩阵的定义和输入可以用矩阵模板来实现。如:A:=matrix(3,4);B:=matrix(3,4,[[1,2,3,4],[5,6,7,8],[9,10,11,12]]);C:=matrix([[1,2,3,4],[5,6,7,8],[9,10,11,12]]);1.零矩阵所有元素全为零的矩阵，称为零矩阵。在Maple中，用函数sparse定义零矩阵。如:a:=array(1..2,1..2,sparse):print(a);2、特殊矩阵2.对称矩阵Maple中用symmetric定义对称矩阵。如：restart;A:=array(symmetric,1..4
3、,1..4):A[1,1]:=1;A[1,2]:=2,A[1,3]:=3,A[1,4]:=4;print(A);3.对角矩阵Maple中用diagonal定义对角矩阵。如：restart;a:=array(1..3,1..3,diagonal):a[2,2]:=1;print(a);4.单位矩阵Maple中用identity定义单位矩阵。如：restart;a:=array(1..3,1..3,identity);print(a);3、矩阵操作1.交换行列with(linalg);a:=array([[1,2,x],[3,4,y]]);swaprow(a,1,2);swapcol(a,2
4、,3);2.行列乘以因子格式：mulcol(A,r,expr)mulrow(A,r,expr);如：restart;A:=array([[1,2],[3,4]]);mulrow(A,2,2);mulcol(A,2,x);3.行(列)相加格式：addrow(A,r1,r2,m)将r1行乘以m加到r2行addcol(A,c1,c2,m).如：restart;A:=array([[1,2,3],[2,3,4],[3,4,5]]);addrow(A,1,2,10);addcol(A,1,2,-x);4.删除行(列)格式：delrows(A,r..s)删除第r到s行delcols(A,r..s).
5、如：restart;A:=matrix(3,3,[1,2,3,4,5,6,7,8,9]);delrows(A,2..3);delcols(A,1..1);4、矩阵运算1.矩阵相等A:=matrix([[x,x+y],[z,3]]):B:=matrix([[2,6],[z*z,w]]):eq1:=A[1,1]=B[1,1]:eq2:=A[1,2]=B[1,2]:eq3:=A[2,1]=B[2,1]:eq4:=A[2,2]=B[2,2]:solve({eq1,eq2,eq3,eq4});2.矩阵加法格式：evalm(A+B)matadd(A,B)(要调用linalg程序包)如：restart
6、;A:=matrix([[x,y],[z,w]]):B:=matrix([[a,b],[c,d]]):C:=evalm(A+B);with(linalg):D1:=matadd(A,B);3.矩阵的数乘格式：evalm(k*A)with(linalg):scalarmul(A,k)如：restart;A:=array([[1,2,3],[2,3,4],[3,4,5]]);evalm(2*A);with(linalg):'4*A'=scalarmul(A,4);4.矩阵的乘法格式：evalm(A&*B)with(linalg):multiply(A,B)如：restart;A:=array
7、([[1,0,3,2],[2,1,0,2]]);B:=matrix([[5,4,2],[-1,1,3],[2,0,2],[0,3,4]]);evalm(A&*B);with(linalg):'A*B'=multiply(A,B);5.矩阵的转置格式：with(linalg):transpose(A)如：A:=array([[1,0,3,2],[2,1,0,2]]);with(linalg);B:=transpose(A);multi
