曲婧佳，周 冉，黃開銀

（1.空軍航空大學(xué) 基礎(chǔ)部，長春130022；2.吉林大學(xué) 數(shù)學(xué)學(xué)院，長春130012；3.吉林大學(xué) 農(nóng)學(xué)部，長春130062）

0 引 言

考慮如下具周期外力的Van der Pol方程［1］：

方程（1）描述了如圖1所示的非線性電子線路，其中：C表示電容；μ表示與電感線圈有關(guān)的量；Vs表示輸入電壓；V表示輸出電壓；ψ（i）表示電子線路的伏安特征曲線，這里選取

在實(shí)際電路中，C和μ通常都非常小，因此方程（1）實(shí)際是帶有兩個(gè)小參數(shù)的奇異攝動(dòng)問題.為討論方便，作時(shí)間尺度變換s＝Ct，則方程（1）變?yōu)?/p>

其中ε＝μ／C，一般假設(shè)ε充分小.

Van der Pol方程是一種經(jīng)典的自激振蕩電路系統(tǒng)，不僅在各種復(fù)雜的動(dòng)力系統(tǒng)建模中廣泛使用，而且在物理學(xué)、生物學(xué)、電子學(xué)和經(jīng)濟(jì)學(xué)等領(lǐng)域應(yīng)用廣泛.文獻(xiàn)［3］給出了具周期外力Van der Pol方程極限環(huán)的存在性和唯一性；文獻(xiàn)［4］研究了在周期外力作用下Van der Pol振蕩系統(tǒng)所產(chǎn)生的同步現(xiàn)象和混沌現(xiàn)象；文獻(xiàn)［5］給出了方程（2）脈沖型解的定義，并探討了脈沖型解的存在性和穩(wěn)定性；文獻(xiàn)［6］利用匹配漸近展開法給出了方程（2）脈沖型解的一致有效漸近逼近，并討論了脈沖型解中脈沖的個(gè)數(shù).

重正化群方法［7－10］是處理微分方程中不規(guī)則現(xiàn)象的主要工具.本文利用建立在包絡(luò)理論上的重正化群方法構(gòu)造方程（2）脈沖型解一致有效的漸近展開式，并給出脈沖型解形成一個(gè)脈沖所用的時(shí)間.

圖1 非線性電子線路Fig.1 Nonlinear electronic circuit

1 主要結(jié)果

為簡單，令

則方程（2）變?yōu)?/p>

如圖2所示，分4段構(gòu)造方程（3）脈沖型解的一致有效逼近.設(shè)解從A點(diǎn)出發(fā)，即初值條件為

其中

1）y＞0，此時(shí)ψ（y）＝K1y.方程（3）變?yōu)?/p>

做時(shí)間伸縮變換τ＝t／ε，則方程（5）變?yōu)?/p>

圖2 構(gòu)造過程Fig.2 Construction process

考慮方程（6）滿足如下初始條件的解：

設(shè)

將式（8）代入式（6），對比ε的同次冪系數(shù)，并分別解出xi（τ），yi（τ）（i＝0，1，2），再代回式（8）得

其中：A，B為任意常數(shù)；τ0為初始時(shí)間.

根據(jù)文獻(xiàn)［7－8］的結(jié)果知，方程（6）的重正化群方程為

即

即

從而

由式（7）得

設(shè)

將式（10）代入式（9）并對比ε同次冪系數(shù)，得

從而

假設(shè)解從A到B所用的時(shí)間為t1，則

設(shè)

（2）規(guī)范化數(shù)據(jù)格式。包括對不規(guī)范數(shù)據(jù)項(xiàng)進(jìn)行字段格式約束定義或數(shù)據(jù)驗(yàn)證，對特殊的數(shù)據(jù)項(xiàng)進(jìn)行合并或拆分處理。如influenza_aa.dat中的第8列是核酸序列的名稱，如 “InfluenzaAvirus (A/swine/Illinois/A01240575/2013(H3N2))”，由于數(shù)據(jù)存儲(chǔ)和檢索分析的需要，需要將該數(shù)據(jù)項(xiàng)拆分為四個(gè)字段來存儲(chǔ)：類型(A)、序列名稱(A/swine/Illinois/A01240575/2013)、H 血清型(3)和 N 血清型(2)。

