高等数学学习笔记 ☞ 微分方程

1.  微分方程的基本概念


1. 微分方程的基本概念:

(1)微分方程:含有未知函数及其导数或微分的方程。

举例说明微分方程:{y}'+y+x=0xdy=ydx

(2)微分方程的阶:指微分方程中未知函数导数最高阶数

举例说明微分方程的阶数:一阶微分方程:{y}'({y}')^{2};二阶微分方程:{y}''{y}''{y}'

(3)微分方程的解:就是把一个函数代入微分方程使其成立,那么这个函数就是微分方程的解。

(4)微分方程的通解:它是微分方程的解,该解具有《解里边含有任意常数的个数等于微分方程的阶数》特性。

举例说明微分方程的通解:{y}'=e^{x}\rightarrow y=e^{x}+C{y}''=e^{x}\rightarrow y=e^{x}+C_{1}x+C_{2}

(5)微分方程的特解:将给定的初始条件代入微分方程的通解,求解任意常数后的通解,称为微分方程的特解。

举例说明微分方程的特解:已知初始条件:y(0)=1,y(1)=2,微分方程{y}''=e^{x}的通解为:y=e^{x}+C_{1}x+C_{2}

                                           将初始条件代入通解,可得:C_{1}=2-e,C_{2}=0,则特解为y=e^{x}+(2-e)x

(6)解微分方程:就是找出未知函数的过程。

2. 微分方程的分类:

(1)根据微分方程中所含有的自变量的个数分类:

  ①:常微分方程:指微分方程中只包含一个自变量的方程,如:y=f(x)

  ②:偏微分方程:指微分方程中包含两个或两个以上自变量的方程,如:z=f(x,y)

备注:偏微分方程里边含有偏导数,如y{z}'_{x}+x{z}'_{y}=0

(2)根据微分方程中未知函数及其导数之间的关系分类:

  ①:线性微分方程:指未知函数及其导数呈线性关系。如:3{y}''+e^{x}{y}'+y=0

  ②:非线性微分方程:指未知函数及其导数不呈线性关系。如3{y}''y+e^{x}{y}'+y=0{y}'+2^{y}=0

备注:微分方程的线性与非线性与线性代数中的线性与非线性是有区别的。


2.  一阶微分方程


2.1 可分离变量的微分方程


1. 定义:就是指自变量和因变量可以分离开来的微分方程。可转化为形如M(x)dx=N(y)dy类型的等式。

说明:①:M(x):与x相关的式子;N(y):与y相关的式子。

           ②:等式两端同时取不定积分,即:\int M(x)dx=\int N(y)dy,再进一步进行求解即可。

2. 举例说明:

(1)求微分方程\frac{1}{x}dx=\frac{1}{y}dy的通解:

解:等式两端同时取不定积分,可得:\int \frac{1}{x}dx=\int \frac{1}{y}dy

       求解等式两端的不定积分可得:\ln |x|=\ln |y|-C_{1}\Rightarrow \ln |y|-\ln |x|=C_{1}\Rightarrow \ln |\frac{y}{x}|=C_{1}

       等式两端同时取e的指数,整理得:y=\pm e^{c_{1}}x

       进一步整理可得:y=Cx

备注:求解微分方程时,根据实际情况适当的对C进行处理,如\ln |y|=\ln |x|+\ln |C|\Rightarrow y=Cx

(2)求微分方程x(y+1)dx=y(x+1)dy的通解:

解:方程两端同时除以(x+1)(y+1),然后等式两端同时取不定积分,可得:\int \frac{x}{x+1}dx=\int \frac{y}{y+1}dy

       求解等式两端的不定积分可得:x-\ln|x+1|=-y+\ln |y+1|+C

       进一步整理可得微分方程的解:x+y-\ln|x+1|-\ln |y+1|=C

注意:当最等式两端做除法运算时,需要考虑相应的定义(验证通解是否包含特殊情况),即:

           当同时除以(x+1)(y+1)时,默认是x\neq -1,y\neq -1的,那么此时就需要验证x= -1,y= -1是否是

           微分方程的解。将x= -1,y= -1代入微分方程可知,等式成立,故x= -1,y= -1也是微分方程的解。

备注:求解微分方程时不要漏解,一般要求最终整理出来的式子不要带分式。


2.2 齐次微分方程


1. 定义:就是指自变量和因变量作为一个整体出现且次数相等的微分方程。可转化为形如\frac{dy}{dx}=f(\frac{y}{x})\frac{dx}{dy}=g(\frac{x}{y})类型的等式。

说明:\frac{dy}{dx}表明:x是自变量和y是因变量;\frac{dx}{dy}表明:y是自变量和x是因变量。

2. 举例说明:

(1)求微分方程(x^{3}+y^{3})dx-3xy^{2}dy=0的通解(固定解法,三步走):

