Upgrade to Pro — share decks privately, control downloads, hide ads and more …

2026年度秋学期 応用数学(解析) 第3部・微分方程式に関する話題 / 第10回 生存時間分...

2026年度秋学期 応用数学(解析) 第3部・微分方程式に関する話題 / 第10回 生存時間分布と半減期 (2026. 11. 20)

関西大学総合情報学部 応用数学(解析)(担当・浅野晃)
http://racco.mikeneko.jp/Kougi/2026a/AMA/

Avatar for Akira Asano

Akira Asano PRO

November 11, 2026

Video

More Decks by Akira Asano

Other Decks in Education

Transcript

  1. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  2. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 時刻 t 以後,時間Δの間に死亡する確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  3. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 時刻 t 以後,時間Δの間に死亡する確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  4. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 時刻 t 以後,時間Δの間に死亡する確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  5. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 時刻 t 以後,時間Δの間に死亡する確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  6. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 単位時間 時刻 t 以後,時間Δの間に死亡する確率 あたり 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  7. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 単位時間 時刻 t 以後,時間Δの間に死亡する確率 あたり 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  8. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 単位時間 次の瞬間 時刻 t 以後,時間Δの間に死亡する確率 あたり 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  9. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 単位時間 次の瞬間 時刻 t 以後,時間Δの間に死亡する確率 あたり l(t) は 時刻tまで生存している人が 次の瞬間に死ぬ危険の度合 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  10. 寿命の確率分布を考える 寿命を表す確率変数 T (時刻0に誕生した人が死亡する時刻) 1 l(t) = lim P (t

    < T < t + ∆|T > t) ∆→0 ∆ 時刻 t までは確かに生存している人が 単位時間 次の瞬間 時刻 t 以後,時間Δの間に死亡する確率 あたり l(t) は 時刻tまで生存している人が 次の瞬間に死ぬ危険の度合 [ハザード関数] 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 4 19
  11. 累積分布関数と「生存関数」 確率変数 T に対して F(t) = P(T ≤ t) [累積分布関数]

    2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 5 19
  12. 累積分布関数と「生存関数」 確率変数 T に対して F(t) = P(T ≤ t) [累積分布関数]

    この場合,寿命が t 以下である確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 5 19
  13. 累積分布関数と「生存関数」 確率変数 T に対して F(t) = P(T ≤ t) [累積分布関数]

    この場合,寿命が t 以下である確率 S(t) = 1 − F (t) = P (T > t) [生存関数] 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 5 19
  14. 累積分布関数と「生存関数」 確率変数 T に対して F(t) = P(T ≤ t) [累積分布関数]

    この場合,寿命が t 以下である確率 S(t) = 1 − F (t) = P (T > t) [生存関数] 時刻 t の時点でまだ生きている確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 5 19
  15. 累積分布関数と「生存関数」 確率変数 T に対して F(t) = P(T ≤ t) [累積分布関数]

    この場合,寿命が t 以下である確率 S(t) = 1 − F (t) = P (T > t) [生存関数] 時刻 t の時点でまだ生きている確率 ハザード関数は「瞬間瞬間の死亡の危険」 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 5 19
  16. 累積分布関数と「生存関数」 確率変数 T に対して F(t) = P(T ≤ t) [累積分布関数]

    この場合,寿命が t 以下である確率 S(t) = 1 − F (t) = P (T > t) [生存関数] 時刻 t の時点でまだ生きている確率 ハザード関数は「瞬間瞬間の死亡の危険」 生存関数は,ある時間がたったとき,まだ生きている確率 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 5 19
  17. (ところで)累積分布関数と確率密度関数 ヒストグラム(だったもの)の「へり」 ヒストグラム 柱の面積が確率を表す [確率密度関数] f(t) 連続型になると t グレーの部分の面積 P(T

    ≤ t) t グレーの部分の面積 = P(T ≤ t) [累積分布関数] F(t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 6 19
  18. (ところで)累積分布関数と確率密度関数 ヒストグラム(だったもの)の「へり」 ヒストグラム 柱の面積が確率を表す [確率密度関数] f(t) 連続型になると t グレーの部分の面積 P(T

    ≤ t) t t グレーの部分の面積 = P(T ≤ t) [累積分布関数] F(t) t + Δt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 6 19
  19. (ところで)累積分布関数と確率密度関数 ヒストグラム(だったもの)の「へり」 ヒストグラム 柱の面積が確率を表す [確率密度関数] f(t) 連続型になると t グレーの部分の面積 P(T

    ≤ t) t グレーの部分の面積 = P(T ≤ t) 青い部分 t [累積分布関数] F(t) の面積 = P(t ≤ T ≤ t + Δt) = F(t + Δt) − F(t) t + Δt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 6 19
  20. (ところで)累積分布関数と確率密度関数 ヒストグラム(だったもの)の「へり」 ヒストグラム 柱の面積が確率を表す [確率密度関数] f(t) 連続型になると t グレーの部分の面積 P(T

    ≤ t) t t グレーの部分の面積 = P(T ≤ t) [累積分布関数] F(t) 青い部分 の面積 = P(t ≤ T ≤ t + Δt) = F(t + Δt) − F(t) 青い部分 F(t + Δt) − F(t) の高さ = Δt t + Δt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 6 19
  21. (ところで)累積分布関数と確率密度関数 ヒストグラム(だったもの)の「へり」 ヒストグラム 柱の面積が確率を表す [確率密度関数] f(t) 連続型になると t グレーの部分の面積 P(T

    ≤ t) t グレーの部分の面積 = P(T ≤ t) [累積分布関数] F(t) 青い部分 の面積 = P(t ≤ T ≤ t + Δt) = F(t + Δt) − F(t) 青い部分 F(t + Δt) − F(t) の高さ = Δt F(t + Δt) − F(t) その Δt → 0 の極限 lim = F′(t) Δt→0 Δt t t + Δt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 6 19
  22. (ところで)累積分布関数と確率密度関数 ヒストグラム(だったもの)の「へり」 ヒストグラム 柱の面積が確率を表す [確率密度関数] f(t) 連続型になると t グレーの部分の面積 P(T

    ≤ t) t グレーの部分の面積 = P(T ≤ t) [累積分布関数] F(t) 青い部分 の面積 = P(t ≤ T ≤ t + Δt) = F(t + Δt) − F(t) 青い部分 F(t + Δt) − F(t) の高さ = Δt F(t + Δt) − F(t) その Δt → 0 の極限 lim = F′(t) Δt→0 Δt t t + Δt すなわち f(t) = F′(t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 6 19
  23. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  24. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  25. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  26. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  27. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  28. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  29. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  30. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) F (t) = P (T ≤ t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  31. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) (累積分布関数の定義) F (t) = P (T ≤ t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  32. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) (累積分布関数の定義) F (t) = P (T ≤ t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  33. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) (累積分布関数の定義) F (t) = P (T ≤ t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  34. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) (累積分布関数の定義) F (t) = P (T ≤ t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  35. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) (累積分布関数の定義) F (t) = P (T ≤ t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  36. 生存関数とハザード関数 寿命 T ハザード関数 l(t) 累積分布関数 F(t) 1 l(t) =

    lim P (t < T < t + ∆|T > t) ∆→0 ∆ 含まれる 1 P {(t < T < t + ∆) and (T > t)} · = lim ∆→0 ∆ P (T > t) (条件付確率の定義) 1 P (t < T < t + ∆) = lim · ∆→0 ∆ P (T > t) (累積分布関数の定義) F (t) = P (T ≤ t) 1 F (t + ∆) − F (t) lim = P (T > t) ∆→0 ∆ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 7 19
  37. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T > t) ∆→0 ∆ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  38. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T > t) ∆→0 ∆ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  39. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T > t) ∆→0 ∆ (微分の定義) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  40. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) (微分の定義) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  41. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) (微分の定義) (確率密度関数) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  42. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) (微分の定義) (確率密度関数) f (t) l(t) = S(t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  43. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) (微分の定義) (確率密度関数) f (t) l(t) = S(t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  44. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) f (t) l(t) = S(t) (微分の定義) (確率密度関数) S(t) = P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  45. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) f (t) l(t) = S(t) (微分の定義) (確率密度関数) (生存関数の定義) S(t) = P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  46. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) f (t) l(t) = S(t) (微分の定義) (確率密度関数) (生存関数の定義) S(t) = P (T > t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  47. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) f (t) l(t) = S(t) (微分の定義) (確率密度関数) (生存関数の定義) S(t) = P (T > t) S ′ (t) = (1 − F (t))′ = −F ′ (t) = −f (t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  48. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) f (t) l(t) = S(t) (微分の定義) (確率密度関数) (生存関数の定義) S(t) = P (T > t) S ′ (t) = (1 − F (t))′ = −F ′ (t) = −f (t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  49. 生存関数とハザード関数 1 F (t + ∆) − F (t) l(t)

    = lim P (T ∆ ( > t)) ∆→0 1 ′ F (t) = P (T > t) ′ f (t) = F (t) f (t) l(t) = S(t) (微分の定義) (確率密度関数) (生存関数の定義) S(t) = P (T > t) S ′ (t) = (1 − F (t))′ = −F ′ (t) = −f (t) S ′ (t) 以上から l(t) = − という微分方程式が得られる S(t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 8 19
  50. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  51. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt (両辺を積分) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  52. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt −  t (両辺を積分) l(u)du = log S(t) + C 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  53. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt −  t (両辺を積分) 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 l(u)du = log S(t) + C 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  54. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt −  t (両辺を積分) l(u)du = log S(t) + C 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 t = 0 のとき S(0) = 1 だから 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  55. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − l(u)du = log S(t) + C 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 t = 0 のとき S(0) = 1 だから 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  56. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − l(u)du = log S(t) + C 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 t = 0 のとき S(0) = 1 だから 0 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  57. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − 0 l(u)du = log S(t) + C 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 t = 0 のとき S(0) = 1 だから 0 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  58. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − 0 l(u)du = log S(t) + C 0 0 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 t = 0 のとき S(0) = 1 だから 1 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  59. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − 0 l(u)du = log S(t) + C 0 0 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 t = 0 のとき S(0) = 1 だから 1 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  60. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 0 l(u)du = log S(t) + C 0 0 1 0 t = 0 のとき S(0) = 1 だから C=0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  61. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 0 l(u)du = log S(t) + C 0 0 1 0 t = 0 のとき S(0) = 1 だから C=0   t  l(u)du という解が得られる よって S(t) = exp − 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  62. 微分方程式を解く l(t) = − ′ S (t) S(t) d =

    − (log S(t)) dt 0 (両辺を積分) t − 時刻0,つまり誕生の瞬間に生存している確率は1 つまり S(0) = 1 0 l(u)du = log S(t) + C 0 0 1 0 t = 0 のとき S(0) = 1 だから C=0   t  l(u)du という解が得られる よって S(t) = exp − 0 ハザード関数と生存関数の関係 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 9 19
  63. ワイブル分布 ハザード関数を l(t) = λp(λt) p−1 と仮定する   t

     S(t) = exp − l(u)du に代入 0 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 11 19
  64. ワイブル分布 ハザード関数を l(t) = λp(λt) p−1 と仮定する   t

     S(t) = exp − l(u)du に代入 0   t  p−1 S(t) = exp − λp(λu) du 0   u=t p p = exp − [(λu) ]u=0 = exp (−(λt) ) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 11 19
  65. ワイブル分布 ハザード関数を l(t) = λp(λt) p−1 と仮定する   t

     S(t) = exp − l(u)du に代入 0   t  p−1 S(t) = exp − λp(λu) du 0   u=t p p = exp − [(λu) ]u=0 = exp (−(λt) ) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 11 19
  66. ワイブル分布 ハザード関数を l(t) = λp(λt) p−1 と仮定する   t

     S(t) = exp − l(u)du に代入 0   t  p−1 S(t) = exp − λp(λu) du 微積分の関係 0   u=t p p = exp − [(λu) ]u=0 = exp (−(λt) ) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 11 19
  67. ワイブル分布 ハザード関数を l(t) = λp(λt) p−1 と仮定する   t

     S(t) = exp − l(u)du に代入 0   t  p−1 S(t) = exp − λp(λu) du 微積分の関係 0   u=t p p = exp − [(λu) ]u=0 = exp (−(λt) ) p F (t) = 1 − S(t) = 1 − exp (−(λt) ) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 11 19
  68. ワイブル分布 ハザード関数を l(t) = λp(λt) p−1 と仮定する   t

     S(t) = exp − l(u)du に代入 0   t  p−1 S(t) = exp − λp(λu) du 微積分の関係 0   u=t p p = exp − [(λu) ]u=0 = exp (−(λt) ) p F (t) = 1 − S(t) = 1 − exp (−(λt) ) この形の累積分布関数をもつ確率分布を[ワイブル分布]とよぶ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 11 19
  69. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p 2026年度秋学期

    応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  70. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  71. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる 死亡・故障する危険が どの時刻でも大きくなる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  72. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる 死亡・故障する危険が どの時刻でも大きくなる p > 1 のときは, 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  73. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) p−1 死亡・故障する危険が どの時刻でも大きくなる の指数が正 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  74. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) p−1 死亡・故障する危険が どの時刻でも大きくなる の指数が正 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  75. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  76. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  77. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる [摩耗故障] 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  78. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる [摩耗故障] 0 < p < 1 のときは, 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  79. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる [摩耗故障] 0 < p < 1 のときは, l(t) = λp(λt) p−1 の指数が負 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  80. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる [摩耗故障] 0 < p < 1 のときは, l(t) = λp(λt) p−1 の指数が負 時間が経つにつれて,死亡・故障する危険が小さくなる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  81. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる [摩耗故障] 0 < p < 1 のときは, l(t) = λp(λt) p−1 の指数が負 時間が経つにつれて,死亡・故障する危険が小さくなる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  82. ワイブル分布のパラメータ l(t) = λp(λt) p−1 パラメータは λ と p λ

    が大きいと,ハザード関数が全体に大きくなる p > 1 のときは, l(t) = λp(λt) 死亡・故障する危険が どの時刻でも大きくなる p−1 の指数が正 時間が経つにつれて,死亡・故障する危険が大きくなる [摩耗故障] 0 < p < 1 のときは, l(t) = λp(λt) p−1 の指数が負 時間が経つにつれて,死亡・故障する危険が小さくなる [初期故障] 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 12 19
  83. ワイブル分布のパラメータ p = 2 の場合と p = 4 の場合 どちらも摩耗故障(時間につれて故障しやすくなる)

    p = 4 のほうが,急激に故障が増える 累積分布関数 (ある時刻までに死亡・ 故障したものの割合) F(t) 2 –t F(t) = 1 – e 4 –t F(t) = 1 – e t 経過時間 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 13 19
  84. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  85. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t)  log log  1 S(t)  両辺の対数を2回とる = log {log (exp ((λt)p ))} 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  86. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t)  log log  1 S(t)  両辺の対数を2回とる = log {log (exp ((λt)p ))} p = log {(λt) } = p(log t + log λ) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  87. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t)  log log Y  1 S(t)  両辺の対数を2回とる = log {log (exp ((λt)p ))} p = log {(λt) } = p(log t + log λ) 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  88. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t)  log log Y  1 S(t)  両辺の対数を2回とる = log {log (exp ((λt)p ))} p = log {(λt) } = p(log t + log λ) X 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  89. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t)  log log Y  1 S(t)  両辺の対数を2回とる = log {log (exp ((λt)p ))} p = log {(λt) } = p(log t + log λ) Y = p(X + log λ) X 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  90. ワイブルプロット 実務では,たくさんの個体で耐久試験を行い, ワイブル分布を仮定して,パラメータを推測する 1 p p = exp ((λt) )

    S(t) = exp (−(λt) ) より S(t)  log log Y  1 S(t)  両辺の対数を2回とる = log {log (exp ((λt)p ))} p = log {(λt) } = p(log t + log λ) Y = p(X + log λ) X 時刻を上の X ,その時刻での生存割合を上の Y に変換してプロット →並びを近似する直線の傾きが p 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 14 19
  91. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  92. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  93. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 死亡・故障する危険が時刻によらず一定 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  94. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 死亡・故障する危険が時刻によらず一定 [偶発故障] 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  95. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 死亡・故障する危険が時刻によらず一定 [偶発故障] 累積分布関数は F (t) = 1 − e 生存関数は S(t) = e −λt −λt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  96. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 死亡・故障する危険が時刻によらず一定 [偶発故障] 累積分布関数は F (t) = 1 − e 生存関数は S(t) = e −λt [指数分布] −λt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  97. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 死亡・故障する危険が時刻によらず一定 [偶発故障] 累積分布関数は F (t) = 1 − e 生存関数は S(t) = e −λt [指数分布] −λt 放射性原子核は,どの時刻においても,その時点で 存在する核のうち一定の割合が崩壊する 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  98. 指数分布 l(t) = λp(λt) p−1 で p = 1 の場合

    ハザード関数は l(t) = λ 死亡・故障する危険が時刻によらず一定 [偶発故障] 累積分布関数は F (t) = 1 − e 生存関数は S(t) = e −λt [指数分布] −λt 放射性原子核は,どの時刻においても,その時点で 存在する核のうち一定の割合が崩壊する ハザード関数が一定で,指数分布にしたがう 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 15 19
  99. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 S(t) = e−λt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  100. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 S(t) = e−λt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  101. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 e −λt′ 1 −λt = e 2 S(t) = e−λt 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  102. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 e −λt′ 1 −λt = e 2 S(t) = e−λt 対数をとる 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  103. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 e 対数をとる −λt′ 1 −λt = e 2 S(t) = e−λt ′ −λt = − log 2 − λt log 2 ′ t −t= λ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  104. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 e 対数をとる −λt′ 1 −λt = e 2 S(t) = e−λt ′ −λt = − log 2 − λt log 2 ′ t −t= 原子核の数が半分になるまでの時間 λ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  105. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 e 対数をとる −λt′ 1 −λt = e 2 S(t) = e−λt ′ −λt = − log 2 − λt log 2 ′ t によらず一定 t −t= 原子核の数が半分になるまでの時間 λ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  106. 半減期 ある時刻に存在する原子核の数が,その半分になるまでの時間は, どの時刻でも一定 時刻 t に存在する原子核の数が半分になる時刻を t′ とする 1 ′

    S(t ) = S(t) 2 指数分布の生存関数 e 対数をとる −λt′ 1 −λt = e 2 S(t) = e−λt ′ −λt = − log 2 − λt log 2 ′ t によらず一定 [半減期] t −t= 原子核の数が半分になるまでの時間 λ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 16 19
  107. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 指数分布の生存関数を S(t) = e −λt 時間 t

    の単位を「年」とする とおくと log 2 半減期 = λ 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 17 19
  108. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 指数分布の生存関数を S(t) = e log 2 半減期

    = λ −λt 時間 t の単位を「年」とする とおくと log 2 半減期は2年なので λ = 2 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 17 19
  109. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 指数分布の生存関数を S(t) = e log 2 半減期

    = λ −λt 時間 t の単位を「年」とする とおくと log 2 半減期は2年なので λ = 2 求めるのは,指数分布の累積分布関数 F(t) = 1 − e −λt において F(1) の値 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 17 19
  110. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 求めるのは,指数分布の累積分布関数 F(t) = 1 − e −λt

    において F(1) の値 log 2 半減期は2年なので λ = 2 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 18 19
  111. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 求めるのは,指数分布の累積分布関数 F(t) = 1 − e −λt

    において F(1) の値 log 2 半減期は2年なので λ = 2 F(1) = 1 − e log 2 − 2 ⋅1 =1−e −1) log(2 2 =1−2 1 −2 =1− 1 2 ≒ 0.293 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 18 19
  112. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 求めるのは,指数分布の累積分布関数 F(t) = 1 − e −λt

    において F(1) の値 log 2 半減期は2年なので λ = 2 F(1) = 1 − e log 2 − 2 ⋅1 =1−e −1) log(2 2 =1−2 1 −2 =1− 1 2 ≒ 0.293 1個の原子が1年以内に崩壊する確率 0.293 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 18 19
  113. 演習問題 ある原子核の半減期が2年であるとする 1年以内に崩壊する原子核の割合は? 求めるのは,指数分布の累積分布関数 F(t) = 1 − e −λt

    において F(1) の値 log 2 半減期は2年なので λ = 2 F(1) = 1 − e log 2 − 2 ⋅1 =1−e −1) log(2 2 =1−2 1 −2 =1− 1 2 ≒ 0.293 1個の原子が1年以内に崩壊する確率 0.293 原子がたくさんあれば,そのうち崩壊する原子の割合が 29.3% 2026年度秋学期 応用数学(解析) / 関西大学総合情報学部 浅野 晃 18 19