林府標(biāo)， 張千宏

（貴州財(cái)經(jīng)大學(xué) 數(shù)統(tǒng)學(xué)院，貴州貴陽(yáng)550025）

群體平衡模型廣泛用于微粒過(guò)程及系統(tǒng)等眾多領(lǐng)域［1-5］.如生物化學(xué)、農(nóng)業(yè)科學(xué)、醫(yī)藥科學(xué)等領(lǐng)域中微生物、細(xì)菌、細(xì)胞的出生、生長(zhǎng)、死亡、傳輸、破損等過(guò)程.細(xì)胞的分裂導(dǎo)致新的子細(xì)胞出生和死亡的微生物群體平衡方程（積分-偏微分方程）［1，3-5］可寫(xiě)成

其中，x表示細(xì)胞質(zhì)量，t表示時(shí)間.進(jìn)一步，f（x，t）表示微生物種群密度，假設(shè)密度函數(shù)f（x，t）足夠光滑且任意階偏導(dǎo)數(shù)均連續(xù)，G（x，t）表示細(xì)胞的生長(zhǎng)率，v（x，t）表示t時(shí)刻細(xì)胞分裂的平均數(shù)量，Γ（x，t）表示細(xì)胞的出生率，p（x，y）表示原質(zhì)量是x的細(xì)胞分裂成質(zhì)量是y的子細(xì)胞的概率，概率函數(shù)滿(mǎn)足正則性條件：

考慮細(xì)胞的生長(zhǎng)率、分裂平均數(shù)量及出生率皆為常數(shù)，對(duì)研究微生物種群細(xì)胞分裂過(guò)程是有價(jià)值和意義的.若選取動(dòng)力學(xué)參變量函數(shù)

其中g(shù)、ν、κ均為正常數(shù)，則在約束限制條件（2）下群體平衡方程（1）可簡(jiǎn)化成

尋找群體平衡方程的精確解及解析解技術(shù)比較困難，實(shí)體工程科學(xué)應(yīng)用領(lǐng)域常采用數(shù)值實(shí)驗(yàn)方案［1-4］進(jìn)行研究，如對(duì)一般簡(jiǎn)單積分-偏微分方程（1）對(duì)應(yīng)的實(shí)體模型，探尋精確解都是比較棘手的問(wèn)題.雖然經(jīng)典李群分析方法［6-10］是計(jì)算常微分方程和純偏微分方程的對(duì)稱(chēng)群的有效方法之一，但不能用于計(jì)算積分-偏微分方程（1）和（3）的對(duì)稱(chēng)群.近年來(lái)，改進(jìn)的李群分析方法［11-12］已被用于研究群體平衡方程的對(duì)稱(chēng)群、約化-積分常微分方程及精確解［13-16］.運(yùn)用改進(jìn)了的李群分析方法探究積分-偏微分方程（3）的對(duì)稱(chēng)群，障礙就是積分-偏微分方程（3）對(duì)應(yīng)的決定方程的求解問(wèn)題.因決定方程仍是積分-偏微分方程及積分類(lèi)型的下限是變量，這些棘手問(wèn)題阻礙了通解的探尋和方法的研究與創(chuàng)新.

本文計(jì)算積分-偏微分方程（3）的細(xì)胞質(zhì)量分布的矩、對(duì)稱(chēng)群、群不變解、約化積分-常微分方程、顯式精確解及分析解的動(dòng)力學(xué)行為性質(zhì)及特征.

1 細(xì)胞質(zhì)量分布的矩

源于統(tǒng)計(jì)學(xué)中矩的概念和方法常用于研究群體平衡模型中的粒子、微?；蚣?xì)胞尺寸、質(zhì)量分布的均值及方差［1-4］.若f＝f（x，t）是方程（3）的任一解，則對(duì)質(zhì)量足夠大的細(xì)胞種群密度分布函數(shù)值必然為零［1-4］，特別地種群密度f(wàn)（x，t）更應(yīng)滿(mǎn)足性質(zhì)細(xì)胞種群密度分布的j階矩定義為

