【微积分】 07 - 微积分的应用

1. 微分的应用

1.1 一元函数的微分

1.1.1 单调性、极值、渐近线

  导数给出了函数的走向,它对我们分析函数的图形性质很有作用,这里就用微分学的知识来了解函数的性质。一阶导数对函数的影响是最直接的,这里先看一阶导数。对于区间上的常值函数f(x)=C,它的导数处处为零,反之由中值定理知,导数恒为零的函数为常值,故函数在区间上f(x)=0的充要条件是f(x)=0。这个结论还说明了导数相同的函数的差函数为常数,这对在证明函数相等很有用,比如可以证明3arccosxarccos3x4x2=π

  利用中值定理容易证明,区间上函数必定单调上升(下降)的充要条件是f(x)0f(x)0)。当等号不成立时,函数还是严格单调上升(下降)的。在部分点等号成立时,利用反证法可知,只要等号不在一个区间恒成立,函数也是严格单调上升(下降)的。

  我们已经知道,对任意函数f(x),如果x0极值点,则有f(x0)=0。反之则不一定(比如x3的零点),使得f(x0)=0的点x0一般称为静止点。如果在x0的领域内可导,则两侧的导数同号时为一般静止点,异号时为极值点(且根据具体情况可判定极大还是极小)。一般地,如果f(x0)=f(x0)==f(n1)(x0)=0,但f(n)0,则由泰勒公式知式(1)成立。进而可以证明,n为奇数时x0是一般静止点,n为偶数时,若f(n)>0x0是极小点,若f(n)<0x0是极大点。

(1)f(x)=f(x0)+1n!f(n)(x0)(xx0)n+o((xx0)n)

  有了以上结论,我们就能画出函数的大致曲线,只要弄清楚零点、单调性即可。另外有时还会有渐近线,它其实就是以下三种情况之一:(1)xx0f(x);(2)xf(x)y0;(3)式(2)分别存在有限极限。前两个分别以x=x0y=y0为渐近线,第三个以y=ax+b为渐近线。

(2)limxf(x)x=a;limx[f(x)ax]=b

1.1.2 凸函数

  最后来看二阶导数在函数图形上的体现,导数可以看做是曲线的切线斜率,那么二阶导数则可刻画了斜率的变化。可以想象,当f(x)>0时函数曲线上凸,而当f(x)<0时函数曲线下凸。如何严格地表述这样的曲线?可以这样说,连接曲线上任意两点形成直线,这两点间的函数值都在直线一侧。受此启发,定义在任意点式(3)都成立的函数为下凸(上凸)函数,等号不成立时也叫严格下凸(上凸)函数

(3)f(tx1+(1t)x2)()tf(x1)+(1t)f(x2),(0<t<1)

  以上凸函数的定义中并未假定函数可导,所以不好描述导数的性质,为此换做观察曲线上变化割线的性质。具体来讲,比如对下凸函数,设x1<x2<x3,利用定义容易证明式(4)成立它们还可以作为凸函数的等价定义。右边的不等式说明任意点x0右侧,f(x)f(x0)xx0随着xx0单调减小,但左边的不等式又说明它是有下界的,从而x0存在右极限(或极限为无穷)。同样可证x0存在左极限(或极限为无穷),当然,如果x0是端点,其中只有一个成立。当x0是区间内点时,显然f(x)x0处连续。

(4)f(x2)f(x1)x2x1f(x3)f(x2)x3x2;f(x2)f(x1)x2x1f(x3)f(x1)x3x1

  其实凸函数不一定可导,比如V字形的f(x)=|x|是下凸函数,但在x0不可导。当f(x)可导时,容易证明f(x)下凸(上凸)的充要条件是,f(x)单调上升(下降)。这个充要条件还等价于:曲线在它任何一条切线的上方。若f(x)二阶可导,还可以证明,f(x)下凸(上凸)的充要条件是f(x)0f(x)0)。这些比较直观,证明也很简单,请自行论证。

  上面的结论说明,如果f(x)连续且f(x0)=0,而在领域内f(x)0,可见f(x)x0左右两侧分别为上、下凸函数,所以曲线在x0左右领域内分别在x0切线的两侧。更一般的,如果f(x)x0处可导,且左右领域的点分别落在切线的两侧,则称x0f(x)拐点

  式(3)对二阶可导的凸函数还有进一步推广,设k=1npk=1,(pk>0),且记X=k=1nxk。对下凸函数f(x)可有式(5)成立,n个式子乘上pk相加便有式(6)成立。凸函数的这个结论,可以用来很容易地证明一些不等式,比如令f(x)=lnx,pk=1n,可以证明xk1nxk