解:方程两端同时除以x^{3},整理可得:\frac{dy}{dx}=\frac{1+(\frac{y}{x})^{3}}{3(\frac{y}{x})^{2}}。①

       三步走:\frac{y}{x}=u,则y=xu,\frac{dy}{dx}=u+x\frac{du}{dx}。②

       ②带入①上式可得:u+x\frac{du}{dx}=\frac{1+u^{3}}{3u^{2}}\Rightarrow\frac{3u^{2}}{1-2u^{3}}du=\frac{1}{x}dx

       对上式方程两端同时取不定积分并整理可得:-\frac{1}{2}\ln |1-2u^{3}|=\ln |x|+C_{1}\Rightarrow x\sqrt{1-2u^{3}}=C_{2}

       进行回代可得:x\sqrt{1-2(\frac{y}{x})^{3}}=C_{2}\Rightarrow 2y^{3}=x^{3}-Cx

(2)求微分方程\frac{dx}{dy}=\frac{2e^{\frac{x}{y}}(\frac{x}{y}-1)}{1+2e^{\frac{x}{y}}}的通解(固定解法,三步走):

解:三步走:\frac{x}{y}=u,则x=yu,\frac{dx}{dy}=u+y\frac{du}{dy}

       将上式带入微分方程可得:u+y\frac{du}{dy}=\frac{2e^{u}(u-1)}{1+2e^{u}}\Rightarrow \frac{1+2e^{u}}{2e^{u}+u}du=-\frac{1}{y}dy

       对上式方程两端同时取不定积分并整理可得:\ln |u+2e^{u}|+\ln|y|=\ln |C|\Rightarrow y(u+2e^{u})=C

       进行回代可得:y(\frac{x}{y}+2e^{\frac{x}{y}})=C\Rightarrow x+2ye^{\frac{x}{y}}=C


2.3 一阶线性微分方程


1. 定义:就是指未知函数的导数的最高阶数为1且未知函数及其导数之间呈线性关系的微分方程,

               可转化为形如\frac{dy}{dx}+P(x)y=Q(x)\frac{dx}{dy}+P(y)x=Q(y)类型的等式。

备注:当xy出现在一起时,默认意识都是y是关于x的函数,但其实x也是关于y的函数。

2. 种类:

(1)一阶齐次线性微分方程:指Q(x)=0Q(y)=0时的一阶线性微分方程。

求解方法:转化为可分离变量的微分方程的形式进行求解。

(2)一阶非齐次线性微分方程:指Q(x)\neq 0Q(y)\neq 0时的一阶线性微分方程。

求解方法:

第一步:将原微分方程转化为\frac{dy}{dx}+P(x)y=Q(x)\frac{dx}{dy}+P(y)x=Q(y)标准形式

              并明确P(x)Q(x)P(y)Q(y)的表达式。

第二步:直接带入通解公式:y=e^{-\int P(x)dx}(\int Q(x)e^{\int P(x)dx}dx+C)

               或x=e^{-\int P(y)dy}(\int Q(y)e^{\int P(y)dy}dy+C)

3. 举例说明:

(1)求微分方程\frac{dy}{dx}+\frac{1}{x}y=1的通解:

解:根据微分方程可知:P(x)=\frac{1}{x},Q(x)=1

       带入通解公式可得:y=e^{-\int\frac{1}{x}dx}(\int 1 \cdot e^{\int \frac{1}{x}dx}dx+C)

                                           =e^{-\ln x}(\int e^{\ln x}dx+C)

                                           =\frac{1}{2}x+\frac{1}{x}C

注意:求解一阶非齐次线性微分方程时,常数C只需要在最后那一项添加即可。

(2)求微分方程y^{3}dx+(2xy^{2}-1)dy=0的通解:

解:根据微分方程可知:P(y)=\frac{2}{y},Q(y)=\frac{1}{y^{3}}

       带入通解公式可得:x=e^{-\int \frac{2}{y}dy}(\int \frac{1}{y^{3}}\cdot e^{\int \frac{2}{y}dy}dy+C)

                                          =e^{\ln\frac{1}{y^{2}}}(\int \frac{1}{y^{3}}e^{\ln y^{2}}dy+C)

                                           =\frac{1}{y^{2}}(\ln |y|+C)

注意:求解一阶非齐次线性微分方程时,若e的指数幂含有\frac{1}{x}\frac{1}{y}形式的不定积分,则积出来的原函数不用加绝对值。


2.4 伯努利方程


1. 定义:形如\frac{dy}{dx}+P(x)y=Q(x)y^{n}(n\neq 0,1)的微分方程。

2. 求解过程:

    第一步:等式两端同时除以y^{n},可得:y^{-n}\frac{dy}{dx}+P(x)y^{1-n}=Q(x)

    第二步:将y^{-n}拿到导数部分,可得:\frac{1}{1-n}\frac{dy^{1-n}}{dx}+P(x)y^{1-n}=Q(x)

    第三步:令y^{1-n}=z,并整理可得:\frac{dz}{dx}+(1-n)P(x)z=(1-n)Q(x)

    第四步:将(1-n)P(x)看成一个整体,(1-n)Q(x)看成一个整体,即转化为一阶线性微分方程。

    第五步:按照一阶线性微分方程的求解方法,求解完微分方程后进行y^{1-n}=z的回代。

