假设某渔场内养殖有某种鱼群(如鲳鱼),分为4个年龄组,分别称为1龄鱼,2龄鱼,3龄鱼和4龄鱼。已知: 各年龄组每条鱼的平均重量分别为5.07,11.55,17.86和22.99克; 各年龄组鱼的自然死亡率均为0.8条/年; 这种鱼为季节性集中产卵繁殖,平均每条4龄鱼的产卵量为1.109×105个,3龄鱼的产卵量为4龄鱼的一半,2龄鱼和1龄鱼不产卵,鱼群产卵和孵化期为每年的最后4个月; 如果顺利则卵孵化并成活为1龄鱼,其成活率(1龄鱼条数与产卵量n之比)为 1.22×10^11/(1.22×10^11+n); 渔业管理部门规定,每年只允许在产卵期前的8个月内进行捕捞作业。 如果每年投入的捕捞能力(如渔船数、下网次数等)固定不变,这时单位时间捕捞量将与各年龄组鱼群条数成正比,比例系数称作捕捞强度系数。通常使用13mm网眼的拉网,这种网只能捕捞3龄鱼和4龄鱼,两类鱼的捕捞强度系数之比为0.42:1。渔业上称这种方式为固定努力量捕捞。用matlab求解(2) 某渔业公司承包这种鱼的捕捞业务5年,合同要求鱼群的生产能力不能受到太大的破坏。已知承包时各年龄组鱼群的数量分别为:122、29.7、10.1、3.29(×109条),如果仍用固定努力量的捕捞方式,该公司采取怎样的策略才能使总收获量最高?
时间: 2023-06-18 08:05:55 浏览: 69
这是一个典型的渔业资源管理问题,可以用动态系统的方法求解。首先,我们需要建立一个数学模型来描述鱼群的数量和捕捞量之间的关系。
设 $N_1, N_2, N_3, N_4$ 分别表示 1 龄鱼、2 龄鱼、3 龄鱼和 4 龄鱼的鱼群数量,$F_3, F_4$ 分别表示 3 龄鱼和 4 龄鱼的年产卵量,$C_3, C_4$ 分别表示 3 龄鱼和 4 龄鱼的捕捞量,$K_3, K_4$ 分别表示 3 龄鱼和 4 龄鱼的捕捞强度系数,则有如下的动态方程组:
$$
\begin{aligned}
\frac{dN_1}{dt} &= \frac{F_3}{2} \times \frac{1.22 \times 10^{11}}{1.22 \times 10^{11} + N_1} \times N_4 - 0.8 N_1 \\
\frac{dN_2}{dt} &= \frac{F_3}{2} \times \frac{1.22 \times 10^{11}}{1.22 \times 10^{11} + N_2} \times N_4 - 0.8 N_2 \\
\frac{dN_3}{dt} &= F_4 \times \frac{1.22 \times 10^{11}}{1.22 \times 10^{11} + N_3} - 0.8 N_3 - K_3 C_3 \\
\frac{dN_4}{dt} &= F_4 \times \frac{1.22 \times 10^{11}}{1.22 \times 10^{11} + N_4} - 0.8 N_4 - K_4 C_4
\end{aligned}
$$
其中,第一个和第二个方程表示 1 龄鱼和 2 龄鱼的数量变化,由于它们不参与产卵,所以只受到自然死亡率的影响和 4 龄鱼产生的后代数量的影响;第三个和第四个方程表示 3 龄鱼和 4 龄鱼的数量变化,它们参与产卵,同时还会被捕捞,所以除了自然死亡率和产生后代的影响外,还要受到捕捞的影响。
由于每年只允许在产卵期前的 8 个月内进行捕捞作业,我们可以将时间离散化,假设每个时间步长为 1 个月,那么可以用欧拉方法来计算鱼群数量的变化:
$$
N_i(t+1) = N_i(t) + \frac{dN_i}{dt} \Delta t
$$
在每个时间步长内,我们需要计算捕捞量和收获量。由于拉网只能捕捞 3 龄鱼和 4 龄鱼,且捕捞强度系数之比为 0.42:1,因此:
$$
\begin{aligned}
C_3 &= K_3 (0.42 N_3 + N_4) \\
C_4 &= K_4 (0.58 N_3 + N_4)
\end{aligned}
$$
同时,每条鱼的平均重量已知,因此可以计算出每个年龄组的总重量,从而计算出总收获量:
$$
\text{总收获量} = 5.07 N_1 + 11.55 N_2 + 17.86 N_3 + 22.99 N_4 - (C_3 + C_4) \times \text{平均重量}
$$
这样,我们就可以用 MATLAB 等工具求解上述动态方程组,得到每个时间步长内鱼群数量的变化以及最终的总收获量。在实际应用中,我们还需要考虑很多其他因素,比如渔业资源的可持续性、经济效益、社会影响等等,需要进行综合考虑和权衡。