(5)f(xk)=f(X)+f(X)(xkX)+12f"(ξ)(xX)2f(X)+f(X)(xkX)

(6)k=1npkf(xk)f(k=1npkxk)

1.2 多元函数的微分

1.2.1 切线、法平面

  现在利用微分的方法复习空间的点线面,请先复习空间解析几何的基本内容。空间曲线的表达式,最简单的就是参数方程x=x(t),y=y(t),z=z(t),一阶连续导数(x(t),y(t),z(t))确定了曲线在(x(t),y(t),z(t))处的切矢量T,切线连续变化的曲线称为光滑曲线。曲线还有可能表示为两个曲面的交集(式(7)左),利用向量值函数隐函数的结论可得到切矢量(1,yx(x),zx(x)),约去分母便得式(7)右。

(7){F(x,y,z)=0G(x,y,z)=0T=((F,G)(y,z),(F,G)(z,x),(F,G)(z,y))

  现在来看空间的曲面F(x,y,z)=0,如果Fx,Fy,Fz都连续,它被称为光滑曲面。考察曲面上经过(x0,y0,z0)的任意曲线,带入曲面方程有F(x(t),y(t),z(t))=0,由曲面的可微性易知曲线光滑,对t求导得式(8)。该式表明所有曲线在(x0,y0,z0)处的切线在同一平面上,这个平面被称为曲面在点(x0,y0,z0)法平面,它的法向量(Fx,Fy,Fz)所示。曲面还可能是用式(9)左边的参数方程表示的,用前两者可以确定隐函数u(x,y),v(x,y)。带入第三个式子就得到曲面表达式,算出法向量后约去分母便可导法向量(式(9)右)。

(8)Fx(x0,y0,z0)xt(t0)+Fy(x0,y0,z0)yt(t0)+Fz(x0,y0,z0)zt(t0)=0

(9){x=x(u,v)y=y(u,v)z=z(u,v)T=((y,z)(u,v),(z,x)(u,v),(x,y)(u,v))

1.2.2 曲率

  不知你有没有注意,平面曲线的二阶导数虽然表示斜率的变化速度,但由于斜率不与角度成正比,二阶导数其实并不能反映曲线的弯曲程度。要准确的度量曲线的弯曲程度,我们必须考察角度本身的变化率,具体讲就是在某点M0切线角度α相比长度s的变化率。如果式(10)的极限存在,则称k为点M0曲率,而1k称为曲率半径

(10)k=|dαds|=|limΔs0ΔαΔs|

  如果曲线以参数方程x(t),y(t)表示,首先有ds=xt2+yt2dt,再由α=arctanytxt也容易得到dα,从而容易有曲率的表达式(11),后者是坐标方程y=y(x)下的结果。对于极坐标方程r=r(θ),可以写成参数方程x=r(θ)cosθ,y=r(θ)sinθ,带入式(11)可得式(12)。特别地,对于圆r=r0,易知其曲率半径就是r0

(11)k=|xtyt2xt2yt|(xt2+yt2)32=|yx2|(1+yx2)32

(12)r=r(θ)k=|r2+2rθ2rrθ2|(r2+rθ2)32

1.2.3 极值

  类似一元函数的结论,对于偏导数处处存在的函数f(x1,,xn),如果fxi(x01,,x0n)=0皆成立,那么x0=(x01,,x0n)称为f静止点。静止点什么时候是极值点?再假设f有连续的二阶偏微分,使用泰勒公式即得式(13)。如果记对称矩阵Q={aij=fxixj(x0)},则式(13)的值取决于Q关于Δx的二次型。二次型正定(负定),则静止点是极大点(极小点),否则就不确定。

(13)Δf(x0)=f(x)f(x0)=12(Δx1x1++Δxnxn)2f(x0+θΔx)

  上面的极值假定变量可以在一个领域内变化,但实际问题中往往还有限制条件。比如已知Gi(x)=0,(i=1,,m),求F(x)的极值,这样的问题被称为条件极值。其实如果在局部(G1,,Gm)(x1,,xm)0,则根据Gi=0可以得到x1,,xm关于xm+1,,xn的隐函数,将它们带入F(x)即将问题转化为无条件极值问题。

  但很多时候,这样的隐函数无法直接写出,或者结果会破坏原本的对称性,从而使计算变得复杂。我们已经有了m个方程Gi=0,现在需要再找(nm)个“好”的方程。我们仍然以xm+1,,xn为自变量考虑问题,由F的极值首先有dF=i=1nFxidxi=0,注意其中dxi,(i=1,,m)为函数。如果能使等式中只有自变量xm+1,,xn的微分,则微分系数都为0,这就得到了另外的nm个方程。

  可以同样对Gi求微分dGi=i=1n(Gi)xidxi=0,由于(G1,,Gm)(x1,,xm)0,则可以选择参数λi,(i=1,,m),使得dF+i=1mλidGidxi,(i=1,,m)的系数为0。这时dxi,(i=m+1,,n)的系数必定是零,它们就是要找的nm个方程。

  现在来总结一下需要解的方程,为方便讨论,把λi也看作是未知数,并记Φ为式(14)左。原先的m个方程Gi=0其实就是Φλj=0,求λim方程其实是Φxi=0,(i=1,,m),而最后的nm个方程便是Φxi=0,(i=m+1,,n)。这个方法称为拉格朗日乘数法,式(14)更便于记忆。但还要注意,我们求得的只是“静止点”,还需根据实际情况确定是否是极值。

(14)Φ(x,λ)=F(x)+i=1mλiGi(x)Φxi=0Φλj=0

2. 积分的应用

2.1 一元函数的积分

2.1.1 平面面积,体积

  之前我们把定积分作为面积的一种定义,现在来看看这个定义的合理性,以及定积分更广泛的应用。首先我们来给出平面图形面积的一个直观定义,对于多边形,它们总可以分割为若干个三角形。对于一般平面图形P,我们总可以构造两个多边形B,ABP围住而AP围住,显然B的面积不小于A的面积。所有满足条件的A的面积有上确界S,所有满足条件的B的面积有下确界S,当S=S时称P可求积,且S=S=S称为P面积

  对于任意图形P,容易证明它可求积的充要条件是,存在多边形序列{Ai},{Bi},它们的面积极限相同。这个条件真好适合定积分的定义,所以对于可积函数,用定积分定义面积是合理的。对于复杂的图形(定义域为[a,b]),记x=x0截得的线段长为g(x)(连续),则图形面积为式(15)左。若xt的参数方程,且x(t)连续,则还可用式(15)右边计算。

(15)SP=abg(x)dx=αβg(x(t))x(t)dt,(x(α)=a,a(β)=b)

  以上定义面积方法其实可以推广开来,如果要求的量Q[a,b]上连续,将它分成若干部分,每一部分使用某个可积分的近似值f(xi)Δxi代替。然后证明误差部分趋于0,这样所求量就等于定积分abf(x)dx。对于每个具体的问题,证明误差部分趋于零是必须的,有时候也是困难的,对f(x)的论证非常必要。但下面的结论,我只打算给出粗略的描述,具体证明请参考教材。

  有些图形用极坐标描述更方便,对定义在[α,β]上的扇形r=r(θ),可以证明其面积为式(16)。类似地,可以用多面体来的近似来定义体积,使用多个棱柱计算更方便。若立体V定义在[a,b]上,且x处的截面面积为S(x),则可以证明其体积为abS(x)dx

(16)r=r(θ)SP=12αβr2(θ)dθ

2.1.2 长度、旋转面、曲线质量

  现在讨论平面里的一条曲线段l,它由参数方程x(t),y(t),t[p,q]给出。在线上取若干个点,用线段按顺序连接它们,然后用这些线段的长度之和的极限定义l长度。若曲线段自身不相交且不封闭,可以证明它的长度为式(17)。当曲线首尾相连时,可以拆成两段计算,容易证明这种情况公式仍然成立。如果把s看成t的函数,s(t)=x2(t)+y2(t)0且连续,从而t1(s)存在。x,y就可以看做s的函数,由ds2=dx2+dy2知公式(18)成立。

(17)sl=pqx2(t)+y2(t)dt=ab1+y2(x)dx=αβr2(θ)+r2(θ)dθ

(18)(dxds)2+(dyds)2=1

  对于一般曲面的面积,在后面会给出一般方法,这里只讨论一类特殊曲面的面积。对于定义在[a,b]上的曲线y(x),将它绕x轴旋转一周,曲线的路径形成旋转面Σ。可以用曲线上的分段线段的旋转面(圆台侧面)作为旋转面面积,已经知道每个线段旋转面的面积是π(yi+yi+1)didi为线段长),从而可以证明旋转面面积为2π0lydsl为曲线长度),整理即得式(19)成立。

