




版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
本文格式為Word版,下載可任意編輯——非齊次微分方程特解二階常系數線性非齊次微分方程特解的求法討
論
幸克堅
(遵義師范學院貴州遵義563002)
摘要:非數學專業《常微分方程》中,“二階常系數線性微分方程〞一般是作為一個單獨的模塊來講授。但在一般非數學專業使用的《高等數學》教材中,特解的介紹往往比較突然和不夠完整,引入不大自然也不易于理解和接受。本文結合非數學專業學生的特點就特解的求法進行了分析和探討。
(二)其次考慮:y?+py?+qy=f(x)中f(x)=asinωx(或f(x)=acosωx)的情形:
即:y?+py?+qy=asinωx(不失一般,取f(x)=asinωx即可)
由于復角“ωx〞的正、余弦函數的微商(或導數)依舊是復角“ωx〞的正、余弦函數,并且y、y?與y?中sinωx與cosωx總是交替出現的。因此,要使方程y?+py?+qy=
**
asinωx的特解y代入原方程后等式成立,y應當形如:*
y=αcosωx+βsinωx
*
將y=αcosωx+βsinωx代入原方程得:
(αcosωx+βsinωx)?+p(αcosωx+βsinωx)?+q(αcosωx+βsinωx)=asinωx求導整理得:
22
—αωcosωx—βωsinωx—pαωsinωx+pβωcosωx+qαcosωx+qβsinωx=asinωx比較系數得:
2(q—ω)α+pωβ=0
2(q—ω)β—pωα=a
2
所以當pω≠0且q—ω≠0時,特解中待定系數α、β可以由原方程的系數p、q、a、ω通過關系式:
?pωaα=β222(q?ω)?(pω)(q?ω2)a=
(q?ω2)2?(pω)2
唯一確定。
22
而pω=0與q—ω=0正好與:(±ωi)+p(±ωi)+q=0等價,即±ωi是特征方程:
2
λ+pλ+q=0的根2
所以,當±ωi不是特征方程λ+pλ+q=0的根時,原方程的特解由原方程的系數p、q、a、ω通過上式唯一確定。
2
但當±ωi是特征方程:λ+pλ+q=0的根時,pω=0
與q—ω=0,此時,若再設:
*
y=αcosωx+βsinωx將無法確定α、β之值。考慮到這一結果正好是由“αcosωx+βsinωx〞及其微商導致的,故可參照f(x)為一次多項式時的情形,考慮設:
*
y=(αcosωx+βsinωx)x
*
則:(y)′=((αcosωx+βsinωx)x)′=(-αωsinωx+βωcosωx)x+(αcosωx+βsinωx)
*
(y)″=((-αωsinωx+βωcosωx)x+(αcosωx+βsinωx))′
22
=(-αωcosωx-βωsinωx)x+2(-
αωsinωx+βωcosωx)代入原方程得:
***y″+py′+qy
22
=(-αωcosωx-βωsinωx)x+2(-α
ωsinωx+βωcosωx)+
+p((-αωsinωx+βωcosωx)x+(α
cosωx+βsinωx))+q(αcosωx+βsinωx)x)
=asinωx比較系數得:
22
(-αωcosωx-βωsinωx)+p(-αωsinω
x+βωcosωx)+q(αcosωx+βsinωx)=0
-2ωα+pβ=apα+2ωβ=0
2
因±ωi是特征方程λ+pλ+q=0的根,所以pω=0與q—ω2
=0,于是前面一式中的α、β可以為任意值,而由后二式解得:α=?22?a2β=2pa2
4??p4??p2
2
所以,當±ωi是特征方程:λ+pλ+q=0的根時,α、β可由原方程的系數通過上式得出。
綜上所述,對于:y?+py?+qy=asinωx形式的非齊次方程,其特解為:
1)當±ωi不是特征方程λ+pλ+q=0的根時,則原
*
方程的特解是:y=αcosωx+βsinωx,
其中α、β由原方程的系數p、q、a、ω通過:
?pωaα=β222(q?ω)?(pω)2
(q?ω2)a=
(q?ω2)2?(pω)2唯一確定。
2
2)當±ωi是特征方程:λ+pλ+q=0的根時,原方程
*
的特解是:y=(αcosωx+βsinωx)x,其中α、β可由原方程的系數通過:
paα=?22?a2β=224??p4??p唯一確定。
而對于:y?+py?+qy=bcosωx形式的方程,可以類似地探討,這里就不贅述了。
bx
(三)最終考慮:y?+py?+qy=f(x)中f(x)=ae即
bx
y?+py?+qy=ae的情形:
bxbx
由于e的各階微商及其原函數均為e的常數倍,所
bx
以,方程y?+py?+qy=ae的特解也應當形如:
*bx
y=αe
2bxbxbx
將其代入原方程求導整理得:αbe+pbαe+qαe=bxae
bx2
由于e≠0,所以由上式得:(b+pb+q)α=a
2
所以,當b+pb+q≠0,即b不是特征方程:
2
λ+pλ+q=0的根時,α=2a
b?pb?q*bx
方程的特解為:y=αe=2ab?pb?q2
2
e
bx
而當b+pb+q=0,即b是特征方程λ+pλ+q=0的根時,由此式不能確定α的值,
*
但若仍參照f(x)為一次多項式時的情形,考慮設:y=bx
αxe代入原方程求導整理得:(b+pb+q)αxe+(2b+p)αe=ae
bx
bx
2
2bxbx
兩邊約去(e≠0)并比較系數得:b+pb+q=0且(2b+p)α=a
2
前一式說明b是特征方程:λ+pλ+q=0的根,
2
于是,在2b+p≠0,即b只是特征方程λ+pλ+q=0的單根,而不是重根的狀況下有:
*
α=a即原方程的特解為:y=α
2b?pxe=
bx
bxaxe2b?p2
在2b+p=0,即b是特征方程λ+pλ+q=0的重根的狀況下,由上式無法確定α之值,但也參照
*2bx
f(x)為多項式時的情形,可考慮設y=αxe代入原方程求導整理得:
***bxbx2y?+py?+qy==2αe+(4b+2p)αxe+(b+2bxbx
pb+q)αxe=ae
即:2α=a,α=a即原方程的特解為:
2y=αxe=axe
*
2
bx
2bx
2綜上所述,在f(x)=ae的狀況下:
2
1)若b不是特征方程λ+pλ+q=0的根,則特解為:
*bxbxy=αe=2ae
b?pb?qbx
2)若b是特征方程λ+pλ+q=0的單根,則特解為:*bxbxy=αxe=axe
2b?p2
3)若b是特征方程λ+pλ+q=0的重根,則特解為:
*2bx
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
評論
0/150
提交評論