其中零階矩M0（t）表示任意時(shí)刻單位體積內(nèi)細(xì)胞總的平均數(shù)量，一階矩M1（t）表示任意時(shí)刻細(xì)胞分裂總的質(zhì)量，j階矩Mj（t）（j≥2）代表的實(shí)體含義可參見(jiàn)文獻(xiàn)［1-4］.運(yùn)用性矩的定義、分部積分法及交換二重積分的積分次序可得矩恒等式

因此，采用這些矩恒等關(guān)系式，一方面對(duì)方程（3）兩邊關(guān)于變量x從0到∞同時(shí)積分.另外，先在方程（3）兩邊皆乘以x，然后關(guān)于變量x從0到∞同時(shí)積分，得到完全封閉的矩方程組滿(mǎn)足的柯西問(wèn)題

若任意t時(shí)刻，微生物種群分裂系統(tǒng)中沒(méi)有細(xì)胞，即可假設(shè)種群細(xì)胞密度f(wàn)（0，t）＝0.在此條件下，常微分方程組（4）給出的零階矩M0（t）和一階矩M1（t）分別為

若f（0，t）≠0，可類(lèi)似地求解矩方程組（4），為了行文簡(jiǎn)潔省略求解過(guò)程.若種群細(xì)胞在分裂或破損過(guò)程中總保持質(zhì)量守恒，則恒有即細(xì)胞的總質(zhì)量密度函數(shù)M1（t）是時(shí)間t的不變量.

2 方程（3）的對(duì)稱(chēng)群

鑒于經(jīng)典李群分析方法［6-10］直接不能用于計(jì)算積分-偏微分方程（3）的對(duì)稱(chēng)群.而改進(jìn)了的李群分析方法［11-12］計(jì)算積分-偏微分方程（3）的對(duì)稱(chēng)群，最大困難是寫(xiě)出積分-偏微分方程（3）的決定方程且求解決定方程.因決定方程仍是積分-偏微分方程，通解的探尋是很棘手的問(wèn)題，如何求解取決于原積分-偏微分方程（3）自身的結(jié)構(gòu)特征［13-16］.基于這些障礙和羈絆對(duì)積分-偏微分方程（3）兩邊同時(shí)關(guān)于變量x求導(dǎo)得

因方程（6）是純偏微分方程，經(jīng)典李群分析方法可直接運(yùn)用.于是假設(shè)偏微分方程（6）接受的無(wú)窮小李對(duì)稱(chēng)算子為

其中X（2）是X的二階延拓算子，系數(shù)函數(shù)ξx、ηxx、ηxt的定義分別為

其中Dx和Dt分別是關(guān)于x和t的全微分算子.通過(guò)計(jì)算系數(shù)函數(shù)ξx、ηxt、ηxx的表達(dá)式為

依據(jù)經(jīng)典李群分析框架及算法，得偏微分方程（6）的決定方程為

其中左下角標(biāo)符號(hào)｜（6）表示決定方程（7）對(duì)偏微分方程（6）的任一解均恒成立.將系數(shù)函數(shù)ξx、ηxx、ηxt的表達(dá)式分別代入決定方程（7）得到

方程（6）可改寫(xiě)成

將f xt的表達(dá)式代入上述決定方程，把決定方程寫(xiě)成關(guān)于f tt、f xx、f t、f x的多項(xiàng)式.再令多項(xiàng)式的各項(xiàng)系數(shù)為零，得到系數(shù)函數(shù)ξ、τ、η滿(mǎn)足的超決定方程組為

運(yùn)用方程ηxf＝0和ηff＝0，函數(shù)η（x，t，f）可假設(shè)為η（x，t，f）＝h（t）f+p（x，t），其中h（t）、p（x，t）是任意待確定的未知函數(shù).將η（x，t，f）的表達(dá)式代入上述超決定方程組的最后一個(gè)方程得到

