加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

二阶常系数非齐次线性方程的通解

(2013-06-11 23:09:40)
标签:

二阶常系数非齐次线性

通解

特解

解方程

高等数学

二阶常系数非齐次线性微分方程的一般形式:y¢¢+py¢+qy=f(x)其中pq均为常数。我们知道,解微分方程的过程就是去导数的过程,所谓微分方程的解就是一个y与x的函数关系即y=g(x),所以解本身也是个函数。对于齐次方程,必有零解,如果有非零解,则解有无穷多个,这无穷多个解是通过方程组列向量的相关关系,乘以某个任意常数得到的,所以这些解并不是杂乱无章的,是可以用某个表达式统一表达出来的,因此我们更青睐于这个统一的表达式,人们命名它为齐次方程的通解;对于非齐次方程,各个列向量需要以某种数量关系,搭配出等号右边的常数列,当然有一种可能解释怎么也搭配不出来,即无解,另外如果非齐次方程有解的话,则必是唯一解。也就是说,通解有很多很多的方程组成,而特解只有一个方程。

由解的结构我们可以知道,这种方程的解的结构是:对应其次方程的通解+非齐次方程的某个特解。那么,我们解方程的步骤有二:首先求齐的通解,再求非齐的特解。

1、y¢¢+py¢+qy=0通解

在求解这类方程的时候有一种“预判”的思想,非常有意思,我还不知道这种方式是否具有严密性,但我想这种方式之所以在方程组求解的过程中行得通,应该是解的结构在支撑,也就是说只要找到原方程众多杂乱无章的解中线性无关的两个,就可以以通解的形式,将所有的解表示出来,当然有无数种表示方法,你怎么找都可以,只要找到两个且线性无关就行。

所以我们知道,如果y1y2是二阶常系数齐次线性微分方程的两个线性无关解, 那么y=C1y1+C2y2就是它的通解那么,首先预判当λ为常数时,指数函数y=eλx和它的各阶导数都只相差一个常数因子。由于指数函数的这个特点,我们可以用y=eλx来尝试,看能否找到适当的常数λ使y=eλx  满足二阶常系数齐次线性微分方程为此将y=eλx代入方程y¢¢+py¢+qy=0(λ 2+pλ+q)eλx =0,因为eλx不为零,由此可见, 只要λ满足代数方程λ2+pλ+q=0, 函数y=eλx就是微分方程的解

虽然方程λ2+pλ+q=0仅仅是原方程的一部分,但它的根也是原方程的根,所以我们称其微分方程y¢¢+py¢+qy=0特征方程”。特征方程的两个根λ1λ2可用公式λ1,2=[-p±(p^2-4q)^(1/2)]/2求 那么我们的关注点也就从原方程转而集中到了其内部的特征方程。

(1)特征方程有两个不相等的实根λ1λ2, 函数eλ1xeλ2x是方程的两个线性无关的解。因此方程的通解为y=c1eλ1x+c2eλ2x 

(2)特征方程有两个相等的实根λ1=λ2, 函数eλx是原方程的一个解,但要求出通解,还需要一个与其线性无关的另外一个解。这个解是谁呢?我们依然使用预判的方法。设y2/y1=u(x),即y2=u(x)eλx ,将y2带回原方程得u¢¢(x)=0,那么取u的最简形式x即可满足条件,即y2=xeλx是原方程的另外一个解,二阶常系数非齐次线性方程的通解不是常数,两解线性无关。因此方程的通解二阶常系数非齐次线性方程的通解

特征方程有一对共轭复根λ1, 2=a±ib, 函数y=e(a+ib)xy=e(a-ib)x是微分方程的两个线性无关的复数形式的解。但为了避免复数形式的出现,我们用欧拉复变公式eix=cosx+isinx将原型是变换, y1=e(a+ib)x=eax(cosbx+isinbx), 

  y2=e(a-ib)x=eax(cosbx-isinbx),

 y1+y2=2eaxcosbx, 二阶常系数非齐次线性方程的通解,

    y1-y2=2ieaxsinbx, 二阶常系数非齐次线性方程的通解. 

eaxcosbxy2=eaxsinbx也是方程解. 可以验证, y1=eaxcosbxy2=eaxsinbx是方程的线性无关解因此方程的通解为y=eax(C1cosbx+C2sinbx ) 

 

 2、求y¢¢+py¢+qy=f(x)的一个特解

 

