当前位置: 首页 > news >正文

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

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

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


相关文章:

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

1. 微分方程的基本概念 1. 微分方程的基本概念: (1)微分方程:含有未知函数及其导数或微分的方程。 举例说明微分方程:;。 (2)微分方程的阶:指微分方程中未知函数的导数…...

【探索 Kali Linux】渗透测试与网络安全的终极操作系统

探索 Kali Linux:渗透测试与网络安全的终极操作系统 在网络安全领域,Kali Linux 无疑是最受欢迎的操作系统之一。无论是专业的渗透测试人员、安全研究人员,还是对网络安全感兴趣的初学者,Kali Linux 都提供了强大的工具和灵活的环…...

四方连续贴图是什么意思

上下左右四个方向都是连续的图案,是由一个纹样或几个纹样组成一个单位,向四周重复地连续和延伸扩展而成的图案形式...

RKNN_C++版本-YOLOV5

1.背景 为了实现低延时,所以开始看看C版本的rknn的使用,确实有不足的地方,请指正(代码借鉴了rk官方的仓库文件)。 2.基本的操作流程 1.读取模型初始化 // 设置基本信息 // 在postprocess.h文件中定义,详见…...

k8s优雅重启

理论上处于terminating状态的pod,k8s 就会把它从service中移除了,只用配置一个优雅停机时长就行了。kubectl get endpoints 验证 因此,优雅重新的核心问题,是怎么让空闲长连接关闭,再等待处理中的请求执行完。 一些底…...

三高“高性能、高并发、高可靠”系统架构设计系列文章

目录 高并发系统的艺术:如何在流量洪峰中游刃有余 《数据密集型应用系统设计》读后感与高并发高性能实践案例 系统稳定性与高可用保障的几种思路 软件系统限流的底层原理解析 技术解决方案调研 延迟队列调研 重试调研 异步回调调研 分库分表调研 分布式事…...

opengrok_使用技巧

Searchhttps://xrefandroid.com/android-15.0.0_r1/https://xrefandroid.com/android-15.0.0_r1/ 选择搜索的目录(工程) 手动在下拉框中选择,或者 使用下面三个快捷按钮进行选择或者取消选择。 输入搜索的条件 搜索域说明 域 fullSearc…...

C++资料

InterviewGuide 首页 - 八股精 Releases halfrost/LeetCode-Go GitHub GitHub - GrindGold/CppGuide: 「C/C学习面试指南」一份涵盖大部分 C 程序员所需要掌握的知识。入门、进阶、深入、校招、社招,准备 C 学习& 面试,首选 CppGuide&#xff0…...

基于模糊PID的孵化箱温度控制系统(论文+源码)

1系统方案设计 本课题为基于模糊PID的孵化箱温度控制系统,其以STM32最小系统与模糊PID控制器为控制核心。系统主要包括数据采集模块、处理器模块、电机控制模块。 数据采集模块由温度传感器构成,通过温度传感器感应温度变化,获得待处理的数据…...

景联文科技加入AIIA联盟数据标注分委会

2025年1月16日,中国人工智能产业发展联盟(简称AIIA)数据委员会数据标注分委会(以下简称“分委会”)正式成立。景联文科技成为第一批AIIA联盟数据标注分委会委员单位。 数据标注分委会的成立旨在搭建数据标注领域产学研…...

1-1 飞机大战项目框架搭建

前言: 基于本人巩固C语言编写,仅供学习参考 1.0 框架搭建 搭建完成后状态 C语言飞机大战框架 使用loadimage时出现波浪线的错误如何解决,这个问题主要是编码格式不对造成的,我们需要修改系统的编码格式,将unicode编码格…...

【C++高并发服务器WebServer】-7:共享内存

本文目录 一、共享内存1.1 shmget函数1.2 shmat1.3 shmdt1.4 shmctl1.5 ftok1.6 共享内存和内存映射的关联1.7 小demo 二、共享内存操作命令 一、共享内存 共享内存允许两个或者多个进程共享物理内存的同一块区域(通常被称为段)。由于一个共享内存段会称…...

RabbitMQ 多种安装模式

文章目录 前言一、Windows 安装 RabbitMq1、版本关系2、Erlang2.1、下载安装 Erlang 23.12.2、配置 Erlang 环境变量 3、RabbitMQ3.1、下载安装 RabbitMQ 3.8.93.2、环境变量3.3、启动RabbitMQ 管理插件3.3、RabbitMQ3.4、注意事项 二、安装docker1、更新系统包:2、…...

C++ 包装器与绑定器的应用之回调函数的实现

