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

第三章 微积分问题的数值实验 通过本章的介绍，加上一些具有代表性的例题，读者可以了解微积分学的一些重要的概念与方法。实际科学与工程研究中，往往只得到一些离散的实验数据，并不知晓函数本身，无法用解析的方法对这些数据进行处理，此时，通过数值的方式进行数值微分与数值积分的运算显得尤为重要。 本章的内容是围绕着单变量与多变量函数微积分、函数极限、级数求和、Taylor幂级数展开Fourier级数展开等问题展开的。为方便读者理解，每节中都附有代表性的例题，并给出了程序代码，读者可以根据自己的需要对程序进行扩充，增补，以实现自己所需要的功能。作者的主要用意是综合运用不同指令解决具体问题，为读者解决实际问题提供一些思路和借鉴。 本章的主要内容包括 ◆ 微积分问题的解析解 ◆ 函数的级数展开与求和3.1 微积分问题的解析解 3.1.1 极限问题的解析解 函数极限问题的定义为：设函数在点x0的某一去心领域内有定义，如果对于任意给定的正数（无论它多么小），总存在正数，使得对于适合不等式的一切x，对应的函数值都满足不等式 那么常数A就叫做函数当时的极限，记作 或(当) 上述定义中的x0可以是某确定的值，也可为无穷大。 当常数A满足：时，即时，可以称函数在点x0连续。函数在某一点连续又可分为左连续和右连续。如果 存在且等于，就说函数在点处左连续。如果 存在且等于，就说函数在点处右连续。 极限问题在符号运算工具箱中提供了很方便的指令limit( )，该通过该函数不同的调用格式，可分别求得函数的极限和左/右极限。其调用格式如下： P = limit(fun, x, x0) 求函数fun关于自变量x在x0处的极限 P = limit(fun, x, x0，left 或 right) 求函数fun关于自变量x在x0处的单边极限 下面通过例子来演示求函数极限的方法。 求解极限问题。 分析：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) 求解单边极限问题 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) 通过看x关于原点邻域内的曲线可以看到，函数在x=0处是连续的，所以无论是求左极限还是右极限，都为1。 另外，上面我们接触的都是函数自变量为单个情形，接下来简单介绍一下多变量函数的极限问题。 多元函数的极限问题同样可以通过嵌套调用matlab指令limit( )来解决。如面对二元函数，若想求得二元函数的极限 嵌套调用limit( )函数： 也可以这样嵌套调用： 值得指出的是，当，不是确定的值，而是另外一个自变量的函数，例如，则极限的求取顺便不能交换。 求二元函数的极限问题 syms a b x y fun=exp(-3/(x^2+y^2))*sin(x)^2/x^2*(1+1/y^2)^(a*x+b^2*y^2) Lim=limit(limit(fun,x,1/sqrt(y)),y,inf) Lim = exp(b^2) 此例题中如某一确定数，那么极限问题中x ,y的顺序可以调换。 .1.2 函数导数的解析解 导数问题的物理意义可以理解为非匀速直线运动的速度和切线的斜率，数学上可以表达为： 更严格的定义为，设函数在某个邻域内有定义，当自变量x在x0处取得增量时，相应地函数y取得增量；如果与之比当时的极限存在，则称函数在点x0处可导，并称这个极限为函数在点x0处的导数，记为，即 中提供的求函数导数的指令为diff( )，它可以解出给定函数的各阶导数，其调用格式为： 求函数fun的关于x的导数 求函数fun的关于x的n阶导数 说明：其中fun为给定的待求导函数，x为自变量，这两个变量都为符号型的，求导的阶次用n来表示，默认为一阶导数。下面举例说明。 给定函数，试求其关于自变量x的四阶导数。 首先看看一阶导数 syms x; f=cos(x)/(x^3+7*x+2); f_1=diff(f,x) pretty(f_1) 2