采用方程ξf＝0，將函數(shù)ξ（x，t，f）按照泰勒級(jí)數(shù)展開(kāi)成，其中ai（t）（i＝0，1，…）是任意待決定的可微函數(shù).因此，將ξ（x，t，f）的表達(dá)式代入方程（8）得到

由此方程可推出

從而ai（t）＝0（i＝0，1，2，…），即ξ（x，t，f）＝0.最后采用τf＝0，τx＝0，將η（x，t，f）的表達(dá)式代入超決定方程組，可求得決定方程（7）的通解為

ξ＝0， τ＝c1， η＝c2f+σ（x，t），

其中c1、c2為任意常數(shù)，σ（x，t）是偏微分方程（6）的任一解.因此，偏微分方程（6）的對(duì)稱(chēng)群的全體生成元構(gòu)成一個(gè)無(wú)窮維李代數(shù)，它包含一個(gè)2維的子李代數(shù)L2＝span｛X1，X2｝，并且有一組基無(wú)窮維的子李代數(shù)L＝span｛X｝，其無(wú)

以及一個(gè)∞σ窮小李對(duì)稱(chēng)算子為

無(wú)窮小李對(duì)稱(chēng)算子X(jué)σ對(duì)應(yīng)的李群為

其中a為群參數(shù).積分-偏微分方程（3）不接受李群Ta.事實(shí)上，采用積分-偏微分方程（3），將李群Ta的表達(dá)式代入積分-偏微分方程

因σ＝σ（x，t）是偏微分方程（6）的任一解，不是積分-偏微分方程（3）的解，對(duì)任意群參數(shù)a恒等式

表明李群Ta不是將積分-偏微分方程（3）的任一解映射變換成同一積分-偏微分方程的解.

無(wú)窮小李對(duì)稱(chēng)算子X(jué)1和X2對(duì)應(yīng)的李群分別為

其中τ0和b分別是李群Tτ0和Tb的群參數(shù).積分-偏微分方程（6）接受李群Tb.事實(shí)上，采用積分-偏微分方程（3）將李群Tb的表達(dá)式代入積分-偏微分方程

對(duì)任意群參數(shù)b恒等式

表明李群Tb將積分-偏微分方程（3）的任一解映射變換成同一積分-偏微分方程的解.同理可證積分-偏微分方程（3）接受李群Tτ0.因此，基于無(wú)窮小李對(duì)稱(chēng)算子（9）被積分-偏微分方程（3）接受的事實(shí)呈現(xiàn)下列定理1和2.

定理1偏微分方程（6）接受的無(wú)窮小李對(duì)稱(chēng)算子，積分-偏微分方程（3）不全部接受.

定理2積分-偏微分方程（3）的對(duì)稱(chēng)群的部分生成元構(gòu)成一個(gè)2維的子李代數(shù)L2＝span｛X1，X2｝，并且有一組基（9）.

針對(duì)定理1激發(fā)的思考，這里給出一個(gè)具有啟發(fā)性的常微分方程例證性的例子.設(shè)y＝y(tǒng)（x），對(duì)方程y′＝x-1兩邊同時(shí)關(guān)于x求導(dǎo)得y″＝1.于是易驗(yàn)證方程y″＝1接受平移算子但方程y′＝x-1卻不接受.

3 方程（3）的群不變解、精確解及解的動(dòng)力學(xué)行為

3.1 方程（3）的群不變解利用子李代數(shù)L2的一維最優(yōu)化子李代數(shù)span｛X2｝，span｛X1+αX2｝，α∈R，計(jì)算積分-偏微分方程（3）的群不變解、約化積分-常微分方程及顯式精確解.

情形span｛X1+αX2｝.李對(duì)稱(chēng)算子X(jué)1+αX2的群不變量為J1＝x，J2＝exp（-αt）f，于是方程（3）的群不變解的表達(dá)式可假設(shè)為

函數(shù)φ（x）滿(mǎn)足約化的積分-常微分方程

