謝英超，程 燕(．陸軍軍官學(xué)院研究生管理大隊，安徽合肥 3003;．陸軍軍官學(xué)院基礎(chǔ)部，安徽合肥 3003)

?

一類非線性時滯微分方程的奇異攝動研究

謝英超1，程 燕2
(1．陸軍軍官學(xué)院研究生管理大隊，安徽合肥 230031;2．陸軍軍官學(xué)院基礎(chǔ)部，安徽合肥 230031)

摘 要:研究了一類依賴于小參數(shù)的小時滯微分方程．首先利用擬合函數(shù)法將雙參數(shù)問題轉(zhuǎn)換為便于分析的單參數(shù)問題，再利用校正函數(shù)法得到了方程的一致有效的漸近解，并利用微分不等式理論給出了證明，最后將其與數(shù)值解進行了精度比較．結(jié)果表明該攝動方法是有效的，從而可以更好地分析這類方程的解的性態(tài)．

關(guān)鍵詞:奇異攝動;時滯微分方程;非線性;漸近解

引用格式:謝英超，程燕．一類非線性時滯微分方程的奇異攝動研究［J］．安徽師范大學(xué)學(xué)報:自然科學(xué)版，2015，38(2) :129－133．

引言

非線性時滯微分方程的研究目前仍然是國內(nèi)外工程界和學(xué)術(shù)界十分關(guān)注的重點和熱點之一［1］，在許多領(lǐng)域都有重要的作用［2］．許多學(xué)者利用奇異攝動的理論和方法研究時滯微分方程，取得了可喜的成果［3－15］．但對于依賴于小參數(shù)的小時滯微分方程研究較少，本文對解這類方程進行了一定的探索．

考慮如下一類依賴于小參數(shù)的小時滯微分方程:

其中:ε是正的小參數(shù)，τ是正的小時滯參數(shù)，Τ為適當(dāng)?shù)恼龜?shù)，φ(t)為初始函數(shù);記［u(t)］= u(t－τ)．問題(1)－(2)是一個具有兩參數(shù)的非線性時滯微分方程的奇異攝動問題．

為研究方便，現(xiàn)作如下假設(shè):

［H1］當(dāng)τ→0時，ε= O(τ) ;

［H2］f和φ在各自的定義區(qū)間內(nèi)為充分光滑函數(shù)，且存在常數(shù)δ1≥0，δ2＞0，使得fu(t，u，［u］)≤－δ1≤0，f［u］(t，u，［u］)≤－δ2＜0，(t，u，［u］)∈［0，T］×R×R; ［H3］退化問題f(t，u(t)，u(t) ) = 0在t∈［0，T］存在唯一的連續(xù)可微解u*(t)．

1　問題的奇異攝動研究

問題(1)－(2)為依賴于小參數(shù)的小時滯微分方程，由于其奇異的雙重性，研究起來較為困難．本文考慮用一個關(guān)于時滯參數(shù)的光滑函數(shù)去擬合這個小參數(shù)，從而轉(zhuǎn)換為只依賴于小時滯參數(shù)的問題．事實上，對于實際模型、工程應(yīng)用等問題，人們更關(guān)心的是給定參數(shù)下的方程的解，所以這種方法在實際應(yīng)用中具有一定的活力．

由假設(shè)［H1］知，當(dāng)ε∈(0，ε0)時，其中ε0為適當(dāng)小的正常數(shù)，ε可以表示為τ的冪級數(shù)之和，即存在足夠小的正常數(shù)τ0，當(dāng)τ∈(0，τ0)時，可選擇適當(dāng)?shù)淖銐蚬饣暮瘮?shù)g(τ)，使得其中g(shù)i= g(i)(0)為適當(dāng)?shù)某?shù)，且g(0) = 0，g1≠0．

對于t∈［0，T］，尋找問題(1)－(2)具有如下形式的解

其中U(t，τ)為外部解，設(shè)其形式為

滿足時滯微分方程(1) ;ξ= t/τ為伸長變量，V(ξ，τ)為初始層校正項，設(shè)其形式為

2　外部解

首先按τ展開u(t－τ)2

將式(3)、(4)和(7)代入(1)式，按τ展開方程兩邊，合并τ的同次冪項．關(guān)于τ0，可得:

f(t，u0(t)，u0(t) ) = 0．

(8)由假設(shè)［H3］知，方程(8)存在唯一解，且

u0(t) = u*(t)．

(9)關(guān)于τi(i = 1，2，…)的系數(shù)，記fu= fu(t，u0(t)，u0(t) )，f［u］= f［u］(t，u0(t)，u0(t) )可分別得到: (fu+ f［u］) ui= Gi－Fi，

(10)

