5.3 补充习题

5.3.1 微分中值定理

菲赫金哥尔茨《微积分学教程》第一卷中有一小节插值法讲了这个内容,这本来是数值分析里讲插值逼近误差分析的内容,但被讲数学分析的人用来构造各种吓人的中值问题。

——WXF

命题 5.3.1 (Birkhoff插值) n次函数g由以下n+1个条件确定:(1)g(ki)(xi)=f(ki)(xi),i=0,1,,n 其满足(记Nn={0,1,,n}):

  • 插值条件互不相同:i,jNnij,都有(ki,xi)(kj,xj)
  • Pólya正则条件(后k阶导的插值条件不超过n+1k个):kNncard{iNnkik}n+1k

则这样的g(以及其k阶导数)由以下行列式唯一确定:(2)|g(k)(x)(1)(k)(x)(k)(xn)(k)f(k0)(x0)(1)(k0)|x=x0(x)(k0)|x=x0(xn)(k0)|x=x0f(k1)(x1)(1)(k1)|x=x1(x)(k1)|x=x1(xn)(k1)|x=x1f(kn)(xn)(1)(kn)|x=xn(x)(kn)|x=xn(xn)(kn)|x=xn|=0

命题的证明是不言自明的,附加条件保证了g(k)(x)对应的余子式可逆。利用这个命题,我们可以构造出各种各样的中值定理。

例 5.3.2 (刘/闫/章·习题4.1.13) 设函数f,g,h[a,b]上连续、在(a,b)内可导,证明:ξ(a,b)使得(3)|f(a)g(a)h(a)f(b)g(b)h(b)f(ξ)g(ξ)h(ξ)|=0

证明(4)F(x):=|f(a)g(a)h(a)f(b)g(b)h(b)f(x)g(x)h(x)|F(a)=F(b)=0 由Rolle定理可知ξ(a,b)使得F(ξ)=0,即(5)F(ξ)=|f(a)g(a)h(a)f(b)g(b)h(b)f(ξ)g(ξ)h(ξ)|=0

例 5.3.3 (刘/闫/章·习题4.1.14) f[a,b]上一阶可导,在(a,b)内二阶可导,且f(a)=f(b)=0f+(a)f(b)>0,证明:

(1)

ξ(a,b),使得f(ξ)+2f(ξ)+f(ξ)=0

(2)

θ(a,b),使得f(θ)2f(θ)+f(θ)=0

(3)

η(a,b),使得f(η)=f(η)

(4)

ζ(a,b),使得f(ζ)=f(ζ)

我们先证明如下引理:

引理 5.3.4 f[a,b]上一阶可导,在(a,b)内二阶可导,且f(a)=f(b)=0f+(a)f(b)>0,则ξ1,ξ2,ξ(a,b),使得f(ξ1)=f(ξ2)=f(ξ)=0

引理的证明 不妨设f+(a)>0,否则可以考虑f。令ε=12f+(a),则δ(0,ba2)使得(6)0<xa<δε<f(x)f(a)xaf+(a)<εf(x)>12f+(a)(xa)>0x1(a,a+δ)(a,a+b2)使得f(x1)=0。同理,x2(bδ,b)(a+b2,b)使得f(x2)=0。由介值定理可知x0(x1,x2)使得f(x0)=0

反复应用Rolle定理可得:ξ1(a,x0)使得f(ξ1)=0ξ2(x0,b)使得f(ξ2)=0ξ(ξ1,ξ2)(a,b)使得f(ξ)=0

证明 (1) 令F(x)=exf(x),则F(a)=F(b)=0,且(7)F+(a)F(b)=ea+bf+(a)f(b)>0 由引理可知ξ(a,b)使得F(ξ)=0,即(8)f(ξ)+2f(ξ)+f(ξ)=0

(2) 令G(x)=exf(x),则G(a)=G(b)=0,且(9)G+(a)G(b)=eabf+(a)f(b)>0 由引理可知θ(a,b)使得G(θ)=0,即(10)f(θ)2f(θ)+f(θ)=0

(3) 令H(x)=exf(x),由引理可知f(ξ1)=f(ξ2)=0,亦即H(ξ1)=H(ξ2)=0。由Rolle定理可知η(ξ1,ξ2)(a,b)使得H(η)=0,即(11)f(η)=f(η)

