ルンゲ・クッタ法を用いた 数値計算 (伝染病モデルより) TA 河合俊典.

Slides:



Advertisements
Similar presentations
1 容積率緩和の便益: 一般均衡論的分析 唐渡 広志 ( 富山大学経済学部 ) 八田 達夫 ( 東京大学空間情報科学研究セン ター ) 2003 年 9 月 19 日.
Advertisements

なぜ貧しい国はなくならないの か 第 1 章 開発経済学とは何か 1. まず定義から始めよう 筆者による定義 「貧しい開発途上国 の貧困削減に貢献する 戦略を研究する学問分 野」 2.
2016 年度 計量経済学 講義内容 担当者: 河田 正樹
中級日本語 第8課 心を伝えるあいさつ 本 文 吉林華橋外国語学院 日本語学部 製作. 第8課 形式段落 第1段落 (ヨーロッパ旅行の経験談1) 外国人の人々は、自国の言葉で声をかけてくれる。 第2段落 (ヨーロッパ旅行の経験談2 ) 日本語で声をかけてくれる人もいて、大変親しみを感じた。 第3段落.
第12回 薬剤師が命を救う 5月11日(金) 午後7:30~7:55 (再)5月12日(土) 午後10:05~10:30
北大における Super-SINET 接続と利用: 2004 年度報告
健康管理上の注意点(デング熱) どんな病気? どんな症状?いつごろ治るの? 診察後, 自宅に帰る場合の注意点
「初夏の水分補給」です。 一日にどれくらい水分をとれば良いのでしょうか?
小樽市における 将来的必要医療・介護病床数 ~行政資料による簡易シミュレーション~
衛生環境課 ノロウイルスによる食中毒について.
口 蹄 疫 侵入防止にご協力をお願いします ◇ 畜産農場へは必要がない限り立ち入らない ◇ 各施設で行われている消毒に協力すること
口蹄疫問題 塩見菜津希.
研究集会 「超大規模行列の数理的諸問題とその高速解法」 2007 年 3 月 7 日 完全パイプライン化シフト QR 法による 実対称三重対角行列の 固有値並列計算 宮田 考史  山本 有作  張 紹良   名古屋大学 大学院工学研究科 計算理工学専攻.
メニューの数学 Lauren Bartlett.
フィードバック制御に基づく 定在波型熱音響エンジンにおける 自励発振条件の特徴付け
情報モラルデザイン A班第2回 DNAのデータベース 吉田.
ソシオン理論における 三者関係のシミュレーション
平成25年度 小平市がん検診受診率向上事業.
時空間データからのオブジェクトベース知識発見
狂牛病(BSE).
社会心理学のStudy -集団を媒介とする適応- (仮)
労務管理の基礎知識 管理職なら 知らないでは済まされない お申込み >> FAX : 0776-63-6202 労務管理セミナー
大阪工業大学 情報科学部 情報システム学科 宇宙物理研究室 B 木村悠哉
薬学物理化学Ⅲ 平成28年 4月15日~.
「真面目さ」と価値観の多様性 ●真面目であること 学生時代: 必ずしも賞賛されない 社会人: 周囲のすべての人から期待される ←「不真面目な人」は運命共同体を危機に陥れる (諸君も、ご両親がまじめでなかったら、大学で学べない。) まじめさこそ社会存立の基盤 ● 「疲れた」と感じたとき → ①身近な人との比較を止める.
4 関数 y=ax 2 1章 関数とグラフ §3 関数 y=ax 2 の値の変化         (5時間)
筑波大学での 勉強、サークル活動、生活について
Mathcad と Excel の比較 PTC ジャパン(株) 部署名 年月日.
Case17 6班 小池・越田・後藤 ・古谷野 3.
誤差の二乗和の一次導関数 偏微分.
症例登録一覧表「円錐切除後HPV感染率調査群」(No. 1)
フィールドセンシング最終課題 Google技術のデータの問題点 総合二年・三崎悠樹・
コンピュータウィルス ~ウィルスへの対処編~
第6章 連立方程式モデル ー 計量経済学 ー.
おさらい 現在、日本では 1日○人のスピードで エイズ感染者が増えている。.
B型・C型肝炎ウイルスについて知ろう 参考資料2 ❏ B型・C型肝炎に関する問い合わせ先 ❏ 広島県健康福祉局薬務課肝炎対策グループ
北大MMCセミナー 第76回 附属社会創造数学センター主催 Date: 2017年10月12日(木) 16:30~18:00
ウイルスについて I98N044 久野耕介 I98N114 藤田和久
経済原論IA 第5回 西村「ミクロ経済学入門」 第5章 消費者需要理論の応用と拡張 京都大学経済学部 依田高典.
第1部 第1篇 第1章 第3節 価値形態または交換価値(A2b)
捕食者と被食者のシミュレーション 望月 陽.
経済学入門7 丹野忠晋 跡見学園女子大学マネジメント学部 2009年6月8日
今後の予定 4日目 10月22日(木) 班編成の確認 講義(2章の続き,3章) 5日目 10月29日(木) 小テスト 4日目までの内容
健康的な職業生活 (保健 最後の授業).
韓国で口蹄疫が再発しました! みなさんの家畜は大丈夫ですか? 牛用 口蹄疫とは・・・
重力レンズ効果による画像の変形と明るさの変化
電機情報工学専門実験 6. 強化学習シミュレーション
~ 獣医師の皆さん、抗菌剤の慎重使用等対策を進め、
シミュレーション論Ⅰ 第14回 シミュレーションの分析と検討.
4. システムの安定性.
シミュレーション物理 大槻東巳.
ベイズ最適化 Bayesian Optimization BO
情報経済システム論:第13回 担当教員 黒田敏史 2019/5/7 情報経済システム論.
カラス撃退装置 ~カラスもたじたじ~ 8班 班長 藤井将之     梅田陽平     大村泰史     藤井秀徳.
バーテックス・ダイナミクス・モデルを用いた上皮組織の創傷治癒における細胞集団運動のシミュレーション
北大MMCセミナー 第72回 附属社会創造数学センター主催 Date: 2017年7月20日(木) 15:00~16:30
1.基本概念 2.母集団比率の区間推定 3.小標本の区間推定 4.標本の大きさの決定
ベイジアンネットワーク概説 第3章 ベイジアンネットワークモデルの 数学的基礎 3.1 ベイジアンネットワークモデルの概要
実験計画法 Design of Experiments (DoE)
大型ホイールのディスク成形における 有限要素シミュレーション 有限要素 シミュレーション 工具と素材形状の最適化 材料の歩留り向上
Mixiの経営戦略 2008年7月2日 吉野 正和.
4) 上級編 ~新型インフルエンザに対する 公衆衛生対応について知る~
北大MMCセミナー 第65回 附属社会創造数学センター主催 Date: 2017年4月20日(木) 16:30~18:00
1.光・音・力.
守谷総合人生大学 (生涯学習の仕組み、学校教育法上の大学ではない) の創設
ATLAS実験におけるSUSY の発見能力
韓国で口蹄疫が再発しました! みなさんの家畜は大丈夫ですか? 豚用 口蹄疫とは・・・
北大MMCセミナー 第23回 Date:2014年3月6日(木) 16:30~18:00 ※通常と曜日が異なります
健康再生療養型リゾート計画.
北大MMCセミナー 第100回 附属社会創造数学センター主催 Date: 2019年7月11日(木) 16:30~18:00
Presentation transcript:

