线性代数试验.docVIP

第4章 线性代数 【教学目的】使学生学会使用MATLAB计算行列式、求逆矩阵、求矩阵的特征值和特征向量；会用MATLAB进行矩阵的初等变换、求线性方程组的解，会求几个特殊矩阵方程的解。 【教学重点、难点】矩阵运算、线性方程组求解 4.1矩阵的生成 4.1.1 4.1.2通过外部数据加载 4.1.3在M文件中创建矩阵 4.1.4通过函数产生矩阵 线性代数中有若干特殊意义的矩阵，在Matlab中可以很容易的通过函数的方式创建它们，见表4.1。其中：matlab函数名必须小写。 表4.1 特殊矩阵函数表 特殊矩阵名称 函数命令格式 功能简介 全零矩阵 zeros(m,n) zeros(n) zeros(size(A)) 创建m行n列的全零矩阵； 创建n行n列的全零矩阵； 产生一个与矩阵A同样大小的零矩阵 全1矩阵 ones(m,n) 创建m行n列的全1矩阵 单位矩阵 eye(m,n) 创建m行n列的对角线为1的矩阵 随机矩阵 rand(m,n) 创建m行n列的(0,1)均匀分布的随机矩阵 正态分布随机矩阵 randn(m,n) 产生均值为0，方差为1的标准正态分布随机矩阵 魔方矩阵 magic(n) 生成一个n阶方阵（其元素由1,2,…n2整数组成，其每行、每列及两条对角线上的元素和都相等） Vandermonde 矩阵 vander(V) 生成以向量V为基础向量的范得蒙矩阵（矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积） Hilbert矩阵 hilb(n) 生成希尔伯特矩阵（第(i,j)个元素为1/(i+j-1)） Hilbert逆矩阵 invhilb(n) 求n阶的希尔伯特矩阵的逆矩阵 Toeplitz矩阵 toeplitz(x,y) 生成一个以向量x为第一列，向量y为第一行的托普利兹矩阵（除第一行第一列外，其他每个元素都与左上角的元素相同） Hankel 矩阵 hankel(c) hankel(c,r) 生成一个以向量c为第一列，其它元素是0； 生成一个以向量c为第一列，最后一行是r向量，其中第一个元素以c的为准，r向量是从第二个元素开始的（汉克尔矩阵是反对角线上元素相等矩阵） 伴随矩阵 compan(p) p是一个多项式的系数向量，高次幂系数排在前，低次幂排在后 帕斯卡矩阵 pascal(n) 生成一个n阶帕斯卡矩阵（由杨辉三角形表组成的矩阵） Hadamard矩阵 hadamard(n) n，n/12或，n/20必须是2的乘方。 （哈达玛矩阵是由+1和-1元素构成的正交方阵） 空矩阵 [] 创建一个空矩阵 示例4.1 求特殊矩阵 (1) 在区间[20,50]内均匀分布的5阶随机矩阵。 x=20+(50-20)*rand(5) (2) 均值为0.6、方差为0.1的5阶正态分布随机矩阵。 y=0.6+sqrt(0.1)*randn(5) (3) 求4阶希尔伯特矩阵及其逆矩阵。 format rat %以有理形式输出H=hilb(4) H=invhilb(4) (4) 已知c = 1:3; r = 7:10;求其Hankel 矩阵。 c = 1:3; r = 7:10; hankel(c,r) (5) 求(x+y)5的展开式。 pascal(6) 矩阵次对角线上的元素1,5,10,10,5,1即为展开式的系数。 4.1.5符号矩阵的生成 表4.2 符号矩阵生成 生成方法 命令格式 功能简介 直接创建 sym([ , , , ]) 矩阵元素可以是任何的符号变量、符号表达式及方程，且元素的长度可以不同。 间接创建 syms 在创建符号矩阵之前，先把符号矩阵的 所有变量定义为符号变量，然后按创建普通矩阵的格式输入矩阵 由数值矩阵转化为符号矩阵 B=sym(A) 将一个数值矩阵A转化为符号矩阵B 注：数值矩阵A转换成符号矩阵B，可以通过whos命令来详细查看这两类矩阵的区别。 4.2矩阵操作 4.2.1元素变换操作 有已知矩阵A，由A的元素构成的各种矩阵分别如表4.3所示。 表4.3 矩阵元素变换操作表 元素操作名称 操作命令格式 功能简介 对角阵 diag(A); diag(A,k); diag(V); diag(V,k) 提取矩阵A主对角线元素产生一个具有min(m,n)个元素的列向量； 提取第k条对角线的元素，主对角线为0，向上为第1，2…条，向下为-1，-2…条； 产生一个m×m对角矩阵，其主对角线元素即为向量V的元素； 产生一个n×n(n=m+|k|)对角阵，其第k条对角线的元素即为向量V的元素 （只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。） 上三角阵 triu(A); triu(A,k)