(4) 令I(x)=ex(f(x)+f(x)),由引理可知ζ1,ζ2(a,b)使得F(ζ1)=F(ζ2)=0,亦即I(ζ1)=I(ζ2)=0。由Rolle定理可知ζ(ζ1,ζ2)(a,b)使得I(ζ)=0,即(12)f(ζ)=f(ζ)

例 5.3.5 (刘/闫/章·习题4.3.9) 设函数fC3[0,1],且f(0)=f(12)=0f(1)=12。证明:ξ(0,1),使得f(ξ)=12

证明f0=f(12),构造三次函数g(x)=a3x3+a2x2+a1x+a0满足(13)g(0)=g(12)=0,g(12)=f0,g(1)=12g(x)=2x32(1+2f0)x2+1+8f02xF(x)=f(x)g(x),则F(0)=F(12)=F(1)=F(12)=0。反复应用Rolle定理可得:

  • ξ1(0,12)使得F(ξ1)=0ξ2(12,1)使得F(ξ2)=0
  • η1(ξ1,12)使得F(η1)=0η2(12,ξ2)使得F(η2)=0
  • ξ(η1,η2)(a,b)使得F(ξ)=0,即f(ξ)=g(ξ)=12

例 5.3.6 (刘/闫/章·习题4.3.11) h>0,函数fC1[x0h,x0+h]。证明:θ(0,1),使得(14)f(x0+h)f(x0h)=[f(x0+θh)+f(x0θh)]h

证明F:[0,1]R满足F(x)=f(x0+xh)f(x0xh),由Lagrange中值定理可得θ(0,1)使得(15)f(x0+h)f(x0h)=F(1)F(0)=F(θ)=[f(x0+θh)+f(x0θh)]h

例 5.3.7 (刘/闫/章·习题4.3.12) 设函数fC1[a,b],在(a,b)内二阶可导,且f(a)=f(b)=0。证明:ξ(a,b),使得(16)|f(ξ)|4(ba)2|f(b)f(a)|

证明 由带Lagrange余项的Taylor展开可得(17)f(a+b2)=f(a)f(a)ba2+f(ξ1)2(ba2)2=f(a)+f(ξ1)8(ba)2=f(b)+f(b)ba2+f(ξ2)2(ba2)2=f(b)+f(ξ2)8(ba)2(18)f(b)f(a)=(ba)28[f(ξ1)f(ξ2)]ξ{ξ1,ξ2}满足|f(ξ)|=max{|f(ξ1)|,|f(ξ2)|},则(19)|f(b)f(a)|(ba)28[|f(ξ1)|+|f(ξ2)|](ba)24|f(ξ)| 得证。

例 5.3.8 (楼红卫·例6.5.4) 设函数fC1[a,b],在(a,b)内三阶可导。证明:ξ(a,b),使得(20)f(b)=f(a)+12(ba)(f(a)+f(b))112(ba)3f(ξ)

证明 构造三次函数g(x)=a3x3+a2x2+a1x+a0满足(21)g(a)=f(a),g(b)=f(b),g(a)=f(a),g(b)=f(b) 这样的g由以下公式确定:(22)|g(x)1xx2x3f(a)1aa2a3f(b)1bb2b3f(a)012a3a2f(b)012b3b2|=0|g(x)0006f(a)1aa2a3f(b)1bb2b3f(a)012a3a2f(b)012b3b2|=0 解得(23)g(x)=12(ba)3[f(b)f(a)]+6(ba)2[f(a)+f(b)]F(x)=f(x)g(x),则F(a)=F(b)=F(a)=F(b)=0。反复应用Rolle定理可得:

  • η(a,b)使得F(η)=0
  • ξ1(a,η)使得F(ξ1)=0ξ2(η,b)使得F(ξ2)=0
  • ξ(ξ1,ξ2)(a,b)使得F(ξ)=0,即(24)f(b)=f(a)+12(ba)(f(a)+f(b))112(ba)3f(ξ)

5.3.2 Taylor展开式(二)

例 5.3.9 阿基米德提出用圆内接正多边形的周长逼近圆的周长,从而计算圆周率的近似值。荷兰的Ludolph van Ceulen用一生的时间计算了262边形,得到圆周率的35位小数。由初等的平面几何知识可得:半径为1的圆内接正32n边形的边长an和周长Ln满足递推关系:(25)an+1=221(an2)2=an12+4an2Ln+1=32n+1an+1=Ln21+1(Ln32n)2 然而,这个公式收敛很慢。试确定常数λ,使得数列{(1λ)Ln+λLn+1}具有最快的收敛速度。

