Post on 07-Jul-2020
個体群動態の数理
• 科目ナンバリングコード:2223011A3 • 開設科目名:個体群動態の数理 • 講義コード:4802000 • 開講期・曜日・時限・教室:前期 水曜日 1・2時限 情報科学講義室(G302) • 対象学生:3回生
1
奈良女子大学理学部 化学生物環境学科 環境科学コース
奈良女子大学理学部・化学生物環境学科 環境科学コース 高須夫悟
連続時間のモデル時間が離散ではなく、連続的に変化するモデル。微分方程式で記述される。
季節変化がない一定条件の環境で、世代が重なっている生物集団に関して個体密度変化が滑らかに変化すると見なせる場合に用いるモデル。
時刻 t における個体密度を n(t) とすると、一般に
という形で書ける
個体密度の時間変化率
時刻 t での n の時間変化率がn の関数 f(n) で与えられるモデル(常微分方程式)
密度に依存した正味の増殖項
一般に f(n) = birth(n) – death(n)
奈良女子大学理学部 化学生物環境学科 環境科学コース
dn
dt= f(n)
<latexit sha1_base64="yMzZuWT81YdCF6wWe++SNmZ5rtU=">AAACeXichVHLSsNAFD2N7/qqj4XgplqU6qJMVVAEpeDGpW2tLWgpSTqtwTQJybRgQ3/AH3DhRgUR9TPc+AMu+gniUsGFLrxNA6JFvcPMnDlzz50zM4qla45grBmQurp7evv6B4KDQ8Mjo6Gx8T3HrNoqz6imbto5RXa4rhk8IzSh85xlc7mi6DyrHG219rM1bjuaaeyKY4vnK3LZ0EqaKguiCqHJg5Itq27RaLhF0QhvhEtRY6EQirAY8yLcCeI+iMCPHTN0jQMUYUJFFRVwGBCEdchwqO0jDgaLuDxc4mxCmrfP0UCQtFXK4pQhE3tEY5lW+z5r0LpV0/HUKp2iU7dJGcYce2Q37IU9sDv2xD5+reV6NVpejmlW2lpuFUZPptJv/6oqNAscfqn+9CxQwprnVSPvlse0bqG29bX66Ut6PTXnzrNL9kz+L1iT3dMNjNqrepXkqTME6QPiP5+7E+wtxeLLsaXkSiSx6X9FP6Yxiyi99yoS2MYOMnRuHee4wW3gXZqRotJiO1UK+JoJfAtp+RMJwJGj</latexit>
指数増加・減少モデル
3
奈良女子大学理学部 化学生物環境学科 環境科学コース
1 個体当たりの増加率 c(n) = f(n) / n を用いると
c(n) が個体密度 n に無関係な定数 r の場合
初期条件 n(0) = n0 を用いると、n(t) = n0 exp[r t]
変数分離法で解ける!
r > 0 の時指数増加、r < 0 の時指数減少のモデル(連続時間)
dn
dt= f(n)
<latexit sha1_base64="yMzZuWT81YdCF6wWe++SNmZ5rtU=">AAACeXichVHLSsNAFD2N7/qqj4XgplqU6qJMVVAEpeDGpW2tLWgpSTqtwTQJybRgQ3/AH3DhRgUR9TPc+AMu+gniUsGFLrxNA6JFvcPMnDlzz50zM4qla45grBmQurp7evv6B4KDQ8Mjo6Gx8T3HrNoqz6imbto5RXa4rhk8IzSh85xlc7mi6DyrHG219rM1bjuaaeyKY4vnK3LZ0EqaKguiCqHJg5Itq27RaLhF0QhvhEtRY6EQirAY8yLcCeI+iMCPHTN0jQMUYUJFFRVwGBCEdchwqO0jDgaLuDxc4mxCmrfP0UCQtFXK4pQhE3tEY5lW+z5r0LpV0/HUKp2iU7dJGcYce2Q37IU9sDv2xD5+reV6NVpejmlW2lpuFUZPptJv/6oqNAscfqn+9CxQwprnVSPvlse0bqG29bX66Ut6PTXnzrNL9kz+L1iT3dMNjNqrepXkqTME6QPiP5+7E+wtxeLLsaXkSiSx6X9FP6Yxiyi99yoS2MYOMnRuHee4wW3gXZqRotJiO1UK+JoJfAtp+RMJwJGj</latexit>
dn
dt= c(n)n
<latexit sha1_base64="g3MTnDyaipDm+BM3Rci2zd+Yu9s=">AAACe3ichVHLSsNAFD2Nr1ofrQoiuAkWRUXK1AeKoAhuXKq1WlCRZDqtwTQJybRQY3/AH3DhSkFE9C/c+AMu/ARxqeBGwds0ICrqDZk5c+49d87M6I5peJKxh4jS1NzS2hZtj3V0dnXHEz29m55ddrnIctu03ZyuecI0LJGVhjRFznGFVtJNsaUfLNfzWxXheoZtbciqI3ZLWtEyCgbXJFF7if6dgqtxP2/V/LysqQsqH7XGVEokWYoFof4E6RAkEcaqnbjEDvKwwVFGCQIWJGETGjz6tpEGg0PcLnziXEJGkBeoIUbaMlUJqtCIPaCxSKvtkLVoXe/pBWpOu5j0u6RUMczu2RV7Znfsmj2yt197+UGPupcqzXpDK5y9+PFA5vVfVYlmif1P1Z+eJQqYC7wa5N0JmPopeENfOTx5zsyvD/sj7Jw9kf8z9sBu6QRW5YVfrIn1U8ToAdLfr/sn2JxMpadSk2vTyaXF8CmiGMQQRum+Z7GEFawiS/se4RzXuIm8K0llXJlolCqRUNOHL6HMfAB5lJJC</latexit>
dn
dt= rn
<latexit sha1_base64="sdZ8+89uznrLUYlfBMeEJwn8t8Q=">AAACeHichVHLSsNAFD2Nr1pf1W4EN8FSH5syrYIiKAU3Lq1aW9AiSTqtwTQJybRQQ3/AH3DhQhRExc9w4w+46CeISwVBXHibBkSLeoeZOXPmnjtnZlTb0F3BWCsk9fT29Q+EByNDwyOjY9HxiV3Xqjkaz2mWYTkFVXG5oZs8J3Rh8ILtcKWqGjyvHq239/N17ri6Ze6Ihs2LVaVi6mVdUwRRB9HYftlRNK9kNr2SaMqrsiMTG2dJ5ofcDVIBiCOITSt6jX2UYEFDDVVwmBCEDShwqe0hBQabuCI84hxCur/P0USEtDXK4pShEHtEY4VWewFr0rpd0/XVGp1iUHdIKSPBHtkte2EP7I49sY9fa3l+jbaXBs1qR8vtg7GTye23f1VVmgUOv1R/ehYoY9n3qpN322fat9A6+vrx6cv2ylbCm2GX7Jn8X7AWu6cbmPVX7SrLt84QoQ9I/XzubrCbTqYWkunsYjyzFnxFGFOYxhy99xIy2MAmcnRuA+e4wW3oXZKlWWm+kyqFAk0M30JKfwKTE5F0</latexit>
dn
n= rdt
<latexit sha1_base64="tFOpAKCzcaseWLf4ry4+QjtZid0=">AAACeHichVHLSsNAFD2N7/qKdiO4CRZfmzKpgiIoghuX1loVVEqSTjWYJiGZFmroD/gDLlyIgmjxM9z4Ay76CeKygiAuvE0DoqLeYWbOnLnnzpkZ3bVMXzDWiEkdnV3dPb198f6BwaFheWR023fKnsFzhmM53q6u+dwybZ4TprD4rutxraRbfEc/Xmvt71S455uOvSWqLj8oaYe2WTQNTRCVlxP7RU8zgoJdC+yasqx4SkHk5SRLsTCUn0CNQBJRbDjyDfZRgAMDZZTAYUMQtqDBp7YHFQwucQcIiPMImeE+Rw1x0pYpi1OGRuwxjYe02otYm9atmn6oNugUi7pHSgWT7JHVWZM9sDv2xN5/rRWENVpeqjTrbS1388OnY9nXf1UlmgWOPlV/ehYoYjH0apJ3N2RatzDa+srJWTO7tDkZTLEr9kz+L1mD3dMN7MqLcZ3hm+eI0weo35/7J9hOp9S5VDozn1xdib6iF+OYwAy99wJWsY4N5OjcKi5wi3rsTVKkaWm2nSrFIk0CX0JKfwCSJ5F0</latexit>
Zdn
n=
Zrdt
<latexit sha1_base64="JCXRA89y38Q4DzGalfUMtRvBE5E=">AAACgnichVHLSgMxFD2O7/qquhFEGCyKuCgZFRRREdy49FUVrJSZaVqD08wwkxZ06MqdP+DClYKIuNNPcOMPuPATxKWCGxfeTgdERb0hycnJPTcnieU5IlCMPTZojU3NLa1t7YmOzq7unmRv32bgln2bZ2zXcf1tywy4IyTPKKEcvu353CxZDt+y9pdq+1sV7gfClRvqwOO7JbMoRUHYpiIqlxzKCqn0bME37TAvq6Gs6vN6xPl6XuWSKZZmUeg/gRGDFOJYcZOXyCIPFzbKKIFDQhF2YCKgtgMDDB5xuwiJ8wmJaJ+jigRpy5TFKcMkdp/GIq12YlbSulYziNQ2neJQ90mpY4Q9sCv2wu7ZNXti77/WCqMaNS8HNFt1LfdyPccD62//qko0K+x9qv70rFDATORVkHcvYmq3sOv6yuHJy/rs2kg4ys7ZM/k/Y4/sjm4gK6/2xSpfO0WCPsD4/tw/weZE2phMT6xOpRYX4q9owyCGMUbvPY1FLGMFGTr3CFe4wa3WpI1rhjZZT9UaYk0/voQ29wFWgZRn</latexit>
log |n| = rt+ C<latexit sha1_base64="mrArtTiX5zq2Gp1Lx1+0vxQBibA=">AAACeXichVHLSsNAFD2Nr1pf9bEQ3ESLIgpl0gqKoBTcuPRVFbSUJE5rME1CMi201R/wB1y4UUFE/Qw3/oALP0FcKrjQhTdpQLSod5iZM2fuuXNmRnNMwxOMPUaklta29o5oZ6yru6e3L94/sOnZZVfnWd02bXdbUz1uGhbPCkOYfNtxuVrSTL6lHSz5+1sV7nqGbW2IqsNzJbVoGQVDVwVR+fjQrmkX5UPZor4gu7KQp+WlfDzBkiwIuRkoIUggjBU7foVd7MGGjjJK4LAgCJtQ4VHbgQIGh7gc6sS5hIxgn+MIMdKWKYtThkrsAY1FWu2ErEVrv6YXqHU6xaTuklLGOHtg1+yF3bNb9sQ+fq1VD2r4Xqo0aw0td/J9x8Prb/+qSjQL7H+p/vQsUMBc4NUg707A+LfQG/pK7eRlfX5tvD7BLtgz+T9nj+yObmBVXvXLVb52ihh9gPLzuZvBZiqppJOp1ZlEZjH8iihGMIZJeu9ZZLCMFWTp3BrOcI2byLs0Kk1KU41UKRJqBvEtpPQnu7+QBA==</latexit>
n(t) = exp[rt+ C]<latexit sha1_base64="l+eFdofMXleBvsFAy8zKNmP2Q94=">AAACeXichVHLSsNAFD2Nr1pf9bEQ3ESLUhXKtBUUQRHcuKyPtkItJYmjBtMkJNNiW/wBf8CFGxVEqp/hxh9w0U8QlwoudOFtGhAV9Q4zc+bMPXfOzKi2obuCsUZAamvv6OwKdod6evv6B8KDQxnXKjkaT2uWYTnbquJyQzd5WujC4Nu2w5WiavCserja3M+WuePqlrklKjbPF5V9U9/TNUUQVQiPmFExLS/JO/zIzsmOkGfl1XwhHGEx5oX8E8R9EIEfKSt8jR3swoKGEorgMCEIG1DgUsshDgabuDxqxDmEdG+f4xgh0pYoi1OGQuwhjfu0yvmsSetmTddTa3SKQd0hpYxJ9sDq7Jnds1v2yN5/rVXzajS9VGhWW1puFwZORjdf/1UVaRY4+FT96VlgDwueV5282x7TvIXW0perp8+bixuTtSl2yZ7I/wVrsDu6gVl+0a7W+cYZQvQB8e/P/RNkErF4MpZYn4usLPtfEcQYJhCl957HCtaQQprOreIcddwE3qRxKSrNtFKlgK8ZxpeQkh8XmJA0</latexit>
指数増加・減少グラフ
4
奈良女子大学理学部 化学生物環境学科 環境科学コース
t t
n(t) n(t)
r = 1.0r = 0.7
r = 0.5
n(t) = n0 exp[r t], n0 = 10
r = –0.5r = –1.0 r = –0.7
1 2 3 4 5t
10
20
30
40
50
n[t]
1 2 3 4 5t
0.2
0.4
0.6
0.8
1.0
n[t]
対数スケール
5
奈良女子大学理学部 化学生物環境学科 環境科学コース
t t
log n(t) log n(t)
r = 1.0
r = 0.7
r = 0.5
r = –0.5
r = –1.0
r = –0.7
n(t) = n0 exp[r t] の両辺の対数をとると
log n(t) = log n0 + r t 傾きが r の直線
1 2 3 4 5t
5
10
50
100
n[t]
1 2 3 4 5t
0.01
0.05
0.10
0.50
1
n[t]
倍加時間と半減時間
6
奈良女子大学理学部 化学生物環境学科 環境科学コース
離散時間の指数増加・減少モデルと同様に、個体密度が 2 倍に増えるのに要する倍加時間と、半減するのに要する半減時間が求められる
倍加時間 Doubling time: Td
n0 exp[r (t + Td)] = 2 n0 exp[r t] より Td = log2 / r = 0.6931 / r
n0 exp[r (t + Th)] = 1/2 n0 exp[r t] より Th = – log2 / r = – 0.6931 / r
半減時間 Half life: Th
r = 0.1 / day の時、Td = 6.93 days
r = –0.2 / week の時、Th = 3.47 weeks
連続時間のロジスティックモデル
7
奈良女子大学理学部 化学生物環境学科 環境科学コース
個体密度 n が高くなると密度効果が働き、1 個体あたりの増加率 c が減少する場合を考える
c(n) が n に比例して減少する場合を考えるn = K の時 c はゼロ。
個体密度 n(t)
1 個体あたりの増加率 c(n)
r
K0
連続時間のロジスティックモデル
r : 内的自然増加率K : 環境収容量
ロジスティックモデルの解
8
奈良女子大学理学部 化学生物環境学科 環境科学コース
非線形微分方程式は一般に解けるとは限らないが、ロジスティックモデルは解ける。
変数を分離
初期条件 n(0) = n0 より
左辺を部分分数の和に分解
ロジスティックモデルのグラフ 1
9
奈良女子大学理学部 化学生物環境学科 環境科学コース
K = 10, r = 1
n0 = 0.1, 1, 5, 10, 12
t
n(t)
十分時間が経つと n(t) は K に収束
K
K は環境収容量
2 4 6 8 10t
2
4
6
8
10
12
n[t]
ロジスティックモデルのグラフ 2
10
奈良女子大学理学部 化学生物環境学科 環境科学コース
K = 10, n0 = 0.1
r = 1, 2, 4
r = 1
r = 2
r = 4
内的自然増加率 r が高いほうが急速に成長
2 4 6 8 10t
2
4
6
8
10
n[t]
t
n(t)
ロジスティックモデルのグラフ 2
11
奈良女子大学理学部 化学生物環境学科 環境科学コース
K = 10, n0 = 0.1
r = 1, 2, 4
r = 1
r = 2
r = 4
内的自然増加率 r が高いほうが急速に成長
t
n(t)
2 4 6 8 10t
0.5
1
5
10
n[t]
縦軸は対数スケール
個体密度 n が小さいとき、n はほぼ指数的に増加(対数スケールで直線)
具体例
12
奈良女子大学理学部 化学生物環境学科 環境科学コース
ゾウリムシ Paramecium caudatum
イースト菌 Yeast cells
r = 0.54 / hour
Cause 1934Bulmer 1994
Carlson 1913Brown and Rothery 1993
微分方程式の Graphical Analysis
13
奈良女子大学理学部 化学生物環境学科 環境科学コース
一般に微分方程式は解析的に解けるとは限らない。
の振る舞いをグラフを用いて調べる方法がある
下準備:横軸に n、縦軸に dn/dt = f(n) のグラフを描く
連続時間のロジスティックモデルの場合、f(n) = r(1–n/K)n
n
f(n)
0
K
dn
dt= f(n)
<latexit sha1_base64="yMzZuWT81YdCF6wWe++SNmZ5rtU=">AAACeXichVHLSsNAFD2N7/qqj4XgplqU6qJMVVAEpeDGpW2tLWgpSTqtwTQJybRgQ3/AH3DhRgUR9TPc+AMu+gniUsGFLrxNA6JFvcPMnDlzz50zM4qla45grBmQurp7evv6B4KDQ8Mjo6Gx8T3HrNoqz6imbto5RXa4rhk8IzSh85xlc7mi6DyrHG219rM1bjuaaeyKY4vnK3LZ0EqaKguiCqHJg5Itq27RaLhF0QhvhEtRY6EQirAY8yLcCeI+iMCPHTN0jQMUYUJFFRVwGBCEdchwqO0jDgaLuDxc4mxCmrfP0UCQtFXK4pQhE3tEY5lW+z5r0LpV0/HUKp2iU7dJGcYce2Q37IU9sDv2xD5+reV6NVpejmlW2lpuFUZPptJv/6oqNAscfqn+9CxQwprnVSPvlse0bqG29bX66Ut6PTXnzrNL9kz+L1iT3dMNjNqrepXkqTME6QPiP5+7E+wtxeLLsaXkSiSx6X9FP6Yxiyi99yoS2MYOMnRuHee4wW3gXZqRotJiO1UK+JoJfAtp+RMJwJGj</latexit>
n
f(n)
n
f(n)
Graphical Analysis
14
奈良女子大学理学部 化学生物環境学科 環境科学コース
連続時間のモデルにおいても、時間的に変化しない状態を平衡点と呼ぶ
平衡点 n* では時間微分がゼロとなる
関数 f(n) と横軸の交点が平衡点
ロジスティックモデルの場合 n* = 0, K の二つ
0 K
0 < n < K の時、dn/dt > 0つまり、n は時間とともに増加
K < n の時、dn/dt < 0つまり、n は時間とともに減少
平衡点 n* = 0 は不安定
平衡点 n* = K は安定
dn
dt= f(n⇤) = 0
<latexit sha1_base64="/vVvbb+yjKhKL+AjkH01di1WBNo=">AAACf3ichVG7SgNBFD1Z3/GRqI1oEww+i3Cjgg9QBBtLX1HBR9jdTHRxs7vsTgK6BKz9AQsrBRWx0H+w8Qcs/ASxVLCx8GazICrqHWbmzJl77pyZ0RzT8CTRY0Spqa2rb2hsija3tLbF4u0dq55ddHWR0W3Tdtc11ROmYYmMNKQp1h1XqAXNFGva3lxlf60kXM+wrRW574itgrpjGXlDVyVT2Xj3Zt5VdT9nlf2cLCemE/lBa3t4iAFl40lKURCJnyAdgiTCWLDjl9hEDjZ0FFGAgAXJ2IQKj9sG0iA4zG3BZ85lZAT7AmVEWVvkLMEZKrN7PO7waiNkLV5XanqBWudTTO4uKxPoowe6ohe6p2t6ovdfa/lBjYqXfZ61qlY42dhR1/Lbv6oCzxK7n6o/PUvkMRF4Ndi7EzCVW+hVfeng+GV5aqnP76czemb/p/RId3wDq/Sqny+KpRNE+QPS35/7J1gdSaVHUyOLY8nZmfArGtGDXgzye49jFvNYQIbPPcQFbnCrRJQBJaVQNVWJhJpOfAll8gNd95IV</latexit>
局所安定性解析(連続時間:微分方程式)
15
奈良女子大学理学部 化学生物環境学科 環境科学コース
の平衡点を n* とする f(n*) = 0
平衡点からの微小なずれを h(t) とする。n(t) = n* + h(t) である
上式に代入してずれ h(t) の式として書き直すと
関数 f(n*+h) を n* の周りでテイラー展開
f(n*) = 0 であり、また、微少量 h の 2 次以上の項を無視して次式を得る
定数項
線型近似
dn
dt= f(n)
<latexit sha1_base64="yMzZuWT81YdCF6wWe++SNmZ5rtU=">AAACeXichVHLSsNAFD2N7/qqj4XgplqU6qJMVVAEpeDGpW2tLWgpSTqtwTQJybRgQ3/AH3DhRgUR9TPc+AMu+gniUsGFLrxNA6JFvcPMnDlzz50zM4qla45grBmQurp7evv6B4KDQ8Mjo6Gx8T3HrNoqz6imbto5RXa4rhk8IzSh85xlc7mi6DyrHG219rM1bjuaaeyKY4vnK3LZ0EqaKguiCqHJg5Itq27RaLhF0QhvhEtRY6EQirAY8yLcCeI+iMCPHTN0jQMUYUJFFRVwGBCEdchwqO0jDgaLuDxc4mxCmrfP0UCQtFXK4pQhE3tEY5lW+z5r0LpV0/HUKp2iU7dJGcYce2Q37IU9sDv2xD5+reV6NVpejmlW2lpuFUZPptJv/6oqNAscfqn+9CxQwprnVSPvlse0bqG29bX66Ut6PTXnzrNL9kz+L1iT3dMNjNqrepXkqTME6QPiP5+7E+wtxeLLsaXkSiSx6X9FP6Yxiyi99yoS2MYOMnRuHee4wW3gXZqRotJiO1UK+JoJfAtp+RMJwJGj</latexit>
d
dt(n⇤ + h) =
dh
dt= f(n⇤ + h)
<latexit sha1_base64="WdQwFEGbLCLk5CFCt/H0crIb1pQ=">AAACl3ichVHLSsNAFD3Gd31V3ShugkXxAWVaBUXwAYK4bNWqYLUk6dQG0yQk04KG/oA/oOBKQUT8Bzcu7A+48BPEpYIbF96mQdGi3mFm7j33njtnZlTb0F3B2GOD1NjU3NLa1h7q6Ozq7gn39m26VtHReEqzDMvZVhWXG7rJU0IXBt+2Ha4UVINvqQfL1fxWiTuubpkb4tDmuwVl39RzuqYIgjLheDrnKJqXLXtZUR4z9ybkSTk/Ls/LAZ73ExTnPpOZcIRFmW9yvRMLnAgCS1jhK6SRhQUNRRTAYUKQb0CBS2MHMTDYhO3CI8whT/fzHGWEiFukKk4VCqEHtO5TtBOgJsXVnq7P1ugUg6ZDTBkj7IFdsxdWYTfsib3/2svze1S1HNKu1rjczvQcD6y//csq0C6Q/2L9qVkgh1lfq07abR+p3kKr8UtHJy/rc2sj3ii7YM+k/5w9sju6gVl61S6TfO0MIfqA2M/nrnc249HYVDSenI4sLQRf0YYhDGOM3nsGS1hFAik69xS3uEdFGpQWpRVptVYqNQScfnwzKfkB0RKa4A==</latexit>
dh
dt= f(n⇤) + f
0(n⇤)h+O(h2)<latexit sha1_base64="b7crU7jFnO6LhLIt57g5vL2QwHY=">AAACkXicSyrIySwuMTC4ycjEzMLKxs7BycXNw8vHLyAoFFacX1qUnBqanJ+TXxSRlFicmpOZlxpaklmSkxpRUJSamJuUkxqelO0Mkg8vSy0qzszPCympLEiNzU1Mz8tMy0xOLAEKxQvoxKQVJSZXp2TUVqeU1CrYKqRp5MVpaSpoK6SpQ1gZQLa/gkZGnJFmvICygZ4BGChgMgyhDGUGKAjIF1jOEMOQwpDPkMxQypDLkMqQx1ACZOcwJDIUA2E0gyGDAUMBUCyWoRooVgRkZYLlUxlqGbiAekuBqlKBKhKBotlAMh3Ii4aK5gH5IDOLwbqTgbbkAHERUKcCg6rBVYOVBp8NThisNnhp8AenWdVgM0BuqQTSSRC9qQXx/F0Swd8J6soF0iUMGQhdeN1cwpDGYAF2aybQ7QVgEZAvkiH6y6qmfw62ClKtVjNYZPAa6P6FBjcNDgN9kFf2JXlpYGrQbAYuYAQYogc3JiPMSM/QWM8o0ETZwQ4aFRwM0gxKDBrA8DZncGDwYAhgCAXaO4lhB8NBhkNMwkyWTA5MThClTIxQPcIMKIDJGwAwAJcA</latexit>
dh
dt⇡ f 0(n⇤)h
<latexit sha1_base64="3xQJGOsWoAhI/LHm5zlWDzVZeYk=">AAAChHichVG7SgNBFD2u7/iK2ggWLob4QsLEB4pFCNhY+koiGA27m4lZ3Owuu5OgLilt/AELKwURsdUvsPEHLPwEsYxgY+HNZkFU1DvMzJkz99w5M6Pahu4Kxp6apOaW1rb2js5QV3dPb1+4fyDtWmVH4ynNMixnS1VcbugmTwldGHzLdrhSUg2eUfeX6/uZCndc3TI3xaHNd0rKnqkXdE0RROXCI9mCo2hevlj18qIqZxXbdqwDuTA+Ye5OTcrFXDjCYswP+SeIByCCIFat8BWyyMOChjJK4DAhCBtQ4FLbRhwMNnE78IhzCOn+PkcVIdKWKYtThkLsPo17tNoOWJPW9Zqur9boFIO6Q0oZUfbIrlmNPbAb9szef63l+TXqXg5pVhtabuf6ToY23v5VlWgWKH6q/vQsUMCi71Un77bP1G+hNfSVo9PaxtJ61BtjF+yF/J+zJ3ZPNzArr9rlGl8/Q4g+IP79uX+C9EwsPhubWZuLJBPBV3RgGKOYoPdeQBIrWEWKzj3GDW5xJ7VJ09KsNN9IlZoCzSC+hJT4AMFalPw=</latexit>
dh
dt= rh r = f 0(n⇤)
<latexit sha1_base64="1kKHU6ouMLbw81TBNbO/lsfavZk=">AAACinichVG7SgNBFD2u7/iK2gg2iyG+ijCJgi8ighaWvqKC0bC7mZjFze66OwnEJT+gH2BhpSAiVrZa2vgDFn6CWCrYWHizWRAV9Q4zc+bMPXfOzKi2obuCscc6qb6hsam5pTXU1t7R2RXu7ll3raKj8ZRmGZazqSouN3STp4QuDL5pO1wpqAbfUPfmq/sbJe64umWuibLNtwvKrqnndE0RRGXC0XTOUTQvm694WVGRk7Ij5+X0flHJEkrKuaFhc2d0JBOOsBjzQ/4J4gGIIIglK3yBNLKwoKGIAjhMCMIGFLjUthAHg03cNjziHEK6v89RQYi0RcrilKEQu0fjLq22AtakdbWm66s1OsWg7pBSRpQ9sEv2wu7ZFXti77/W8vwaVS9lmtWaltuZrsO+1bd/VQWaBfKfqj89C+Qw6XvVybvtM9VbaDV96eD4ZXV6JeoNsjP2TP5P2SO7oxuYpVftfJmvnCBEHxD//tw/wXoiFh+LJZbHI3OzwVe0oB8DGKb3nsAcFrGEFJ17hGvc4FZqlxLSlDRTS5XqAk0vvoS08AG7z5Yn</latexit>
n
グラフで見る局所安定性
16
奈良女子大学理学部 化学生物環境学科 環境科学コース
平衡点 n* からの微小なずれ h(t) の時間変化は指数的に変化
f '(n*) > 0 の時、ずれは指数的に増加:平衡点 n* は不安定
f '(n*) < 0 の時、ずれは指数的に減少:平衡点 n* は局所的に安定
f(n) = r(1–n/K)n f'(n) = r(1–2n/K)
f'(0) = r > 0
f'(K) = –r < 0
n* = 0 は不安定
n* = K は安定
h(t) / exp[rt]<latexit sha1_base64="4L3vAveyXoBJR+tSu7EPCuzD+B4=">AAACe3ichVHLSgMxFD0dX7W+qoIIbgaLoiIlrYriQgQ3Llu1KrSlzIypDk5nwkxa1OoP+AMuXFkQkfoXbvwBF36CuFRwo+DtdEBU1BuSnJzcc3OS6MIyPcnYQ0hpaW1r7wh3Rrq6e3r7ov0Dm55Tdg2eMRzLcbd1zeOWafOMNKXFt4XLtZJu8S19f6Wxv1Xhrmc69oY8FDxf0nZts2gamiSqEB3am5CTak64jpCOmuMHIuuqMl+Ixlic+aH+BIkAxBBEyoleIYcdODBQRgkcNiRhCxo8alkkwCCIy6NKnEvI9Pc5ThAhbZmyOGVoxO7TuEurbMDatG7U9Hy1QadY1F1Sqhhj9+yaPbM7VmeP7O3XWlW/RsPLIc16U8tFoe90eP31X1WJZom9T9WfniWKWPC9muRd+EzjFkZTXzk6e15fXBurjrMaeyL/F+yB3dIN7MqLcZnma+eI0Ackvj/3T7CZjCdm4sn0bGx5KfiKMEYwigl673ksYxUpZOjcY9RQx03oXYkpU8p0M1UJBZpBfAll7gOY2ZJX</latexit>
r = f 0(n⇤)<latexit sha1_base64="P2/YiVn93fI4Dwt41lhDwQGRDDI=">AAACbnichVHLSsNAFD2Nr1ofrQoiiFgs1eqiTKugCErBjcs+rAq1ShKnNTRNQpIWavEH3IsLQVEQET/DjT/gop8gboQKblx4mwZERb1hMmfO3HPnzB3JUBXLZqzhETo6u7p7vL2+vv6BQX9gaHjT0iumzLOyrurmtiRaXFU0nrUVW+XbhsnFsqTyLam01trfqnLTUnRtw64ZPF8Wi5pSUGTRJipnBleChZmItjs3uxcIsShzIvgTxFwQghtJPXCDHexDh4wKyuDQYBNWIcKiL4cYGAzi8qgTZxJSnH2OI/hIW6EsThkisSX6F2mVc1mN1q2alqOW6RSVhknKIMLskd2yJntgd+yJvf9aq+7UaHmp0Sy1tdzY8x+PZd7+VZVptnHwqfrTs40ClhyvCnk3HKZ1C7mtrx6eNjPL6XB9ml2xZ/J/yRrsnm6gVV/l6xRPn8FHDxD73u6fYDMejc1H46mFUGLVfQovxjGFCPV7EQmsI4ms07ETnOPC8yKMChPCZDtV8LiaEXwJIfIBETyMpw==</latexit>
n
f(n)
ロジスティックモデルの場合 n* = 0, K の二つ
0
K
微分方程式の振る舞いのまとめ
17
奈良女子大学理学部 化学生物環境学科 環境科学コース
微分方程式一般 平衡点 n* での傾き f ‘(n*)で 局所安定性が決まる
dn
dt= f(n)
<latexit sha1_base64="yMzZuWT81YdCF6wWe++SNmZ5rtU=">AAACeXichVHLSsNAFD2N7/qqj4XgplqU6qJMVVAEpeDGpW2tLWgpSTqtwTQJybRgQ3/AH3DhRgUR9TPc+AMu+gniUsGFLrxNA6JFvcPMnDlzz50zM4qla45grBmQurp7evv6B4KDQ8Mjo6Gx8T3HrNoqz6imbto5RXa4rhk8IzSh85xlc7mi6DyrHG219rM1bjuaaeyKY4vnK3LZ0EqaKguiCqHJg5Itq27RaLhF0QhvhEtRY6EQirAY8yLcCeI+iMCPHTN0jQMUYUJFFRVwGBCEdchwqO0jDgaLuDxc4mxCmrfP0UCQtFXK4pQhE3tEY5lW+z5r0LpV0/HUKp2iU7dJGcYce2Q37IU9sDv2xD5+reV6NVpejmlW2lpuFUZPptJv/6oqNAscfqn+9CxQwprnVSPvlse0bqG29bX66Ut6PTXnzrNL9kz+L1iT3dMNjNqrepXkqTME6QPiP5+7E+wtxeLLsaXkSiSx6X9FP6Yxiyi99yoS2MYOMnRuHee4wW3gXZqRotJiO1UK+JoJfAtp+RMJwJGj</latexit>
f(n)
n
局所的に安定な平衡点
不安定な平衡点
n
dn/dt
局所安定と大域安定
18
奈良女子大学理学部 化学生物環境学科 環境科学コース
初期条件に関わらず収束する平衡点を大域安定な平衡点と呼ぶ
平衡点から少しずらしても元に戻る平衡点は局所安定
大域的に安定な平衡点は局所的にも安定。しかし逆は必ずしも成立しない。
連続時間のロジスティックモデルの K (環境収容量)は大域的に安定な平衡点
初期値に依存して収束先は必ずしも1 つとは限らない
n
f(n)
0
K
n
dn/dt
モデルの拡張
19
奈良女子大学理学部 化学生物環境学科 環境科学コース
ロジスティック成長をしている集団へ外部からの移入・移出がある場合
個体密度:n(t)内的自然増加率:r
環境収容量:K 外部からの移入が個体密度 n に比例する場合
移入 a n
a > 0 :移入a < 0 :移出
内的自然増加率が r + a環境収容量が K(r+a)/r のロジスティック増加と同じ
a :移入率
漁獲のモデル
20
奈良女子大学理学部 化学生物環境学科 環境科学コース
ロジスティック成長をする生物資源を考える
常に一定量の割合で収穫を行う場合、資源量 n の時間変化は
h :収穫率(定数)
h = 0
0 < h < r K2/4 の時、局所的に安定な平衡点 n*が存在
n* は次式の解
生物資源の枯渇を招かずに将来にわたって持続可能な収穫を行いたい(でもたくさん獲りたい)。どうしたらよいか?
n
dn/dt
最大収穫量
h = r K2/4
21
奈良女子大学理学部 化学生物環境学科 環境科学コース
収穫率 h の最大値は r K2/4
しかし、このとき資源量 n が平衡点 K / 2 よりも少しでも減少すると資源量は一気に減少してしまう。
一定の収穫率 h で生物資源を限界まで収奪する方法は賢明ではない
K/2
n
dn/dt
h = 0
持続可能な収穫
22
奈良女子大学理学部 化学生物環境学科 環境科学コース
ロジスティック成長をする生物資源を考える
現存する資源量に比例した収穫を行う場合、資源量 n の時間変化は
a :収穫努力係数(定数)収穫量
これは、内的自然増加率 r–a、環境収容量 K(r–a)/r のロジスティック成長と同じ
K(r–a)/r
資源量は n* = K(r–a)/r に収束(大域安定)
この時の収穫量 H = a n* = K(r–a)a/r
H を最大化する努力係数 a は r/2。この時、収穫量は r K/4 で最大かつ、持続可能n
f(n)
dn
dt= r
⇣1� n
K
⌘n� an = (r � a)
✓1� n
K(r � a)/r
◆n
<latexit sha1_base64="X6GRS7GhjgYIww1rzieUbLWcucY=">AAACy3ichVFNSxtBGH7cak3TD9N6KfQyNFjiwTirBUWwCL0IRfCjUcFImN3MxsHN7jI7CcTtHnvpH/DgqUIppT+jlx68evAnlN5qoYd66JvNYrHS9l125p3neZ93nplxIl/FhvOzIevG8MjN0cKt4u07d++Nle4/2IzDjnZlzQ39UG87Ipa+CmTNKOPL7UhL0XZ8ueXsP+/zW12pYxUGL00vkrtt0QqUp1xhCGqUvLqnhZs0gzRpmpQtMs3qvvRMxWZTbMAR9SJlda1ae2aSBYQLGhdZRU+JybyaXS3PqGmdXooapTKv8izY9cTOkzLyWA1L71FHEyFcdNCGRABDuQ+BmL4d2OCICNtFQpimTGW8RIoiaTtUJalCELpPY4tWOzka0LrfM87ULu3i069JyTDBT/kHfs4/84/8C7/4a68k69H30qPZGWhl1Bh783Djx39VbZoN9n6r/unZwMN85lWR9yhD+qdwB/ruweH5xsL6RPKEH/Ov5P8tP+Of6ARB97v7bk2uH6FID2D/ed3Xk82Zqj1bnVl7Wl56lj9FAY/wGBW67zksYRmrqNG+J/iGn7iwVqzYOrBeDUqtoVwzjithvf4FzYSt/w==</latexit>
資源管理と数理モデル
23
奈良女子大学理学部 化学生物環境学科 環境科学コース
持続可能な最大収穫量はどれだけか?
資源量変動に関する数理モデルの組立
モデル解析を通じた最適資源利用法の考察
資源管理への提言
イワシ、サバ、マグロ等の資源量の変動捕鯨問題、等
検索エンジンを用いて、資源管理、数理モデルのキーワードで検索
時間遅れの効果
24
奈良女子大学理学部 化学生物環境学科 環境科学コース
これまで扱ってきた微分方程式は、個体密度の時間変化率(微分)が、その時刻の個体密度に依存して決まるモデルであった。
時刻 t での変化率 時刻 t での個体密度
実際の生物では、例えば、誕生後一定の期間を経ないと繁殖しないという時間遅れ Time delay の効果が働く
また、今現在たくさん餌を食べたとしても、獲得したエネルギーが子孫の繁殖に変換されるのには時間がかかる。
時刻 T だけ過去に遡った個体密度時刻 t での変化率
時間遅れを伴う微分方程式
時間遅れを伴ったロジスティックモデル
25
奈良女子大学理学部 化学生物環境学科 環境科学コース
密度効果は時刻 T だけ過去にさかのぼった時点の密度に依存(T :時間遅れ)
時間遅れを導入すると、n(t) が周期的に変動する解がありうる 時間遅れが無ければ(T = 0 )周期解はあり得ない
t1
K
t1+T t2 t2+T
n(t1) = K であれば、時刻 t1+T で個体数は最大となる(微分がゼロ)
周期解の予感
n(t2) = K であれば、時刻 t2+T で個体数は最小となる(微分がゼロ)
t
n
r T < π / 2 の時、n(t) は K に収束
r T > π / 2 の時、K は不安定。振動する周期解、カオス
時間遅れと振動解
26
奈良女子大学理学部 化学生物環境学科 環境科学コース
平衡解 n(t) = K の周りで線型近似をすることにより
• r T = 0.3(K に収束)• r T = 1.2(K に収束)• r T = 1.8(振動解)
K = 1
具体例
27
奈良女子大学理学部 化学生物環境学科 環境科学コース
レミングの個体数変動
黒の実線:個体密度の観察値赤の点線:時間遅れロジスティックモデル
r = 3.33 / year, T = 0.72 year
レミングは個体数が増えると互いに攻撃的になり、雌は子供をあまり産まなくなる。個体数が減っても攻撃性がしばらく残るため、時間遅れが生じる、と考えられている。
From Case 2000
具体例2
28
奈良女子大学理学部 化学生物環境学科 環境科学コース
r T = 2.1赤線:時間遅れのロジスティックモデル
羊の害虫 Blowfly の個体数変動データ
時間遅れ T は、ハエの幼虫が成虫になるまでの時間
From Case 2000
血液濃度の減少は、数日の時間遅れを伴って骨髄での血液生成を促進させる。
図 (a) 時間遅れモデル T = 6 日
図 (b) 時間遅れモデル T = 20 日
図 (c) 白血病患者の白血球濃度
c(t): 血液濃度λ : 血液生産関数g : 血液消費率T: 時間遅れ
血液生成のモデル
29
奈良女子大学理学部 化学生物環境学科 環境科学コース
From Murray 1993
dc
dt= �(c(t� T ))� gc(t)
<latexit sha1_base64="LrPHYfj9x7c5RdEBFPAEa6Ykldw=">AAACj3icSyrIySwuMTC4ycjEzMLKxs7BycXNw8vHLyAoFFacX1qUnBqanJ+TXxSRlFicmpOZlxpaklmSkxpRUJSamJuUkxqelO0Mkg8vSy0qzszPCympLEiNzU1Mz8tMy0xOLAEKxQtoxaQVJSZXpyTXVqeU1CrYKsTkADWnJGooJGuU6IZoKmgq6Cqkgzia8QLKBnoGYKCAyTCEMpQZoCAgX2A5QwxDCkM+QzJDKUMuQypDHkMJkJ3DkMhQDITRDIYMBgwFQLFYhmqgWBGQlQmWT2WoZeAC6i0FqkoFqkgEimYDyXQgLxoqmgfkg8wsButOBtqSA8RFQJ0KDKoGVw1WGnw2OGGw2uClwR+cZlWDzQC5pRJIJ0H0phbE83dJBH8nqCsXSJcwZCB04XVzCUMagwXYrZlAtxeARUC+SIboL6ua/jnYKki1Ws1gkcFroPsXGtw0OAz0QV7Zl+SlgalBsxm4gBFgiB7cmIwwIz1DYz2jQBNlBztoVHAwSDMoMWgAw9ucwYHBgyGAIRRo7wSGbQz7GQ4wCTKZM9kxOUCUMjFC9QgzoAAmTwDLz5dj</latexit>
1 種系連続時間ダイナミクスのまとめ
30
奈良女子大学理学部 化学生物環境学科 環境科学コース
時間遅れが無いモデルでは、周期解はあり得ない
発散するか、安定平衡点に収束するかのどちらか
発散するか、安定平衡点に収束するか、もしくは周期振動解、カオス
時間遅れがあるモデルでは、周期振動解がありうる
離散時間のモデル Nt+1 = f(Nt) は、1 時刻ステップ前の状態に依存して個体密度が
決まるから本質的に時間遅れを含む離散時間ロジスティックモデルで周期解、カオスが生じることは既出
微分方程式モデル
差分方程式モデル
d
dtn(t) = f(n(t))
<latexit sha1_base64="0wg5WZz8FeJjdEV1DGz4duJiCYY=">AAACgXichVG7SgNBFD2u7/iK2ghaLAYlacJNFBRFEWwsfUUFlbC7mejiZnfZnQR0SWPpD1hYKYiKlf6CjT9g4SeIpYKNhTebBVFR7zAzZ87cc+fMjO5api+JHhuUxqbmlta29lhHZ1d3T7y3b813yp4hcoZjOd6GrvnCMm2Rk6a0xIbrCa2kW2Jd35uv7a9XhOebjr0q912xXdJ2bLNoGppkKh8f2ip6mhEUqkFBVu2kTKkzajGphiiVjycoTWGoP0EmAglEsejEL7CFAhwYKKMEARuSsQUNPrdNZEBwmdtGwJzHyAz3BaqIsbbMWYIzNGb3eNzh1WbE2ryu1fRDtcGnWNw9VqoYoQe6ohe6p2t6ovdfawVhjZqXfZ71ula4+Z6jgZW3f1UlniV2P1V/epYoYjL0arJ3N2RqtzDq+srB8cvK1PJIMEpn9Mz+T+mR7vgGduXVOF8SyyeI8Qdkvj/3T7CWTWfG0tml8cTcbPQVbRjEMJL83hOYwwIWkeNzD3GJG9wqjUpKISVbT1UaIk0/voQy/QHlWJK+</latexit>
d
dtn(t) = f(n(t), n(t� T ))
<latexit sha1_base64="QcE3yelyW/bNQRRZqt3rZRr7gVc=">AAACinichVHLSsNAFD3GV62vqhvBTbBUWtAyrYIvFEEXLn1VBZWSpFMNTZOQTAsa+gP6AS5cKYiIK7e6dOMPuPATxKWCGxfepAFRUW/I3DNn7rlzZka1Dd0VjD02SI1NzS2tkbZoe0dnV3esp3fdtSqOxnOaZVjOpqq43NBNnhO6MPim7XClrBp8Qy3N++sbVe64umWuiX2b75SVXVMv6poiiMrHEttFR9G8Qs0riJqZFCl5Ri4m5QAN+2lkLSWn8rE4S7Mg5J8gE4I4wliyYhfYRgEWNFRQBocJQdiAApe+LWTAYBO3A484h5AerHPUECVthao4VSjElmjcpdlWyJo093u6gVqjXQz6HVLKSLAHdsle2D27Yk/s/ddeXtDD97JPWa1ruZ3vPuxffftXVaYssPep+tOzQBETgVedvNsB459Cq+urB8cvq1MrCW+InbFn8n/KHtkdncCsvmrny3zlBFF6gMz36/4J1rPpzGg6uzwWn5sNnyKCAQwiSfc9jjksYgk52vcI17jBrdQhZaVJabpeKjWEmj58CWnhA9O4lTg=</latexit>
問題 1
31
奈良女子大学理学部 化学生物環境学科 環境科学コース
離散時間の指数増加モデル Nt+1 = R Nt と、
連続時間の指数増加モデル とが等しく増加する場合、
R と r の関係式を求めよ
等しく増加するとは、数列 {N0, N1, N2, ..., Nt, ...}と
{n(0), n(1), n(2), ..., n(t), ...} とが等しいことを意味する。
r = 0.01 の時、R の値はいくらになるか?
dn
dt= rn
<latexit sha1_base64="sdZ8+89uznrLUYlfBMeEJwn8t8Q=">AAACeHichVHLSsNAFD2Nr1pf1W4EN8FSH5syrYIiKAU3Lq1aW9AiSTqtwTQJybRQQ3/AH3DhQhRExc9w4w+46CeISwVBXHibBkSLeoeZOXPmnjtnZlTb0F3BWCsk9fT29Q+EByNDwyOjY9HxiV3Xqjkaz2mWYTkFVXG5oZs8J3Rh8ILtcKWqGjyvHq239/N17ri6Ze6Ihs2LVaVi6mVdUwRRB9HYftlRNK9kNr2SaMqrsiMTG2dJ5ofcDVIBiCOITSt6jX2UYEFDDVVwmBCEDShwqe0hBQabuCI84hxCur/P0USEtDXK4pShEHtEY4VWewFr0rpd0/XVGp1iUHdIKSPBHtkte2EP7I49sY9fa3l+jbaXBs1qR8vtg7GTye23f1VVmgUOv1R/ehYoY9n3qpN322fat9A6+vrx6cv2ylbCm2GX7Jn8X7AWu6cbmPVX7SrLt84QoQ9I/XzubrCbTqYWkunsYjyzFnxFGFOYxhy99xIy2MAmcnRuA+e4wW3oXZKlWWm+kyqFAk0M30JKfwKTE5F0</latexit>
問題 2
32
奈良女子大学理学部 化学生物環境学科 環境科学コース
体長 L が次式にしたがって成長する魚を考える。
ただし、k, Lmax > 0
1)誕生時、魚の体長は十分小さく L(0) =0 と見なせる場合、体長 L(t) を解け
2)体長 L(t) の時間変化をグラフに描け
3)k = 0.1 / day の時、魚の体長が最大体長の半分(Lmax/2)に成長するのに要する時間を求めよ。
dL
dt= k(Lmax � L)
<latexit sha1_base64="ikqb3JID78IkD8dj0/36ggXN84Y=">AAAChHichVG7SgNBFD1Z3/EVtREsXAyKgoZJoiiCIthYpPAVFVTC7maiS/bF7iSoS0obf8DCSkFE0uoX2PgDFn6CWCrYWHizWRAV9Q4zc+bMPXfOzKiOoXuCsceI1NDY1NzS2hZt7+js6o719G54dsnVeFazDdvdUhWPG7rFs0IXBt9yXK6YqsE31eJibX+zzF1Pt611cejwXVPZs/SCrimCqFxscKfgKpqfz1T8vKjIc3JRHs3kfFM5qMgTcmYsF4uzBAtC/gmSIYgjjGU7doUd5GFDQwkmOCwIwgYUeNS2kQSDQ9wufOJcQnqwz1FBlLQlyuKUoRBbpHGPVtsha9G6VtML1BqdYlB3SSljmD2wa/bC7lmVPbH3X2v5QY2al0Oa1bqWO7nuk/61t39VJs0C+5+qPz0LFDATeNXJuxMwtVtodX356PRlbXZ12B9hF+yZ/J+zR3ZHN7DKr9rlCl89Q5Q+IPn9uX+CjVQimU6kVibjC/PhV7RiAEMYpfeexgKWsIwsnXuMKm5wKzVL41JamqqnSpFQ04cvIc1/AE8WlEk=</latexit>
問題 3
33
奈良女子大学理学部 化学生物環境学科 環境科学コース
下記の式に従う動態を考える
1)グラフを用いた解析により、上式はパラメータの値に依存して複数個の平衡点を持つことを示せ。
2)適当なパラメータと初期値を設定して、数値計算により N の振る舞いを確認せよ。