3. 举例说明:求微分方程\frac{dy}{dx}+\frac{y}{x}=(a\ln x)y^{2}的通解:

解:根据微分方程可知:n=2,则\frac{1}{y}=z\Rightarrow y=\frac{1}{z}

       将y=\frac{1}{z}带入原微分方程可得:-\frac{1}{z^{2}}\frac{dz}{dx}+\frac{1}{xz}=\frac{a\ln x}{z^{2}}\Rightarrow \frac{dz}{dx}-\frac{1}{x}z=-a\ln x

       由上式可得:P(x)=-\frac{1}{x},Q(x)=-a\ln x

       所以:z=e^{-\int-\frac{1}{x}dx}(\int( -a\ln x)e^{\int -\frac{1}{x}dx}dx+C)\Rightarrow x(-a\cdot \frac{1}{2}\cdot \ln ^{2}x+C)

       将y=\frac{1}{z}回代可得:xy(C- \frac{a}{2} \ln ^{2}x)=1


3.  高阶微分方程


3.1 可降阶的高阶微分方程


1. 形如:y^{(n)}=f(x)型的微分方程

    例:求微分方程{y}''=x-\sin x的通解。

    解:对微分方程{y}''=x-\sin x连续积分两次可得:{y}'=\frac{1}{2}x^{2}+\cos x+C_{1}y=\frac{1}{6}x^{3}+\sin x+C_{1}x+C_{2}

备注:y^{(n)}=f(x)型微分方程的通解含有n个任意常数(微分方程的阶数)

2. 形如:y^{(n)}=f(x,y^{(n-1)})型的微分方程

    例:求微分方程xy^{(4)}-y^{(3)}=0的通解。

    解:y^{(3)}=p,则y^{(4)}={p}'代入原微分方程可得:x{p}'-p=0\Rightarrow x\frac{dp}{dx}=p\Rightarrow \frac{dp}{p}=\frac{dx}{x}

           对\frac{dp}{p}=\frac{dx}{x}两端同时取积分可得:p=C_{1}x

           又知p=y^{(3)},则有y^{(3)}=C_{1}x

           对y^{(3)}=C_{1}x连续积分3次可得:y=C_{1}x^{4}+C_{2}x^{2}+C_{3}x+C_{4}

备注:替换后得到的一阶微分方程的形式不同,根据相应的形式去求解。

