基于有限元方法求解帶有年齡結構的種群模型

郝永樂,左 平,朱 青

(1.周口師范學院 數學與統計學院,河南 周口 466001;2.空軍航空大學 基礎部,長春 130022)

0 引 言

在生物種群中,年齡不同將導致種群內部例如生育率、死亡率、捕食能力等方面存在明顯差異.因此,將種群按照年齡結構進行劃分的模型能更好地反映個體的生理特征以及種群擴散程度.目前,對帶有年齡結構種群模型的研究已取得了豐富的成果[1-6].偏微分方程求解數值解是近年來的研究熱點,目前對具有年齡結構數學模型的研究主要集中在解的定性理論上,而對復雜偏微分生物模型的數值求解問題研究報道較少.本文考慮帶有隨機擴散項的種群模型[7]:

1 模型處理

為了對上述模型進行數值求解,先要對模型進行適當的處理和簡化,步驟如下:

2) 對方程(1)和方程(4)關于年齡a進行積分;

3) 針對反應擴散系數進行處理[8-10],得到如下簡化系統:

(8)

其中:

(9)

表示種群在t時刻的出生率;

(10)

表示種群在t時刻的死亡率;

(11)

表示種群在t時刻、位置x處的種群總量;

(12)

2 有限元方法

下面對所提出的偏微分模型進行數值解求解.首先給出方程(8)的變分形式:

(13)

將空間區域Ω剖分成有限個單元的集合Th,形成有限元網格.在一維空間中常用等分剖分,二維空間中常用三角剖分或四邊形剖分.通常情況下,有限元方法要求網格是正規的.其次,構造有限元子空間,其由Th上所有屬于C(Ω)的、多項式次數不大于p-1的分片多項式組成的子空間Sh形成.本文的數值實驗使用等分剖分作為有限區間的網格.