MATLAB语言与程序设计第2章 数值运算.pptVIP

MATLAB采用find函数进行矩阵元素查找。一般格式为： ind=find(X) 查找满足条件X的元素，返回元素单下标。 [row,col]=find(X) 查找满足条件X的元素，返回元素二维坐标。 A=[1 3 5;7 9 11;13 15 16]; ind=find(A10) [row,col]=find(A10) MATLAB使用函数sort对矩阵元素进行排序，默认升序排列，返回排序后的矩阵。一般格式为： Y=sort(A) 对矩阵按照升序进行排列。 [Y,I]=sort(A,DIM,MODE) 当dim=1时，按列排序，当dim=2时，按行排列。MODE指定排序方式，默认为ascend，按照升序排列，参数为descend，按照降序排列。输出参数I中的元素表示B中对应元素在输入参数中的位置。 A=[3 7 5;0 4 2]; sort(A,1) sort(A,2) sort(A,descend) [Y,I]=sort(A,2) MATLAB中提供sum函数对矩阵元素求和，一般格式为： S=sum(X) 对矩阵各列元素求和。 S=sum(X,DIM) 当dim=1时，计算矩阵各列元素的和。当dim=2时，计算矩阵各行元素的和。 a=[1 2 3; 4 5 6; 7 8 9] a = 1 2 3 4 5 6 7 8 9 sum(a) ans = 12 15 18 sum(a,2) ans = 6 15 24 MATLAB中提供prod函数用于对矩阵中元素的求积，一般格式为： B=prod(A)：对矩阵各列元素求积。 B=prod(A,dim)：当dim=1时，计算矩阵各列元素的积。当dim=2时，计算矩阵各行元素的积。 A=[1 4 7 8;2 5 8 6;3 7 9 11]; p1=prod(A) p1 = 6 140 504 528 p2=prod(A,2) p2 = 224 480 2079 2.4矩阵分析 矩阵信息量获取 1.行列式 把一个n阶矩阵看作一个行列式，根据行列式计算规则求值，得到的值称为矩阵对应的行列式的值。采用det函数用来求方阵行列式值，一般格式为： B=det(A) 2.矩阵的逆 对于一个n阶矩阵A，如果存在一个与其同阶的矩阵B，使得 A*B=B*A=E 其中E为与A同阶的单位矩阵，则称B和A互为逆矩阵。MATLAB采用inv函数求一个矩阵的逆，一般格式为： B=inv(A) 3.矩阵的秩 矩阵的秩是对矩阵行或列线性无关数的评估。所谓线性无关，是指某行或列不能由其他行或列以线性组合表达，已知矩阵A，则矩阵的秩可用rank(A)求得。 4.矩阵的迹 矩阵的迹就是矩阵的主对角线上所有元素之和，用trace函数表示。一般形式为： trace(A) 5.矩阵特征值和特征向量 对于n阶矩阵A，如果存在一个以特征值组成的对角矩阵d和n阶矩阵v，使Av=vd，矩阵v的每一个列向量对应于特征值的特征向量，用[v,d]=eig(A)表示。 6.矩阵范数 矩阵范数是从向量范数引申出来的，同样具有长度的意义。一般格式为： n=norm(A) n=norm(A,p) 结果返回一个标量，就是矩阵A的最大奇异值，p值决定返回范数的类型。 矩阵分解 矩阵分解是通过把一个复杂矩阵分解为比较简单的若干矩阵乘积的组合，便于理论分析或数值计算。矩阵分解方法有特征分解﹑奇异值分解﹑cholesky分解﹑LU分解和QR分解等。 特征分解是将矩阵分解为由其特征值和特征向量表示的乘积形式。一般格式为： [V,D]=eig(A) 奇异值分解，返回一个与矩阵大小相同的对角矩阵s和两个酉矩阵u和v，且满足A=u*s*v，若A为m*n阵，则u为m*m,v为n*n，奇异值在s主对角线上，且为非负降序排列。所谓酉矩阵是它的逆矩阵等于其共轭转置矩阵。 [u,s,v]=svd(A) 第二章 数值运算 1.数据类型 基本数据类型 字符串 单元数组 结构体 基本数据类型 整数型 MATLAB支持1﹑2﹑4和8字节的有符号整数型和无符号整数型。无符号整型可以表示0和正数，有符号整型可以表示正数和负数。 clear clc a=240 b1=int8(a) %将a转换为有符号1字节整数 b2=int16(a) %将a转换为有符号2字节整数 b3=uint8(a) %将a转换为无符号1字节整数 b4=uint16(a) %将a转换为无符号2字节整数 浮点型 MATLAB有单精度和双精度两种浮点数，都是依据I