MATLAB语言及应用培训讲义.pptVIP

第3章 MATLAB 在现代科学计算中的应用 数据处理与多项式计算 数据统计与分析 数值插值 曲线拟合 多项式计算 函数的最大值与最小值 数据统计与分析 1. 求矩阵最大和最小元素 (1)求向量的最大最小元素 ①y=max(X) 返回向量X的最大元素存入y。 ②[y,I]=max(X) 返回向量X的最大元素存入y，最大元素的序号存入I。 (2)求矩阵的最大和最小元素 ①max(A) 返回一个行向量，向量的第i个元素是A矩阵的第i列上的最大元素。 ②[Y,U]=max(A) 返回两个行向量，Y向量记录A的每列的最大元素，U向量记录每列最大元素的行号。 ③max(A,[],dim) dim取1或2。dim取1时，该函数和max(A)完全相同。dim取2时，该函数返回一个列向量，其第i个元素是A矩阵的第i行上的最大元素。 (3)两个向量或矩阵对应元素的比较 ①U=max(A,B) A,B是两个同型的向量或矩阵。结果U是与A,B同型的向量或矩阵，U的每个元素等于A,B对应元素的较大者。 ②U=max(A,n) n是一个标量。结果U是与A同型的向量或矩阵，U的每个元素等于A对应元素和n中的较大者。 min函数的用法和max完全相同。 例1求矩阵A的每行及每列的最大和最小元素，并求整个矩阵的最大和最小元素。 命令如下： A=[13,-56,78;25,63,-235;78,25,563;1,0,-1]; max(A,[],2) %求每行最大元素 min(A,[],2) %求每行最小元素 max(A) %求每列最大元素 min(A) %求每列最小元素 max(max(A)) %求整个矩阵的最大元素 min(min(A)) %求整个矩阵的最小元素 2. 求矩阵的平均值和中值 求矩阵和向量元素的平均值的函数是mean，求中值的函数是median。它们的调用方法和max函数完全相同。 3. 矩阵元素求和与求积 矩阵和向量求和与求积的基本函数是sum和prod，其使用方法和max类似。 例2 求矩阵A的每行元素的乘积和全部元素的乘积。 命令如下： A=[1,2,3,4;5,6,7,8;9,10,11,12]; S=prod(A,2) prod(S) %求A的全部元素的乘积 4. 矩阵元素累加和与累乘积 MATLAB中，使用cumsum和cumprod函数能方便地求得向量和矩阵元素的累加和与累乘积向量，函数的用法和sum及prod相同 例3 求向量X=(1！，2！，3！，…，10！)。 命令如下： X=cumprod(1:10) 5. 标准方差 MATLAB中，提供了计算数据序列的标准方差的函数std。对于向量X，std(X)返回一个标准方差。对于矩阵A，std(A)返回一个行向量，它的各个元素便是矩阵A各列的标准方差。std函数的一般调用格式为： std(A) 例3 产生8×6阶的正态分布随机数矩阵R，求其各列的平均值和标准方差，并求全体的平均值和标准方差。 6. 元素排序 MATLAB中对向量X是排序函数是sort(X)，函数返回一个对X中的元素按升序排列的新向量。 sort函数也可以对矩阵A的各列(或行)重新排序，其调用格式为： [Y,I]=sort(A,dim) 其中dim指明对A的列还是行进行排序，若dim=1，则按列排，若dim=2，则按行排。Y是排序后的矩阵，而I记录Y中的元素在A中位置。 例4 对矩阵做各种排序。 命令如下： A=[1,-8,5;4,12,6;13,7,-13]; sort(A) %对A的每列按升序排序 -sort(-A,2) %对A的每行按降序排序 [X,I]=sort(A) %对A按列排序，并将每个元素所在行号送矩阵I 数值插值 1. 一维数值插值 interp1函数调用格式为： Y1=interp1(X,Y,X1,method) 函数根据X、Y的值，计算函数在X1处的值。X、Y是两个等长的已知向量，分别描述采样点和样本值，X1是一个向量或标量，描述欲插值的点，Y1是一个与X1等长的插值结果。method是插值方法，允许的取值有linear(线性插值)、nearest(最近插值)、spline(三次样条插值)、cubic（三次多项式插值），缺省值是linear。 例5 用不同的插值方法计算sin(x)在π/2点的值。 这是一维插值问题。在MATLAB命令窗口，输入命令： X=0:0.2:pi;Y=sin(X); %给出X、Y interp1(X,Y,pi/2) %用缺省方法计算sin(π/2