MATLAB7.x基础教程第4章数学运算基础.ppt

4.1.1 矩阵分析　　MATLAB提供的矩阵分析函数如表4-1所示。 　　1．向量和矩阵的范数　　对于一个n维向量，可以用一个数?||x||?来度量该n维向量的大小。如果?||x||?满足以下三个性质，则称?||x||?为向量x的范数。　　● 非负性：对于一切x都有　　　||x||??0，且?||x||?=?0的充要条件是x=0 (4-1)　　● 正齐性：对任何实数α和向量x，有　　　　　　　||αx||?=?|α|?||x|| (4-2)　　● 三角不等式：对任何向量x和y，有　　　　　　　||x+y||≤||x||?+?||y|| ? 4-3) 　　MATLAB中，常用norm函数计算向量或矩阵的1-范数、2-范数和?∞-范数，其调用格式为norm(x,p)或norm(A,p)，作用是计算向量x或矩阵A的p-范数。其中，p为1、2或inf。命令norm(x)或norm(A)相当于命令norm(x,2)或norm(A,2)，而normest(A)实现矩阵A的2-范数的快速估算。　　【例】 求向量的1-范数、2-范数和∞-范数。　　在命令窗输入：　　　 x = [1 2 3];　　　 [norm(x,1),norm(x,2),norm(x,inf)] 　　运行结果：　　　ans =　　　 6.0000 3.7417 3.0000　　【例】 求矩阵的1-范数、2-范数和∞-范数。　　在命令窗输入：　　　 A = [1 2 3　　　　　 4 5 6　　　　　 7 8 9];　　　 [norm(A,1), norm(A,2), norm(A,inf)]　　运行结果：　　ans =　　 18.0000 16.8481 24.0000 　　2．矩阵的秩　　对于一个矩阵A，如果A=0，则A的秩为零；如果A≠0，则称A中非零子式的最高阶数为A的秩。MATLAB中常用命令rank(A)来以默认允许误差计算矩阵A的秩，而用命令rank(A,tol)来以给定容许误差tol计算矩阵的秩。　　【例】 求矩阵的秩。　　在命令窗输入：　　　 A = [1 2 3　　　 4 5 6　　　 7 8 9];　　　 rank(A) 　　运行结果：　　ans =　　 2　　3．矩阵的行列式　　当矩阵A为方阵时，可以利用命令det(A)来计算A的行列式。MATLAB中行列式的计算是通过高斯消去法先实现矩阵的LU分解，然后再计算下三角矩阵L和上三角矩阵U的行列式之积(对角线元素之积)，即A的行列式。 　　【例】 求矩阵的行列式。　　在命令窗输入：　　　 A = [1 3 2　　　　　 6 5 4　　　　　　 7 8 9];　　　 det(A)　　运行结果：　　ans =　　 -39 　　4．矩阵的迹　　矩阵的迹定义为主对角线的元素之和。无论矩阵是否为方阵，MATLAB中都可以用命令trace(A)来计算矩阵A的迹。　　【例】 求矩阵的迹。　　在命令窗输入：　　 A = [1 3 2　　　　 6 5 4];　　 trace(A)　　运行结果：　　　ans =　　　 6 　　5．化零矩阵　　矩阵A的化零矩阵Z满足A*Z的元素近似为零，并且Z‘*Z=I。MATLAB中可以用命令Z = null(A)求矩阵A的化零矩阵，用命令Z = null(A,’r‘)求矩阵A的有理形式的化零矩阵，有理形式的化零矩阵Z’*Z≠I。　　【例】 求矩阵的化零矩阵。　　在命令窗输入：　　 A = [1 2 3　　　　 3 2 1　　　　 1 2 3];　　 Z = null(A);　　 A*Z 　　运行结果：　　ans =　　 1.0e-015 *　　　 0　　 -0.6661　　　 0 　　【例】 求矩阵的有理形式的化零矩阵。　　在命令窗输入：　　 A = [1 2 3　　　　 3 2 1