第5章 数值计算功.pptVIP

5.2.2转换成满矩阵full(S)函数可以将稀疏矩阵S转换成一个满矩阵。【例5-16】将例5-15中的结果再转换成满矩阵。full(as)ans=100010001第29页，共78页，星期日，2025年，2月5日5.2.3对非零元素进行操作nnz(A)函数可以求矩阵A中非零元素的个数，它既可求满矩阵也可求稀疏矩阵。spy(A)函数可以画出稀疏矩阵A中非零元素的分布。也可用在满矩阵中，在这种情况下，只给出非零元素的分布。spy(A，cstr，size)函数可以用指定的颜色cstr在size规定的范围内画出稀疏矩阵A中非零元素的分布。nonzeros(A)函数可以按照列的顺序找出矩阵A中非零的元素。spones(A)函数可以把矩阵A中的非零元素全换为1。nzmax(A)函数可以给出为矩阵A中非零元素分配的内存数.不一定和nnz(A)得到的数相同。issparse(A)函数在矩阵A是稀疏矩阵的情况下返回1;否则返回0。spfun(fcn，A)函数可以用A中所有非零元素对函数fcn求值，如果函数不是对稀疏矩阵定义的，同样也可以求值。spfun(A)函数可以求稀疏矩阵A的结构秩。对于所有的矩阵来说，都有sprank(A≥rank(A)。第30页，共78页，星期日，2025年，2月5日【例5-17】求3阶单位矩阵的非零元素的个数，并画出非零元素的分布。nnz(eye(3))ans=3spy(eye(3))结果如图5.1所示。图5.1矩阵A的稀疏表示第31页，共78页，星期日，2025年，2月5日5.2.4稀疏矩阵的特例MATLAB中有四个基本稀疏矩阵，它们是单位矩阵，随机矩阵，对称随机矩阵和对角矩阵。下面介绍各自的生成函数。speye(n)函数可以生成n×n的单位稀疏矩阵。speye(m，n)函数可以生成m×n的单位稀疏矩阵。sprand(A)函数可以生成与A有相同结构的随机稀疏矩阵，且元素服从均匀分布。sprand(m，n，dens)函数可以生成一个m×n的服从均匀分布的随机稀疏矩阵，有dens×m×n个非零元素，0≤dens≤1.参数dens是非零元素的分布密度。sprandn(A)函数可以生成与A有相同结构的随机稀疏矩阵，且元素服从正态分布。sprandsym(S)函数可以生成一个随机对称稀疏矩阵。它的下三角及主对角线部分与S的结构相同，矩阵元素服从正态分布。sprandsym(n，dens)函数可以生成一个m×n的随机对称稀疏矩阵.矩阵元素服从正态分布，分布密度为dens。[B，d]=spdiags(A)函数可以求出A中所有的对角元，对角元保存在矩阵B中，它们的下标保存在向量d中。spdiags(A，d)函数可以生成一个矩阵，这个矩阵包含有矩阵A中向量d规定的对角元。spdiags(B，d，A)函数可以生成矩阵A，用矩阵B中的列替换d定义的对角元。A=spdiags(B，d，m，n)用保存在由d定义的B中的对角元创建稀疏矩阵A。第32页，共78页，星期日，2025年，2月5日表5.6给出了处理稀疏矩阵的常用函数。表5.6常用处理稀疏矩阵的函数函数说明full将稀疏矩阵转换为完全矩阵nnz返回矩阵的非零元素数nzmax返回为非零元素分配的储存量sparse建立稀疏矩阵或将完全矩阵转换为稀疏矩阵sprand建立稀疏均匀的随机分布矩阵issparse测试矩阵是否为稀疏矩阵nonzeros返回矩阵的非零元素spalloc为稀疏矩阵分配存储空间speye建立一个稀疏等同矩阵第33页，共78页，星期日，2025年，2月5日5.3MATLAB与线性代数5.3.1矩阵的特征参数1求矩阵的行列式行列式对于查明一个方程组是否有解很有用。行列式是一个特殊的方形阵列，并且还可以简化为一个数。用户使用竖条来表示一个行列式，而使用方括号来代表矩阵。将一个22的行列式简化为一个数的规则如下所示：数学上通过人笔算来得出行列式的值，现在可以通过MATLAB利用det函数轻松算出行列式的值。第34页，共78页，星期日，2025年，2月5日【例5-18】计算如下的行列式：D=[345;123;369];det(D)ans=0第35页，共78页，星期日，2025年，2月5日2矩阵的秩矩阵的秩是矩阵行或列的数值线性独立的度量。如果一个向量线