回调函数的实现 在消息队列和网络库的框架中,当接收到消息(报文)时,回调用户自定义的函数对象,把消息(报文)参数传给它,由它决定如何处理。 queue参考文章:C queue(STL queue&…...

Baichuan大模型Base、Chat、Instruct等版本的区别

Baichuan大模型Base与Instruct等版本的区别解析 Baichuan大模型作为国内领先的开源语言模型,其不同版本(如Base、Chat、Instruct等)在训练目标、应用场景和性能特点上存在显著差异。以下是基于公开技术文档和行业分析的详细对比:…...

3.DrawCall的概念

DrawCall是渲染管线中的一个重要概念,指的是CPU向GPU发送的一个绘制命令,告诉GPU:“请根据我提供的数据,画一个物体(或一部分物体)。” 通俗易懂讲解:DrawCall就像给画师下订单 想象你是一个老…...

ubuntu电脑调用摄像头拍摄照片

一、 1、先装环境 conda create -n text python3.8 -y conda activate text 2、 pip install opencv-python -i https://pypi.tuna.tsinghua.edu.cn/simple 1、连接摄像头拍摄收集数据集 capture_image5.py import cv2 as cv import os import datetime import n…...

PyQt4 的图片切割编辑器

一、 编辑器功能明确 允许用户加载图片、选择切割模式、对切割后的图片片段进行操作(如移动、复制、粘贴、删除等),并支持撤销和重做操作。 环境:Py2.7 PyQt 4.11 二、导入模块介绍 sys: 用于访问与 Python 解释器强相关的变…...

mac 电脑上安装adb命令

在Mac下配置android adb命令环境,配置方式如下: 1、下载并安装IDE (android studio) Android Studio官网下载链接 详细的安装连接请参考 Mac 安装Android studio 2、配置环境 在安装完成之后,将android的adb工具所在…...

Webrtc (1) - Windows 编译

最近项目上遇到webrtc wgc 的几个test case无法通过,与webrtc人员沟通后决定要自行修复一下(因为他们不想管…) 参考文档 https://webrtc.org/support/contributinghttps://chromium.googlesource.com/chromium/src//main/docs/#checking-out-and-building 以上两…...

学习数据结构(1)算法复杂度

1.数据结构和算法 (1)数据结构是计算机存储、组织数据的方式,指相互之间存在⼀种或多种特定关系的数据元素的集合 (2)算法就是定义良好的计算过程,取一个或一组的值为输入,并产生出一个或一组…...

GCC之编译(8)AR打包命令

GCC之(8)AR二进制打包命令 Author: Once Day Date: 2025年1月23日 一位热衷于Linux学习和开发的菜鸟,试图谱写一场冒险之旅,也许终点只是一场白日梦… 漫漫长路,有人对你微笑过嘛… 全系列文章请查看专栏: Linux实践记录_Once-Day的博客-C…...

RocketMQ原理—4.消息读写的性能优化

大纲 1.Producer基于队列的消息分发机制 2.Producer基于Hash的有序消息分发 3.Broker如何实现高并发消息数据写入 4.RocketMQ读写队列的运作原理分析 5.Consumer拉取消息的流程原理分析 6.ConsumeQueue的随机位置读取需求分析 7.ConsumeQueue的物理存储结构设计 8.Cons…...

(Halcon)轮廓等分切割(项目分析)

目标:获取绿色圆所在位置(可用于点焊/点胶引导) 实现思路 一,相机标定板标定(如果实战用于点焊/点胶引导需要做图像畸变校正以减小误差) 相机标定 如何做一个C#仿Halcon Calibration插件-CSDN博客 二&…...

NIO 和 Netty 在 Spring Boot 中的集成与使用

Netty到底是个啥,有啥子作用 1. Netty 的本质:对 NIO 的封装 NIO 的原生问题: Java 的 NIO 提供了非阻塞 I/O 和多路复用机制,但其使用较为复杂(如 Selector、Channel、Buffer 的配置和管理)。开发者需要自…...

【更正版】梯级水光互补系统最大化可消纳电量期望短期优化调度模型

目录 1 主要内容 目标函数: 约束条件: 线性化处理: 流程示意: 2 部分代码 3 程序结果 4 下载链接 1 主要内容 该程序参考文献《梯级水光互补系统最大化可消纳电量期望短期优化调度模型》,构建了以最大化整体可…...

基于AnolisOS 8.6安装GmSSL 3.1.1及easy_gmssl库测试国密算法

测试环境 Virtual Box,AnolisOS-8.6-x86_64-minimal.iso,4 vCPU, 8G RAM, 60 vDisk。最小化安装。需联网。 系统环境 关闭防火墙 systemctl stop firewalld systemctl disable firewalld systemctl status firewalld selinux关闭 cat /etc/selinux/co…...

vue3 实际应用 将一个日期使用 moment.js 实现星期 今天 明天 ...

数据源 ["2025-01-23","2025-01-24","2025-01-25","2025-01-28","2025-01-26","2025-01-27" ] 后端给返回了一个这样的数据 日期数据 实际应用中实现的是这样的显示效果 日期需要这样显示的tabs 栏 我们需要…...

LLM幻觉(Hallucination)缓解技术综述与展望

LLMs 中的幻觉问题(LLM 幻觉:现象剖析、影响与应对策略)对其可靠性与实用性构成了严重威胁。幻觉现象表现为模型生成的内容与事实严重不符,在医疗、金融、法律等对准确性要求极高的关键领域,可能引发误导性后果&#x…...

Unity入门2 背景叠层 瓦片规则

切割场景 瓦片调色盘 放在Assets里面新建瓦片地图,palettes tile 瓦片 palettes调色板 上下窗口是分开的 拖进这个格子窗 瓦片太碎,要封装 装好之后,只是把瓦片放上去了,但是还没有画布,显示是这样的 no valid target 新建“…...