，…，故Gi為逐次已知的函數(shù)，可以證明Fi也為逐次已知的函數(shù)，特別地

其中fuu= fuu(t，u0(t)，u0(t) )，fu［u］= fu［u］(t，u0(t)，u0(t) )，f［u］［u］(t，u0(t)，u0(t) )．

方程(10)退化為代數(shù)方程，由假設(shè)［H2］可知，(fu+ f［u］)≠0，故可依次求得其唯一解為ui(t) (t∈［0，T］)，i = 1，2，…．至此，外部解已經(jīng)構(gòu)造完成，但它未必滿足初始條件(2)，故還需要構(gòu)造初始層校正項．

3　初始層校正項和解的漸近表達式

令t = 0代入式(4)，并由條件(2)可得

因為完全展開式(4)和外部展開式(5)滿足時滯微分方程(1)，所以初始層校正項V(ξ，τ)必須滿足

對于初始層校正項vi(ξ)可以通過逐步積分法，在區(qū)間s≤ξ≤s + 1(s = 0，1，2，…)上逐次求得．由式(3)、(5)、(6)、(11)、(12)和(13)，關(guān)于τ0，v0(ξ)，必為

在初始條件v0(0) =φ(0)－u0(0)的連續(xù)解，通過逐段積分可求得v0(ξ)．在滿足假設(shè)［H1］和［H2］的前提下，可以證明v0(ξ)是當(dāng)ξ→∞時指數(shù)型地趨于零的函數(shù)．關(guān)于τ1，v1(ξ)必為

在初始條件v1(0) =－u1(0)的連續(xù)解，通過逐段積分可求得v1(ξ)．在滿足假設(shè)［H1］和［H2］的前提下，可以證明v1(ξ)也是當(dāng)ξ→∞時指數(shù)型地趨于零的函數(shù)．

關(guān)于τi(i = 2，3，…)，vi(ξ)必為

在初始條件vi(0) =－ui(0)的連續(xù)解，其中，為逐次已知的函數(shù)，可以證明i和也為逐次已知的函數(shù)．在滿足假設(shè)［H］1

和［H2］的前提下，可以證明vi(ξ)也是當(dāng)ξ→∞時指數(shù)型地趨于零的函數(shù)．’’至此，外部解和初始層校正項構(gòu)造完成，得到問題(1)－(2)解的形式漸近展開式為

4　漸近解的一致有效性

我們有如下定理:

定理 在滿足假設(shè)［H1］－［H2］的前提下，當(dāng)ε∈(0，ε0)，τ∈(0，τ0)時，依賴于小參數(shù)的小時滯微分方程初值問題(1)－(2)存在一個解u，并具有形如式(17)的一致有效的漸近展開式．

證明 首先構(gòu)造如下形式的輔助函數(shù)α(t，τ)和β(t，τ)

(18) (19)其中為r足夠大的正數(shù)，在證明中將給定，

顯然有

對足夠大r的和τ∈(0，τ0)，有

下面證

事實上，將式(3)、(18)和(20)代入(23)，對于0≤ξ≤1和τ∈(0，τ0)，存在一個常數(shù)R＞0，有

由式(21)、(22)、(23)和(24)，并利用微分不等式理論［16］，可得

由此可得，依賴于小參數(shù)的小時滯微分方程初值問題(1)－(2)存在一個解u，并具有形如式(17)的一致有效的漸近展開式．定理證畢．

5　漸近解與數(shù)值解的精度比較

為比較所求漸近解與數(shù)值解的精度，舉一個特殊的例子來說明，取定問題(1)－(2)的一組參數(shù)為:ε= 0．01，τ= 0．01，并且f(t，u(t)，［u(t)］) = 2t－u(t)－［u(t)］，φ(t) = ei，選取適當(dāng)函數(shù)g(τ) =τ，則g1= 1，gk= 0(k = 2，3，…)．用上述攝動方法，可得問題(1)－(2)的外部解為

u0(t) = t，ui(t) = 0(i = 1，2，…)初始層校正項v0(ξ)為

在初始條件v0(0) = 1的連續(xù)解，通過逐段積分可求得v0(ξ)．故解的一次漸近展開式為

問題(1)－(2)的攝動解與數(shù)值解在ε=τ= 0．01，φ(t) = ei時的變化曲線如圖1所示．

由圖1可見，攝動解與數(shù)值解曲線非常接近，表明所述攝動方法是一個簡單而有效的方法．

6　結(jié)束語