易见(26)Ln=32n+1sin2π32n+1=32n+1[2π32n+116(2π32n+1)3+O(125n)]=2ππ32722n+O(124n) 于是(27)(1λ)Ln+λLn+1=2π(43λ)π32722n+2+O(124n)λ=43可得(28)L~n=13Ln+43Ln+1=2π+O(124n)

例 5.3.10

(1)

证明:x[1,1](29)arctanx=xx33+x55+(1)nx2n+12n+1+

(2)

利用arctan12+arctan13=π4,求π的近似值。

证明 (1) 注意到(30)(arctanx)=11+x2=1x2+x4+(1)nx2n+o(x2n) 由Taylor展开式的唯一性可知ξ(1,1)使得(31)arctanx=xx33+x55+(1)nx2n+12n+1+arctan(2n+2)(ξ)(2n+2)!x2n+2 形式上有(32)(arctanx)=11+x2=12i(1x+i1xi) 求导可得(33)arctan(n+1)(x)=(1)nn!2i(1(xi)n+11(x+i)n+1) 故有(34)|arctan(n+1)(x)(n+1)!|22(n+1)max{1|xi|n+1,1|x+i|n+1}1n+1 因此(35)|arctan(2n+2)(ξ)(2n+2)!x2n+2||x|2n+22n+212n+20,n+,ξ(1,1)(36)arctanx=xx33+x55+(1)nx2n+12n+1+

(2) 由题可得(37)π4=arctan12+arctan13=n=1+(1)n12n1(122n1+132n1) 其前22项的和就能得到π的近似值为3.141592653589790,它的前14位小数是精确的。

5.3.3 函数的凹凸性

定理 5.3.11 f:IR,则凸函数的以下5个定义等价:

(1)

x1,x2It[0,1],有(38)f(tx1+(1t)x2)tf(x1)+(1t)f(x2)

(2)

x1,x2,x3Ix1<x2<x3,有(39)f(x2)f(x1)x2x1f(x3)f(x1)x3x1f(x3)f(x2)x3x2

(3)

fI上可微,则fI上单调不减。

(4)

fI上可微,则x1,x2I,有(40)f(x2)f(x1)+f(x1)(x2x1)

(5)

fI上二阶可微,则xI,有f(x)0

证明 (1)(2):当t=0,1时,不等式显然成立。当t(0,1)时,令x2=tx1+(1t)x3,则(41)f(tx1+(1t)x3)tf(x1)+(1t)f(x3)f(x2)f(x1)x2x1f(x3)f(x1)x3x1f(x3)f(x2)x3x2

(2)(3)x1,x3Ix1<x3,分别令x2x1+x2x3,由(2)可得(42)f(x1)=limx2x1+f(x2)f(x1)x2x1f(x3)f(x1)x3x1limx2x3f(x3)f(x2)x3x2=f(x3)

(3)(4):由Lagrange中值定理可知ξ(x1,x2)使得(43)f(x2)f(x1)x2x1=f(ξ)f(x1)f(x2)f(x1)+f(x1)(x2x1)

(4)(2):由(4)可得(44){f(x1)f(x2)+f(x2)(x1x2)f(x3)f(x2)+f(x2)(x3x2)f(x2)f(x1)x2x1f(x2)f(x3)f(x2)x3x2 再注意到(45)f(x3)f(x1)x3x1=x2x1x3x1f(x2)f(x1)x2x1+x3x2x3x1f(x3)f(x2)x3x2x2x1x3x1f(x2)f(x1)x2x1+x3x2x3x1f(x2)f(x1)x2x1=f(x2)f(x1)x2x1f(x3)f(x1)x3x1x2x1x3x1f(x3)f(x2)x3x2+x3x2x3x1f(x3)f(x2)x3x2=f(x3)f(x2)x3x2

(3)(5):由f单调不减可知xIhR满足x+hI,都有f(x+h)f(x)h同号,故(46)f(x+h)f(x)h0f(x)=limh0f(x+h)f(x)h0

(5)(3)x1,x2Ix1<x2,由Lagrange中值定理可知ξ(x1,x2)使得(47)f(x2)f(x1)x2x1=f(ξ)0f(x2)f(x1)

fI上严格凸,则(1)(2)(4)中的不等号需改为严格不等号,(3)中f需改为严格增;(5)不再等价,并且仅有(5)(3)f(x)>0f严格增)成立。其余的证明过程基本相近,但需要注意(2)(3)x1,x5Ix1<x5,取x1<x3<x5,分别令x2x1+x4x5,由(2)可得(48)f(x1)=limx2x1+f(x2)f(x1)x2x1f(x3)f(x1)x3x1<f(x5)f(x3)x5x3limx4x5f(x5)f(x4)x5x4=f(x5)

