MATLAB解方程与函数极值b.pptVIP

第3章 MATLAB数据分析与多项式计算 3.1 数据统计处理 3.2 数据插值 3.3 曲线拟合 3.4 离散傅立叶变换 3.5 多项式计算;;;;;3．两个向量或矩阵对应元素的比较 函数max和min还能对两个同型的向量或矩阵进行比较，调用格式为： (1) U=max(A,B)：A,B是两个同型的向量或矩阵，结果U是与A,B同型的向量或矩阵，U的每个元素等于A,B对应元素的较大者。 (2) U=max(A,n)：n是一个标量，结果U是与A同型的向量或矩阵，U的每个元素等于A对应元素和n中的较大者。 min函数的用法和max完全相同。 例3-3 求两个2×3矩阵x, y所有同一位置上的较大元素构成的新矩阵p。;;;3.1.3 平均值和中值 求数据序列平均值的函数是mean，求数据序列中值的函数是median。两个函数的调用格式为： mean(X)：返回向量X的算术平均值。 median(X)：返回向量X的中值。 mean(A)：返回一个行向量，其第i个元素是A的第i列的算术平均值。 median(A)：返回一个行向量，其第i个元素是A的第i列的中值。 mean(A,dim)：当dim为1时，该函数等同于mean(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的算术平均值。 median(A,dim)：当dim为1时，该函数等同于median(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的中值。 例3-5 分别求向量x与y的平均值和中值。;3.1.4 累加和与累乘积 在MATLAB中，使用cumsum和cumprod函数能方便地求得向量和矩阵元素的累加和与累乘积向量，函数的调用格式为： cumsum(X)：返回向量X累加和向量。 cumprod(X)：返回向量X累乘积向量。 cumsum(A)：返回一个矩阵，其第i列是A的第i列的累加和向量。 cumprod(A)：返回一个矩阵，其第i列是A的第i列的累乘积向量。 cumsum(A,dim)：当dim为1时，该函数等同于cumsum(A)；当dim为2时，返回一个矩阵，其第i行是A的第i行的累加和向量。 cumprod(A,dim)：当dim为1时，该函数等同于cumprod(A)；当dim为2时，返回一个向量，其第i行是A的第i行的累乘积向量。 ;3.1.5 标准方差与相关系数 1．求标准方差 在MATLAB中，提供了计算数据序列的标准方差的函数std。对于向量X，std(X)返回一个标准方差。对于矩阵A，std(A)返回一个行向量，它的各个元素便是矩阵A各列或各行的标准方差。std函数的一般调用格式为： Y=std(A,flag,dim) 其中dim取1或2。当dim=1时，求各列元素的标准方差；当dim=2时，则求各行元素的标准方差。flag取0或1，当flag=0时，按σ1所列公式计算标准方差，当flag=1时，按σ2所列公式计算标准方差。缺省flag=0，dim=1。 例3-7 对二维矩阵x，从不同维方向求出其标准方差。;;;;例3-9 对二维矩阵做各种排序。 3.2 数据插值 3.2.1 一维数据插值 在MATLAB中，实现这些插值的函数是interp1，其调用格式为： Y1=interp1(X,Y,X1,method) 函数根据X,Y的值，计算函数在X1处的值。X,Y是两个等长的已知向量，分别描述采样点和样本值，X1是一个向量或标量，描述欲插值的点，Y1是一个与X1等长的插值结果。method是插值方法，允许的取值有‘linear’、‘nearest’、‘cubic’、‘spline’。;;例3-11 某观测站测得某日6:00时至18:00时之间每隔2小时的室内外温度(℃)，用3次样条插值分别求得该日室内外6:30至17:30时之间每隔2小时各点的近似温度(℃)。 设时间变量h为一行向量，温度变量t为一个两列矩阵，其中第一列存放室内温度，第二列储存室外温度。命令如下： h =6:2:18; t=[18,20,22,25,30,28,24;15,19,24,28,34,32,30]; XI =6.5:2:17.5 YI=interp1(h,t,XI,‘spline’) %用3次样条插值计算;;例3-12 设z=x2+y2，对z函数在[0,1]×[0,2]区域内进行插值。 例3-13 某实验对一根长10米的钢轨进行热源的温度传播测试。用x表示测量点0:2.5:10(米)，用h表示测量时间0:30:60(秒)，用T表示测试所得各点的温度(℃)。试用线性插值求出在一分钟内每隔20秒、钢轨每隔1米处的温度TI。 命令如下： x=0:2.5:10; h=[0:30:60]; T=[95,14