本文研究了一類依賴于小參數(shù)的小時滯微分方程，在一定假設(shè)條件下對這類方程的解進行了探索，利用擬合函數(shù)法和校正函數(shù)法得到了方程的一致有效的漸近解．結(jié)果表明本文所述方法是一個簡單而有效的近似解析方法，不同于數(shù)值解法，它還能繼續(xù)進行解析運算．所以，能夠利用得到的漸近展開式更好地對這類方程解的性態(tài)進行更深層次的分析和研究，從而能更好地應(yīng)用到物理、機械和工程技術(shù)等領(lǐng)域中．

參考文獻:

［1］ 徐鑒，裴利軍．時滯系統(tǒng)動力學(xué)近期研究進展與展望［J］．力學(xué)進展，2006，36(1) :17－29．

［2］魏俊杰，王洪濱，蔣衛(wèi)華．時滯微分方程的分支理論及應(yīng)用［M］．北京:科學(xué)出版社，2012．

［3］劉樹德，魯世平，姚靜蓀，等．奇異攝動邊界層和內(nèi)層理論［M］．北京:科學(xué)出版社，2012．

［4］ 歐陽成，林萬濤，程榮軍，等．一類厄爾尼諾海－氣時滯振子的漸近解［J］．物理學(xué)報，2013，62(6) :060201．

［5］ 莫嘉琪，林萬濤，林一驊．厄爾尼諾/拉尼娜－南方濤動機制時滯海－氣振子的漸近解［J］．物理學(xué)報，2011，60(8) :080202．

［6］ 朱敏，林萬濤，林一驊，等．一類厄爾尼諾時滯海－氣振子攝動解［J］．物理學(xué)報，2011，60(3) :030204．

［7］ 莫嘉琪，溫朝暉．奇攝動時滯反應(yīng)擴散方程［J］．應(yīng)用數(shù)學(xué)和力學(xué)，2010，31(6) :739－744．

［8］ 汪娜，倪明康．經(jīng)典物理中的擾動時滯模型解［J］．物理學(xué)報，2011，60(5) :050203．

［9］ 趙永祥，肖愛國．兩步W－方法關(guān)于時滯奇異攝動初值問題的誤差分析［J］．?dāng)?shù)學(xué)物理學(xué)報，2011，31A(5) :1239－1252．

［10］ 任景莉，葛渭高．具非線性邊界條件的半線性時滯微分方程邊值問題奇攝動［J］．應(yīng)用數(shù)學(xué)和力學(xué)，2003，24(12) :1285－1290．

［11］ 任景莉，葛渭高．具非線性邊界條件的Volterra型時滯微分方程邊值問題奇攝動［J］．?dāng)?shù)學(xué)物理學(xué)報，2003，23A(4) :504－512．

［12］ 魯世平．具非線性邊界條件的Volterra型泛函微分方程邊值問題奇攝動［J］．應(yīng)用數(shù)學(xué)和力學(xué)，2003，24(12) :1276－1284．

［13］ 魯世平．奇攝動非線性時滯微分方程邊值問題［J］．?dāng)?shù)學(xué)研究與評論，2003，23(2) :304－308．

［14］ 歐陽成．具有小延遲的微分－差分方程漸近解［J］．吉林大學(xué)學(xué)報:理學(xué)版，2008，46(4) :628－632．

［15］ AGARWAL R P，CAPASSO，Ch G，PHILOS P，TSAMATOS Ch．Global solutions of a singular initial value problem to second order nonlinear delay differential equations［J］．Mathematical and Computer Modelling，2006(43) :854－869．

［16］周明儒，杜增吉，王廣瓦．奇異攝動中的微分不等式理論［M］．北京:科學(xué)出版社，2012．

Singularly Perturbed Study on a Class of Nonlinear Delay Differential Equation

XIE Ying-chao1，CHENG Yan2
(1．Graduate Management Unit，Army Officer Academy，Hefei 230031，China;2．Basic department，Army Officer Academy，Hefei 230031，China)

Abstract:In this paper，a class of small delay differential equation depending on a small parameter is considered．Firstly using the fitting function method，the two parameters problem is converted to a single parameter problem for facilitating analysis．Then using the correction function method，the uniformly valid asymptotic solution of the equation is obtained，which is proved by using the theory of differential inequalities．Finally it’s compared with the numerical solution．The result shows that the perturbation method is effective，so we can analyze the behavior of solution for the class of delay differential equation much better．

Key words:singular perturbation; delay differential equation; nonlinear; asymptotic solution

作者簡介:謝英超(1989－)，男，碩士研究生，研究方向:應(yīng)用數(shù)學(xué)．

基金項目:國家自然科學(xué)基金項目(11202106)，安徽省自然科學(xué)基金項目(1408085MA06)．

收稿日期:2013－12－12

DOI:10．14182/J．cnki．1001－2443．2015．02．005

文章編號:1001－2443(2015) 02－0129－05

文獻標志碼:A

中圖分類號:O175．7