　精密工学科の学生の皆さん．慣れない自宅での履修苦労されているかと思います．いま社会で起こっている新型コロナウイルス感染拡大の最前線の現場は病院等の医療現場です．しかし，私たちエンジニアもけっして別の世界に居るわけではありません．たとえば，マシンダイナミクスで学習するバネ・マス・ダッシュポット系の運動は微分方程式で表すことができます．実は同じように感染症の流行も微分方程式で表現をすることができます．

　学科開講科目のマシンダイナミクスは機械系の振動問題についての講義が主ですが，この機会に，その中で番外編として少しだけ紹介している感染症流行のシミュレーションについて再度紹介をしたいと思います．すでにマシンダイナミクスを受講した学生も受講していない学生もいることでしょう．受講済の学生は思い出してみてください．これから受講する学生はお楽しみに！

　よいエンジニアになるためには，目の前で起こっていることを批判的に（ホントかどうか？），定量的に（何秒？何％？）評価できるようになることが大切です．ずっと家にいて窮屈な日々を過ごしているかと思いますが，エンジニアとしてできることは何か？，ニュースで流れている情報は本当か？，ノートと鉛筆そして，エンジニアだから電卓（パソコン）を使ってじっくり考えてみましょう．

１．意外と単純！？微分方程式を使ったモデル化

　さて本題に入ります．古典的な感染症の流行モデルにSIRモデルというのがあります．Sは未感染者，Iは感染者，Rは免疫保持者のことです．この３つの集団を考えて，その振る舞いを数式で表現することにしましょう．即ち変数として，

　S：未感染者数，I：感染者数，R：回復者数（免疫保持者数）

を考えます．S + I + R は考えている系（たとえば日本）の全体の人数になります．ここでは，一度感染して回復した人には免疫ができて二度と感染しなくなると仮定します．

　まず， I（感染者数）がどのように増えるのかを考えてみます．感染するシチュエーションを想像してください．感染は感染者と未感染者が接触することで起こります．したがって，「感染者という集団に属する人数(I)」が多いほど，かつ「未感染者という集団に属する人数(S)」が多いほど接触の機会が多くなります．したがって，I × Sに比例して感染者が増えやすくなると考えます．また，長い時間が経過するほど接触の機会が増えますので，いま考えている時間経過の長さΔtにも比例すると考えます．ここまでの議論から，感染者数Iは S × I × Δt に比例して増えそうだ！と，とても単純化して考えます．でもちょっとまってください．接触したら確実に感染するわけではないですよね．そこで，感染率βを導入します．感染率βとは，単位時間あたりに一人の感染者が未感染者の中から新たに感染者を作り出す割合のことです．以上より，いま考えている時間Δtの間に新たに感染する人数は β × I × S × Δt となります．

　どうでしょうか？「ちょっと待てよ」と思った人もいるでしょう．たしかに，新たに感染する人数は上記のβ × I × S × Δt ですが，時間が経過するにしたがって，回復して，感染者（I）の集団から離脱する人も居るはずです．その人たちも考慮して感染者数（I）の変化を考える必要があります．その回復者はどのように考えればよいでしょうか．そもそも感染している人が居なければ回復する人も居ないので，考えている時間Δtの間に回復する人の数は感染者の数(I)に比例すると考えることにします．したがって，Δtの間に新たに回復する人の数はI × Δtに比例すると考えます．ここでも，全員がその時間の間に回復するわけではありませんから，回復率γを導入します．回復率γは単位時間あたりに感染者の中から回復する割合です．これを掛けると，Δtの間に新たに回復する人の数はγ × I × Δtとなります．

　以上では新たに感染する人数はβ × I × S × Δt であり，新たに回復する人数はγ × I × Δtであることを説明しました．このことから，いま考えている時間Δtの間の感染者の変化ΔIを求めると，単純に引き算をして，

　　　　　　ΔI = （　β × S × I × Δt　） -（ γ × I × Δt　）

となります．

　これを，Δtは限りなく小さいと考えて極限をとると，

　　　　　dI/dt =　β × S × I -　γ × I

が得られます．これが，感染者数（I）の振舞いを表す微分方程式です．また，上記の考察から落ち着いて考えると，未感染者数については，

　　　　　dS/dt = - β × S × I

回復者（免疫保持者）については，

　　　　　dR/dt = γ × I

となります．ここは，ゆっくり落ち着いて考えてみたら理解できると思います．

2020.04