数模论坛

 找回密码
 注-册-帐-号
搜索
热搜: 活动 交友 discuz
查看: 3062|回复: 3

求助!一道传染病建模(SIR)及拟合的题。万分感谢!

[复制链接]
发表于 2005-1-5 22:44:38 | 显示全部楼层 |阅读模式
<>某种传染病的潜伏期为数小时至21天,其每天的发病与死亡数据如下,请高手帮我建一个SIR微分方程并拟合。SIR模型是Reed和Frost在1920年的一篇未发表论文中首先提出的.它把人群分为三组:未患病但可感染(suscepti ble)、病人且可传染别人(infected)、治愈且终生免疫(re covered).这就是SIR模型,多谢!</P>
<>数据为:</P>
<>
<TABLE  cellSpacing=0 cellPadding=0 width=360 border=0 x:str>
<COLGROUP>
<COL style="WIDTH: 54pt" width=72>
<COL style="WIDTH: 54pt" span=4 width=72>

<TR  height=28>
<TD class=xl24  width=72 height=28><FONT face=Tahoma size=1>天数</FONT></TD>
<TD class=xl24  width=72><FONT face=Tahoma size=1>日增发病数</FONT></TD>
<TD class=xl24  width=72><FONT face=Tahoma size=1>日增死亡数</FONT></TD>
<TD class=xl24  width=72><FONT face=Tahoma size=1>日累积发病数</FONT></TD>
<TD class=xl24  width=72><FONT face=Tahoma size=1>日累积死亡数</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>1</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>280</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>280</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>280</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>280</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>2</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>280</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>280</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>3</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2000</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>1300</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2280</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>1580</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>4</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>1400</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>1400</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>3680</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2980</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>5</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>3680</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2980</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>6</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>3680</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2980</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>7</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>3680</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2980</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>8</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>147</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>80</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>3827</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>3060</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>9</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>4820</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>4820</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>8647</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>7880</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>10</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>33211</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>33211</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>41858</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>41091</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>11</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>41858</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>41091</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>12</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2560</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>1390</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>44418</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>42481</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>13</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>7509</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>5266</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>51927</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>47747</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>14</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>51927</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>47747</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>15</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>57801</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>57545</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>109728</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>105292</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>16</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>109728</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>105292</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>17</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>109728</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>105292</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>18</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>17566</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>13631</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127294</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>118923</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>19</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2810</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>2233</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>130104</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>121156</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>20</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>130104</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>121156</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>21</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>10794</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>4708</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>140898</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>125864</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>22</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>1161</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>753</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142059</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>126617</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>23</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142059</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>126617</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>24</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>449</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>449</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142508</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127066</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>25</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>432</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>410</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142940</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127476</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>26</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142940</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127476</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>27</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142940</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127476</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>28</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>0</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>142940</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127476</FONT></TD></TR>
<TR  height=19>
<TD class=xl24  width=72 height=19 x:num><FONT face=Tahoma size=1>29</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>219</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>219</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>143159</FONT></TD>
<TD class=xl24  width=72 x:num><FONT face=Tahoma size=1>127695</FONT></TD></TR></TABLE></P>
 楼主| 发表于 2005-1-5 22:47:35 | 显示全部楼层
<>我的QQ是77554044,E-MAIL <a href="mailtlijing3719@163.com" target="_blank" >lijing3719@163.com</A>,</P><>我是新手,肯请各位帮助或指点应该学哪方面的知识来解决这个模型问题?谢谢!</P>
发表于 2005-9-15 01:15:40 | 显示全部楼层

<>[em05]</P>
<>[em07]</P>
发表于 2005-9-14 04:08:19 | 显示全部楼层
[em06]
您需要登录后才可以回帖 登录 | 注-册-帐-号

本版积分规则

小黑屋|手机版|Archiver|数学建模网 ( 湘ICP备11011602号 )

GMT+8, 2024-11-30 11:38 , Processed in 0.050169 second(s), 19 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表