重庆大学数学实验实验四.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1421 总 成 绩 教师签名 数 理 学 院 制 开课学院、实验室： 实验时间 ： 年 月 日 课程 名称 实验项目 名 称 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 实验目的 [1] 归纳和学习求解常微分方程(组)的基本原理和方法； [2] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； [3] 熟悉MATLAB软件关于微分方程求解的各种命令； [4] 通过范例学习建立微分方程方面的数学模型以及求解全过程； 通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是十分必要的。 应用实验（或综合实验） 一、实验内容 最优捕鱼策略 为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源）的开发必须适度。一种合理、简化的策略是，在实现可持续收获的前提下，追求最大产量或最佳效益。 ????考虑对某种鱼（鲳鱼）的最优捕捞策略： ????假设这种鱼分４个年龄组：称１龄鱼，……，４龄鱼。各年龄组每条鱼的平均重量分别为5.07,11.55,17.86,22.99（克）；各年龄组鱼的自然死亡率均为0.8（１/年）；这种鱼为季节性集中产卵繁殖，平均每条４龄鱼的产卵量为1.109×105(个）；３龄鱼的产卵量为这个数的一半，２龄鱼和１龄鱼不产卵，产卵和孵化期为每年的最后４个月；卵孵化并成活为１龄鱼，成活率（１龄鱼条数与产卵总是n之比）为1.22×1011/(1.22×1011+n). ????渔业管理部门规定，每年只允许在产卵孵化期前的８个月内进行捕捞作业。如果每年投入的捕捞能力（如渔船数、下网次数等）固定不变，这时单位时间捕捞量将与各年龄组鱼群条数成正比。比例系数不妨称捕捞强度系数。通常使用１３mm网眼的拉网，这种网只能捕捞３龄鱼和４龄鱼，其两个捕捞强度系数之比为0.42:1。渔业上称这种方式为固定努力量捕捞。 ????１）建立数学模型分析如何可持续捕获（即每年开始捕捞时渔场中各年龄组鱼群条数不变），并且在此前提下得到最高的年收获量（捕捞总重量）。 ????２）某渔业公司承包这种鱼的捕捞业务５年，合同要求5年后鱼群的生产能力不能受到太大的破坏。已知承包时各年龄组鱼群的数量分别为：122,29.7,10.1,3.29(×109条），如果仍用固定努力量的捕捞方式，该公司采取怎样的策略才能使总收获量最高。 二、问题分析 对自然死亡率的理解： 本题中给出的鱼的自然死亡率是指平均死亡率，即单位时间鱼群死亡数量与现有鱼群数量的比例系数，它与环境等其它外在因素无关；这是一个有量纲的量，它既不是简单的百分率又不是简单的变化速率，实际上它是百分比率的变化率。它应该理解为以每年死亡80%的速率减少，并不是在一年内恰好死亡80%。另一方面，鱼群的数量是连续变化的，且1，2龄鱼在全年及3，4龄鱼在后4个月的数量只与死亡率有关。由此可知，各龄鱼的变化满足： 对捕捞强度系数的理解： 捕捞强度系数是单位时间内捕捞量与各年龄组鱼群条数的比例系数，单位时间4龄鱼捕捞量与4龄鱼群总数成正比，捕捞强度系数是一定的，且只在捕捞期内（即每年的前8个月）捕捞3，4龄鱼。所以，捕捞强度系数k影响了3，4龄鱼在捕捞期内的数量变化： 则有 对卵的成活率的理解： 1，2龄鱼不产卵，3，4龄鱼在每年的后四个月产卵，我们假设了在9月初一次产卵，因此可将每年的产卵量n表示为： 对于问题一，要实现可持续捕获，即每年开始捕捞时渔场中各年龄组鱼群条数不变，因此我们要算出每年初各龄鱼组的数量。 前面已对自然死亡率，捕捞强度系数和卵成活率作出了解释，即1，2龄鱼仅受自然死亡率的影响；而3，4龄鱼不仅受自然死亡率的影响，还受捕捞强度系数的影响；因为该种鱼的最高寿命为4，所以在后四月中4龄鱼都不存活