数学实验讲义4.ppt

* * 5.18 * Volterra-Lotoka方程中的参数估计 由于Volterra-Lotoka方程的解析解和x,y的显式解很难求出,所以就不能直接用MATLAB函数拟和估计出参数?、?了。下面将给出一种估计方法。 将Volterra-Lotoka方程变形为 在[ti-1,ti]上对两式积分,得 其中 于是得到参数的方程组 其中 故参数的最小二乘解为 由梯形法可近似求出 * 麦蚜与天敌种群数量的模拟 麦蚜种群数量记录(间隔两天调查一次) 29.7 33.1 32.5 69.1 134.2 236.0 269.6 162.3 69.6 39.8 34.0 20.7 21.7 37.6 57.6 124.6 215.8 272.7 195.7 95.0 41.9 25.7 10.9 22.6 33.6 48.1 92.5 183.3 268.5 230.6 111.1 天敌种群数量记录(间隔两天调查一次) 1.6 1.3 1.1 1.2 1.1 1.3 1.8 2.2 2.4 2.1 1.9 1.5 1.5 1.2 0.9 1.1 1.3 1.6 2.3 2.4 2.1 1.7 1.8 1.4 1.2 1.0 0.9 1.1 1.3 1.9 2.3 * 麦蚜与天敌模拟的MATLAB程序 * 估计Volterra-Lotoka方程参数的另一方法 Volterra-Lotoka方程组的两个式子相除可得一个变量可分离方程： 可以求出解析解，再通过拟合确定参数。 * * 生命表组建实验 生命表是表示特定年龄的同年群个体生存和死亡概率的表格,它是种群活力和动态的有效总结,显示出种群死亡数生存数以及未来发展趋势,反映出种群从出生到死亡的数量动态,在种群动态模拟中具有重要意义。 组建一张完整的生命表,一般要通过已知的各时间段(x,x+n)内同年群个体存活数lx来计算:同时段内死亡数dx、死亡概率qx、群体总的生存时间Llx、以后群体总的生存时间Tx、观测时每个个体的平均存活时间(期望寿命)ex。 * 符号 意义 计算公式 T(k)=(x,x+n) 第k段年龄区间 已知,不用算 lx(k) 第k段对应的存活数 已知,不用算 dx(k) 第k段对应的死亡数 qx(k) 第k段的死亡概率 Llx(k) 第k段区间内群体存活时间 Tx(k) 第n(k-1)天后群体总存活时间 ex(k) 第n(k-1)天时观察的个体期望寿命 生命表模型建立与寿命分析 * 果蝇的寿命分析-建立生命表 雄性果蝇存活数据 年龄区间/天 (x,x+n) 0~5 5~10 10~15 15~20 20~25 25~30 30~35 x天时存活数 lx 270 268 264 261 254 251 248 年龄区间/天 (x,x+n) 35~40 40~45 45~50 50~55 55~60 60+ x天时存活数 lx 232 166 130 76 34 13 雄性果蝇寿命分析及生命表组建见EXCEL文件 “生命表.xls”。 雌性果蝇生命表也见文件“生命表.xls”。 实验：用MATLAB程序来实现组建生命表; 并分别画出雌雄果蝇的存活曲线与死亡曲线 (年龄-存活数,年龄-死亡数)。 * 人口预报模型 目前常用的人口预报模型有两个： 模型1： 用指令N=dsolve(DN=r*N*(1-N/K))求解并整理得 即Logistic模型。 模型2： 用指令N=dsolve(DN=-a*N*log(N/K)) 求解并整理得 * 两模型预报效果比较 通过预测误差比较来确认预报效果。 比较指标：相对误差的绝对值之和。 1-16期 人口数据作为训练拟合的样本，17-19期数据作为预测的检验。 * 最优捕鱼策略问题实验 背景：某种鱼分4个年龄组,称1龄鱼,…,4龄鱼。各年龄鱼的平均体重分别为5.07,11.55, 17.86, 22.99(克),它们的自然死亡率均为0.8(/年); 这种鱼是季节性集中产卵平均每条4龄鱼产卵1.109e+9个,3龄鱼减半,1、2龄鱼不产卵,产卵、孵化期为每年最后4个月;卵成活率为1.22e11/(1.22e11+n)。渔业部门规定,每年1-8月可用13mm网眼的拉网捕3、4龄鱼,两种鱼的捕捞强度系数(单位时间内捕捞量与各龄鱼群条数的比值)之比为0.42：1。这种方式称固定量捕捞。 问题：1.如何在每年开始时各龄鱼量基本不变前提下,获取全年最大捕捞总量？2.某渔业公司承包此鱼的捕捞业务5年,承包初各龄鱼数量分别为1.22,29.7,10.1,3.29(?109条);如用固定量捕捞方式,且保证5年后鱼群生产能力不受较大破坏,应采用怎样策略保证总收获量最大。 * 捕鱼问题的分析 假设及符号说明：xi(t)—在t时刻i龄鱼的条数,i=1,2,3,4