則由式（11），（12）知，p1，q1和r1滿足如下方程：

于是

其中：

考慮方程（18）從B點(diǎn)出發(fā)的解，即滿足初值條件

的解.設(shè)

假設(shè)解從A到C所用的時(shí)間為t2，則

設(shè)

則由式（21），（22）知，p2，q2滿足如下方程：

于是

其中：

3）y＜0，此時(shí)ψ（y）＝K1y＋（K2－K1）i0.令a＝（K1－K2）i0，用類似于1）的方法可得

假設(shè)解從A到D所用的時(shí)間為t3，則

設(shè)

則由式（26），（27）知，p3，q3和r3滿足如下方程：

于是

其中：

4）y＜0，此時(shí)ψ（y）＝K2y.計(jì)算方法同1），解得

假設(shè)解從A到E所用的時(shí)間為t4，則

設(shè)t4＝t3＋p4εlnε＋q4ε＋…，則由式（32），（33）知，p4，q4滿足如下方程：

于是xⅣ（t4）＝α4＋β4εlnε＋γ4ε＋…，其中

若該解是周期解，點(diǎn)E必須與點(diǎn)A重合，從而有β0＝β4，γ0＝γ4.

綜上，可得：

定理1 假設(shè)對周期輸入函數(shù)f（t），系統(tǒng)（2）至少存在具有一個(gè)脈沖的脈沖型解，則該解形成一個(gè)脈沖所用的時(shí)間為t0＝p＋qεlnε＋rε＋…，其中：

p1，p2，p3，p4，q1，q2，q3，q4，r1，r3分別由式（14），（23），（28），（34），（15），（24），（29），（35），（16），（30）給出.

［1］Van der Pol B.Forced Osillations in a Circuit with Non－linear Resistance［J］.Phil Magaz Ser 7，1927，3（13）：65－80.

［2］唐志凱，劉隆和，姜永華，等.Van der Pol振蕩器的模擬仿真 ［J］.現(xiàn)代雷達(dá)，2007，29（5）：89－92.（TANG Zhikai，LIU Longhe，JIANG Yonghua，et al.Simulation of Van der Pol Oscillators［J］.Modern Radar，2007，29（5）：89－92.）

［3］Levinson N，Smith O K.A General Equation for Relaxation Oscillations［J］.Duke Math J，1942，9（2）：382－403.

［4］Elabbasy E M，El－Dessoky M M.Synchronization of Van der Pol Oscillator and Chen Chaotic Dynamical System［J］.Chaos，Solitons and Fractals，2008，36（5）：1425－1435.

［5］Chow S N，HUANG Wenzhang.On the Number of Spikes of Solution for a Forced Singularly Perturbed Differential Equation［J］.Ann Mat Pura Appl，2006，185（5）：129－169.

［6］Chow S N，LIN Ping，SHI Shaoyun.Spike Solutions of a Nonlinear Electric Circuit with a Periodic Input［J］.Taiwanese Journal of Mathematics，2005，9（4）：551－581.

［7］Kunihiro T.The Renormalization－Group Method Applied to Asymptotic Analysis of Vector Fields ［J］.Prog Theor Phys，1997，97（2）：179－200.

［8］Kunihiro T.A Geometrical Formulation of the Renormalization Group Method for Global Analysis ［J］.Prog Theor Phys，1995，94（4）：503－514.

［9］Kirkinis E.The Renormalization Group：A Perturbation Method for the Graduate Curriculum ［J］.SIAM Rev，2012，54（2）：374－388.

［10］Delamotte B.An Introduction to the Non－perturbative Renormalization Group ［M／OL］.2007－02－15.http：／／arxiv.org／abs／cond－mat／0702365vl.