第3章节 MATLAB数值计算1.pptVIP

第三章 MATLAB数值计算 3.1矩阵运算 3.2稀疏矩阵 3.3 数据处理与多项式运算 3.1矩阵运算 3.1.1矩阵变换 3.1.2矩阵的逆与秩 3.1.3矩阵的特征值与特征向量 3.1.1矩阵变换 MATLAB提供了几个矩阵变换函数及调用方式: B=flipud(A);A对于水平轴沿上下方向进行列维翻转 B=fliplr(A); A对于垂直轴沿左右方向进行列维翻转 B=flipdim(A,dim); A沿着特定维进行翻转. B=rot90(A);矩阵A逆时针方向旋转90度 B=rot90(A,k);矩阵A逆时针方向旋转k*90度(k为整数) 例 A=[1,2;3,4;5,6]; B=fliplr(A); B = 2 1 4 3 6 5 B=flipud(A); B = 5 6 3 4 1 2 B=rot90(A); B = 2 4 6 1 3 5 B=rot90(A,2); B = 6 5 4 3 2 1 3.1.2矩阵的逆与秩 3x1+ x2 - x3 = 3.6 x1+2x2+4x3 = 2.1 -x1+4x2+5x3 = -1.4 A=[3 1 -1;1 2 4;-1 4 5];b=[3.6;2.1;-1.4]; Ax=b 解为:x=A-1b x=inv(A)*b x = 1.4818 -0.4606 0.3848 矩阵的秩 函数:rank() c=rank(A); c= 3 3.1.3矩阵的特征值与特征向量 矩阵的特征值 函数:eig 格式: E= eig(A) 求A的全部特征值构成向量E [V,D]= eig(A) 先对A作相似变换,然后求矩阵A的全部特征值,构成对角阵D,并求它的特征向量以构成V的列向量 [V,D]= eig(A,’nobalance’) 直接求A的特征向量和特征值. 例: A=[1 2 3;4 5 6;7 8 9] E=eig(A) E = 16.1168 -1.1168 -0.0000 [V,D]= eig(A,’nobalance’) V = -0.2833 -1.0000 0.5000 -0.6417 -0.1104 -1.0000 -1.0000 0.7792 0.5000 D = 16.1168 0 0 0 -1.1168 0 0 0 -0.0000 求根函数 例： a=[4,9,0,-5,3,-20]; x=roots(a); x = -1.8403 + 0.6250i -1.8403 - 0.6250i 1.1422 0.1442 + 1.0668i 0.1442 - 1.0668i 3.2稀疏矩阵 定义：矩阵中含有大量的零元素 存储方式：全元素存储和稀疏存储 全元素存储：每个元素需要相同的存储空间 稀疏存储：MATLAB仅存储那些非零元素及下标,对于有较多零元素的大矩阵这种存储方式可以减少存储空间. 3.2.2创建稀疏矩阵 全元素存储与稀疏存储的转换 例: A=[0 0 1;0 0 0;0 0 2]; s=sparse(A); f=full(s); s = f= 0 0 1 (1,3) 1 0 0 0 (3,3) 2 0 0 2 直接创建稀疏矩阵 格式: s=sparse(i,j,s,m,n); 说明: i和j分别是矩阵非零元素的行和列的下标向量,s是对应行列处的非零元素值,m和n分别是矩阵的行数和列数 例: s=sparse([1,2,3],[2,3,2],[5,3,6],3,3); s = (1,2) 5 (3,2) 6 (2,3) 3 3.2.3稀疏矩阵的查看与运算 函数nnz: 返回稀疏矩阵非零元素个数 例: nnz(s) ans = 3 函数nonzeros: 返回稀疏矩阵所有非零元素值 nonzeros(s) ans = 5