魏美華，常金勇，馬崛

1.榆林學(xué)院數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，陜西榆林 719000

2.中國(guó)科學(xué)院信息工程研究所，北京 100093

一類(lèi)活化基質(zhì)模型非常數(shù)正平衡解的全局結(jié)構(gòu)

魏美華1，常金勇2，馬崛1

1.榆林學(xué)院數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，陜西榆林 719000

2.中國(guó)科學(xué)院信息工程研究所，北京 100093

1 引言

對(duì)于該活化基質(zhì)模型的研究工作主要集中在從常微分方程角度討論極限環(huán)的存在性和唯一性[8-10]。關(guān)于已有常微分方程結(jié)果，文獻(xiàn)[8-9]運(yùn)用數(shù)值和實(shí)驗(yàn)方法說(shuō)明了極限環(huán)的存在性，我國(guó)陳蘭蓀先生[10]運(yùn)用常微分方程理論數(shù)學(xué)上嚴(yán)格地證明了極限環(huán)的存在性和唯一性?？紤]到反應(yīng)物分子擴(kuò)散的情況，這涉及到偏微分方程。而關(guān)于已有的偏微分方程結(jié)果，文獻(xiàn)[11]討論了常數(shù)平衡解的存在性和穩(wěn)定性，而運(yùn)用拓?fù)涠壤碚摵徒獾南闰?yàn)估計(jì)研究了非常數(shù)平衡解的存在性和不存在性[12]。

在此基礎(chǔ)上，本文在一維空間下運(yùn)用分歧理論和度理論研究該模型非常數(shù)平衡解的局部和全局結(jié)構(gòu)，運(yùn)用數(shù)值模擬給予解釋和驗(yàn)證，并闡明該模型具有活化基質(zhì)模型所具有的本質(zhì)特征。所得理論分析結(jié)果推廣了文獻(xiàn)[13]。

2 預(yù)備知識(shí)

引理2表明若d1比較小或d2比較大時(shí)，方程（1）沒(méi)有非常數(shù)正解。

假設(shè)條件：

成立。它是該糖酵解模型為活化基質(zhì)模型以及常數(shù)平衡解無(wú)擴(kuò)散穩(wěn)定的充分條件。

3 局部分歧

本章在一維空間Ω=(0，l)，l＞0下，以d1作為分歧參數(shù)運(yùn)用分歧理論分析式（1）的非常數(shù)平衡解的局部結(jié)構(gòu)。

在一維情況下，系統(tǒng)式（1）的平衡解相應(yīng)于系統(tǒng)：的特征值λj=(πj/l)2，j=0，1，…，相應(yīng)的正規(guī)化的特征函數(shù)為：

則{φj(x):j=0，1，…}構(gòu)成L2(0，l)的一組標(biāo)準(zhǔn)正交基。設(shè)Hilbert空間Y=L2(0，l)×L2(0，l)的內(nèi)積為：

考慮L的伴隨算子：

4 全局分歧

圖1 模型式（1）的非常數(shù)平衡解

設(shè)Uˉ是系統(tǒng)式（6）的解，則可利用反射和周期延拓構(gòu)造系統(tǒng)式（3）的解：

易知(，(0，0))是系統(tǒng)式（6）的分歧點(diǎn)。同理可證系統(tǒng)式（6）的分歧曲線要么延伸到無(wú)窮，要么連接((0，0))，k′＞k。若第二種情況發(fā)生，則與k的定義矛盾，所以?！延伸到無(wú)窮。綜上可知，Γj延伸到無(wú)窮。而由引理1知u是有界的，因此Γj隨延向無(wú)窮。

注2定理2說(shuō)明Γj不會(huì)最終連接到某個(gè)分歧點(diǎn)，但是有可能Γj連接某個(gè)分歧點(diǎn)并最終延向無(wú)窮。

注3圖1刻畫(huà)了非常數(shù)平衡解的存在性，證實(shí)了定理1和定理2。圖2體現(xiàn)了活化基質(zhì)模型的本質(zhì)特征，即活化劑越多則基質(zhì)越少，這是因?yàn)榛罨瘎┰蕉鄷?huì)導(dǎo)致基質(zhì)的消耗越多。

5 結(jié)束語(yǔ)

圖2 基質(zhì)u和活化劑v的對(duì)比

基于已有的工作，參見(jiàn)文獻(xiàn)[12]，為了進(jìn)一步討論帶Neumann邊界條件的活化基質(zhì)模型（即糖酵解模型）正平衡解的存在性，本文在一維空間上以擴(kuò)散系數(shù)d1為分歧參數(shù)，運(yùn)用局部和全局分歧理論以及度理論研究了該模型常數(shù)平衡解產(chǎn)生的局部分歧解結(jié)構(gòu)和全局分歧解走向，并利用數(shù)值模擬得以驗(yàn)證，所得結(jié)果證實(shí)了圖靈模式的形成。