ルンゲ・クッタ法を用いた 数値計算 (伝染病モデルより) TA 河合俊典

伝染病モデルとは dx/dt = -ixy +ry dy/dt = ixy -ry 伝染病の伝わる様子をシミュレート 次式で表される(最もシンプルな場合) dx/dt = -ixy +ry (← 非感染者数増減の式) 感染 治癒 dy/dt = ixy -ry (← 感染者数増減の式) x:非感染者数 y:感染者数 i :感染率 r :治癒率

結果 これを踏まえて、次は応用編へ・・・ 非感染者 感染者 感染率 i = 0.005 (一定) t=12の時にワクチン開発 (t<12の時 r=0、 t≧12の時 r=0.6) ワクチン投与開始      ↓ 非感染者 感染者 これを踏まえて、次は応用編へ・・・

映画「バイオハザードⅢ」公開記念(?) 某工大バイオハザード シミュレーション

映画見たことない人のために ・「Tウイルス」なる恐ろしいウイルスが存在 ・Tウイルスに感染するとゾンビに・・・ ・ただし、感染後数時間は理性が保たれており、 この間に抗ウイルス剤を投与すれば治る。 ・投与せずに放置すればそのままゾンビに。 理性を失ってさまよい歩き、生きている人を襲う。 ・ゾンビにかじられたり、引っかかれたりすると感染。 ・ゾンビになってしまうと、もう抗ウイルス剤は効かない。

