MATLAB语言基础 ch6数值计算.pptVIP

ch6 数值计算 5.1 数值微积分 5.1.1 近似数值极限和导数 5.1.2 数值求和与近似数值积分 5.2 矩阵和代数方程 5.2.1 矩阵运算和特征参数 5.2.2 矩阵的变换和特征值分解 5.2.3 线性方程的解 5.2.4 一般代数方程的解 5.1.1 近似数值极限和导数 （1）有限精度浮点的离散本质 （2）自变量的增量选取一定要大于原数据相对精度的10倍以上 （3）尽量使用Matlab提供的现成指令 Dx=diff(X) 求差分 FX=gradient(F) 求一元函数梯度 [FX,FY]=gradient(F) 求二元函数梯度 5.1.2 数值求和与近似数值积分 Sx=sum(X) 沿列方向求和 Scs=cumsum(X) 沿列方向求累计和 St=trapz(x,y) 采用梯形法沿列方向求函数y关于x的积分 Sct=cumtrapz(x,y) 采用梯形法沿列方向求函数y关于x的累计积分 5.2.1 矩阵运算和特征参数 1. 矩阵运算 2. 矩阵的标量特征参数 秩 rank(A) 迹 trace(A) 行列式 det(A) 5.2.2 矩阵的变换和特征值分解 [R,ci]=rref(A) 借助初等变换把A变换成行阶梯矩阵R X= null(A) A矩阵零空间的全部正交基，满足AX=0 Z=orth(A) A矩阵空间的全部正交基，满足span(Z)=span(A) [V,D]=eig(A) A矩阵的特征值、特征向量分解，使AV=AD 5.2.3 线性方程的解 1 直接解法 对于线性方程组Ax=b，可以利用左除运算符“\”求解： x=A\b 例1 用直接解法求解下列线性方程组。 命令如下： A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; b=[13,-9,6,0]; x=A\b 数据分析 1、查取最大值 MAX函数的命令格式有： [Y,I]= max (X)：将max(X)返回矩阵X的各列中的最大元素值及其该元素的位置赋予行向量Y与I；当X为向量时，则Y与I为单变量。 [Y,I]=max(X,[],DIM)：按数组X的第DIM维的方向查取其最大的元素值及其该元素的位置赋予向量Y与I。 例查找下面数列x的最大值。 x=[3 5 9 6 1 8] % 产生数列x x = 3 5 9 6 1 8 y=max(x) % 查出数列x中的最大值赋予y y = 9 [y,l]=max(x) % 查出数列x中的最大值及其该元素的位置赋予y,l y = 9 l = 3 例分别查找下面3×4的二维数组x中各列和各行元素中的最大值。 x=[1 8 4 2;9 6 2 5;3 6 7 1] % 产生二维数组x x = 1 8 4 2 9 6 2 5 3 6 7 1 y=max(x) % 查出二维数组x中各列元素的最大值产生赋予行向量y y = 9 8 7 5 [y,l]=max(x) % 查出二维数组x中各列元素的最大值及其这些 % 元素的行下标赋予y,l y = 9 8 7 5 l = 2 1 3 2 [y,l]=max(x,[ ],1) % 本命令的执行结果与上面命令完全相同 y = 9 8 7 5 l = 2 1 3 2 [y,l]=max(x,[ ],2) % 由于本命令中DIM=2,故查找操作在各行中进行 y = 8 9 7 l = 2 1 3 例试取下面两个2×3的二维数组x、y所有同一位置上的元素值大者构成一个新矩阵p。 x=[4 5 6;1 4 8] % 产生二维数组x x = 4