鄱阳湖国家级自然保护区东方白鹳的种群数量的变化(数学建模).docxVIP

鄱阳湖国家级自然保护区东方白鹳的种群数量的变化 及越冬气候对其影响的影响 问题重述 东方白鹳（ Ciconia boyciana）全球种群现存仅~3 000 只个体，被 IUCN 红皮书列为濒危物种，也是我国国家 I 级重点保护野生动物。东方白鹳繁殖地主要位于中国和俄罗斯交界的黑龙江和乌苏里江流域，越冬地为我国长江中、下游地区（ Birdlife International， 2013）。其中，鄱阳湖区是东方白鹳在长江中、下游地区的主要越冬地。近年来，该区越冬东方白鹳种群数量表现出较大波动。 研究表明，气候变化会对鸟类的种群数量变化产生较强的影响。例如，极端气候事件可能会直接导致死亡率的迅速提高，甚至导致局部种群的灭绝。天气条件也可能通过改变食物资源的数量和质量影响成体存活率和繁殖成功率，进而影响鸟类的种群数量。该影响具有一定的时滞性，即在生活史某一阶段的气候条件对种群的影响需经过一定的时间方才显现。某些留鸟成体的存活率与非繁殖期的天气条件存在显著相关性，越冬期的不利天气将导致迁徙鸟类繁殖成功率下降，第二年种群数量下降。大多数迁徙鸟类种群数量的年际变化主要受其繁殖地生存条件的影响，但也有研究发现某些迁徙鸟类的种群数量受其越冬地生存条件的影响。 通过构建基于越冬地气候条件的种群数量波动模型，来分析鄱阳湖国家级自然保护区（以下简称保护区） 1982—2011 年越冬东方白鹳种群数量的年际变化与越冬地气候变化的相关性及时滞效应。 问题分析 （一）种群数目 表1 1982-2011年鄱阳湖保护区东方白鹳的种群数量变化（单位：只） 年份 1982 1983 1984 1985 1986 1987 1988 1989 1990 种群数目 200 110 403 319 205 758 539 425 500 年份 1991 1992 1993 1994 1995 1996 1997 1998 1999 种群数目 507 321 199 621 684 257 137 178 968 年份 2000 2001 2002 2003 2004 2005 2006 2007 2008 种群数目 1780 500 1920 1873 928 2544 2376 1000 1865 年份 2009 2010 2011 种群数目 2167 3446 4052 温度变化 表3 1982-2011年鄱阳湖保护区温度变化（单位：℃） 年份 1982 1983 1984 1985 1986 1987 1988 1989 1990 温度 17.5 17.9 16.8 17.4 17.6 17.7 17.5 17.4 18.0 年份 1991 1992 1993 1994 1995 1996 1997 1998 1999 温度 17.7 17.8 17.4 18.2 17.8 17.6 17.8 18.8 18.2 年份 2000 2001 2002 2003 2004 2005 2006 2007 2008 温度 17.9 18.2 18.4 18.6 18.7 18.2 18.6 19.3 18.4 年份 2009 2010 2011 温度 18.8 18.5 18.4 （来源：南昌气象站） 符号说明 t～年份（年）； xt～第k年东方白鹳种群数目（只）； rt～用前差公式计算的东方白鹳在第t年的增长率； r～东方白鹳的固有增长率； N～东方白鹳的最大容量。 模型建立求解和模型检验 模型A—自然增长模型（Malthus） 脚本： t=1:1:30;//此处将年份进行修改，将1982对应为1,1983对应为2，以此类推// x(t)=[200,110,403,319,205,758,539,425,500,507,321,199,621,684,257,137,178,968,1780,500,1920,1873,928,2544,2376,1000,1865,2167,3446,4052]; y=log(x(t));a=polyfit(t,y,1); r=a(1);x0=exp(a(2)); x1=x0*exp(r*t); plot(t,x(t),r,t,x1,b); 即得：x(t)=200e0.0900t （二）模型B—离散阻滞模型（logistic） 即 Curvefit_fun2.m脚本为： function f=curvefit_fun2 (a,t) f=a(1)./(1+(a(1)/200-1)*exp(-a(2)*(t-1982))); End G2脚本为： x=1982:1:2011; y=[200,110,403,319,205,758,539,425,500,507,321,199,621,684,257,137,178,968,1780