伝染病の方程式
最近、豚インフルエンザの話を聞いていて大昔に読んだ伝染病の方程式のことを思い出した。元ネタはこの数学・物理100の方程式です。
全員でA人いるとして、病人の数をy、病気にかかっていない人の数をxとする。
病気の伝染速度はまあy’みたいだと思うと、だいたい速度は病人がどれだけかかっていない人に接触するかだからxとyの積に比例するとする。でとりあえず自然に治っていくとすると、
y' = αxy - βy = α(A-y)y - βy
となる。βが治癒を表して、αがかかりやすさを表す。でこれって単なるロジスティック方程式なので普通に(数値解析しなくても)解けてしまう。
y(t) = (A-β/α) / ( 1 + ( (A-β/α)/y0 - 1 ) exp ( -α (A-β/α)t)
となる。
A-β/α>0のときは
のようにどこかで飽和する。
A-β/α<0のときは
のように駆逐されてしまう。
まあ実際の伝染はこんな単純じゃないので単なる遊びですが、面白いなと。
« 御香宮へ行ってきた。 | トップページ | AudiのCMのピアノ曲はAsturias(アストゥリアス)ですね。 »
「学問・資格」カテゴリの記事
- Visual Basic (VB.NET)でC#用の数値計算ライブラリMath.NET Numericsを使う(1)複素行列を定義して一次方程式や逆行列、行列式などを計算する。(2023.03.24)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース(3/16) IEEE Microwave Magazineは女性マイクロ波研究者特集、Microwave Journalで285GHz帯で30GHz帯域のOTA測定!QorvoがUWB室内ナビのデモ、STMとonsemiのBluetooth新商品、IDTechExの6Gレポート、など。(2023.03.16)
- KeysightのADSで位相を±180°に限らずに連続にする関数をよく聞かれるがいつも忘れる…unwrap()だ。PythonのNumPyでもあるので(matlabにもある)いい加減に覚えたい。とりあえずPythonでやってみて記憶する。(2023.03.02)
- 50万人が毎年受ける試験で採用、“謎”のプログラミング言語「DNCL」という記事を見た。大学入試センターで2022年1月に出している仕様(日本語プログラムぽい)と、令和7年度試験の問題作成の方向性,試作問題等で出している仕様(Pythonを日本語にしたっぽい)違うのか…(2023.03.01)
- LC共振回路のモンテカルロシミュレーションでL,Cを一様乱数で振って共振周波数を見る、、、と全然一様じゃない。そりゃそうだ。何を勘違いしてたんだろう…(2023.02.27)
トラックバック
この記事へのトラックバック一覧です: 伝染病の方程式:
» Schweinegrippe [ごんざぶログ]
こちらのニュースでも出ています。
Erste Schweinegrippe-Flle in Ja [続きを読む]
« 御香宮へ行ってきた。 | トップページ | AudiのCMのピアノ曲はAsturias(アストゥリアス)ですね。 »
コメント