这里f(x)有两种特殊形式 

   2.1  f(x)=Pm(x)elx 

    f(x)=Pm(x)elx,f(x)m多项式与指数函数的乘积,我们进行预判,多项式与指数函数乘积的导数仍然是多项式与指数函数的乘积,所以方程的特解也应具有这种形式因此, 设特解形式为y*=Q(x)elx, 将其代入方程, 因为eλx不为零,得等式Q¢¢(x)+(2l+p)Q¢(x)+(l2+pl+q)Q(x)=Pm(x)。依然是对这个特征方程进行求解。非齐是建立在齐的基础之上的,我们两次的预判解和带入后的方程之间是有联系的。两次的预判都含有lelx,那么这两个l什么关系?

(1)如果非齐l不是齐的特征方程λ2+pλ+q=的根, l2+pl+q¹0. 要使上式成立, Q(x)应设为m 次多项式: 

        Qm(x)=b0xm+b1xm-1+ × × × +bm-1x+bm , 

通过比较等式两边同次项系数, 可确定b0, b1, × × × , bm, 并得所求特解 

        y*=Qm(x)elx. 

    (2)如果非齐l是齐的特征方程 λ2+pλ+q=的单根, l2+pl+q=0, 2l+p¹0, 要使等式

        Q¢¢(x)+(2l+p)Q¢(x)+(l2+pl+q)Q(x)=Pm(x). 

成立, Q(x)应设为m+次多项式: 

        Q(x)=xQm(x), 

        Qm(x)=b0xm +b1xm-1+ × × ×  +bm-1x+bm , 

通过比较等式两边同次项系数, 可确定b0, b1, × × ×  , bm, 并得所求特解  

        y*=xQm(x)elx. 

    (3)如果非齐l是齐的特征方程 λ2+pλ+q=0的二重根, l2+pl+q=0, 2l+p=0, 要使等式

        Q¢¢(x)+(2l+p)Q¢(x)+(l2+pl+q)Q(x)=Pm(x). 

成立, Q(x)应设为m+2次多项式: 

        Q(x)=x2Qm(x), 

        Qm(x)=b0xm+b1xm-1+ × × × +bm-1x+bm , 

通过比较等式两边同次项系数, 可确定b0, b1, × × × , bm , 并得所求特解

        y*=x2Qm(x)elx. 

    综上所述, 我们有如下结论: 如果f(x)=Pm(x)elx, 则二阶常系数非齐次线性微分方程y¢¢+py¢+qy =f(x)有形如

        y*=xk Qm(x)elx

的特解, 其中Qm(x)是与Pm(x)同次的多项式, k l不是特征方程的根、是特征方程的单根或是特征方程的的重根依次取为012. 

2.2方程y¢¢+py¢+qy=elx[Pl (x)coswx+Pn(x)sinwx]的特解形式

    应用欧拉公式可得

    elx[Pl(x)coswx+Pn(x)sinwx]

    二阶常系数非齐次线性方程的通解

    二阶常系数非齐次线性方程的通解 

    二阶常系数非齐次线性方程的通解, 

其中二阶常系数非齐次线性方程的通解, 二阶常系数非齐次线性方程的通解. m=max{l, n}. 

    设方程y¢¢+py¢+qy=P(x)e(l+iw)x的特解为y1*=xkQm(x)e(l+iw)x, 

二阶常系数非齐次线性方程的通解必是方程二阶常系数非齐次线性方程的通解的特解, 

其中kl±iw不是特征方程的根或是特征方程的根依次取01. 

    于是方程y¢¢+py¢+qy=elx[Pl(x)coswx+Pn(x)sinwx]的特解为

        二阶常系数非齐次线性方程的通解

          二阶常系数非齐次线性方程的通解

          =xk elx[λ(1)m(x)coswx+λ(2)m(x)sinwx]. 

    综上所述, 我们有如下结论: 

    如果f(x)=elx [Pl(x)coswx+Pn(x)sinwx], 则二阶常系数非齐次线性微分方程

        y¢¢+py¢+qy=f(x)

的特解可设为

        y*=xk elx[λ(1)m(x)coswx+λ(2)m(x)sinwx], 

其中λ(1)m(x)λ(2)m(x)m次多项式, m=max{l, n}, k l+iw (l-iw)不是特征方程的根或是特征方程的单根依次取01. 

 

 标准化方法如下:

1、求齐次微分方程的通解

   列出特征方程 ...=0

   若b^2-4ac>0,则通解为

   若b^2-4ac>0,则通解为

   

 

 

 

 


0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有