Matlab70Simulink60应用实例仿真与高效算法开发(第三章).docVIP

第三章 微积分问题的数值实验 通过本章的介绍，加上一些具有代表性的例题，读者可以了解微积分学的一些重要的概念与方法。实际科学与工程研究中，往往只得到一些离散的实验数据，并不知晓函数本身，无法用解析的方法对这些数据进行处理，此时，通过数值的方式进行数值微分与数值积分的运算显得尤为重要。 本章的内容是围绕着单变量与多变量函数微积分、函数极限、级数求和、Taylor幂级数展开Fourier级数展开等问题展开的。为方便读者理解，每节中都附有代表性的例题，并给出了程序代码，读者可以根据自己的需要对程序进行扩充，增补，以实现自己所需要的功能。作者的主要用意是综合运用不同指令解决具体问题，为读者解决实际问题提供一些思路和借鉴。 本章的主要内容包括 ◆ 微积分问题的解析解 ◆ 函数的级数展开与求和 ◆ 数值微分问题 ◆ 数值积分问题 ◆ 曲线积分与曲面积分的计算 3.1 微积分问题的解析解 3.1.1 函数极限问题的定义为：设函数 SKIPIF 1 0 在点x0的某一去心领域内有定义，如果对于任意给定的正数 SKIPIF 1 0 （无论它多么小），总存在正数 SKIPIF 1 0 ，使得对于适合不等式 SKIPIF 1 0 的一切x，对应的函数值 SKIPIF 1 0 都满足不等式 SKIPIF 1 0 那么常数A就叫做函数 SKIPIF 1 0 当 SKIPIF 1 0 时的极限，记作 SKIPIF 1 0 或 SKIPIF 1 0 (当 SKIPIF 1 0 ) 上述定义中的x0可以是某确定的值，也可为无穷大。 当常数A满足： SKIPIF 1 0 时，即 SKIPIF 1 0 时，可以称函数 SKIPIF 1 0 在点x0连续。函数在某一点连续又可分为左连续和右连续。如果 SKIPIF 1 0 存在且等于 SKIPIF 1 0 ，就说函数 SKIPIF 1 0 在点 SKIPIF 1 0 处左连续。如果 SKIPIF 1 0 存在且等于 SKIPIF 1 0 ，就说函数 SKIPIF 1 0 在点 SKIPIF 1 0 处右连续。 极限问题在MATLAB符号运算工具箱中提供了很方便的指令limit( )，该通过该函数不同的调用格式，可分别求得函数的极限和左/右极限。其调用格式如下： P = limit(fun, x, x0) 求函数fun关于自变量x在x0处的极限 P = limit(fun, x, x0，left 或 right) 求函数fun关于自变量x在x0处的单边极限 下面通过例子来演示MATLAB求函数极限的方法。 【例3-1】求解极限问题 SKIPIF 1 0 。 分析：Matlab中求解此类问题首先要定义符号变量，然后定义极限式子，接着才调用函数球给定函数的极限。 syms a b x; fun=log(a+b*x^2)/(sec(x)-cos(x)); Lim = limit(fun, x, 1) Lim = -log(a+b)*cos(1)/(-1+cos(1)^2) 【例3-2】 求解单边极限问题 SKIPIF 1 0 。 syms x; fun=x^(sin(x)); lim=limit(fun,x,0,right) lim = 1 我们还可以绘制出(-0.1, 0.1)区间的函数曲线，如下图所示 x=-0.1:0.0001:0.1; f_x=x.^(sin(x)); plot(x,f_x,black-,[0],[1],blacko) 图 3-1 x=0 附近的曲线 通过看x关于原点邻域内的曲线可以看到，函数 SKIPIF 1 0 在x=0处是连续的，所以无论是求左极限还是右极限，都为1。 另外，上面我们接触的都是函数自变量为单个情形，接下来简单介绍一下多变量函数的极限问题。 多元函数的极限问题同样可以通过嵌套调用matlab指令limit( )来解决。如面对二元函数 SKIPIF 1 0 ，若想求得二元函数的极限 SKIPIF 1 0 嵌套调用limit( )函数： SKIPIF 1 0