指数関数的増加・減少エクスプローラー

N(t) = N0 · e^(kt) で増加と減少をモデル化する

インタラクティブシミュレーションを読み込んでいます...

レッスン

理論 — 指数関数的増加・減少エクスプローラー

指数関数的変化とは、変化率が存在する量に比例するときに生じる現象である — 預金利息、原子核の崩壊、増加の余地がある個体群など。成長(増加)と減衰で別々の公式が存在するわけではなく、1つの公式と符号の違いがあるのみである。

各記号の意味

N₀
t = 0 における初期値である。曲線全体を拡大縮小するのみで形状自体は変化させないため、2つの速度関連の行(「倍加時間」と「半減期」)のいずれにも現れない。
k
速度定数であり、その符号がすべてを決定する。正であれば量は無限に増加し、負であればゼロに達することなくゼロに向かって減少する。正確にゼロであれば何も変化しない。
t
k で指定された単位に基づく時間である。この2つの単位は一致していなければならない — 年単位の速度定数に対して月単位の時間を指定することは、誤った計算結果を引き起こす最もよくある原因である。
N(t)
指定時間後の値である。N₀ = 1000k = 0.08t = 10 のとき、値は 2225.5409 と表示される — 単発の増加ではなく連続的に複利成長したため、2倍以上に増加している。

公式の導き方

  1. 「変化率は存在する量に比例する」を数式で表すと dN/dt = kN となる。これが唯一の前提であり、以下のすべてはこの式から導かれる。
  2. 変数分離を行って積分すると N(t) = N₀·e^(kt) が得られる。底 e は便宜上の選択ではなく — 段階的ではなく連続的に増加が適用されるときに定数が到達する形式である。
  3. 倍加時間を求めるには、N2N₀ に達する時間を考える。e^(kt) = 2 より、kt = ln 2 となり、t = ln 2 / k となる。初期値は相殺されるため、倍加時間は速度定数のみによって決まる特性である。
  4. 半減期は 2 の代わりに ½ を用いた同様の計算であり、t = ln 2 / |k| となる — 速度定数 −k で半減するのにかかる時間は、速度定数 +k で倍加する時間と全く同じであるため、同じ式になる。k = 0.08 のとき「倍加時間」の行は 8.6643 を示し、その4分の1の速度定数である k = 0.02 では 34.6574 と正確に4倍の長さを示す。

表示の読み方

2つの変化率の行は互いに排他的であり、ダッシュ「—」は重要な役割を果たしています。k = 0.08 では倍加時間が 8.6643 で半減期は と表示され、k = −0.2 に切り替えると半減期が 3.4657 となって倍加時間のほうが空欄になります。増加する量には半減期がなく、減少する量には倍加時間がないため、この空欄の行は単なる抜け漏れではなく、一つの主張なのです。

前提
全期間にわたって速度定数が一定であり、量に対する制約が存在しないことを前提としている。現実にはいずれの前提も長くは維持されない。個体群は環境収容力に直面し、金利は再交渉されるからである。
成り立たない場合
k を正確に 0 に設定すると、両方の速度関連の行(「倍加時間」と「半減期」)に同時に と表示される。これは誠実な回答であり、熟考に値する。計算が壊れたわけではない。速度定数がゼロの場合、量は永久に 1000 のままであり、倍加することも半減することもなく、ln 2 / 0 は無限大ではなく未定義となる。公式は水平な直線へと退化し、2つの空欄表示は、計算の失敗ではなくゼロ除算を正しく回避した結果である。

なぜ70をパーセンテージで割ると倍加時間が求められるのか 🖖

この簡便法が成り立つのは ln 2 = 0.693147 だからであり、率を小数ではなくパーセントで数えると、これは姿を変えた70となる。年率7%のとき、このツールは倍加時間を 9.9021 と報告するのに対し、70/7 では 10 となり — 誤差は約100分の1(およそ1%)にすぎない。8%では 70/8 = 8.75 に対して 8.6643 を示し、率が大きくなるほど乖離は広がる。なぜなら、この近似は k が小さい極限においてのみ厳密だからだ。70の代わりに72を引き合いに出す銀行家たちも失敗したわけではない。72は 2, 3, 4, 6, 8, 9 で割り切れるため、彼らは口頭でできる計算のためにわずかな精度を犠牲にしているのだ。