一維子李代數(shù)span｛X2｝對(duì)應(yīng)的群不變量沒(méi)有找到，故相應(yīng)方程（3）的群不變解和約化方程皆沒(méi)有找到.

3.2 方程（3）的顯式精確解利用平移變換群Tτ0的平移作用，方程（3）的精確解可寫(xiě)成采用觀察試湊函數(shù)方法［13-16］研究約化積分-常微分方程（10），可找到許多精確解，再結(jié)合相應(yīng)群不變解的表達(dá)式可獲得方程（3）的精確解.事實(shí)上，假設(shè)積分-常微分方程（10）的精確解的表達(dá)式可寫(xiě)成

其中常數(shù)γ＞0，di（i＝1，…，n）為待決定的未知常數(shù).將表達(dá)式（11）代入積分-常微分方程（10），通過(guò)計(jì)算整理成關(guān)于xj（j＝0，1，…）的多項(xiàng)式，分別令xj（j＝0，1，…）的各項(xiàng)系數(shù)為零，得到關(guān)于未知參數(shù)di（i＝1，…，n）、α、ν、κ、γ、g的方程組.利用吳消元法［17］結(jié)合數(shù)學(xué)軟件Reduce［18］計(jì)算解得di（i＝1，…，n）、α、ν、κ、γ、g.采用積分-常微分方程（10）對(duì)應(yīng)的群不變解的表達(dá)式以及平移變換李群Tτ0的平移作用，得到積分-偏微分方程（3）的精確解為

當(dāng)取n＝1，…，7時(shí)，積分-偏微分方程（3）的精確解（12）的計(jì)算結(jié)果列于表1.當(dāng)n取其余值的情況可根據(jù)微生物種群細(xì)胞分裂實(shí)體模型的需要類(lèi)似計(jì)算.表1中解的動(dòng)力學(xué)行為分析僅考慮n＝1的情形，其它的解可類(lèi)似分析.

表1 積分-偏微分方程（3）的顯式精確解Tab.1 Explicit exact solutions of the integro-partial differential equation（3）

3.3 方程（3）解的動(dòng)力學(xué)行為選取n＝1時(shí)，精確解（12）的表達(dá)式為精確解（13）滿(mǎn)足動(dòng)力學(xué)的相關(guān)性質(zhì).首先當(dāng)x→∞時(shí)，f（x，t）→0，這表明對(duì)質(zhì)量足夠大的細(xì)胞，細(xì)胞分裂或破損進(jìn)化分布函數(shù)值必然為零.另外f（x，t）→0，當(dāng)ν＜1且t→∞時(shí)，這表明精確解是漸進(jìn)穩(wěn)定的.而f（x，t）→∞，當(dāng)ν＞1且t→∞時(shí)，這表明精確解是不穩(wěn)定的.若取ν＝1，精確解（13）是平凡的.精確解（13）對(duì)應(yīng)的邊值條件和柯西問(wèn)題的初值條件分別為

精確解（13）在矩形區(qū)域［0，L］×［0，T］上對(duì)應(yīng)的邊值條件和柯西問(wèn)題的初值條件分別為

通過(guò)計(jì)算得精確解（13）對(duì)應(yīng)的零階矩M0（t）和一階矩M1（t）分別為

3.4 動(dòng)力學(xué)參數(shù)與種群密度分布實(shí)際微生物化學(xué)工程科學(xué)系統(tǒng)中，任意時(shí)刻細(xì)胞分裂的平均數(shù)量v（x，t）的最小值為2，如微生物種群一個(gè)細(xì)胞分裂成兩個(gè)，此時(shí)恒有v（x，t）＝2.又因v（x，t）用的是平均值，故v（x，t）的值可不要求一定為正整數(shù).在多重細(xì)胞分裂過(guò)程中，v（x，t）的取值常依賴(lài)實(shí)驗(yàn)數(shù)據(jù)而決定.增長(zhǎng)率G（x，t）代表質(zhì)量是x的細(xì)胞的變化率表示每單位時(shí)間內(nèi)細(xì)胞的質(zhì)量及沿著x軸方向細(xì)胞的傳輸速度.對(duì)任意給定的t＝t0時(shí)刻，將種群密度分布函數(shù)（13）看作是關(guān)于細(xì)胞的質(zhì)量x、分裂平均數(shù)量ν、增長(zhǎng)率g及出生率κ的四元函數(shù)

