MATLABch4_数值计算(阅读).pptVIP

第4章MATLAB的数值计算功能 — matlab 具有出色的数值计算能力，占据世界上数值计算软件的主导地位 数值微积分 矩阵与代数方程 数据统计分析 多项式与插值 数值求和与数值积分 sum（X）按列求和 cumsum （X） 按列求元素累计和 trapz(x , y)给定数据点x和y，计算y=f(x)下的梯形面积积分。 cumtrapz(x , y) 给定数据点x和y，计算y=f(x)下的梯形面积累计积分。 函数极值的数值求解 常微分方程的数值求解 单个高阶常微分方程处理方法 4.2 矩阵与线性代数-矩阵运算 加(+)、减(-) 乘：(*） 矩阵之间的乘、向量与矩阵相乘、标量与矩阵相乘 除：右除（/）、左除（\） 数学上没有矩阵除法的定义。 幂 ^ 转置 ‘ 其他： 指数运算 Y = expm(A) 还有:expm1、expm2、expm3 开方 Y = sqrtm(A) 对数 Y = logm(A) 行列式和逆 行列式：d = det(A) 计算方阵的行列式值 矩阵的逆：C= inv(A) AC=CA=I ，则：C=A-1 其中A是非奇异方阵 实际数值计算中较少使用 矩阵特征值 特征方程：Ax=λx由|A-λI|=0 求特征值λi由Ax=λ i x 求特征向量x i eig求解Ax=λx特征值问题 d = eig(A) 只返回特征值向量d [V,D] = eig(A) 返回特征向量阵V和特征值对角阵D [V,D] = eig(A,‘nobalance’)矩阵中有与截断误差（eps)相当的元素时,计算精度更高。 [V,D] = eig(A,B) 计算广义特征向量阵V和广义特征值阵D，使AV=BVD成立。 【例】求下列对称矩阵A的特征值/特征向量 求解线性方程组 矩阵方程AX=b的解 X=A\b 矩阵方程XA=b的解 X=A/b 三种方程组情况的求解（对m×n矩阵A） m=n, 恰定方程组，可以精确求解 mn, 超定方程组，可以求得最小二乘解 mn, 欠定方程组，求得一个基解 恰定方程组的解 方程ax＝b(a为非奇异)x=a-1 b矩阵逆 两种解: x=inv(a)?b — 采用求逆运算解方程 x=a\b — 采用左除运算解方程方程ax=b a=[1 2;2 3];b=[8;13]; ?x=inv(a)*bx=a\b x =x =2.002.003.003.00 一般代数方程求解 4.3 概率分布和统计数据分析—数据集 单变量的统计数据存储在1×n或n×1的向量里。 多变量的数据用矩阵来表示，不同的变量放在不同的列。 数据预处理 MATLAB用NaN表示非数 用NaN处理缺少的数据为数据分析带来方便 TF = isnan(A) 返回与A同维的矩阵TF，A中为NaN的得到逻辑1，否则为逻辑0 a = [-2 -1 0 1 2] isnan(0./a) Warning: Divide by zero. ans =00100 删除NaN i = find(~isnan(x));x = x(i) x = x(find(~isnan(x))) x = x(~isnan(x)); x(isnan(x)) = []; X(any(isnan(X)‘),:) = []; ％删除任意含有NaN的行 function X = excise(X) X(any(isnan(X)),:) = []; 样本分布的频数直方图描述 基本数据分析函数 max(X) :求各列最大值 min(X) :求各列最小值 mean(X):求各列平均值 median(X) :各列求中值 std(X) :求各列标准差 var (X) :求各列方差 Cov(X) :求X各列之间的协方差阵 Corrcoef(x):求X各列之间的相关系数 1、求取最大值(max) [Y,I]= 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 =359618 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 =184296253671 y=max(x) %查出二维数组x中各列元素的最大值产生赋予行向量y y =98