第10章matlab高级数值计算.ppt

《第10章matlab高级数值计算.ppt》由会员分享，可在线阅读，更多相关《第10章matlab高级数值计算.ppt（33页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第10章matlab高级数值计算技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.1 多项式l多项式在数学理论分析、数值计算等方面具有很好的性质，这使得多项式在插值、回归分析、曲线拟合、微分方程求解等众多领域都有重要应用。本节将介绍多项式表示、多项式求值、多项式求根、多项式微积分、有理分式展开，为后续各节内容的展开奠定基础。技术凝聚实力技术凝聚实力专业创新出版专业创新出版技术凝聚实力技术凝聚实力专业创新出版专业创新出版技术凝聚实力技术凝聚实力专业创新出版专业创新出版技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.1.4 求多项式的根l多项式的根即是使的。MATLAB提供专门的函数ro

2、ots用于求多项式的根，函数roots的调用格式如下：ls=roots(p)l其中p为多项式表示，返回值s为 解向量，N为 多项式的阶数。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.1.5 多项式卷积和反卷积技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.1.6 多项式微积分技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.1.7 有理式的部分展开技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.2 插值技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.2.1 一维插值技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.2.2 二维插值l被插值函数为二元函数时，插

3、值过程为二维插值，依次类推，有三维插值、高维插值，这些内容将在下面的章节中介绍。图是二维插值的简单示意。lMATLAB利用函数interp2实现二维插值，其一般的调用格式为：lZI=interp2(X,Y,Z,XI,YI,method)技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.2.3 高维插值lMATLAB支持三维及三维以上的高维插值，分别由函数interp3和interpn实现。高维插值与三维插值类似，这里仅介绍三维插值，至于高维插值，读者可以参考MATLAB帮助文档和三维插值的例子。三维插值函数interp3的一般调用格式为：lVI=interp3(X,Y,Z,V,XI,YI,

4、ZI,method)l其中X、Y、Z、V 是具有相同大小的三维数组，X,Y,Z为三维数据网格，V是数据网格上的函数值；XI、YI、ZI、VI是具有相同大小的三维数组，返回值VI是三维插值网格XI,YI,ZI上的函数值估计；method为字符串，表示不同的插值方法，主要有以下四种：lmethod=nearest，最近邻插值。lmethod=linear，三次线性插值。lmethod=cubic，三次立方插值。lmethod=spline，样条插值。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.2.4 样条插值l利用分段多项式逼近函数可以降低插值多项式的阶数，使曲线连接处更加光滑，这种插值

5、方法称为样条插值，分段插值多项式称为样条函数，采样点称为节点。样条插值广泛地应用于各种制造业的计算机辅助设计（CAD）、各种图形的绘制工作、地理信息系统、实验数据的拟合、以及现在“热门”的计算机动画制作等。在样条函数中，应用最广的是三次样条函数。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.2.5 插值方法比较l不同的插值方法本质上是对插值函数的约束条件不同，相应地，插值的效果及效率也有很大的差别，这里对MATLAB中常用的四种插值方法总结如下：l最近邻插值法，利用阶梯函数作插值，速度快，内存消耗少，但是得到的插值数据光滑性能差。l线性插值法，利用分段线性函数作插值，速度快，内存消耗少

6、，但是在采样点处的光滑性能较差。l立方插值法，利用三次多项式函数作插值，在采样点处的光滑性能好，但是效率低，内存消耗大。l样条插值法，利用分段三次多项式函数作插值，速度较快，得到的插值数据光滑性能好。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.3 回归分析l回归分析和下节将要介绍的曲线拟合都是统计学中非常重要的数据分析方法，在信号处理、经济学等众多领域中都有广泛的应用。为此MATLAB对回归分析和曲线拟合提供了强大的支持，并专门提供了一个工具箱。l假设得到了以下实验观测数据：x=0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1，y=0.99567,0.9

7、9334,1.0413,1.0929,1.1485,1.2619,1.3719,1.4896,1.6433,1.8117,.9981，试找出x、y的约束关系。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.3.1 线性回归分析技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.3.2 多项式回归分析l当 为多项式，且 时，l此时线性回归分析称为多项式回归分析。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.3.3 多分量回归分析l前面介绍了单个变量的回归分析，实际中的问题需要考虑多个方面因素的影响，研究多个变量之间的联系，就是本小节将要研究的多个变量回归分析。对多个变量的回归分析

8、类似于单变量回归分析，首先对函数建模，要求函数模型与未知参数的关系是线性的，代入观测数据得Fa=y，由最小二乘拟合a=F/y求得线性参数，从而得到y=f(x)，这里的x为自变量组成的向量。下面通过一个简单的例子介绍多变量回归分析。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.4 曲线拟合l曲线拟合是要找到一条光滑曲线，使其最佳地拟合数据，该曲线不必经过数据点，这正是曲线拟合及回归分析与插值的不同之处。曲线拟合涉及两个基本问题：光滑曲线的形式是怎样的，以及最佳拟合的概念。第一个问题即是函数建模的问题，多项式模型，或是指数模型，以及别的更复杂的函数模型。最佳拟合是在某个误差准则下的最佳，常

