第三章 ATLAB线性方程组.pptVIP

第三章 线性代数方程组及矩阵特征值;§ 3.0预备知识; ;二、 用于专门学科中的矩阵 (1) 魔方矩阵 魔方矩阵是每行、每列及两条对角线上的元素和都相等的矩阵。对于n阶魔方阵，其元素由1,2,3,…,n2共n2 个整数组成. magic(n) 生成一个n阶魔方阵,各行各列及两条 对角线和为(1+2+3+...+n2 )/n ;例 magic(5) ans = 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 例：将101~125等25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。命令如下： M=100+magic(5);(2) 范得蒙矩阵 范得蒙(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。 vander(V) 生成以向量V为基础向量的范得蒙矩阵。 例 A=vander([1;2;3;5]) A = 1 1 1 1 8 4 2 1 27 9 3 1 125 25 5 1; (3) 希尔伯特矩阵（Hilbert）  Hilbert矩阵的每个元素hij=1/(i+j-1) hilb(n) 生成n阶希尔伯特矩阵 invhilb(n) 专求n阶的希尔伯特矩阵的逆矩阵 注意1：高阶Hilbert矩阵一般为病态矩阵，所以直接求逆可能出现错误结论，故应该采用invhilb(n) 注意2：由于Hilbert矩阵本身接近奇异，所以建议处理该矩阵时建议尽量采用符号运算工具箱，若采用数值解时应该验证结果的正确性; (4) 托普利兹矩阵 (toeplitz) toeplitz矩阵除第一行第一列外，其他每个元素都与左上角的元素相同。 toeplitz(x,y) 生成一个以x为第一列，y为第一行的托普利兹矩阵。这里x, y均为向量，两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。 例 T=toeplitz([1:5],[1:7]) T = 1 2 3 4 5 6 7 2 1 2 3 4 5 6 3 2 1 2 3 4 5 4 3 2 1 2 3 4 5 4 3 2 1 2 3; (5) 帕斯卡矩阵 二次项(x+y)n展开后的系数随n的增大组成一个三角形表，称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。 pascal(n) 生成一个n阶帕斯卡矩阵。; ; ;常用的向量范数:;范数意义下的单位向量: X=[x1, x2]T;常用的矩阵范数:; ; ;例：Hilbert矩阵（非常有名的病态矩阵）：;;说明1:对于Hilbert求逆时，不建议用inv,可用 invhilb直接产生逆矩阵 说明2：符号工具箱中也对符号矩阵定义了inv( )函数,即使对更高阶的奇异矩阵也可以精确的求解出逆矩阵;说明3：对于奇异矩阵用数值方法的inv( )函数，会给出错误的结果，但符号工具箱中的inv( )会给出正确的结论;但用符号工具箱的inv:; ; ;DD=det(D); x1=det(D1)/DD; x2=det(D2)/DD; x3=det(D3