3. 形如:{y}''=f(y,{y}')型的微分方程

    例:求微分方程y{y}''-({y}')^{2}=0的通解。

    解:{y}'=p,则{y}''=p\frac{dp}{dy}代入原微分方程可得:yp\frac{dp}{dy}-p^{2}=0\Rightarrow \frac{dp}{p}=\frac{dy}{y}

           对\frac{dp}{p}=\frac{dy}{y}两端同时取积分可得:p=C_{1}y

           又知p={y}',则有{y}'=C_{1}y\Rightarrow \frac{dy}{dx}=C_{1}y\Rightarrow \frac{dy}{y}=C_{1}dx

           对\frac{dy}{y}=C_{1}dx同时取积分可得:\ln |y|=C_{1}x+C_{2}\Rightarrow y=C_{2}e^{C_{1}x}

核心思想:设{y}'=p,即:\frac{dy}{dx}=p,则{y}''={p}',即:{y}''=\frac{d(\frac{dy}{dx})}{dx}=\frac{dp}{dx}=\frac{dp}{dy}\cdot \frac{dy}{dx}=p\frac{dp}{dy}。 


3.2 二阶常系数齐次线性微分方程


1. 二阶齐次线性微分方程:形如{y}''+P(x){y}'+Q(x)y=0的微分方程。

2. 二阶齐次线性微分方程解的结构:

(1)y_{1}(x)y_{2}(x){y}''+P(x){y}'+Q(x)y=0的解,则y=C_{1}y_{1}(x)+C_{2}y_{2}(x)也是{y}''+P(x){y}'+Q(x)y=0的解。

备注:按照微分方程的解的定义即可求证。

         但y=C_{1}y_{1}(x)+C_{2}y_{2}(x)不一定是{y}''+P(x){y}'+Q(x)y=0的通解。eg:y_{1}(x)=x,y_{2}(x)=5x

备注:按照微分方程的通解的定义即可求证。

         那么如何保证y=C_{1}y_{1}(x)+C_{2}y_{2}(x){y}''+P(x){y}'+Q(x)y=0的通解呢?

         想要保证y=C_{1}y_{1}(x)+C_{2}y_{2}(x)是通解,那么通解必须含有两个任意常数,即这两个常数不会被整合成一个常数。

         此时,就需要保证y_{1}(x)y_{2}(x)不能成比例,也就是保证y_{1}(x)y_{2}(x)是线性无关的。

所以:y_{1}(x)y_{2}(x){y}''+P(x){y}'+Q(x)y=0的两个线性无关的解,那么y=C_{1}y_{1}(x)+C_{2}y_{2}(x)就是

           {y}''+p{y}'+qy=0的通解。(定理)

(2)若y_{1}(x)y_{2}(x)...y_{n}(x)n阶齐次线性微分方程的n个线性无关的解,那么n阶齐次线性微分方程的通解为

          y=C_{1}y_{1}(x)+C_{2}y_{2}(x)+...+C_{n}y_{n}(x)

3. 二阶常系数齐次线性微分方程:形如{y}''+p{y}'+qy=0的微分方程。

4. 二阶常系数齐次线性微分方程求解:

    若求{y}''+p{y}'+qy=0的通解,只需要求解出该微分方程的两个线性无关的解,那么根据相关定理,通解就可以表示为:

    y=C_{1}y_{1}(x)+C_{2}y_{2}(x)。故当前目标就是求解出该微分方程的两个线性无关的解。

    假设y=e^{rx}(r为常数)为{y}''+p{y}'+qy=0的解,但到底是不是该微分方程的解谁知道呢,需要进一步验证。

    根据y=e^{rx}可得:{y}'=re^{rx}{y}''=r^{2}e^{rx}。将其代入{y}''+p{y}'+qy=0,整理可得:e^{rx}(r^{2}+pr+q)=0

    因为e^{rx}\neq 0,故若使得e^{rx}(r^{2}+pr+q)=0恒成立,则有:r^{2}+pr+q= 0成立。

    由此可知,只需要解算出r的值,并带入y=e^{rx},就可以得到{y}''+p{y}'+qy=0的解y=e^{rx},故y=e^{rx}是微分方程

    {y}''+p{y}'+qy=0的解。

    此时,我们将r^{2}+pr+q= 0称为微分方程{y}''+p{y}'+qy=0的特征方程。并通过求解特征方程r^{2}+pr+q= 0的特征根

    获取微分方程 {y}''+p{y}'+qy=0的两个线性无关的解。分析如下:

    特征方程r^{2}+pr+q= 0是一个一元二次方程,它的根有以下三种情形:

 (1)特征方程r^{2}+pr+q= 0有两个不相等的实根,即r_{1}\neq r_{2}

          经过上述分析可知,y=e^{rx}是微分方程的解,所以我们就直接找到了微分方程的两个解,即y_{1}=e^{r_{1}x}y_{2}=e^{r_{2}x}

          因为\frac{y_{2}}{y_{1}}=\frac{e^{r_{2}x}}{e^{r_{1}x}}=e^{(r_{2}-r_{1})x}不是成比例关系的,故y_{1}=e^{r_{1}x}y_{2}=e^{r_{2}x}线性无关。

          所以微分方程{y}''+p{y}'+qy=0的通解为:y=C_{1}e^{r_{1}x}+C_{2}e^{r_{2}x}

举例说明,求微分方程{y}''-2{y}'-3y=0的通解:

解:根据所给的微分方程可得:特征方程r^{2}-2r-3=0

       解得:r_{1}=-1,r_{2}=3

       所以通解为:y=C_{1}e^{-x}+C_{2}e^{3x}

 (2)特征方程r^{2}+pr+q= 0有两个相等的实根,即r_{1}= r_{2}

          经过上述分析可知,y=e^{rx}是微分方程的解,所以我们就直接找到了微分方程的第一个解,即y_{1}=e^{r_{1}x}

          接下来寻找第二个解,第二个解只需要与第一个解线性无关即可。

          因为\frac{y_{2}}{y_{1}}不能成比例关系,所以设\frac{y_{2}}{y_{1}}=u(x),即:y_{2}=u(x)y_{1}=u(x)e^{r_{1}x}

          假设y_{2}=u(x)e^{r_{1}x}{y}''+p{y}'+qy=0的解,但到底是不是该微分方程的解谁知道呢,需要进一步验证。

          根据y_{2}=u(x)e^{r_{1}x}可得:{​{y}'_{2}}=e^{r_{1}x}({u}'(x)+r_{1}u(x)){y}''_{2}=e^{r_{1}x}({u}''(x)+2r_{1}{u}'(x)+r_{1}^{2}u(x))

          将其代入{y}''+p{y}'+qy=0,整理可得:{u}''(x)+(2r_{1}+p){u}'(x)+(r_{1}^{2}+pr_{1}+q)u(x)=0

          因为r_{1}是特征方程r^{2}+pr+q= 0的根,则有r_{1}^{2}+pr_{1}+q= 0

          又因为特征方程r^{2}+pr+q= 0有两个相等的实根,则特征方程可表示为(r-r_{1})^{2}=r^{2}-2r_{1}r+r_{1}^{2}=0

          通过等式对比可知:p=-2r_{1}

          故{u}''(x)+(2r_{1}+p){u}'(x)+(r_{1}^{2}+pr_{1}+q)u(x)=0\Rightarrow {u}''(x)=0

          我们的目的是寻找一个u(x),使得\frac{y_{2}}{y_{1}}不成比例,所以根据{u}''(x)=0,可选取u(x)=x

          所以微分方程的第二个解为y_{2}=xe^{r_{1}x}

          因为\frac{y_{2}}{y_{1}}=x不是成比例关系的,故y_{1}=e^{r_{1}x}y_{2}=xe^{r_{1}x}线性无关。

          所以微分方程{y}''+p{y}'+qy=0的通解为:y=C_{1}e^{r_{1}x}+C_{2}xe^{r_{1}x}=(C_{1}+C_{2}x)e^{r_{1}x}

举例说明,求微分方程{y}''+2{y}'+y=0的通解:

解:根据所给的微分方程可得:特征方程r^{2}+2r+1=0

       解得:r_{1}=r_{2}=-1

       所以通解为:y=C_{1}e^{-x}+C_{2}xe^{-x}=(C_{1}+C_{2}x)e^{-x}

 (3)特征方程r^{2}+pr+q= 0有一对共轭复根,即r_{1}=\alpha + i\betar_{2}=\alpha - i\beta

          经过上述分析可知,y=e^{rx}是微分方程的解,所以我们就直接找到了微分方程的两个解,即y_{1}=e^{r_{1}x}y_{2}=e^{r_{2}x}

          代入整理可得: y_{1}=e^{(\alpha +i\beta )x}y_{2}=e^{(\alpha -i\beta )x},此时的解为复数形式,为得到实数形式的解,为此进一步将其转化。

          利用欧拉公式:e^{i\theta }=\cos \theta +i\sin \theta可以把 y_{1}=e^{(\alpha +i\beta )x}y_{2}=e^{(\alpha -i\beta )x}改写为:

          y_{1}=e^{(\alpha +i\beta )x}=e^{\alpha x}(\cos \beta x+i\sin \beta x), y_{2}=e^{(\alpha -i\beta )x}=e^{\alpha x}(\cos \beta x-i\sin \beta x)

          因为y_{1},y_{2}微分方程的两个解,则\bar{y}_{1}=\frac{1}{2}(y_{1}+y_{2})\bar{y}_{2}=\frac{1}{2i}(y_{1}-y_{2})依然是微分方程的两个解(代入证明可得结果),

          即:\bar{y}_{1}=\frac{1}{2}(y_{1}+y_{2})=e^{\alpha x}\cos \beta x\bar{y}_{2}=\frac{1}{2i}(y_{1}-y_{2})=e^{\alpha x}\sin \beta x

          因为\frac{y_{2}}{y_{1}}=\tan \beta x不是成比例关系的,故\bar{y}_{1}=e^{\alpha x}\cos \beta x\bar{y}_{2}=e^{\alpha x}\sin \beta x线性无关。

          所以微分方程{y}''+p{y}'+qy=0的通解为:y=e^{\alpha x}(C_{1}\cos \beta x + C_{2}\sin \beta x)

举例说明,求微分方程{y}''-2{y}'+5y=0的通解:

解:根据所给的微分方程可得:特征方程r^{2}-2r+5=0

       解得:r_{1}=1+2i,r_{2}=1-2i。故:\alpha =1,\beta =2

       所以通解为:y=e^{ x}(C_{1}\cos 2 x + C_{2}\sin 2 x)

5. n阶常系数齐次线性微分方程求解:

对于n阶常系数齐次线性微分方程,根据特征方程的特征根,可以写出其对应的微分方程的解如下所示:
特征方程的根微分方程通解中的对应项
单实根r给出一项:Ce^{rx}
一对单复根r_{1,2}=\alpha \pm i\beta给出两项:e^{\alpha x}(C_{1}\cos \beta x + C_{2}\sin \beta x)
k重实根r给出k项:e^{rx}(C_{1}+C_{2}x+C_{3}x^{2}+...+C_{k}x^{k-1})
一对k重复根r_{1,2}=\alpha \pm i\beta给出2k项:e^{\alpha x}[(C_{1}+C_{2}x+C_{3}x^{2}+...+C_{k}x^{k-1})\cos \beta x + (D_{1}+D_{2}x+D_{3}x^{2}+                        ...+D_{k}x^{k-1})\sin \beta x](C,D\, is\, constant)

举例说明:若微分方程的特征根别r_{1}=0r_{2}=r_{3}=1r_{4,5}=1\pm 2ir_{6,7}=r_{8,9}=1\pm 3i

则通解为:y=C_{1}+e^{ x}(C_{2}+C_{3}x)+e^{ x}(C_{4}\cos 2 x + C_{5}\sin 2 x)+ e^{ x}[(C_{6}+C_{7}x)\cos 3 x + (C_{8}+C_{9}x)\sin 3 x]


3.3 二阶常系数非齐次线性微分方程


1. 二阶非齐次线性微分方程:形如{y}''+P(x){y}'+Q(x)y=f(x)的微分方程。

2. 二阶非齐次线性微分方程解的结构:

(1)y^{*}{y}''+P(x){y}'+Q(x)y=f(x)的特解,Y{y}''+P(x){y}'+Q(x)y=0的通解,

         则:{y}''+P(x){y}'+Q(x)y=f(x)的通解为y^{*}+Y

(2)若y^{*}_{1}y^{*}_{2}{y}''+P(x){y}'+Q(x)y=f(x)的特解,则:y^{*}_{1}-y^{*}_{2}{y}''+P(x){y}'+Q(x)y=0的特解。

(3)若y^{*}_{1}{y}''+P(x){y}'+Q(x)y=f(x)的特解,y^{*}_{2}{y}''+P(x){y}'+Q(x)y=g(x)的特解,

         则:y^{*}_{1}+y^{*}_{2}{y}''+P(x){y}'+Q(x)y=f(x)+g(x)的特解。

3. 二阶常系数非齐次线性微分方程:形如{y}''+p{y}'+qy=f(x)的微分方程。

4. 二阶常系数非齐次线性微分方程求解:

    由二阶非齐次线性微分方程解的结构可知:

    二阶常系数非齐次线性微分方程的通解等于二阶常系数非齐次线性微分方程的特解 对应的二阶常系数齐次线性微分方程的通解

    上一节已解决二阶常系数齐次线性微分方程的通解的求法,此时需要解决的问题就是求二阶常系数非齐次线性微分方程的特解。

    以下是微分方程{y}''+p{y}'+qy=f(x)中的f(x)取两种常见形式时特解的求法。

(1){y}''+p{y}'+qy=e^{\lambda x}P_{m}(x)

备注:\lambda是常数;P_{m}(x)是关于x的一个m次多项式。

     观察{y}''+p{y}'+qy=e^{\lambda x}P_{m}(x)可知,等号右边是多项式与指数函数的乘积。而多项式与指数函数的乘积的导数依然是多项式

     与指数函数的乘积,所以我们有理由怀疑y^{*}=e^{\lambda x}R(x)可能是微分方程的特解。其中R(x)是某个多项式。

     假设y^{*}=e^{\lambda x}R(x)是微分方程的特解,则{y^{*}}'=e^{\lambda x}[\lambda R(x)+{R}'(x)]{y^{*}}''=e^{\lambda x}[\lambda^{2} R(x)+2\lambda {R}'(x)+{R}''(x)]

     将{y^{*}}'{y^{*}}''带入{y}''+p{y}'+qy=e^{\lambda x}P_{m}(x)可得如下等式:

                                                             {R}''(x)+(2\lambda +p){R}'(x)+(\lambda ^{2}+p\lambda +q)R(x)=P_{m}(x)\blacksquare

     根据{y}''+p{y}'+qy=e^{\lambda x}P_{m}(x)可得对应的{y}''+p{y}'+qy=0的特征方程为:r^{2}+pr+q= 0

     ①:若\lambda不是特征方程r^{2}+pr+q= 0的根,则有\lambda ^{2}+p\lambda +q\neq 0

            因为P_{m}(x)是关于x的一个m次多项式,若使得\blacksquare式恒成立,则等号左端也应该是关于x的一个m次多项式,由于R(x)求导

            之后,次数降低,所以此时只能R(x)是关于x的一个m次多项式。

            设R(x)=R_{m}(x)=b_{0}x^{m}+b_{1}x^{m-1}+...+b_{m},将R(x),{R}'(x),{R}''(x)带入\blacksquare式,利用待定系数法,即可得到

            b_{i}(i=0,1,...m)的值。所以此时微分方程的特解为:y^{*}=R_{m}(x)e^{\lambda x}

     ②:若\lambda是特征方程r^{2}+pr+q= 0的单根,则有\lambda ^{2}+p\lambda +q= 02\lambda +p\neq 0

            因为P_{m}(x)是关于x的一个m次多项式,若使得\blacksquare式恒成立,则等号左端也应该是关于x的一个m次多项式,由于R(x)的系数

            为零,所以此时只能{R}'(x)是关于x的一个m次多项式,那么R(x)则是关于x的一个m+1次多项式。

            设{R}'(x)=R_{m}(x)=b_{0}x^{m}+b_{1}x^{m-1}+...+b_{m},那么R(x)可用x{R}'(x)=xR_{m}(x)来表示,将R(x),{R}'(x),{R}''(x)

            带入\blacksquare式,利用待定系数法,即可得到b_{i}(i=0,1,...m)的值。所以此时微分方程的特解为:y^{*}=xR_{m}(x)e^{\lambda x}

     ③:若\lambda是特征方程r^{2}+pr+q= 0的重根,则有\lambda ^{2}+p\lambda +q= 02\lambda +p= 0

            因为P_{m}(x)是关于x的一个m次多项式,若使得\blacksquare式恒成立,则等号左端也应该是关于x的一个m次多项式,由于R(x){R}'(x)

            的系数为零,所以此时只能{R}''(x)是关于x的一个m次多项式,则R(x)是关于x的一个m+2次多项式。

            设{R}''(x)=R_{m}(x)=b_{0}x^{m}+b_{1}x^{m-1}+...+b_{m},那么R(x)可用x^{2}{R}''(x)=x^{2}R_{m}(x)来表示,将R(x),{R}'(x),{R}''(x)

            带入\blacksquare式,利用待定系数法,即可得到b_{i}(i=0,1,...m)的值。所以此时微分方程的特解为:y^{*}=x^{2}R_{m}(x)e^{\lambda x}

综上所述,总结如下:

如果f(x)=e^{\lambda x}P_{m}(x),那么二阶常系数非齐次线性微分方程的特解可以表示为:y^{*}=x^{k}R_{m}(x)e^{\lambda x}

其中:R_{m}(x)P_{m}(x)是同次(m次)的多项式。

           \lambda不是特征方程的根,则k=0\lambda是特征方程的单根,则k=1\lambda是特征方程的重根,则k=2

(2){y}''+p{y}'+qy=e^{\lambda x}[P_{l}(x)\cos wx+P_{n}(x)\sin wx]

备注:\lambda,\omega是常数;P_{l}(x),P_{n}(x)分别是关于x的一个l次、n次多项式,且有一个多项式可以为零。

     如果f(x)=e^{\lambda x}[P_{l}(x)\cos wx+P_{n}(x)\sin wx],那么二阶常系数非齐次线性微分方程的特解可以表示为:

     y^{*}=x^{k}e^{\lambda x}[R_{m}^{(1)}(x)\cos wx+R_{m}^{(2)}(x)\sin wx]

     其中:R_{m}^{(1)}(x)R_{m}^{(2)}(x)m次的多项式(eg:R_{m}^{(1)}(x)=ax+bR_{m}^{(2)}(x)=cx+d)m= \max \left \{ l,n \right \}

                \lambda+\omega i(或\lambda-\omega i)不是特征方程的根,则k=0\lambda+\omega i(或\lambda-\omega i)是特征方程的单根,则k=1

备注:根据题设写出特解,将其代入原微分方程,利用待定系数法即可求解多项式系数。


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/6891.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【优选算法】9----长度最小的子数组

----------------------------------------begin-------------------------------------- 铁子们,前面的双指针算法篇就算告一段落啦~ 接下来是我们的滑动窗口篇,不过有一说一,算法题就跟数学题一样,只要掌握方法,多做…

计算机网络之链路层

本文章目录结构出自于《王道计算机考研 计算机网络_哔哩哔哩_bilibili》 02 数据链路层 在网上看到其他人做了详细的笔记,就不再多余写了,直接参考着学习吧。 1 详解数据链路层-数据链路层的功能【王道计算机网络笔记】_wx63088f6683f8f的技术博客_51C…

Kubernetes v1.28.0安装dashboard v2.6.1(k8s图形化操作界面)

准备工作 Kubernetes v1.28.0搭建教程请参考:Kubernetes v1.28.0集群快速搭建教程-CSDN博客 查看当前集群nodes都是ready状态 查看当前pods都是running状态 下载并修改配置文件 下载 recommended.yaml ,下载好之后,进入文件编辑 下载地址…

Spring 框架:配置缓存管理器、注解参数与过期时间

在 Spring 框架中,可通过多种方式配置缓存具体行为,常见配置方法如下。 1. 缓存管理器(CacheManager)配置 基于内存的缓存管理器配置(以SimpleCacheManager为例) SimpleCacheManager 是 Spring 提供的简单…

2025.1.21——六、BUU XSS COURSE 1 XSS漏洞|XSS平台搭建

题目来源:buuctf BUU XSS COURSE 1 目录 一、打开靶机,整理信息 二、解题思路 step 1:输入框尝试一下 step 2:开始xss注入 step 3:搭建平台 step 4:利用管理员cookie访问地址 三、小结 二编&#…

基于Python机器学习的双色球数据分析与预测

python统计分析2003-2024年所有的中奖记录,通过人工智能机器学习预测双色球,个人意见,仅供参考. 声明:双色球具有随机性,任何工具无法预测。本文章仅作为技术交流,提供学习参考。本文所涉及的代码均为python之机器学习的代码。双色球为公益事…

警企联动齐发力、共筑反诈“防护墙”

2025年1月10日是第五个中国人民警察节,南通移动联合南通公安反诈中心,深入社区商圈,开展防范电信网络诈骗宣传活动,进一步增强广大人民群众的反诈意识和能力,全力守护好群众的“钱袋子”。 当日,活动现场一大早就呈现出一片忙碌景象,工作人员支起摊位,将各类精心制作的反诈宣传…

mysql 学习3 SQL语句--整体概述。SQL通用语法;DDL创建数据库,查看当前数据库是那个,删除数据库,使用数据库;查看当前数据库有哪些表

SQL通用语法 SQL语句分类 DDL data definition language : 用来创建数据库,创建表,创建表中的字段,创建索引。因此成为 数据定义语言 DML data manipulation language 有了数据库和表以及字段后,那么我们就需要给这个表中 添加数…

九、CSS工程化方案

一、PostCSS介绍 二、PostCSS插件的使用 项目安装 - npm install postcss-cli 全局安装 - npm install postcss-cli -g postcss-cli地址:GitHub - postcss/postcss-cli: CLI for postcss postcss地址:GitHub - postcss/postcss: Transforming styles…

Baklib助力企业如何搭建内容中台的全新策略解析

内容概要 在当今快速变化的商业环境中,企业面临着数据爆炸和信息管理的巨大挑战。内容中台作为一种新兴的管理理念,专注于内容的集中管理与高效利用,旨在帮助企业优化信息流动,提高工作效率。通过构建一个灵活、可扩展的内容中台…

AR智慧点巡检系统探究和技术方案设计

一、项目背景 随着工业生产规模的不断扩大和设备复杂度的提升,传统的人工点巡检方式效率低下、易出错,难以满足现代化企业对设备运行可靠性和安全性的要求。AR(增强现实)技术的发展为点巡检工作带来了新的解决方案,通…

LabVIEW项目中的工控机与普通电脑选择

工控机(Industrial PC)与普通电脑在硬件设计、性能要求、稳定性、环境适应性等方面存在显著差异。了解这些区别对于在LabVIEW项目中选择合适的硬件至关重要。下面将详细分析这两种设备的主要差异,并为LabVIEW项目中的选择提供指导。 ​ 硬件设…

计算机图形学:实验四 带纹理的OBJ文件读取和显示

一、程序功能设计 在程序中读取带纹理的obj文件,载入相应的纹理图片文件,将带纹理的模型显示在程序窗口中。实现带纹理的OBJ文件读取与显示功能,具体设计如下: OBJ文件解析与数据存储 通过实现TriMesh类中的readObj函数&#x…

C# OpenCV机器视觉:红外体温检测

在一个骄阳似火的夏日,全球却被一场突如其来的疫情阴霾笼罩。阿强所在的小镇,平日里熙熙攘攘的街道变得冷冷清清,人们戴着口罩,行色匆匆,眼神中满是对病毒的恐惧。阿强作为镇上小有名气的科技达人,看着这一…

PAT甲级-1020 Tree Traversals

题目 题目大意 给出一棵树的后序遍历和中序遍历,要求输出该树的层序遍历。 思路 非常典型的树的构建与遍历问题。后序遍历和中序遍历可以得出一个树的结构,用递归锁定根节点,然后再遍历左右子树,我之前发过类似题目的博客&…

C语言进阶——3字符函数和字符串函数(2)

8 strsrt char * strstr ( const char *str1, const char * str2);查找子字符串 返回指向 str1 中第一次出现的 str2 的指针,如果 str2 不是 str1 的一部分,则返回 null 指针。匹配过程不包括终止 null 字符,但会在此处停止。 8.1 库函数s…

python学opencv|读取图像(四十二)使用cv2.add()函数实现多图像叠加

【1】引言 前序学习过程中,掌握了灰度图像和彩色图像的掩模操作: python学opencv|读取图像(九)用numpy创建黑白相间灰度图_numpy生成全黑图片-CSDN博客 python学opencv|读取图像(四十)掩模:三…

基于C语言的数组从入门到精通

简介:本篇文章主要介绍了一维数组,二维数组,字符数组的定义,数组的应用,数组的核心代码解析,适用于0基础的初学者. C语言数组 1.一维数组 1.1定义 1.1.1声明 语法:数据类型 数组名[数组大小];示例:int arr[5]; 1.1.2初始化 a.静态初始化 完全初始化:int arr[5] {1…

【kong gateway】5分钟快速上手kong gateway

kong gateway的请求响应示意图 安装 下载对应的docker 镜像 可以直接使用docker pull命令拉取,也可以从以下地址下载:kong gateway 3.9.0.0 docker 镜像 https://download.csdn.net/download/zhangshenglu1/90307400, postgres-13.tar http…

python 安装插件 requests 下载免费简历(自学7)

安装 requests 库: 他们三个 按一个就行 pip install requests 或者 pip3 install requests 或者 conda install requests 代码 每次只可以下载一页的 简历模板 需要手动修改 id ### import requests from lxml import etree import osif __name__ "__…