专题1-MATLAB求解方程.pptVIP

MATLAB 求解方程 制作: 陈学明 基本概念 稀疏矩阵(sparse matrix)：矩阵中非零元素的个数远远小于矩阵元素的总数），并且非零元素的分布没有规律。与之相区别的是，如果非零元素的分布存在规律（如上三角矩阵、下三角矩阵、对称矩阵），则称该矩阵为特殊矩阵。 稠密矩阵:非0元素占所有元素比例较大的矩阵。 若n阶矩阵A的行列式不为零，即 |A|≠0，则称A为非奇异矩阵，否则称A为奇异矩阵。 把矩阵A的行换成相应的列，得到的新矩阵称为A的转置矩阵，记作A‘。 AA‘=E（E为单位矩阵，）或A’A=E，则n阶实矩阵 A称为正交矩阵 。 基本概念 求矩阵A的秩 rank(A) 求矩阵A的迹 trace(A) 求矩阵A的行列式 det(A) 求矩阵V 的1范数 norm(V,1) 求矩阵V 的2范数 norm(V)或norm(V,2) 求矩阵V 的∞范数 norm(V,inf) 魔方矩阵 魔方矩阵是每行、每列及两条对角线上的元素和都相等的矩阵。对于n阶魔方阵，其元素由1,2,3,…,n2共n2 个整数组成. magic(n) ：生成一个n阶魔方阵,各行各列及两条 对角线和为(1+2+3+...+n2 )/n 范得蒙矩阵 范得蒙(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 希尔伯特矩阵 Hilbert矩阵的每个元素hij=1/(i+j-1) hilb(n) 生成n阶希尔伯特矩阵 invhilb(n) 求n阶的希尔伯特矩阵的逆矩阵 注意1：高阶Hilbert矩阵一般为病态矩阵，所以直接求逆可能出现错误结论，故应该采用invhilb(n) 注意2：由于Hilbert矩阵本身接近奇异，所以建议处理该矩阵时建议尽量采用符号运算工具箱，若采用数值解时应该验证结果的正确性 托普利兹矩阵 (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 帕斯卡矩阵 二次项(x+y)n展开后的系数随n的增大组成一个三角形表，称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。 pascal(n) 生成一个n阶帕斯卡矩阵。 矩阵分解 矩阵分解通过将复杂矩阵表示成形式简单或具有良好数学性质（统称为简单矩阵）的组合，以便于理论分析或数值计算。通常矩阵分解将复杂矩阵分解为几个简单矩阵的乘积。求解线性方程组不可避免地要用到矩阵分解的概念。 MATLAB中，线性方程组的求解主要用到三种基本的矩阵分解，即对称正定矩阵的cholesky分解、一般方程的gaussian消去法和矩阵的正交分解。这三种分解由函数chol、lu和qr完成。 矩阵的逆 求方阵A的逆可用 inv(A) 说明1:对于Hilbert求逆时，不建议用inv,可用 invhilb直接产生逆矩阵 说明2：符号工具箱中也对符号矩阵定义了inv( )函数,即使对更高阶的奇异矩阵也可以精确的求解出逆矩阵 说明3：对于奇异矩阵用数值方法的inv( )函数，会给出错误的结果，但符号工具箱中的inv( )会给出正确的结论