定理 5.3.12 I为开区间、f:IR为凸函数,则xIf的左右导数都存在,从而f连续,并且f(x)f+(x)

证明 由于x为内点,故δ>0使得[xδ,x+δ]I。令g:[δ,0)(0,δ]R满足g(h):=f(x+h)f(x)h,由凸函数的等价定义可得(49)h1<h2<0f(x)f(x+h1)h1f(x)f(x+h2)h2g(h1)g(h2)h1<0<h2f(x)f(x+h1)h1f(x+h2)f(x)h2g(h1)g(h2)0<h1<h2f(x+h1)f(x)h1f(x+h2)f(x)h2g(h1)g(h2)g在定义域内单调不减。当h<0时,g(h)有上界g(δ),因此limh0g(h)=f(x)存在;当h>0时,g(h)有下界g(δ),因此limh0+g(h)=f+(x)存在。由g单调不减可知f(x)f+(x)

例 5.3.13 证明:x1,,xnR+,都有(50)x1x1xnxn(x1++xnn)x1++xn

证明 考虑函数f(x):=xlnxf(x)=1x>0x>0),因此f严格凸,从而由Jensen不等式可得(51)f(x1++xnn)f(x1)++f(xn)n(52)x1++xnnlnx1++xnnx1lnx1++xnlnxnn 不等式两端同取exp即可得到待证不等式。

5.3.4 曲线的弯曲性质与渐近线

例 5.3.14 f在区间[a,+)上是凸函数,y=kx+by=f(x)x+时的一条渐近线。证明:

(1)

f可微,则limx+f(x)=k

(2)

f可微,则f(x)kx+bxa

(3)

f严格凸,则f(x)>kx+bx>a

证明 (1) 由f凸知f单调不减,于是limx+f(x)=AR{+}存在,AR当且仅当f有上界。由渐近线的定义和L’Hôpital法则可得(53)k=limx+f(x)x=limx+f(x)=A

(2) 考虑函数g(x):=f(x)kxb,则g(x)=f(x)k,由g凸知g单调不减,故由(1)知g(x)g(+)=0,因此g单调不增,从而g(x)g(+)=0,即f(x)kx+b

(3) 考虑函数g(x):=f(x)kxb,则g严格凸,g(+)=0

假设x1a使得g(x1)0,若存在x2>x1使得g(x2)>g(x1),则x>x2,都有(54)g(x)g(x2)xx2>g(x2)g(x1)x2x1 从而(55)g(x)>g(x2)+g(x2)g(x1)x2x1(xx2)+,x+ 这与g(+)=0矛盾。因此x>x1,都有g(x)g(x1),于是0=g(+)g(x1),从而g(x1)=0

x>x1g(x)g(x1)=0,仿照上述可证明g(x)=0,亦即g(x)0xx1),这与g严格凸矛盾。因此x>a,都有g(x)>0,即f(x)>kx+b

例 5.3.15 讨论函数f(x)=2x2x+1的凹凸性和渐近线。

计算可得(56)f(x)=2(x1)+2x+1,f(x)=4(x+1)3 所以当x<1时,f(x)<0f严格凹;当x>1时,f(x)>0f严格凸。

注意到limx1f(x)=,所以x=1y=f(x)的竖直渐近线。又f(x)=2x2+o(1)x),所以y=2x2y=f(x)x的渐近线;当x<1时,曲线y=f(x)位于这条渐近线的下方;当x>1时,曲线y=f(x)位于这条渐近线的上方。

例 5.3.16 讨论平面曲线x3+y3=3xy的渐近线和曲线的位置关系。

引入参数t=yx,得到曲线的参数方程(57)x(t)=3tt3+1,y(t)=3t2t3+1,tR{1}t1时,x(t),此时(58)y(t)x(t)=t1,y(t)+x(t)t=3t(1+t)1+t31,t1 所以y=x1是曲线在无穷远处的渐近线。注意到(59)y(t)+x(t)+1=t3+3t2+3t+11+t3=(t+1)21t+t2>0 所以曲线位于渐近线的上方。