第五讲矩阵计算及其应用.docVIP

第五讲 矩阵计算及其应用 〖实验目的〗 矩阵计算是线性代数中最基本的运算形式，矩阵的有些计算很很枯燥也很容易出错.这一讲将用Mathematic软件完成矩阵的一些基本计算并通过试验提高对矩阵应用的理解,这样既可以把重复繁杂的工作交给程序进行,同时提高发现问题和解决问题的能力. 〖实验内容〗 这一讲包括两个部分的内容:矩阵计算的基本函数；四个数学试验. 一、矩阵计算的有关函数 (一)矩阵的生成函数 1. Table可以生成多层表,一层表是行向量(行矩阵),二层表是矩阵. 输入: 输出: {1,4,9,16,25,36,49,64,81,100} 输入: 输出: {-27,-8,-1,0,1,8,27} 输入: TableForm[Table[Binomial[i,j],{i,0,8},{j,0,8}]] 输出: 输入: TableForm[Table[Binomial[i,j],{i,0,8},{j,0,i}]] 输出: 2. Array具有和Table相似的用途. 输入: Array[a,3] 输出: {a[1],a[2],a[3]} 输入: Array[a,{2,3},0] {{a[0,0],a[0,1],a[0,2]},{a[1,0],a[1,1],a[1,2]}} 输入: Array[a,{2,3}] 输出: {{a[1,1],a[1,2],a[1,3]},{a[2,1],a[2,2],a[2,3]}} 输入: Array[a,{2,3},1,Plus] 输出: a[1,1]+a[1,2]+a[1,3]+a[2,1]+a[2,2]+a[2,3] 输入: Array[a,{2,3},{10,20}] 输出: {{a[10,20],a[10,21],a[10,22]},{a[11,20],a[11,21],a[11,22]}} 3. DiagonalMatrix生成对角矩阵 输入: IdentityMatrix[3] 输出: {{1,0,0},{0,1,0},{0,0,1}} 输入: L1={1,-2,3,4};A=DiagonalMatrix[L1] {MatrixQ[L1],VectorQ[L1]} 输出: {{1,0,0,0},{0,-2,0,0},{0,0,3,0},{0,0,0,4}} {False,True} 输入: ma=Table[m[i,j],{i,3},{j,4}] Dimensions[ma] Tr[ma] MatrixPower[A,n]//MatrixForm 输出: {{m[1,1],m[1,2],m[1,3],m[1,4]},{m[2,1],m[2,2],m[2,3],m[2,4]},{m[3,1],m[3,2],m[3,3],m[3,4]}} {3,4} m[1,1]+m[2,2]+m[3,3] 4. SparseArray生成稀疏矩阵 输入: 输出: SparseArray[ 6 ,{5,5}] SparseArray[ 6 ,{5,5}] True SparseArray[ 36 ,{10,10}] 5. RandomSquare函数可以自动生成一个随机方阵 输入: Unprotect[RandomSquare]; Clear[RandomSquare]; RandomSquare[dim_Integer:3,max_:5,detmax_Integer:3,del_List:{0}]:=Module[{t,tt,i,j},tt=Table[m,{m,-detmax,detmax}]~Complement~del; t=Table[Random[Integer,{-max,max}],{dim},{dim}]; While[FreeQ[tt,Det[t]],i=Random[Integer,{1,dim}]; j=Random[Integer,{1,dim}]; t[[i,j]]=Random[Integer,{-max,max}]]; t]/;Positive[dim]NumberQ[max]Positive[max]Positive[detmax] Protect[RandomSquare]; RandomSquare[3] 输出(具有随机性): {{5,4,2},{2,5,4},{3,-2,-3}} 思考题： （1）编写实现随机矩阵（不一定是方阵）的函数. （2）编写产生随机向量的函数. （二） 扩展程序包 LinearAlgebra扩展程序中包括: MatrixManipulation,Orthogonalization,Tridiagonal等. 输入