(19)SΣ=2πpqy(t)x2(t)+y2(t)dt=2πaby(x)1+y2(x)dx

  更一般地,曲线l每一点的密度为f(x,y)(也可能是其它意义),那么曲线的质量是多少呢?同样的方法,将曲线分割为若干小段Δs,用所有段的质量和的极限作为l重量的定义。这样极限也被称为第一型曲线积分,记作lf(x,y)ds。类似长度的分析,可以用关于t参数方程表示x,y,s,并将第一型曲线积分转化为一元积分(式(20))。

(20)lf(x,y)ds=pqf(x(t),y(t))x2(t)+y2(t)dt

2.2 多元函数的积分

2.2.1 平面面积、体积

  重积分本身就是对面积(体积)的积分,因此将积分函数设为1便可求平面面积体积(式(21)),然后可以通过累次积分或换元法求得重积分。

(21)S=Ddxdy;V=Ωdxdydz

2.2.2 曲面面积、曲面质量

  现在来看一般空间曲面Γ的面积,先介绍一个基本结论:设平面π1,π2之间的夹角为θ,则容易证明π1上任何图形在π2上的垂直投影的面积是原图形的cosθ。为了定义曲面面积,我们将曲面f(x,y)分割为多个小区域Γ1,,Γn,每个区域在xy平面上的垂直投影是Di。对于每个区域Γi,可以用它上面的任一点ξi,ηi的法平面被投影分割的部分Ti来近似,为此还要假设f(x,y)有连续偏导数。

  设Ti,Di的面积分别为Δτi,Δσi,由于Ti的法向量为(fx(ξi,ηi),fy(ξi,ηi),1),故Ti,Di的夹角满足式(22)左。所以曲面的近似面积如式(22)右所示,它其实就是D上的一个积分和,因此曲面面积为式(23)的重积分。如果x,y,z由参数u,v给出,重新计算便得式(24)的重积分。

(22)cosθ=11+fx2(ξ,η)+fy2(ξ,η);i=1nΔτi=i=1nΔσicosθi

(23)S=D1+fx2(x,y)+fy2(x,y)dxdy

(24)S=DA2+B2+C2dudv,(A=(y,z)(u,v),B=(z,x)(u,v),C=(x,y)(u,v))

  类似于第一型曲线积分,如果曲面Γ上的密度为f(x,y,z)(或其它意义),则曲面质量被称作第一型曲面积分。该积分记作Γf(x,y,z)dτ,上面的曲面面积其实就是Γdτ。可以将微分dτ展开,从而将第一型曲面积分转化为二重积分(比如式(25),也可以写成关于参变量u,v的重积分)。

(25)Γf(x,y,z)dτ=Df(x,y,z(x,y))1+zx2+zy2dxdy

posted on   卞爱华  阅读(1781)  评论(0编辑  收藏  举报

编辑推荐:
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
阅读排行:
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 提示词工程——AI应用必不可少的技术
· Open-Sora 2.0 重磅开源!
· 周边上新:园子的第一款马克杯温暖上架

导航

点击右上角即可分享
微信分享提示