!!! サイト改修中のため表示が乱れる場合があります(1月末頃まで) !!!
数値解法 / 数値シミュレーション

5-22. 摂取カロリーと運動カロリーによる体重変化の数理モデルとシミュレーション

やること

「数理モデルって実際どんなもんなの?」という質問を受けたので(受けてない定期)、毎日の摂取カロリー・運動カロリーによって体重がどのように増減するかを数理モデル化・シミュレーションしてみましょう。

数理モデル化の参考になれば幸いです。コードは付きません。

数理モデル

適当に作ります。xは体重 [kg] です。※テキトーなので実際の健康管理には用いないでください

     \begin{equation*} \begin{split} \frac{dx}{dt} =  k_1 (a - (b_1 + b_2)) + k_2 (60 - x) + k_3 max(60 - x, 0) \end{split} \end{equation*}

パラメータはこちら

#標準体重 [kg]
60(モデルにベタ書き)

#標準摂取カロリー [kcal/day]
a = 1800

#基礎代謝 [kcal/day]
b1 = 1500

#標準運動カロリー [kcal/day]
b2 = 300

#係数
k1 = 0.0015
k2 = 0.0120
k3 = k2 * 0.8

シミュレーション条件はこちら

#体重初期値 [kg]
x_init = 60

#シミュレーション時間 [day]
t_end = 365

#時間きざみ幅 [day]
dt = 0.01

シミュレーション結果

まず、摂取カロリーの違いによる体重変化です。毎日1800 kcal 摂取すれば、基礎代謝1500 kcal と標準運動カロリー300 kcal でちょうど打ち消されるので標準体重60 kg を維持できます。たくさん食べると太ります。

次に、運動カロリーの違いによる体重変化です。標準運動カロリーは300 kcal ですが、毎日100 kcal しか運動しないと1年後には84 kg になります。

次に、摂取カロリーと運動カロリーを変えたときの平衡体重(終端体重)をヒートマップにしました。「目指す体重」と「通勤で発生する運動カロリー」から摂取カロリーを逆算したりできます。

半年間はあまり運動せず、その後半年間はちょっとだけ運動を頑張った場合の体重推移です。頑張ったと言っても標準運動カロリーには及んでいませんが、それでも平衡体重である72 kg に向かって痩せていきます。

最後に、偶数月だけチートデーして奇数月は標準的なご飯にした場合です。体重は80 kg 前後で振動します。

さいごに

数理モデルの構築はお任せください ᐠ(  ᐢ ᵕ ᐢ )ᐟ

リアクションのお願い

「参考になった!」「刺激された!」と思ったらぜひリアクションをしましょう。エンジニアの世界はGive and Takeによって成り立っています。これからも無料で良質な情報にアクセスできるよう、Giveする人への感謝をリアクションで示しましょう!

この記事をシェアする

自身のブログ等で使用する場合は引用を忘れずに!

また、寄付も受け付けています。コーヒー1杯でとても喜びます(*˘︶˘*)

 Amazonでギフト券(アマギフ)を贈る

こちらのリンク から金額を指定してお贈りください。(デフォルトで10000円になっているのでご変更ください)

配送:Eメール
受取人:staffあっとvigne-cla.com
贈り主:あなたのお名前やニックネーム
メッセージ:◯◯の記事が参考になりました。など

のようにご入力ください。見返りはありませんのでご了承ください。

 Amazonで食事券(すかいらーく優待券)を贈る

500円 1000円 2000円 5000円 からお贈りください。

配送:Eメール
受取人:staffあっとvigne-cla.com
贈り主:あなたのお名前やニックネーム
メッセージ:◯◯の記事が参考になりました。など

のようにご入力ください。見返りはありませんのでご了承ください。

 その他、ギフト券やクーポン券をメールで贈る

デジタルのギフト券/クーポン券はメールアドレス(staffあっとvigne-cla.com)までお送りください。受領の返信をいたします。
紙のギフト券/クーポン券は 「郵便物はこちらへ」の住所 まで送付してください。名刺やメールアドレスを同封していただければ受領の連絡をいたします。
余った株主優待券等の処理におすすめです。
いずれも見返りはありませんのでご了承ください。

不明点はSNSでお気軽にご連絡ください

ビネクラのTwitter・Youtubeでコメントをください!


Slack・Discordの場合はこちらの公開グループに参加してShoya YasudaまでDMをください!


※当ブログに関することは何でもご相談・ご依頼可能です。

この記事を書いた人
Yasuda

博士(理学)。専門は免疫細胞、数理モデル、シミュレーション。米国、中国で研究に携わった。遺伝的アルゴリズム信者。物価上昇のため半額弁当とともに絶滅寸前。

タイトルとURLをコピーしました