第五.MATLAB数值计算.pptVIP

第5章　MATLAB数值计算;5.1 数值线性代数问题及求解;(1)提取矩阵的对角线元素 设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素产生一个具有min(m,n)个元素的列向量。 diag(A)函数还有更进一步的形式diag(A,k)，其功能是提取第k条对角线的元素。 (2)构造对角矩阵 设V为具有m个元素的向量，diag(V)将产生一个m×m对角矩阵，其主对角线元素即为向量V的元素。 diag(V)函数也有更进一步的形式diag(V,k)，其功能是产生一个n×n(n=m+)对角阵，其第k条对角线的元素即为向量V的元素。; ; ; ; ; ;5.1.2 矩阵的简单变换; ; ; ; ; ;;;;5.1.9 矩阵分解; ;;5.2 数据处理与多项式计算;(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完全相同。;例5.25 求矩阵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类似。;;5. 标准方差 MATLAB中，提供了计算数据序列的标准方差的函数std。对于向量X，std(X)返回一个标准方差。对于矩阵A，std(A)返回一个行向量，它的各个元素便是矩阵A各列或各行的标准方差。std函数的一般调用格式为： std(A,FLAG,dim) 其中dim取1或2。当dim=1时，求各列元素的标准方差；当dim=2时，则求各行元素的标准方差。FLAG取0或1。;6. 元素排序 MATLAB中对向量X是排序函数是sort(X)，函数返回一个对X中的元素按升序排列的新向量。 sort函数也可以对矩阵A的各列(或行)重新排序，其调用格式为： [Y,I]=sort(A,dim) 其中dim指明对A的列还是行进行排序，若dim=1，则按列排，若dim=2，则按行排。Y是排序后的矩阵，而I记录Y中的元素在A中位置。;;5.2.2 数值插值 1. 一维数值插值 Interp1函数调用格式为： Y1=interp1(X,Y,X1,method) （1代表1维） 函数根据X、Y的值，计算函数在X1处的值。X、Y是两个等长的已知向量，分别描述采样点和样本值，X1是一个向量或标量，描述欲插值的点，Y1是一个与X1等长的插值结果。method是插值方法，允许的取值有linear(线性插值)、nearest(最近插值)、spline(三次样条插值)、cubic（三次多项式插值），缺省值是linear。;; ;2. 二维数值插值 MATLAB中，提供了解决二维插值问题的函数。其调用格式为： Z1=interp2(X,Y,Z,X1,Y1,method) 其中X、Y是两个向量，分别描述两个参数 的采样点，Z是与参数采样点对应的采样变 量的样本值，X1、Y1是两个向量或标量， 描述欲插值的点。method的取值与一维插 值函数相同。;;3. 三维数值插值 对三维函数插值的函数是interp3，其使用方法和interp2相同。其调用格式为： W1=interp3(X,Y,Z,W,X1,Y1,Z1,method) 函数返回三维插值结果。其中X、Y、Z是三个向量，分别描述三个参数的采样点，W是与参数采样点对应的采样变量的样本值，X1、Y1、Z1是三个向量或标量，描述欲插值的点。method是插值方法，可选，其缺省值是 ‘line。method的取值与一、二维插值函数相同。;;;POLYVAL Evaluate polynomial. Y = POLYVAL(P,X), when P is a vector of length N+1 whose elements