可分別選取坐標(biāo)平面（x，t），（x，ν），（x，g），（x，κ）研究函數(shù)值H的變化特征與性質(zhì)，為生物化學(xué)中細(xì)胞分裂過(guò)程選取合適的真實(shí)動(dòng)力學(xué)參數(shù)值ν、κ、g提供參考.

選取坐標(biāo)平面（x，t）和動(dòng)力學(xué)參數(shù)ν＝2，κ＝0.02，g＝1，τ0＝0和ν＝4，κ＝0.1，g＝1.5，τ0＝0，細(xì)胞分裂進(jìn)化動(dòng)力學(xué)行為分布函數(shù)（14）在矩形區(qū)域［0，150］×［0，4］和［0，20］×［0，4］上的空間圖像分別見(jiàn)圖1的（a）和（b）.選取坐標(biāo)平面（x，ν）及動(dòng)力學(xué)參數(shù)τ0＝0，t＝1，κ＝0.2，g＝1.2，細(xì)胞分裂進(jìn)化動(dòng)力學(xué)行為分布函數(shù)（14）在矩形區(qū)域［0，8］×［2，8］上的空間圖像見(jiàn)圖2的（c）.選取坐標(biāo)平面（x，g）及動(dòng)力學(xué)參數(shù)τ0＝0，t＝1，κ＝1.2，ν＝2，細(xì)胞分裂進(jìn)化動(dòng)力學(xué)行為分布函數(shù)（14）在矩形區(qū)域［0，8］×［0.1，8］上的空間圖像見(jiàn)圖2的（d）.分布函數(shù)（14）在其余坐標(biāo)平面上的情形可類(lèi)似討論，為了行文簡(jiǎn)潔，省略了相關(guān)動(dòng)力學(xué)行為分布函數(shù)在空間的圖像性質(zhì)及特征分析.

＝圖1 圖（a）和（b）分別是ν2，4時(shí)細(xì)胞分裂進(jìn)化行為分布函數(shù)（13）的空間圖像Fig.1 The figures（a）and（b）are spatial images of evolution behaviour of cells fission distribution function（13）withν＝2，4 respectively

＝圖2 圖（c）和（d）分別是平面（x，ν），（x，g）上t1時(shí)細(xì)胞分裂進(jìn)化行為函數(shù)（14）的空間圖像Fig.2 The figures（c）and（d）are spatial images of evolution behaviour of cells fission distribution function（14）with t＝1 on（x，ν），（x，g）plane respectively

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

經(jīng)典李群分析方法［6-10］不能直接計(jì)算積分-偏微分方程（3）的對(duì)稱(chēng)群.本文獲得了積分-偏微分方程（3）的對(duì)稱(chēng)群、約化積分-常微分方程、群不變解及顯式精確解.分析了部分精確解的動(dòng)力學(xué)行為性質(zhì)和特征.如何直接利用改進(jìn)了的李群分析方法［11-12］計(jì)算積分-偏微分方程（3）的對(duì)稱(chēng)群，值得在今后的工作研究中嘗試和探索.

致謝貴州省教育廳創(chuàng)新群體項(xiàng)目（黔教合KY字［2021］015）、貴州省教育廳青年科技人才成長(zhǎng)項(xiàng)目（黔教合KY字［2017］150）、2018年度貴州財(cái)經(jīng)大學(xué)校級(jí)科研基金（2018XYB04）和貴州財(cái)經(jīng)大學(xué)創(chuàng)新探索及學(xué)術(shù)新苗項(xiàng)目（黔科合平臺(tái)人才［2017］5736-020）對(duì)本文給予了資助，謹(jǐn)致謝意.