張宇1000題高等數學 第十五章 微分方程

目錄

AA

10.微分方程(1+x2)y2xy=0(1+x^2)y''-2xy'=0的通解爲______。

  令y=py'=p,則方程(1+x2)y2xy=0(1+x^2)y''-2xy'=0化爲p=2x1+x2pp'=\cfrac{2x}{1+x^2}p,兩端分別積分dpp=2x1+x2dx\displaystyle\int\cfrac{\mathrm{d}p}{p}=\displaystyle\int\cfrac{2x}{1+x^2}\mathrm{d}x,得lnp=ln(1+x2)+C0\ln|p|=\ln(1+x^2)+C_0,即y=C1(1+x2)y'=C_1(1+x^2)。積分可得y=C1(x+x33)+C2y=C_1\left(x+\cfrac{x^3}{3}\right)+C_2,其中C1C_1爲非零常數,C2C_2爲任意常數。(這道題主要利用了換元法求解

BB

2.微分方程dydx=2xyx2+y\cfrac{\mathrm{d}y}{\mathrm{d}x}=\cfrac{2xy}{x^2+y}的通解爲______。

  原方程改寫爲dxdy=x2+y2xy=x2y+12x\cfrac{\mathrm{d}x}{\mathrm{d}y}=\cfrac{x^2+y}{2xy}=\cfrac{x}{2y}+\cfrac{1}{2x},這是一個伯努利方程。令z=x2z=x^2,有dzdyzy=1\cfrac{\mathrm{d}z}{\mathrm{d}y}-\cfrac{z}{y}=1,得x2=z=e1ydy(e1ydydy+C)=y(lny+C)x^2=z=e^{\int\frac{1}{y}\mathrm{d}y}\left(\displaystyle\int e^{-\int\frac{1}{y}\mathrm{d}y}\mathrm{d}y+C\right)=y(\ln|y|+C),其中CC爲任意常數。(這道題主要利用了伯努利方程求解

3.微分方程ydxxdy=x2ydyy\mathrm{d}x-x\mathrm{d}y=x^2y\mathrm{d}y的通解爲______。

  將方程改寫爲ydyydxxdyx2=0y\mathrm{d}y-\cfrac{y\mathrm{d}x-x\mathrm{d}y}{x^2}=0,此爲全微分方程,即d(y22+yx)=0\mathrm{d}\left(\cfrac{y^2}{2}+\cfrac{y}{x}\right)=0。通解爲y22+yx=C\cfrac{y^2}{2}+\cfrac{y}{x}=C,其中CC爲任意常數。

10.設a>0a>0,函數f(x)f(x)[0,+)[0,+\infty)內連續有界,證明:微分方程y+ay=f(x)y'+ay=f(x)的解在[0,+)[0,+\infty)內有界。

  原方程的通解爲y(x)=eax[C+0xf(t)eatdt]y(x)=e^{-ax}\left[C+\displaystyle\int^x_0f(t)e^{at}\mathrm{d}t\right],其中CC爲任意常數。
  因f(x)f(x)[0,+)[0,+\infty)內有界,設f(x)M|f(x)|\leqslant M,則當x0x\geqslant0時,有
y(x)=eax[C+0xf(t)eatdt]Ceax+eax0xf(t)eatdtC+Meax0xeatdt=C+Ma(1eax)C+Ma. \begin{aligned} |y(x)|&=\left|e^{-ax}\left[C+\displaystyle\int^x_0f(t)e^{at}\mathrm{d}t\right]\right|\leqslant\left|Ce^{-ax}\right|+e^{-ax}\left|\displaystyle\int^x_0f(t)e^{at}\mathrm{d}t\right|\\ &\leqslant|C|+Me^{-ax}\displaystyle\int^x_0e^{at}\mathrm{d}t=|C|+\cfrac{M}{a}(1-e^{-ax})\leqslant|C|+\cfrac{M}{a}. \end{aligned}
  即y(x)y(x)[0,+)[0,+\infty)內有界。(這道題主要利用了放縮法求解

18.求(5x2y32x)y+y=0(5x^2y^3-2x)y'+y=0的通解。

  交換x,yx,y地位,認定xxyy的函數,則原方程化爲ydxdy+2x=5x2y3-y\cfrac{\mathrm{d}x}{\mathrm{d}y}+2x=5x^2y^3,是n=2n=2的伯努利方程。令z=x12=1xz=x^{1-2}=\cfrac{1}{x},則上述方程可化爲ydzdy+2z=5y3y\cfrac{\mathrm{d}z}{\mathrm{d}y}+2z=5y^3,進一步化爲一階線性微分方程的標準形式dzdy+2yz=5y2\cfrac{\mathrm{d}z}{\mathrm{d}y}+\cfrac{2}{y}z=5y^2。由通解公式得z=e2ydy(5y2e2ydydy+C)=1y2(5y2y2dy+C)=y3+Cy2z=e^{-\int\frac{2}{y}\mathrm{d}y}\left(\displaystyle\int5y^2e^{\int\frac{2}{y}\mathrm{d}y}\mathrm{d}y+C\right)=\cfrac{1}{y^2}\left(\displaystyle\int5y^2\cdot y^2\mathrm{d}y+C\right)=y^3+\cfrac{C}{y^2},故原方程的通解爲1x=y3+Cy2\cfrac{1}{x}=y^3+\cfrac{C}{y^2},其中CC爲任意常數。(這道題主要利用了伯努利方程求解

25.設f(x)f(x)在區間[0,+)[0,+\infty)上具有連續的一階導數,且滿足f(0)=1f(0)=1f(x)+f(x)11+x0xf(t)dt=0f'(x)+f(x)-\cfrac{1}{1+x}\displaystyle\int^x_0f(t)\mathrm{d}t=0

(1)求導函數f(x)f'(x)

  由題設可知f(0)+f(0)=0f'(0)+f(0)=0,即f(0)=f(0)=1f'(0)=-f(0)=-1,且(x+1)[f(x)+f(x)]=0xf(t)dt(x0)(x+1)[f'(x)+f(x)]=\displaystyle\int^x_0f(t)\mathrm{d}t(x\geqslant0)
  將上式兩端對xx求導並整理,得(x+1)f(x)+(x+2)f(x)=0(x+1)f''(x)+(x+2)f'(x)=0,這是以f(x)f'(x)爲未知函數的變量可分離的方程,故可求其通解爲f(x)=Cex1+xf'(x)=C\cfrac{e^{-x}}{1+x}
  代入初值條件f(0)=1f'(0)=-1,得C=1C=-1,所以f(x)=ex1+x(x0)f'(x)=-\cfrac{e^{-x}}{1+x}(x\geqslant0)

(2)證明:當x>0x>0時,有ex<f(x)<1e^{-x}<f(x)<1

  對f(x)f(x)在區間[0,x][0,x]上應用拉格朗日中值定理,得f(x)f(0)=xf(ξ)=xeξ1+ξ(0<ξ<x)f(x)-f(0)=xf'(\xi)=-x\cfrac{e^{-\xi}}{1+\xi}(0<\xi<x),所以f(x)<f(0)=1f(x)<f(0)=1
  爲了證明ex<f(x)e^{-x}<f(x),構造輔助函數F(x)=f(x)exF(x)=f(x)-e^{-x}。因爲F(x)=f(x)+ex=ex1+x+ex=xex1+xF'(x)=f'(x)+e^{-x}=-\cfrac{e^{-x}}{1+x}+e^{-x}=\cfrac{xe^{-x}}{1+x},所以函數F(x)F(x)在區間[0,+)[0,+\infty)上單調遞增,故當x>0x>0時,F(x)>F(0)=0F(x)>F(0)=0,即f(x)>exf(x)>e^{-x}。(這道題主要利用了構造函數求解

26.求微分方程{y+y=x,xπ2,y+4y=0,x>π4\begin{cases}y''+y=x,&x\leqslant\cfrac{\pi}{2},\\y''+4y=0,&x>\cfrac{\pi}{4}\end{cases}滿足條件yx=0=0,yx=0=0y\biggm\vert_{x=0}=0,y'\biggm\vert_{x=0}=0且在x=π2x=\cfrac{\pi}{2}處可導的特解。

  先求解當xπ2x\leqslant\cfrac{\pi}{2}時的初值問題{y+y=x,y(0)=y(0)=0.\begin{cases}y''+y=x,\\y(0)=y'(0)=0.\end{cases}
  易知,方程y+y=xy''+y=x的通解爲y=C1cosx+C2sinx+xy=C_1\cos x+C_2\sin x+x。根據條件y(0)=y(0)=0y(0)=y'(0)=0可解得C1=0,C2=1C_1=0,C_2=-1,所以相應的特解爲y=xsinx(xπ2)y=x-\sin x\left(x\leqslant\cfrac{\pi}{2}\right),此時,有yx=π2=π21,yx=π2=1y\biggm\vert_{x=\frac{\pi}{2}}=\cfrac{\pi}{2}-1,y'\biggm\vert_{x=\frac{\pi}{2}}=1
  進一步,當x>π2x>\cfrac{\pi}{2}時,欲使所求的解在x=π2x=\cfrac{\pi}{2}處可導(因而必連續),這就歸結爲求解新的初值問題{y+4y=0,yx=π2=π21,yx=π2=1\begin{cases}y''+4y=0,\\y\biggm\vert_{x=\frac{\pi}{2}}=\cfrac{\pi}{2}-1,y'\biggm\vert_{x=\frac{\pi}{2}}=1\end{cases}。易知,方程y+4y=0y''+4y=0的通解爲y=C3cos2x+C4sin2xy=C_3\cos2x+C_4\sin2x。再由初始條件yx=π2=π21,yx=π2=1y\biggm\vert_{x=\frac{\pi}{2}}=\cfrac{\pi}{2}-1,y'\biggm\vert_{x=\frac{\pi}{2}}=1可解得C3=1π2,C4=π2C_3=1-\cfrac{\pi}{2},C_4=-\cfrac{\pi}{2}。所以相應的特解爲y=(1π2)cos2x12sin2x(x>π2)y=\left(1-\cfrac{\pi}{2}\right)\cos2x-\cfrac{1}{2}\sin2x\left(x>\cfrac{\pi}{2}\right)
  因此,原方程滿足所給條件的特解爲y={xsinx,xπ2,(1π2)cos2x12sin2x,x>π2.y=\begin{cases}x-\sin x,&x\leqslant\cfrac{\pi}{2},\\\left(1-\cfrac{\pi}{2}\right)\cos2x-\cfrac{1}{2}\sin2x,&x>\cfrac{\pi}{2}.\end{cases}
  可以驗證所求函數yyx=π2x=\cfrac{\pi}{2}處可導。(這道題主要利用了分段函數求解

27.已知平面上的曲線y=f(x)y=f(x)與曲線0x2+yet2dt=2yxcosx\displaystyle\int^{x^2+y}_0e^{-t^2}\mathrm{d}t=2y-x\cos x相切於點(0,0)(0,0),且f(x)f(x)滿足微分方程y2y3y=3x1y''-2y'-3y=3x-1。求函數f(x)f(x)

  記y=g(x)y=g(x)是由方程0x2+yet2dt=2yxcosx\displaystyle\int^{x^2+y}_0e^{-t^2}\mathrm{d}t=2y-x\cos x相確定的隱函數。對方程兩邊關於xx求導,並利用隱函數求導法則,及積分上限函數的求導法則,有e(x2+y)2(2x+y)=2ycosx+xsinxe^{-(x^2+y)^2}(2x+y')=2y'-\cos x+x\sin x
  將x=0,y=0x=0,y=0代入上式,可得y(0)=1y'(0)=1,即g(0)=1g'(0)=1
  因爲兩曲線都經過點(0,0)(0,0)且在該點處具有公切線,所以f(0)=0,f(0)=g(0)=1f(0)=0,f'(0)=g'(0)=1。因此,函數f(x)f(x)是初值問題{y2y3y=3x1,yx=0=0,yx=0=1\begin{cases}y''-2y'-3y=3x-1,\\y\biggm|_{x=0}=0,y'\biggm|_{x=0}=1\end{cases}的解。下面先求方程y2y3y=3x1y''-2y'-3y=3x-1的通解。
  對應齊次方程的特徵方程r22r3=0r^2-2r-3=0有互異根r1=1,r2=3r_1=-1,r_2=3,則齊次方程的通解爲Y=C1ex+C2e3xY=C_1e^{-x}+C_2e^{3x}。非齊次方程的自由項3x13x-1,寫成Pm(x)eλxP_m(x)e^{\lambda x}的形式,應有P1(x)=3x1,λ=0P_1(x)=3x-1,\lambda=0。而λ=0\lambda=0不是特徵方程的根,故應設特解爲y=ax+by^*=ax+b。代入原方程,解得a=1,b=1a=-1,b=1,則y=x+1y^*=-x+1。因此,方程的通解爲f(x)=Y+y=C1ex+C2e3xx+1f(x)=Y+y^*=C_1e^{-x}+C_2e^{3x}-x+1
  最後,由f(0)=0,f(0)=1f(0)=0,f'(0)=1解得C1=54,C2=14C_1=-\cfrac{5}{4},C_2=\cfrac{1}{4},所以f(x)=54ex+14e3xx+1f(x)=-\cfrac{5}{4}e^{-x}+\cfrac{1}{4}e^{3x}-x+1。(這道題主要利用了積分上限函數的求導法則求解

CC

1.微分方程d2ydx2+(x+siny)(dydx)3=0\cfrac{\mathrm{d}^2y}{\mathrm{d}x^2}+(x+\sin y)\left(\cfrac{\mathrm{d}y}{\mathrm{d}x}\right)^3=0滿足初值條件y(0)=0,y(0)=23y(0)=0,y'(0)=\cfrac{2}{3}的特解是______。

  由反函數的求導法則dydx=1dxdy,d2ydx2=d2xdy2(dxdy)3\cfrac{\mathrm{d}y}{\mathrm{d}x}=\cfrac{1}{\cfrac{\mathrm{d}x}{\mathrm{d}y}},\cfrac{\mathrm{d}^2y}{\mathrm{d}x^2}=-\cfrac{\cfrac{\mathrm{d}^2x}{\mathrm{d}y^2}}{\left(\cfrac{\mathrm{d}x}{\mathrm{d}y}\right)^3}可知,原方程可化爲xx關於yy的二階常係數線性方程。將兩式代入原方程,原方程化爲d2xdy2x=siny\cfrac{\mathrm{d}^2x}{\mathrm{d}y^2}-x=\sin y,解得xx關於yy的通解爲x=C1ey+C2ey12sinyx=C_1e^y+C_2e^{-y}-\cfrac{1}{2}\sin y,當x=0x=0時,y=0y=0,代入得C1+C2=0C_1+C_2=0。將上式對yy兩邊求導,有dxdy=C1eyC2ey12cosy\cfrac{\mathrm{d}x}{\mathrm{d}y}=C_1e^y-C_2e^{-y}-\cfrac{1}{2}\cos y,當x=0x=0時,dxdy=1dydx=32\cfrac{\mathrm{d}x}{\mathrm{d}y}=\cfrac{1}{\cfrac{\mathrm{d}y}{\mathrm{d}x}}=\cfrac{3}{2},代入上式,有32=C1C212\cfrac{3}{2}=C_1-C_2-\cfrac{1}{2},解得C1=1,C2=1C_1=1,C_2=-1
  於是得特解x=eyey12sinyx=e^y-e^{-y}-\cfrac{1}{2}\sin y。(這道題主要利用了反函數的求導法則求解

4.已知y=f(x)y=f(x)是微分方程xyy=2xx2xy'-y=\sqrt{2x-x^2}滿足初值條件f(1)=0f(1)=0的特解。則 01f(x)dx=\displaystyle\int^1_0f(x)\mathrm{d}x=______。


01f(x)dx=xf(x)0101xf(x)dx=01[f(x)+2xx2]dx=01f(x)dx012xx2dx,01f(x)dx=12012xx2dx=12011(x1)2dx=x1=sint1201cos2tdt=π8. \begin{aligned} \displaystyle\int^1_0f(x)\mathrm{d}x&=xf(x)\biggm\vert^1_0-\displaystyle\int^1_0xf'(x)\mathrm{d}x=-\displaystyle\int^1_0[f(x)+\sqrt{2x-x^2}]\mathrm{d}x\\ &=-\displaystyle\int^1_0f(x)\mathrm{d}x-\displaystyle\int^1_0\sqrt{2x-x^2}\mathrm{d}x, \end{aligned}\\ \begin{aligned} \displaystyle\int^1_0f(x)\mathrm{d}x&=-\cfrac{1}{2}\displaystyle\int^1_0\sqrt{2x-x^2}\mathrm{d}x=-\cfrac{1}{2}\displaystyle\int^1_0\sqrt{1-(x-1)^2}\mathrm{d}x\\ &\xlongequal{x-1=\sin t}-\cfrac{1}{2}\displaystyle\int^1_0\cos^2t\mathrm{d}t=-\cfrac{\pi}{8}. \end{aligned}
這道題主要利用了分部積分法求解

8.設y(x)y(x)是方程y(4)y=0y^{(4)}-y''=0的解,且當x0x\to0時,y(x)y(x)xx33階無窮小,求y(x)y(x)

  由泰勒公式y(x)=y(0)+y(0)x+12y(0)x2+13!y(0)x3+ο(x3)(x0)y(x)=y(0)+y'(0)x+\cfrac{1}{2}y''(0)x^2+\cfrac{1}{3!}y'''(0)x^3+\omicron(x^3)(x\to0)。當x0x\to0時,y(x)y(x)x3x^3同階,則y(0)=0,y(0)=0,y(0)=0,y(0)=Cy(0)=0,y'(0)=0,y''(0)=0,y'''(0)=C,其中CC爲非零常數,由這些初值條件,現將方程y(4)y=0y^{(4)}-y''=0兩邊積分得0xy(4)(t)dt0xy(t)dt=0\displaystyle\int^x_0y^{(4)}(t)\mathrm{d}t-\displaystyle\int^x_0y''(t)\mathrm{d}t=0,即y(x)Cy(x)=0y'''(x)-C-y'(x)=0,兩邊再積分得y(x)y(x)=Cxy''(x)-y(x)=Cx
  易知,上述方程有特解y=Cxy^*=-Cx,因此它的通解是y=C1ex+C2exCxy=C_1e^x+C_2e^{-x}-Cx
  由初值y(0)=0,y(0)=0y(0)=0,y'(0)=0C1+C2=0,C1C2=CC_1+C_2=0,C_1-C_2=C,即C1=C2,C2=C2C_1=\cfrac{C}{2},C_2=-\cfrac{C}{2}
  故y[12(exex)x]Cy\left[\cfrac{1}{2}(e^x-e^{-x})-x\right]C,其中CC爲非零常數。(這道題主要利用了泰勒公式求解

11.求方程x2d2ydx22y=x2x^2\cfrac{\mathrm{d}^2y}{\mathrm{d}x^2}-2y=x^2的通解。

  當x>0x>0時,令x=etx=e^t,有t=lnxt=\ln x,經計算化原方程爲d2ydt2dydt2y=e2t\cfrac{\mathrm{d}^2y}{\mathrm{d}t^2}-\cfrac{\mathrm{d}y}{\mathrm{d}t}-2y=e^{2t},其特徵方程爲r2r2=0r^2-r-2=0,特徵根r1=2,r2=1r_1=2,r_2=-1,則其齊次方程的通解爲Y=C1e2t+C2etY=C_1e^{2t}+C_2e^{-t}。設特解y=Ate2ty^*=Ate^{2t},可得A=13A=\cfrac{1}{3}。從而得通解爲y=C1e2t+C2et+13te2t=C1x2+C2x+13x2lnxy=C_1e^{2t}+C_2e^{-t}+\cfrac{1}{3}te^{2t}=C_1x^2+\cfrac{C_2}{x}+\cfrac{1}{3}x^2\ln x
  當x<0x<0時,令x=ux=-u,原方程化爲yy關於uu的方程u2d2ydu22y=u2u^2\cfrac{\mathrm{d}^2y}{\mathrm{d}u^2}-2y=u^2,得通解y=C1u2+C2u+13u2lnu=C1x2C2x+13x2ln(x)y=C_1u^2+\cfrac{C_2}{u}+\cfrac{1}{3}u^2\ln u=C_1x^2-\cfrac{C_2}{x}+\cfrac{1}{3}x^2\ln(-x)
  合併兩種情形得原方程的通解爲y=C1x2+C2x+13x2lnx=C1x2C2x+13x2ln(x)y=C_1x^2+\cfrac{C_2}{|x|}+\cfrac{1}{3}x^2\ln|x|=C_1x^2-\cfrac{C_2}{x}+\cfrac{1}{3}x^2\ln(-x),其中C1,C2C_1,C_2爲任意常數。(這道題主要利用了換元法求解

14.求一條凹曲線,已知其上任意一點處的曲率k=12y2cosαk=\cfrac{1}{2y^2|\cos\alpha|},其中α\alpha爲該曲線在相應點處的切線的傾斜角,且該曲線在(1,1)(1,1)點處的切線爲水平方向。

  由曲率計算公式及曲線爲凹知,k=y[1+(y)2]32k=\cfrac{y''}{[1+(y')^2]^{\frac{3}{2}}}
  因爲α\alpha爲曲線在相應點處的切線的傾斜角,所以cosα=1secα=11+(y)2|\cos\alpha|=\cfrac{1}{|\sec\alpha|}=\cfrac{1}{\sqrt{1+(y')^2}}
  於是,由條件k=12y2cosαk=\cfrac{1}{2y^2|\cos\alpha|}推知,y[1+(y)2]32=1+(y)22y2\cfrac{y''}{[1+(y')^2]^{\frac{3}{2}}}=\cfrac{\sqrt{1+(y')^2}}{2y^2},整理得微分方程2y2y=[1+(y)2]22y^2y''=[1+(y')^2]^2
  此爲缺xx的可降階二階方程。令p=dydx,d2ydx2=dpdx=dydydydx=pdpdyp=\cfrac{\mathrm{d}y}{\mathrm{d}x},\cfrac{\mathrm{d}^2y}{\mathrm{d}x^2}=\cfrac{\mathrm{d}p}{\mathrm{d}x}=\cfrac{\mathrm{d}y}{\mathrm{d}y}\cdot\cfrac{\mathrm{d}y}{\mathrm{d}x}=p\cfrac{\mathrm{d}p}{\mathrm{d}y},代入上述微分方程,化簡爲2y2pdpdy=(1+p2)22y^2p\cfrac{\mathrm{d}p}{\mathrm{d}y}=(1+p^2)^2。分離變量得2pdp(1+p2)2=dyy2\cfrac{2p\mathrm{d}p}{(1+p^2)^2}=\cfrac{\mathrm{d}y}{y^2},解得y=(p2+1)+y(p2+1)C1y=(p^2+1)+y(p^2+1)C_1。由於曲線在點(1,1)(1,1)處切線水平,故y(1)=1,y(1)=0y(1)=1,y'(1)=0。於是有1=1+C1,C1=01=1+C_1,C_1=0。故得y=p2+1y=p^2+1,即dydx=±y1\cfrac{\mathrm{d}y}{\mathrm{d}x}=\pm\sqrt{y-1}。由於曲線是凹的,y=1y=1不是解,再將dydx=±y1\cfrac{\mathrm{d}y}{\mathrm{d}x}=\pm\sqrt{y-1}分離變量後積分得±y1=x+C2\pm\sqrt{y-1}=x+C_2
  由y(1)=1y(1)=1,所以C2=1C_2=-1,得±y1=x1\pm\sqrt{y-1}=x-1,化簡得4(y1)=(x1)24(y-1)=(x-1)^2。(這道題主要利用了曲率求解

寫在最後

  如果覺得文章不錯就點個贊吧。另外,如果有不同的觀點,歡迎留言或私信。
   歡迎非商業轉載,轉載請註明出處。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章