第2章 数值计算功能.pptVIP

spones(A)函数可以把矩阵A中的非零元素全换为1。nzmax(A)函数可以给出为矩阵A中非零元素分配的内存数.不一定和nnz(A)得到的数相同。issparse(A)函数在矩阵A是稀疏矩阵的情况下返回1;否则返回0。spfun(fcn，A)函数可以用A中所有非零元素对函数fcn求值，如果函数不是对稀疏矩阵定义的，同样也可以求值。sprank(A)函数可以求稀疏矩阵A的结构秩。对于所有的矩阵来说，都有sprank(A）≥rank(A)。第29页，共66页，星期日，2025年，2月5日【例2-16】求3阶单位矩阵的非零元素的个数，并画出非零元素的分布。nnz(eye(3))ans=3spy(eye(3))结果如图2.1所示。图2.1矩阵A的稀疏表示第30页，共66页，星期日，2025年，2月5日2.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有相同结构的随机稀疏矩阵，且元素服从正态分布。第31页，共66页，星期日，2025年，2月5日sprandsym(S)函数可以生成一个随机对称稀疏矩阵。它的下三角及主对角线部分与S的结构相同，矩阵元素服从正态分布。sprandsym(n，dens)函数可以生成一个n×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页，共66页，星期日，2025年，2月5日2.3MATLAB与线性代数2.3.1矩阵的特征参数1.求矩阵的行列式行列式对于查明一个方程组是否有解很有用。行列式是一个特殊的方形阵列，并且还可以简化为一个数。用户使用竖条来表示一个行列式，而使用方括号来代表矩阵。将一个2×2的行列式简化为一个数的规则如下所示：数学上通过人笔算来得出行列式的值，现在可以通过MATLAB利用det函数轻松算出行列式的值。第33页，共66页，星期日，2025年，2月5日【例2-17】计算如下的行列式：D=[345;123;369];det(D)ans=0第34页，共66页，星期日，2025年，2月5日2.矩阵的秩矩阵的秩是矩阵行或列的数值线性独立的度量。如果一个向量线性独立于另外一些向量组，那意味着这一个向量不能写成它们的线性组合。rank函数可以实现求矩阵的秩。对于带有n个未知量的m个线性系统方程：Ax=b把b连结A上构成了增广矩阵：[Ab]当且仅当rank(A)=rank(Ab)时系统有解。如果秩等于n，那么系统有唯一解，但如果秩小于n，那么系统有无数解。如果秩条件吻合并且秩与未知数个数相等，解就可以用左除计算得到。下面举例具体说明以下。第35页，共66页，星期日，2025年，2月5日【例2-18】求下列方程的解。A=[1-21;345;-217];b=[12;20;11];C=[Ab]%创建增广矩阵C=1-211234520-21711第36页，共66页，星期日，2025年，2月5日现在我们检查一下A的秩：rank(A)ans=3增广矩阵的秩为：rank(C)ans=3由于秩相同，因此解存在。这里有三个未知量，我们也注意到秩r满足r=n。这意味着解唯一。我们用左除求得解：x=A\bx=4.3958-2.22923.1458第37页，共66页，星期日，2025年，2月5日3.矩阵的特征值和特征向量矩阵特征值与特征向量的定义：设A是n阶矩阵，如果存在