全プロセスの詳細解説

  1. 1000個の細胞から始まり1時間あたり8%で増殖する培養 5 ステップ

    ある培養細胞群が 1000 個からスタートし、1 時間あたり 8% の割合で 10 時間成長する。3 人に聞けば 3 通りの異なる答えが返り、その全員が自身の答えを正当化できる。指数モデルが導き出す答えはどれか、そしてそれらの答えがどれほど離れているかを求めよ。

    1. 第 1 の答えは、1 時間あたり 8% を固定量として扱う。すなわち 1 時間に 80 個、10 時間で 800 個である。これは給与における「時給」と同じ意味であり、多くの人が思い浮かべる計算である。しかし、これは前提条件自体に照らしても誤っている唯一の答えである。なぜなら、第 1 時間目に生成された細胞は第 2 時間目にも分裂を続けているからである。

    2. 第 2 の答えは、1 時間ごとに 1 回複利計算を行うことでそれを修正する。すなわち 1.08 を 10 回掛けるのである。これで新しい細胞も成長に寄与する。これは十分に正当化できるモデルであり、個体群がちょうど毎時間の区切りで段階的に増殖する場合には完全に正しい。

    3. 第 3 の答えは極限をとる。細胞は時計を待たないため、成長率は連続的に適用される。これが微分方程式 dN/dt = kN である。その解は指数関数であり、第 2 の答えとの差のすべては e^0.8 に集約される。

    4. 倍加時間も同じ解から導かれ、その過程で N₀ は相殺される。これこそが指数成長が直感に反する理由である。2 倍になるまでの時間は現在の個体数の大きさには依存しないため、培養細胞が 100 万個から 200 万個になるのにかかる時間は、1000 個から 2000 個になるのにかかる時間と同じである。

    5. 答えの開かれ方(ばらつき)こそが重要な点である。同じ提示レートで同じ 10 時間経過したとき、連続成長は単純な読み方を 24% 上回り、毎時の複利計算を 3% 上回る。

    解答

    ツールは 2225.5409、および倍加時間 8.6643 時間 を出力する。3 つの候補となる答えは 1800、2158.92、および 2225.54 であった。最初と最後の答えの差は 425 個の細胞、すなわち初期個体数の 40% 以上 であり、この差は率についての見解の相違から生じたものではない。成長がどの時点で計上されるかについての見解の相違のみから生じている。これは「70 の法則」の由来でもある。真の値 8.6643 に対し 70/8 = 8.75 と約 1% 高くなるのは、ln 2 = 0.693 であり、約数が多く扱いやすい分子が好まれたためである。

  2. 細胞 1 個はどれだけ生きるのか、そして曲線が個数を数えなくなるのはいつか 6 ステップ

    「一定の減少」プリセットでは、当初 1000 個あった細胞が 10 時間後に 135.3353 個となり、3.4657 ごとに半分になっていく。細胞 1 個の平均寿命を求めよ。そして、この曲線がどの時点から培養細胞群を記述しなくなるのかを判断せよ。

    1. N(t) = 1000e^(−0.2t) であるから、t = 10 では指数部が −2 となり、パネルには 135.3353 と表示される。

    2. 半減期は e^(−0.2T) = ½ から求まり、T = ln 2 / 0.2 = 3.4657 となって、倍加時間が空けた行に収まる。10 時間はその 2.8854 個分にあたり、1000 × 2^(−2.8854) を計算すればふたたび 135.3353 が得られる。別の道をたどって同じ値に着いたわけである。

    3. 半減期は指数の中にある定数ではない。同じ曲線を N = N₀e^(−t/τ) と書き直せば τ = 1/|k| = 5 時間となる。5 時間の時点での個数は 1000/e = 367.8794 である。

    4. この τ が細胞 1 個の平均寿命である。ある細胞が時刻 t にまだ残っている確率は e^(−t/τ) であり、これを 0 から無限大まで積分すればちょうど τ になる。したがって τ/T = 1/ln 2 = 1.4427 であり、この比はあらゆる指数関数的減衰に共通する。平均寿命は半減期より常に 44% ほど長い。半減期は中央値であって、長く生き延びるわずかな個体が平均をその先へ引き延ばすからである。

    5. 次はフラスコが空になる話である。曲線がゼロに達することはないが、フラスコのほうは空になる。生き残る個数の期待値はちょうど N(t) であるから、最後の 1 個が消えると期待されるのは曲線が 1 を横切る時点、すなわち t = τ ln N₀ = 5 ln 1000 = 34.5388 時間である。

    6. その途中のどこかで、表示された数値は個数であることをやめる。t = 30 では曲線が 2.4788 を与えるが、2.4788 個の細胞を抱えた培養など存在したためしがない。この数値は実験を何度も繰り返したときの平均であって、1 つのフラスコが抱えているのは 0 個、1 個、2 個、あるいは 3 個である。

    解答

    半減期 3.4657 に対し、細胞 1 個の平均寿命は 5 時間である。そして 1/ln 2 というこの比は、炭素 14 でもセシウム 137 でも変わらない。フラスコが空になると期待されるのは 34.5388 時間の時点であり、これは新たな計算ではなく、曲線自身が N = 1 を横切る点にほかならない。初期個数を 2 倍にしても、その時点が延びるのは 5 ln 2 = 3.4657 時間、すなわち半減期 1 回分だけであり、それは初期個数がいくつであっても変わらない。もう一方のプリセットで見た倍加時間の議論を、符号を反転させて眺めたものである。N(t) は、値が大きいうちは個数として読み、小さくなったら期待値として読むこと。公式は変わらない。変わるのは、それが何を語っているのかである。

参考文献 (2)

例題

  • 一定の増加 - k = 0.08 は 1000 を t = 10 までに 2225.5409 にし、倍加時間は 8.6643。半減期の行は — を示す
  • 一定の減少 - k = −0.2 では 1000 のうち 135.3353 が残り、半減期は 3.4657。今度は倍加時間のほうが空欄になる
  • 緩やかな増加 - k = 0.02 の倍加には 34.6574 かかり、k = 0.08 の 8.6643 のちょうど4倍。倍加時間が ln(2)/k だからである