9、用的误差准则是误差平方和最小，也可以是最大误差最小准则，误差的绝对值和最小等。以下的讨论均是以误差平方和最小作为误差准则。l根据函数模型的不同，曲线拟合可以分为多项式拟合、指数函数拟合等，用户甚至可以自己定义新的拟合形式。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.4.1 多项式拟合l顾名思义，多项式拟合是利用多项式最佳地拟合观测数据，使得在观测数据点处的误差平方和最小。MATLAB中利用函数polyfit和polyval进行多项式拟合。函数polyfit 根据观测数据及用户指定的多项式阶数得到光滑曲线的多项式表示，polyfit的一般调用格式为：lP=polyfit(x,y,n)l

10、其中x为自变量，y为应变量，n为多项式阶数。至于polyval的用法已经在10.1.3小节介绍过，这里不再赘述。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.4.2 指数函数拟合l指数函数拟合是利用指数函数 ，对观测数据进行拟合，使误差平方和最小。MATLAB对指数函数拟合没有提供专门的函数支持，通常利用一阶多项式拟合来解决指数函数拟合问题。对上式两边取对数，得 ，因此x、lny构成一阶多项式拟合问题。这里以对 的采样数据拟合为例。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.4.3 交互式曲线拟合工具lMATLAB为用户提供了一个交互式曲线拟合工具，即Basic Fittin

11、g interface。通过该工具，用户无须编写代码就可以完成一些常用的曲线拟合。下面以MATLAB自带的census data数据拟合为例介绍Basic Fitting interface的使用方法。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.5 傅立叶分析l傅立叶分析在信号处理、图像处理等众多的领域中都有极其重要的应用。傅立叶变换（Fourier Transform）将函数表示为不同频率的正弦、余弦函数之和，对于离散数据，傅立叶变换对应地为DFT，即离散傅立叶变换。l对长度为的输入序列，其DFT是长度为的向量，且，MATLAB利用快速傅立叶变换函数fft实现离散傅立叶变换。另外，

12、由的DFT也可以得到，称为的IDFT（逆傅立叶变换），MATLAB利用逆快速傅立叶变换函数ifft实现IDFT。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.5.1 快速傅立叶变换（FFT）、逆快速傅立叶变换（IFFT）技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.5.2 FFT的幅度和相位技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.5.3 傅立叶分析的应用例子l太阳黑子的活动具有一定的周期性，MATLAB提供了过去300年内太阳黑子活动的数据，下面通过对该数据的傅立叶分析研究太阳黑子活动的周期性。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.6 常微分方程l

13、解微分方程是数学中一个重要的问题，在物理学、控制理论、信号处理等方面都有很重要的应用。本节主要讨论解微分方程中一个基础而又重要的主题，即常微分方程的求解。本节首先介绍一阶常微分方程的求解，主要是ODE函数的用法，随后讨论如何将高阶常微分方程转换为一阶常微分方程，进而方便高微分方程的求解。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.6.1 一阶常微分方程l考虑下面的运动学问题：某物体初始位置为0，从时刻开始，以速度沿轴正向运动，求物体在时刻的位移。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.6.2 ODE函数的选择lMATLAB提供了8种ODE函数用于常微分方程的求解，它们分

14、别是ode45、ode23、ode113、ode15s、ode23s、ode23t、ode23tb、ode15i，前面介绍的ode45是比较常用的一种。这8个ODE函数分别用于不同类型常微分方程的求解，使用的算法也各不相同，不同ode函数的比较ODE Solver针对的常微分方程类型算 法ode45非刚性微分方程Kutta Rungeode23非刚性微分方程Runge-Kuttaode113非刚性微分方程Adamsode15s刚性微分方程和 DAEsNDFs(BDFs)ode23s刚性微分方程Rosenbrockode23t中等刚性微分方程DAEsTrapezoidal ruleode23tb

15、刚性微分方程TR-BDF2ode15i全隐式微分方程BDFs技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.6.3 高阶常微分方程lMATLAB对高阶常微分方程的求解是基于一阶常微分方程的，用户需要将高阶常微分方程转换为一阶常微分方程。技术凝聚实力技术凝聚实力专业创新出版专业创新出版10.7 小结l本章讨论数值计算的一些高级主题，如插值、拟合、微分方程求解等，通过本章的学习，读者应熟练掌握以下的内容。lMATLAB的多项式运算l插值l回归分析和曲线拟合l傅立叶分析l微分方程求解此此课件下件下载可自行可自行编辑修改，修改，仅供参考！供参考！感感谢您的支持，我您的支持，我们努力做得更好！努力做得更好！谢谢!