第二章数值计算课件.pptVIP

MATLAB 的数值计算 制作: 陈学明 矩阵的拆分 1 利用冒号表达式获得子矩阵 ① A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列的元素。 ② A(i:i+m,:)表示取A矩阵第i～i+m行的全部元素；A(:,k:k+m)表示取A矩阵第k～k+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第i～i+m行内，并在第k～k+m列中的所有元素。此外，还可利用一般向量和end运算符来表示矩阵下标，从而获得子矩阵。end表示某一维的末尾元素下标。 矩阵的拆分 2 利用空矩阵删除矩阵的元素 在MATLAB中，定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。注意，X=[]与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为0。 特殊矩阵 1．通用的特殊矩阵常用的产生通用特殊矩阵的函数有：zeros：产生全0矩阵(零矩阵)。ones：产生全1矩阵(幺矩阵)。eye：产生单位矩阵。rand：产生0～1间均匀分布的随机矩阵。randn：产生均值为0，方差为1的标准正态分布随机矩阵。 特殊矩阵 1．通用的特殊矩阵常用的产生通用特殊矩阵的函数有：zeros：产生全0矩阵(零矩阵)。ones：产生全1矩阵(幺矩阵)。eye：产生单位矩阵。rand：产生0～1间均匀分布的随机矩阵。randn：产生均值为0，方差为1的标准正态分布随机矩阵。 特殊矩阵 2．用于专门学科的特殊矩阵 (1) 魔方矩阵魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵，其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n)，其功能是生成一个n阶魔方阵。 特殊矩阵 2．用于专门学科的特殊矩阵 (2) 范得蒙矩阵范得蒙(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中，函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如，A=vander([1;2;3;5])即可得到上述范得蒙矩阵。 特殊矩阵 2．用于专门学科的特殊矩阵 (3) 希尔伯特矩阵在MATLAB中，生成希尔伯特矩阵的函数是hilb(n)。使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数invhilb(n)，其功能是求n阶的希尔伯特矩阵的逆矩阵。 矩阵的其它运算 inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 diag —— 对角矩阵 ’ —— 矩阵转置 sqrt —— 矩阵开方 矩阵的一些特殊操作 矩阵的变向 rot90:矩阵整体反时针旋转90度; rot90(a) fliplr:矩阵左右翻转; flipud:矩阵上下翻转; 矩阵的抽取 diag:抽取主对角线; tril: 抽取左下三角; triu:抽取右上三角; 稀疏矩阵 存储稀疏矩阵时，往往只存放其中的非零元。稀疏矩阵的三元组表法是顺序存储方法的一种。采用这种方法时，线性表中的每个结点对应稀疏矩阵的一个非零元素，其中包括3个字段，分别为该元素的行下标、列下标和值，结点间的次序按矩阵的行优先顺序排列。另外，用第0行的三个元素分别存储矩阵的行数、列数和非零元数目。例如，矩阵A： 5??0 0 7 可以用三元组表示为 0??2 0 0 3 4 4 0??0 0 0 1 1 5 1 4 7 2 2 2 在具体编程过程中，往往可以用一个简单的n*3二数组来表示此三元组。即稀疏矩阵A可以用数组a[5][3]={{3,4,4},{1,1,5},{1,4,7},{2,2,2}}?来表示。 稀疏矩阵的生成 在MATLAB 7中，生成稀疏矩阵用特殊的函数来进行，这些函数有speye、spones、spdiags、sparse、find、full、spalloc、sprand和sprandn 等。 稀疏矩阵与全元素矩阵的相互转换 用来将稀疏矩阵和全元素矩阵相互转换的函数有sparse、full和find等3个函数 。 注：find 函数生成非中非零元