[1]Turing A M.The chemical basis of morphogenesis[J].Phil Trans R Soc London Ser B，1952，237：37-72.

[2]Murray J D.Mathematical biology II：spatial models and biomedical applications[M].3rd ed.[S.l.]：Springer，2003.

[3]葉其孝，李正元.反應(yīng)擴(kuò)散方程引論[M].北京：科學(xué)出版社，1994.

[4]Ni W M，Kanako S，Izumi T.The dynamics of a kinetic activator-inhibitor system[J].J Differential Equations，2006，229：426-465.

[5]Wei J，Winter M.Existence，classification and stability analysis of multiple-peaked solutions for the Gierer-Meinhardt system inR1[J].Methods Appl Anal，2007，14：119-163.

[6]Hale J K，Peletier L A，Troy W C.Exact homoclinic and heteroclinic solutions of the Gray-Scott model for autocatalysis[J].SIAM J Appl Math，2000，61：102-130.

[7]Li Y.Steady-state solution for a general Schnakenberg model[J].Nonlinear Anal RWA，2011，12（4）：1985-1990.

[8]Tyson J，Kauffman S.Control of mitosis by a continuous biochemical oscillation[J].J Math Biol，1975，1：289-310.

[9]Segel L A.Mathematical models in molecular and cellular biology[M].Cambridge：Cambridge University Press，1980.

[10]陳蘭蓀，王東達(dá).一個(gè)生物化學(xué)反應(yīng)的振動(dòng)現(xiàn)象[J].數(shù)學(xué)物理學(xué)報(bào)，1985（3）：261-266.

[11]Ashkenazi M，Othmer H G.Spatial patterns in coupled biochemical oscillators[J].J Math Biol，1978，5：305-350.

[12]魏美華，吳建華.一類(lèi)糖酵解模型正平衡解的存在性分析[J].數(shù)學(xué)學(xué)報(bào)，2011（4）：553-560.

[13]Jang J，Ni W M，Tang M.Global bifurcation and structure of turing patterns in 1-D Lengyel-Epstein model[J]. J Dynam Differential Equations，2005，16：297-320.

[14]Crandall M，Rabinowitz P.Bifurcation from simple eigenvalues[J].J Functional Anal，1971，8：321-340.

[15]Wei M H，Wu J H，Guo G H.Turing structures and stability for the 1-D Lengyel-Epstein system[J].J Math Chem，2012，50：2374-2396.

[16]Rabinowitz P H.Some global results for nonlinear eigenvalue problems[J].J Functional Anal，1971，7：487-513.

WEI Meihua1,CHANG Jinyong2,MA Jue1

1.School of Mathematics and Statistics,Yulin University,Yulin,Shaanxi 719000,China
2.Institute of Information Engineering,Chinese Academy of Sciences,Beijing 100093,China

An activator-substrate system under Neumann boundary condition is considered in one-dimensional space.Taking the diffusion coefficientd1as bifurcation parameter,the local and global bifurcation of constant steady-state solution are studied by bifurcation theory and degree theory.Moreover,the theoretical results are confirmed by numerical simulations, and also continue the previous work.It is shown that the nonconstant steady-state branches join up with infinity.

activator-substrate model;bifurcation;steady-state solutions;global structure

在一維空間上，研究一類(lèi)帶Neumann邊界條件的活化基質(zhì)模型。以擴(kuò)散系數(shù)d1為分歧參數(shù)，運(yùn)用分歧理論和度理論研究該模型常數(shù)平衡解的局部分歧和全局分歧，利用數(shù)值模擬得以證實(shí)。理論結(jié)果表明非常數(shù)平衡解分支延伸向無(wú)窮，這是已有工作的繼續(xù)。

活化基質(zhì)模型；分歧；平衡解；全局結(jié)構(gòu)

A

O175.26

10.3778/j.issn.1002-8331.1402-0295

WEI Meihua,CHANG Jinyong,MA Jue.Global structure of nonconstant steady-state solutions for activator-substrate system.Computer Engineering and Applications,2014,50（18）：50-53.

國(guó)家自然科學(xué)基金（No.11271236）；陜西省教育廳科研計(jì)劃資助項(xiàng)目（No.14JK1862）；榆林學(xué)院博士科研啟動(dòng)基金項(xiàng)目（No.13GK04）。

魏美華（1981—），女，博士，講師，主要研究方向?yàn)槠⒎址匠碳捌淇梢暬怀＝鹩拢?982—），男，博士生，講師；馬崛（1981—），女，講師。

2014-02-26

2014-05-30

1002-8331（2014）18-0050-04

CNKI網(wǎng)絡(luò)優(yōu)先出版：2014-07-25，http://www.cnki.net/kcms/detail/10.3778/j.issn.1002-8331.1402-0295.html