14.2 习题课讲解

14.2.1 有关一阶微分方程的证明题

例 14.2.1 (2023秋期末考试 · 17)设f∈C[0,1],且满足(1)|f(x)|≤1+∫0x|f(t)|dt 证明:|f(x)|≤ex,∀x∈[0,1]。

证明 令u(x):=e−x∫0x|f(t)|dt,则(2)u′(x)=e−x[|f(x)|−∫0x|f(t)|dt]≤e−x 故∀x∈[0,1],都有(3)u(x)=u(0)+∫0xu′(t)dt≤∫0xe−tdt=1−e−x 因此(4)|f(x)|≤1+∫0x|f(t)|dt=1+exu(x)≤ex ◻

另证 令g(x):=∫0x|f(t)|dt,则g(0)=0且(5)g′(x)=|f(x)|≤1+g(x)⟹g′(x)1+g(x)≤1 两边从0到x积分可得(6)ln⁡[1+g(x)]≤x⟹g(x)≤ex−1 因此(7)|f(x)|≤1+g(x)≤ex ◻

例 14.2.2 (Grönwall不等式的积分形式)设x0∈[a,b],f,g,h∈C[a,b]满足g非负,h在[a,x0]单调不增、在[x0,b]单调不减,且(8)f(x)≤∫x0xf(t)g(t)|dt|+h(x) 证明:(9)f(x)≤e∫x0xg(t)|dt|h(x),∀x∈[a,b] 其中(10)∫x0xφ(t)|dt|={∫x0xφ(t)dt,x≥x0∫xx0φ(t)dt,x≤x0

证明 为了突出强调x≥x0和x≤x0两种情况下符号和不等式方向的变化,这里将两部分证明合写在一起。方便起见,证明过程中的双重不等号均为不严格不等号,即≷表示一一对应的≥或≤,例如:(11)f≶g±h⟹(1)  f≤g+h(2)  f≥g−h 如果读者阅读此证明较为困难,可以将其手动分为两部分阅读。

先考虑一个简单的情形:h可微。设x≷x0,记F(x)=∫x0xf(t)g(t)|dt|+h(x),则(12)F′(x)=±g(x)⏟≷0f(x)⏟≤F(x)+h′(x)≶±g(x)F(x)+h′(x)⟹F′(x)∓g(x)F(x)≶h′(x) 其中h′(x)≷0。利用积分因子e−∫x0xg(t)|dt|可得(13)[F(x)e−∫x0xg(t)|dt|]′=e−∫x0xg(t)|dt|⏟>0[F′(x)∓F(x)g(x)]⏟≶h′(x)≶e−∫x0xg(t)|dt|⏟≤1h′(x)⏟≷0≶h′(x) 故有(14)F(x)e−∫x0xg(t)|dt|=F(x0)+∫x0x[F(ξ)e−∫x0ξg(t)|dt|]′⏟≶h′(ξ)dξ⏟≷0≤h(x0)+∫x0xh′(ξ)dξ=h(x) 综上所述,我们有(15)f(x)≤F(x)≤e∫x0xg(t)|dt|h(x)

对于一般的h和x,我们需要修改以上证明。从下式出发(红色表示待转化的式子,即需要h′可微的式子),利用分部积分法可得:(16)[F(x)e−∫x0xg(t)|dt|]′dξ≤?e−∫x0xg(t)|dt|h′(x)dξF(x)e−∫x0xg(t)|dt|≤?F(x0)+∫x0xe−∫x0ξg(t)|dt|h′(ξ)dξ=h(x0)+e−∫x0ξg(t)|dt|h(ξ)|ξ=ax−∫x0xh(ξ)[e−∫x0ξg(t)|dt|]′dξ=e−∫x0xg(t)|dt|h(x)±∫x0xe−∫x0ξg(t)|dt|g(ξ)h(ξ)dξ=e−∫x0xg(t)|dt|h(x)+∫x0xe−∫x0ξg(t)|dt|g(ξ)h(ξ)|dξ| 这等价于证明(17)e−∫x0xg(t)|dt|∫x0xf(t)g(t)|dt|≤?∫x0xe−∫x0ξg(t)|dt|g(ξ)h(ξ)|dξ| 此时不等式左侧可微,故可利用(18)e−∫x0xg(t)|dt|∫x0xf(t)g(t)|dt|=∫x0x[e−∫x0ξg(t)|dt|∫x0ξf(t)g(t)|dt|]′dξ=±∫x0xe−∫x0ξg(t)|dt|g(ξ)[f(ξ)−∫x0ξf(t)g(t)|dt|]dξ=∫x0xe−∫x0ξg(t)|dt|g(ξ)[f(ξ)−∫x0ξf(t)g(t)|dt|]|dξ|≤∫x0xe−∫x0ξg(t)|dt|g(ξ)h(ξ)|dξ| 因此(19)F(x)e−∫x0xg(t)|dt|≤e−∫x0xg(t)|dt|h(x)+∫x0xe−∫x0ξg(t)|dt|g(ξ)⏟≥0h(ξ)⏟≤h(x)|dξ|≤e−∫x0xg(t)|dt|h(x)+h(x)∫x0xe−∫x0ξg(t)|dt|g(ξ)|dξ|=h(x)[e−∫x0xg(t)|dt|−e−∫x0ξg(t)|dt||ξ=ax]=h(x) 亦即(20)f(x)≤F(x)≤e∫x0xg(t)dth(x) ◻

注 由于我们拓展了定积分的定义,即不必要求积分下限小于积分上限,故在证明时需要特别注意下式仅在b≥a时成立:(21)f(x)≤g(x)⟹∫abf(x)dx≤∫abg(x)dx,b≥a 或者可以一般地写成(22)f(x)≤g(x)⟹∫abf(x)|dx|≤∫abg(x)|dx|,∀a,b 三角(绝对值)不等式也可以一般地写成(23)|∫abf(x)dx|≤∫ab|f(x)||dx|

例 14.2.3 (Grönwall不等式的微分形式)设φ,ψ是非负的连续函数,η是可微函数,η′是Riemann可积函数,且(24)η′(t)≤φ(t)η(t)+ψ(t),∀t≥t0 证明:(25)η(t)≤e∫t0tφ(s)ds[η(t0)+∫t0tψ(s)ds],∀t≥t0

证明 考虑函数(26)F(t):=e−∫t0tφ(s)dsη(t) 则(27)F′(t)=e−∫t0tφ(s)ds[η′(t)−φ(t)η(t)]≤e−∫t0tφ(s)dsψ(t) 所以(28)F(t)=F(t0)+∫t0tF′(s)ds≤η(t0)+∫t0te−∫t0sφ(u)duψ(s)ds≤η(t0)+∫t0tψ(s)ds 因此(29)η(t)=e∫t0tφ(s)dsF(t)≤e∫t0tφ(s)ds[η(t0)+∫t0tψ(s)ds] ◻

注 Grönwall不等式的微分形式是积分形式的推论,积分形式放宽了对函数的可微性要求。

例 14.2.4 设x0,y0∈R,区间I∋x0,利用Grönwall不等式证明:以下初值问题的解在I上唯一存在。(30){y′=2xy2y(x0)=y0

证明 设y1,y2都是该初值问题的解,则(31)y1(x)=y0+∫x0xy1′(s)ds=y0+∫x0x2s(y1(s))2dsy2(x)=y0+∫x0xy2′(s)ds=y0+∫x0x2s(y2(s))2ds 令y=y1−y2,则有y(x0)=0,且(32)y(x)=y1(x)−y2(x)=∫x0x2s(y1(s)+y2(s))y(s)ds 由绝对值不等式可得(33)|y(x)|≤∫x0x|2s(y1(s)+y2(s))|⋅|y(s)|⋅|ds| 对|y(x)|使用Grönwall不等式可得|y(x)|≤0,因此y=0,即y1=y2。故解唯一存在。 ◻

注 本题可扩展至:设x0,y0∈R,区间I∋x0,若初值问题(34){y′=f(x,y)y(x0)=y0 满足f关于y Lipschitz连续,即存在M>0使得∀y1,y2∈R,有|f(x,y1)−f(x,y2)|≤M|y1−y2|,则解唯一存在。证明方法同上。

例 14.2.5 证明一阶线性方程(35)xy′−(2x2+1)y=x2,x>0 有且仅有一个解y∗(x)当x→+∞存在有限极限。写出y∗(x)的表达式,并求这个极限。

解 原方程的通解为(36)y(x)=xex2[C0+∫1xe−t2dt] 注意到∫1xe−t2dt收敛,而xex2→+∞,所以y有界仅当(37)C0=−∫1+∞e−t2dt 因此(38)y∗(x)=xex2[∫1xe−t2dt−∫1+∞e−t2dt]=−xex2∫x+∞e−t2dt 此时(39)limx→+∞y∗(x)=limx→+∞−∫x+∞e−t2dt1xe−x2=limx→+∞e−x2(−1x2−2)e−x2=−12 ◻

例 14.2.6 (2020秋期末考试 · 17)设f∈C(R)且为有界函数。

(1)

证明:常微分方程y′+y=f(x)的每个解y=y(x)都是[0,+∞)上的有界函数。

(2)

当x≤0时,常微分方程y′+y=f(x)是否存在有界解?若存在,有几个?

解 (1) 常微分方程的通解为(40)y(x)=e−x[∫0xf(t)etdt+y(0)] 当x≥0时,有(41)|y(x)|≤|y(0)|e−x+e−x|∫0xf(t)etdt|≤|y(0)|+e−x∫0x|f(t)|etdt≤|y(0)|+Me−x∫0xetdt≤|y(0)|+Me−x⋅ex=M1

(2) 当x≤0时,e−x→+∞,故常微分方程存在有界解仅当(42)y(0)=∫−∞0f(x)exdx 此时(43)y(x)=e−x∫−∞xf(t)etdt 验证可得(44)|y(x)|≤e−x|∫−∞xf(t)etdt|≤e−x∫−∞x|f(t)|etdt≤Me−x∫−∞xetdt≤Me−x⋅ex=M2 因此常微分方程存在有界解当且仅当y(0)=∫−∞0f(x)exdx,此时有界解的个数为1。 ◻

例 14.2.7 (2023秋期末考试 · 16)考虑一阶线性常微分方程(45)dydx+a(x)y=b(x) 其中a,b∈C(R),且

  • ∃c>0使得a(x)≥c,∀x≥0;
  • limx→+∞b(x)=0。

证明:该方程的每个解y=y(x)均满足limx→+∞y(x)=0。

证明 利用常数变易法或积分因子e∫0xa(s)ds可将方程的解y表示为(46)y(x)=e−∫0xa(s)ds[y(0)+∫0xb(t)e∫0ta(s)dsdt]=y0+∫0xb(t)e∫0ta(s)dsdte∫0xa(s)ds 由题设a(x)≥c>0知分母e∫0xa(s)dx在(0,+∞)严格增且(47)e∫0xa(s)ds≥ecx→+∞ 由L’Hôpital法则可得(48)(y0+∫0xb(t)e∫0ta(s)dsdt)′(e∫0xa(s)ds)′=b(x)e∫0xa(s)dsa(x)e∫0xa(s)ds=b(x)a(x) 由夹挤定理知(49)0≤|b(x)a(x)|≤|b(x)|c→0⟹limx→+∞b(x)a(x)=0 因此(50)limx→+∞y(x)=limx→+∞y0+∫0xb(t)e∫0ta(s)dsdte∫0xa(s)ds=0 ◻

另证 令A(x):=∫0xa(s)ds,同理可知eA(x)严格增且趋于+∞,由L’Hôpital法则可得(51)limx→+∞y(x)=limx→+∞eA(x)y(x)eA(x)=limx→+∞eA(x)[y′(x)+a(x)y(x)]eA(x)a(x)=limx→+∞b(x)a(x)=0 ◻

例 14.2.8 设函数f∈C(R)且恒非负,又存在常数A和a(a>0)使得(52)f(x)+a∫x−1xf(t)dt=A 证明:

(1)

F(x):=eaxf(x)单调增;

(2)

若infx∈Rf(x)=0,则f≡0。

证明 (1) 易证f∈C1(R),且有(53)f′(x)+af(x)=af(x−1) 所以(54)F′(x)=eax[f′(x)+af(x)]=eaxaf(x−1)≥0 因此F单调增。

(2) 由于F单调增且非负,注意到(55)f(x)=A−a∫x−1xf(t)dtF(x)=Aeax−a∫x−1xeaxf(t)dt≥Aeax−a∫x−1xea(t+1)f(t)dt=Aeax−aea∫x−1xF(t)dt≥Aeax−aeaF(x) 所以(56)f(x)≥A−aeaf(x)⟹A≤(1+aea)infx∈Rf(x)=0 因此(57)0≥A=f(x)+a∫x−1xf(t)dt≥f(x)≥0⟹f≡0 ◻

例 14.2.9 设f是连续的周期函数,周期为T>0,对方程(58)y′−λy=f(x) 讨论:(1) 有界解的个数;(2) 以T为周期的解的个数。

解 (1) 方程的解为(59)y(x)=eλx[y(x0)e−λx0+∫x0xe−λtf(t)dt] 因此(60)y(x0+T)=eλ(x0+T)[y(x0)e−λx0+∫x0x0+Te−λtf(t)dt]=eλTy(x0)+∫0Te−λsf(x0+s)ds 取an=y(nT),则有递推关系式(61)an+1=eλTan+∫0Te−λsf(s)ds

当λ=0时,{an}有界当且仅当∫0Tf(s)ds=0,这与y(0)无关;因此要么没有有界解,要么有无穷多个有界解。

当λ≠0时,令β满足(62)an+1+β=eλT(an+β)⟹β=11−eλT∫0Te−λsf(s)ds 从而有(63)an=enλT(a0+β)−β=enλT(y(0)+β)−β 故{an}有界当且仅当(64)y(0)=−β=1eλT−1∫0Te−λsf(s)ds 此时有界解的个数为1。

(2) 由(1)的讨论可知,(1)中所有的有界解都是以T为周期的解。 ◻

14.2.2 有关高阶微分方程的证明题

例 14.2.10 设f∈C[a,+∞),limx→+∞f(x)=0。证明:微分方程(65)y″+3y′+2y=f(x) 的所有解y=y(x)均满足limx→+∞y(x)=0。

证明 左侧的微分算子可因式分解得(66)(d2dx2+3ddx+2)y=(ddx+2)(ddx+1)y 利用两次积分因子法和∗∞型L’Hôpital法则可得(67)limx→+∞(y+y′)=limx→+∞e2x(y+y′)e2x=L′Hlimx→+∞e2x(y″+3y′+2y)2e2x=limx→+∞f(x)2=0limx→+∞y(x)=limx→+∞exyex=L′Hlimx→+∞ex(y′+y)ex=limx→+∞(y+y′)=0 ◻

注 对于任意y″+py′+qy=f(x),若特征方程λ2+pλ+q=0有两个负实根,则结论limx→+∞y(x)=0成立。证明方法同上。

例 14.2.11 对于微分方程(68){y″+p(x)y′+q(x)y=r(x),x∈(a,b)y(a)=A,y(b)=B 其中q(x)<0、A,B为常数,证明:若其在[a,b]上有连续的解,则解必定唯一。

证明 设y1,y2均为满足给定边界条件的微分方程的解,令y=y1−y2,则y满足(69){y″+p(x)y′+q(x)y=0,x∈(a,b)y(a)=y(b)=0 考虑Sturm-Liouville标准式(70)(p~y′)′+q~y=0⟹y″+p~′p~y′+q~p~y=0 所以(71)p~(x)=e∫axp(t)dt>0,q~(x)=p~(x)q(x)<0 因此(72)0≤∫ab−q~y2dx=∫aby(p~y′)′dx=yp~y′|x=ab−∫abp~(y′)2dx=−∫abp~(y′)2dx≤0 由于y∈C(a,b),故y≡0,亦即y1=y2。 ◻

注 本题提供的微分方程定解条件并非初始条件(初值问题),而是边界条件(边值问题)。本题是Sturm-Liouville定理的推论。

例 14.2.12 (2023秋期末考试 · 18,Sturm零点分离定理)设y1,y2为二阶线性齐次常微分方程y″+p(x)y′+q(x)y=0的两个线性无关解,其中p,q为开区间J上的连续函数。证明y1,y2的零点相互分离,即在y1的任意两个零点之间,必存在y2的一个零点,反之亦然。

证明 由于y1,y2线性无关,则W(x):=W[y1,y2](x)≠0,∀x∈J。由于W连续,不妨设W(x)<0,∀x∈J。设x0,x1为y1的相邻零点,则W(x0)=−y1′(x0)y2(x0),故y1′(x0),y2(x0)同号,不妨设它们均为正数。同时在x=x1处,有W(x1)=−y1′(x1)y2(x1)。由于x0,x1为y1的相邻零点,故y1′(x1)<0。设y1′(x1)>0,则∃x2∈(x0,x0+δ)使得y1(x2)>y1(x0)=0,∃x3∈(x1−δ,x1)使得y1(x3)<y1(x1)=0,由介值定理可得∃x4∈(x2,x3)使得y1(x4)=0,这与x0,x1为y1的相邻零点矛盾。

因而y2(x1)<0,故∃x2∈(x0,x1)使得y2(x2)=0。假设∃x3∈(x0,x1)且x3≠x2使得y2(x3)=0,则对(x2,x3)或(x3,x2)重复上述操作可得∃x4∈(x2,x3)使得y1(x4)=0,这与x0,x1为y1的相邻零点矛盾。故y2在(x0,x1)上有唯一零点x2。

综上所述,命题得证。 ◻