[工学]第4章 MATLAB数值计算 1.pptVIP

MATLAB 7.X程序设计 第4章 数值计算 (1) 主要内容： ①数据分析； ②矩阵分析； ③多项式运算； ④数值插值和拟合； ⑤数值积分和微分 。 4.1 概述 数值计算的特点： 1）数值计算是定义在数值数组的基础上，在计算之前必须定义变量并赋值。 2）数值计算有误差。 3）数值计算的计算速度较快。 4）数值计算除定义了算术运算、关系运算和逻辑运算三种，还提供了许多运算函数。 4.2 数据分析 4.2.1统计分析 1 最大值和最小值 1）求向量的最大和最小值 例4-1 已知x=[-43,72,9,16,23,47]， 求向量x的最大值和最小值。 y=max(x) %求向量x中的最大值 [y,l]=max(x) %求向量x中的最大值及其该元素的位置 z=min(x) %求向量x中的最小值 [z,m]=min(x) %求向量x中的最小值及其该元素的位置 2）求矩阵A的最大值和最小值 max(A) 返回一个行向量，向量的第i个元素是矩阵A的第i列上的最大值。 [Y,U]=max(A) 返回行向量Y和U，Y向量记录A的每列的最大值，U向量记录每列最大值的行号。 max(A,[],dim) dim取1或2。dim取1时，该函数和max(A)完全相同；dim取2时，该函数返回一个列向量，其第i个元素是A矩阵的第i行上的最大值。 2 平均值和中值 求数据序列平均值的函数是mean，求数据序列中值的函数是median。 median(A) 返回一个行向量，其第i个元素是A的第i列的中值。 例4-4 已知x=[1 8 4 2;9 6 2 5;3 6 7 1]，从 不同维方向求出其平均值和中值。 median(x) median(x,1) %按列方向，求数组的中值 3 求和与求积 sum(X) 返回向量X各元素的和。 sum(A,dim) 当dim为1时，该函数等同于sum(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的各元素之和。 例4-5 已知x=[4 5 6;1 4 8]，分析矩阵x 的每行、每列元素的乘积和全部元 素的乘积。 sum(x) sum(x,1) %求数组各列元素的和 sum(x,2) %求数组各行元素的和 sum(sum(x)) %求数组所有元素的和 prod(x) prod(x,1) %求数组各列元素的乘积 prod(x,2) %求数组各行元素的乘积 prod(prod(x)) %求数组所有元素的乘积 4 累加和与累乘积 cumsum(X) 返回向量X累加和向量。 cumprod(X) 返回向量X累乘积向量。 cumsum(A) 返回一个矩阵，其第i列是A的第i列的累加和向量。 cumprod(A) 返回一个矩阵，其第i列是A的第i列的累乘积向量。 例4-6 已知a = [1 2 3; 3 9 6; 4 10 8;  4 0 7]，求矩阵a的每行、每列元 素的乘积和全部元素的累加和。 cumsum(a) %求矩阵的各列元素的累加和 cumsum(a,2) %求矩阵的各行元素的累加和 cumprod(a) cumprod(a,1) %求矩阵的各列元素的累乘积 cumprod(a,2) %求矩阵的各行元素的累乘积  对于向量x，std(x，flag)返回一个标准方差。flag取0或1。 对于矩阵A，std(A)返回一个行向量，它的各个元素便是矩阵A各列或各行的标准方差。 Y=std(A,flag,dim) 其中dim取1或2。当dim=1时，求各列元素的标准方差；当dim=2时，则求各行元素的标准方差。flag取0或1。缺省flag=0，dim=1。 6 排序 sort(X) 返回一个对向量X中的元素按升序排列的新向量。 按升序排列行函数sortrows的调用格式： B = sortrows(A) 按升序排列A 的列，其中A为矩阵或行向量。 B = sortrows(A,column) 以指定的列向量排序矩阵 [B,index] = sortrows(A) 排序时返回向量的序号。 如果A 为列向量，B = A(index)； 如果A 为m-by-n 矩阵，B = A(index,:)。 例4-8 已知a = [1 2 3; 3 9 6; 4  10 8; 4 0 7]，对矩阵排序。 so