常秀芳

（山西大同大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，山西大同037009）

變系數(shù)的線性微分方程的求解一般沒有固定方法可循，然而現(xiàn)實(shí)生活中常常遇到諸如熱傳導(dǎo)、振動(dòng)諧波、電磁波、等變系數(shù)的線性微分方程問題，它們正好是特殊的變系數(shù)的線性微分方程——?dú)W拉方程[1-4]。

1 概念

定義形如

的方程，稱為歐拉方程。其中：p1,p2,…,pn-1,pn是常 數(shù)且p1≠0, 當(dāng)f(x)≡0 時(shí) ，方 程xny(n)+p1xn-1xy(n-1)+… +pn-1xy′+pny=0 稱為n階線性齊次歐拉方程。

當(dāng)f(x)≠0 時(shí)，方程xny(n)+p1xn-1y(n-1)x+…+pn-1xy′+pny=f(x)

稱為n階線性非齊次歐拉方程。

歐拉方程用微分算子表示為

由解的結(jié)構(gòu)知：只要能求得n階線性齊次歐拉的通解以及n階線性非齊次歐拉方程一個(gè)特解，則歐拉方程的通解即可寫出。至于特解的形式可參閱n階常系數(shù)線性微分方程的特解形式進(jìn)行假設(shè)，用代入法求得，或知齊次方程的通解用常數(shù)變易法解之，因此，關(guān)鍵是求n階線性齊次歐拉的通解。

2 特點(diǎn)

歐拉方程是一個(gè)n階變系數(shù)線性微分方程，其特點(diǎn)：（1）它的左邊的每一項(xiàng)是都由冪函數(shù)與未知函數(shù)導(dǎo)數(shù)的乘積組成；（2）每一項(xiàng)冪函數(shù)的指數(shù)與未知函數(shù)導(dǎo)數(shù)的階數(shù)相等；（3）k階導(dǎo)數(shù)Dky的系數(shù)是k次的冪函數(shù)pn-kxk(k=0,1,2,…,n，且P0=1)。

3 一般解法

當(dāng)n=1時(shí)，是一階歐拉方程

變形為

則此方程為一階線性微分方程，其通解是

當(dāng)n≠1時(shí)，設(shè)x=eu，因

將上述各項(xiàng)代入方程

得一個(gè)關(guān)于未知函數(shù)是y，自變量為u的常系數(shù)線性方程

Pn(D)y=0，

求得此常系數(shù)線性齊次方程后，因x=eu，所以用u=lnx回代，則得齊次歐拉的解。

例1求方程x3y?+3x2y″+xy′-y=xlnx的解。

解設(shè)x=eu，因原方程為

[D(D-1)(D-2)+3D(D-1)+D-1]y=ueu，即

(D3-1)y=ueu，

其特征方程為r3-1=0，特征根為

則(D3-1)y=0 的通解為

設(shè)方程 (D3-1)y=ueu的特解為y?=u(Au+B)eu，代入方程得

(6Au+6A+3B)eu=ueu，

因此方程(D3-1)y=ueu的通解為

則所求方程的通解為

4 獨(dú)特解法

由冪函數(shù)導(dǎo)數(shù)仍為冪函數(shù)的特點(diǎn)，不妨設(shè)歐拉方程

代入原方程為

由于xλ≠0，則得一個(gè)關(guān)于λ的n次一元方程

不妨規(guī)定此方程為歐拉方程的特征方程。

4.1 特征方程有n個(gè)不同的特征根

設(shè)歐拉方程的特征方程有n個(gè)不同的特征根為λ1,λ2,…,λn。

因λ=λ1,λ2,…,λn而 且λ1,λ2,…,λn互 不 相等，于是xλ1,xλ2,…,xλn線性無關(guān)。則歐拉方程的通解為

解特征方程λ(λ-1)+λ-1=0 ，特征根為λ=±1，所以，所求方程的通解為

4.2 特征方程有兩個(gè)同的特征根

設(shè)歐拉方程的特征方程有兩n個(gè)相同的特征根為λ1、λ2。

因λ=λ1=λ2，所以y1=xλ1是齊次歐拉方程的解。

設(shè)y2=xλ1lnx，因

由于λ1是特征方程的重根，將代入齊次歐拉方程的左端，化簡得

因此，y2=xλ1lnx也是齊次歐拉方程的解。又y=xλ1與y=xλ1

lnx線性無關(guān)，所以歐拉方程的通解含有 (C1+C2lnx)xλ1項(xiàng)。

例3求方程x2y″-xy′+y=0 的通解。

綜合考慮研究區(qū)6個(gè)時(shí)相遙感影像，決定在NDVI概率累計(jì)表上取概率為99.5%的值為NDVImax，取概率為0.5%的為NDVImin，并利用ERDAS IMAGINE 2013 軟件中的Modeler實(shí)現(xiàn)植被覆蓋度定量轉(zhuǎn)換模型，得到1989—2015年6期植被覆蓋度專題圖。

解特征方程λ(λ-1)-λ+1=0，即(λ-1)2=0特征根為λ=1 ，所以，所求方程的通解為y=(C1+C2lnx)x。

4.3 特征方程有共軛復(fù)數(shù)的特征根

設(shè)歐拉方程的特征方程有一對(duì)的共軛復(fù)數(shù)的特征根為α±iβ。

因λ=α±iβ，則是歐拉方程的解。

由歐拉公式知

由線性齊次微分方程解的結(jié)構(gòu)知：

又y1與y2線性無關(guān)，所以齊次歐拉方程的通解含有xα[C1cos(βlnx)+C2sin(βlnx)]項(xiàng)。

例4求方程x2y″+xy′+y=0 的通解。

解特征方程λ(λ-1)+λ+1=0，即λ2+1=0，

特征根為λ=±i，所以，所求方程的通解為y=C1cos(lnx)+C2sin(lnx)。

4.4 特征方程有k重共扼復(fù)數(shù)的特征根

若λ=α±iβ是歐拉方程的k重共軛復(fù)數(shù)的特征根，則歐拉方程的通解一定含有

5 求齊次歐拉方程通解的具體步驟

第一步：依據(jù)方程寫出特征方程

第二步：求出特征根；

第三步：依特征根寫出方程的通解如表1。

表1 齊次歐拉方程特征根與通解關(guān)系的對(duì)照表