バイオハザード発生!! ここは、東京都内の、とある国立大学、 某京工業大学、通称「某工大」。 そこで、事件は起こった・・・。 なんと、本館の地下研究施設で極秘に     Tウイルスの研究を行われていたのだ・・・ しかし、実験中のミスにより、ウイルスが蔓延。 実験中だった研究員2人が知らぬ間に感染し、 そのままゾンビに・・・ 生きている人間を見かけては、襲いかかる! バイオハザード発生!!

しかし、どんどん広がる感染! 果たして某工大の運命やいかに・・・ それでも、感染は広がるばかり。 事件発生から5分後 保健管理センターに一報が入る。 感染拡大防止のために全ての門が閉鎖。 キャンパス内にいた2000人が閉じこめられる。 事件発生から15分後 保健管理センターのおばちゃん、 事件現場に到着。 感染者に対し、抗ウイルス剤の投与を始める。 それでも、感染は広がるばかり。 始めは2人だったゾンビの数も次第に増え始める。 事件発生から60分後 自衛隊が到着 - ゾンビ撃退へ しかし、どんどん広がる感染! 果たして某工大の運命やいかに・・・

この様子をシミュレーション dx/dt = -ixz +ry dy/dt = ixz -ry -dy dz/dt = dy -kz (← 非感染者数増減の式) 感染 治癒 dy/dt = ixz -ry -dy (← 感染者数増減の式) ゾンビ化 dz/dt = dy -kz (← ゾンビ数増減の式) ゾンビ撃退 本シミュレーションでは次のように設定  x(0) = 2000  y(0) = 0  z(0) = 2  i = 0.005  r = 0.0 (t<15) r = 0.3 (t≧15)  d= 0.01  k= 0.0 (t<60) k = 0.3 (t≧60) x:非感染者数 i:感染率 y:感染者数 r:治癒率 z:ゾンビ数 d:ゾンビ化率 k:ゾンビ撃退率

某工大の運命は・・・ 非感染者 感染者 ゾンビ 事件発生から約120分後にはほぼ事態は収束。 どうやら1200人ほどは助かったみたいです。 おばちゃん、頑張る    ↓ 自衛隊到着     ↓ 非感染者 感染者 犠牲者のご冥福をお祈り致します ゾンビ 事件発生から約120分後にはほぼ事態は収束。 どうやら1200人ほどは助かったみたいです。

もしも・・・ ・自衛隊の到着が遅れていたら (t≧90 で k = 0.3) ・保健管理センターのおばちゃんの        処置がなかったら (ずっと治癒率 r = 0) ・某工大生が不真面目だったら (x(0) 減少) ・偶然にも、学園祭の日だったら (x(0) 上昇) ・のら猫にも感染していたら (感染率 i 上昇) いろいろパラメータをいじってみよう! 某工大の運命は果たしてどう変化するのか・・・?