MATLAB教 第5章 MATLAB数值计算.pptVIP

MATLAB基础及应用第2版 主编：于润伟 朱晓慧 第5章 MATLAB数值计算 本章要点 数据的分析与统计 数值插值 曲线拟合 求解常微分方程 5.1 数据分析 一、 数据统计 二、离差和相关 5.2 数值运算 一、 多项式 二、插值与拟合 1．数值插值 （1）一维插值 函数：Z1=interp1（X，Y，X1，‘参数’） 说明：X是向量，表示采样点；Y是采样点上的样本值，与X等长；X1可是向量或标量，表示欲插值的点；Z1是与X1等长的插值结果。 （2）二维插值 函数：Z1=interp2(X，Y，Z，X1，Y1，‘参数’) 说明：X是长度为M的向量、Y是长度为N的向量，表示采样点；Z是与采样点对应的样本值，长度为N×M；X1、Y1是向量或标量，长度可以不等，表示欲插值的点；Z1是插值结果。二维插值是对双变量函数同时做插值。 线性插值主要参数 2．曲线拟合 MATLAB的曲线拟合函数如下： 函数：[P，S]=polyfit(X，Y，N) 说明：X、Y是两个等长的向量，X是采样点，Y是采样点函数值，N是多项式的次数；P是一个长度为N+1的向量，代表N次多项式；S是采样点的误差向量。 三、函数的极值和零点 5.3 常微分方程的数值求解 一、常微分方程的解法 数值解法的基本思想是：先取一系列离散的点， 通常取等步长h，使， 再求每个点对应的y(xn)，用y(xn)近似（n=1,2,…）的值。主要有欧拉法、线性多步法、预估校正法、龙格—库塔法等，其中以龙格—库塔法使用最多， 二、 龙格—库塔法的实现 基于龙格—库塔法，MATLAB提供了求常微分方程数值解的函数，其函数格式如下： 函数：[X，Y]=ode23(‘f ’，[x0,xn]，y0) [X，Y]=ode45(‘f’，[x0,xn]，y0) 说明：X，Y是两个向量。X对应自变量x在求解区间[x0，xn]的一组采样点，其采样密度是自适应的，无需指定；Y是与X对应的一组解。f是一个M函数文件，代表待求解方程。[x0，xn]代表自变量的求解区间。y0=y(x0)，由方程的初值给定。函数在求解区间[x0，xn]内，自动设定采样点向量X，并求出解函数y在采样点X处的样本值。 renshaobin523@ 第四章 计算机网络 计算机网络技术应用 * 第5章 MATLAB数值计算 * 求x各列元素之积 prod(x) 使x的各列元素按递增排序 sort(x) 找x各列的中间值元素 median(x) 求x各列元素之和 sum(x) 求x各列的平均值 mean(x) 找x各列的最小元素 min(x) 找x各列的最大元素 max(x) 功能 函数名称 功能 函数名称 两个矩阵x和y的相关系数 corr2(x，y) x的协方差阵 cov(x) 两个矩阵x和y的互相关系数，结果为方阵 corrcoef(x，y) x各列的极差 range(x) x的自相关阵 corrcoef(x) x各列的标准差 std(x) 两个矩阵x和y的协方差 cov(x，y) x各列的方差 var(x) 功能 函数名称 功能 函数名称 P1多项式与P2多项式相除后的导函数，导函数的分子放入q，分母放入r [q,r]=polyder(P1,P2) P1多项式与P2多项式乘积的导函数 P=polyder(P1,P2) 多项式P的导函数 p=polyder(P) 多项式求导 P1多项式与P2多项式相除。q为商，r为余数 [q,r]=deconv P1,P2) 多项式除法 P1多项式与P2多项式相乘 conv(P1,P2) 多项式乘法 m为n×n阶方阵。求x分别等于m中每一个元素时，多项式的值（结果为n×n阶方阵） polyvalm(P,m) 当A为标量时，求多项式P在自变量x=A时的值；当A为向量时，求x分别等于A中每个元素时，多项式的值 polyval(P，A) 求值 求该多项式的根，以列向量的形式给出 roots(P) 求根 A为向量。创建以向量A中元素为根的多项式 P=poly(A) P为多项式（以下各函数中P均为多项式），a0 a1 a2 … an-1 an为按降幂顺序排列的多项式系数 P=[ a0 a1 a2 …an-1 an] 创建多项式 说明 函数格式 名称 二维插值函数独有。插值点处的值和该点值的导数都连续 双立方插值法。利用已知的数据点拟合一个双立方曲面，然后根据插值点的坐标插值，每个插值点的值由该点附近的六个点的坐标确定 bicubic 需要较多的内存和运算时间，平滑性很好 立方插值法，也称三次