最佳捕鱼方案。
夏天，郑巨耀，王芽。
摘要：本文解决的是一个最佳捕鱼方案设计的单目标线性规划问题，目的是制定每天的捕鱼策略，使得总收益最大。根据题设条件，结合实际情况，我们设计了成本与损失率随天数的增加成反比变化的函数曲线（见图三所示），并导出总收益的表达式：
。由于**是关于**量的分段函数（见图一所示），我们引入“0－1”变量法编写程序（程序见附录一），并用数学软件lingo求解，得到最大收益(w)为441291.4元，分21天捕捞完毕。其中第1～16天，日捕捞量在1030～1070公斤之间，第17～21天的日捕捞量为1610～1670公斤之间（具体数值见正文）。
由结果分析，我们对模型提出了优化方向，例如人工放水来降低成本。
关键词：“0-1”整数规划，单目标线性规划，离散型分布。
一． 问题重述。
一个水库，由个人承包，为了提高经济效益，保证优质鱼类有良好的生活环境，必须对水库里的杂鱼做一次彻底清理，因此放水清库。水库现有水位平均为15米，自然放水每天水位降低0.5米，经与当地协商水库水位最低降至5米，这样预计需要二十天时间，水位可达到目标。
据估计水库内尚有草鱼二万五千余公斤，鲜活草鱼在当地市场上，若日**量在500公斤以下，其**为30元/公斤；日**量在500—1000公斤，其**降至25元/公斤，日**量超过1000公斤时，**降至20元/公斤以下，日**量到1500公斤处于饱和。捕捞草鱼的成本水位于15米时，每公斤6元；当水位降至5米时，为3元/公斤。同时随着水位的下降草鱼死亡和捕捞造成损失增加，至最低水位5米时损失率为10%。
承包人提出了这样一个问题：如何捕捞鲜活草鱼投放市场，效益最佳？
二． 模型假设。
1． 池塘中草鱼的生长处于稳定状态，不考虑种群繁殖以及其体重增减，即在捕捞过程中草鱼总量保持在25，000公斤不变。
2． 第一天捕捞时水位为15m，每天都在当天的初始水位捕捞草鱼，水库水位每天按自然放水0.5m逐渐降低，20天后刚好达到最低要求水位5m。
3． 在水库自然放水的21内将草鱼捕完。
4． 在草鱼日**量未达饱和的之前，市场**量等于销售量。
5． 每天草鱼的捕捞成本随着每天水位的降低呈等差数列递增分布。
6． 随着水库水位的下降，草鱼的种群密度逐渐变大，存在着对空间、食物、氧气的竞争，种群死亡率逐渐升高。题设中给定草鱼死亡及捕捞损失率随着水位的降低而升高，在这里我们假设草鱼损失率是一个统计学概念，即已经综合了因自然死亡和捕捞等其他原因共同造成的损失。
7． 草鱼损失率与水库水位成反比关系，每天捕捞量的损失率与当天池塘总鱼量的损失率是一致的，以每次捕捞时池塘总鱼数为当次基数。
8． 捕捞上的草鱼中的死鱼将另行处理，不会放回水库也不会与活鱼一起**。
9． 日**量在1000---1500公斤时，我们假定草鱼**为20元每公斤这一常数。总体**随**量变化关系，如图五所示：
图一。三． 问题的分析。
1. 在符合题意并且与实际情况较吻合的情况下，我们应寻求对最优解的精确求解以及依据草鱼捕捞的可行性方案来捕捞使得承包人获益(w)最大。
2. 我们在追求收益最大的同时，需要求出草鱼捕捞的天数以及每天的捕鱼量，这是一个单目标线性规划问题，原题中给定的草鱼日**量不同的情况下草鱼的单价也不一样，这样每天的草鱼****均取决于当日的草鱼**量，于是在模型求解过程中我们采用“0-1整数规划”来解决这个问题，并运用数学软件lingo来求解，最后对所得的解进行讨论和分析。
四． 模型的建立及求解。
1. 符号的说明：
—第i天水库水位（米）； 第i天**量（公斤）；
—第i天草鱼的损失率第i天草鱼捕捞成本(元/公斤）；
—第i天草鱼销售收益（元)； 捕捞期内草鱼销售总收益（元）；
—第i天的捕捞量（公斤）； 第i天的售价(元）；
2. 模型的建立。
根据假设5，随着水位自然地下降，草鱼的捕捞成本呈等差数列递减分布，第一天捕捞时水位仍然维持在15m，共需21天。故每公斤草鱼捕捞成本为：
根据假设7，损失率与水位成反比：
第i天草鱼的损失率为（如图三散点图所示）：
图二。图三。
第i天的捕捞量与**量之间的关系式是：
收益=销售额-成本。
即我们的目标函数：
又根据已知条件可得以下5个约束条件：
a) 总的捕鱼量不大于总鱼量。
b) 由日**量到1500公斤达到饱和则。
c) 再根据题意的三个在不同的**量之间的**的不同可得到如下关系式：
3. 模型的求解。
根据上述解题思路，我们用数学软件lingo对模型进行求解，根据模型的约束条件和软件的特点，我们采用引入0－1变量的方法编写解题程序：
—0,1变量引入，使得 .
上述式子需改动为：
约束条件：由程序运行结果得到，21天内草鱼捕捞总收益w=441291.4元。
每天草鱼的具体捕捞策略见表一：
表一。五． 结果分析和检验。
1. （公斤），且大部分捕捞量集中在1000－1100 之间，说明每天的捕捞量基本相对稳定。
图四。2. 由附录一得，除了外，其他变量的reduced cost全部为零，即变量的微小变化对结果没有影响，而只能取0，1整数，本身没有微小变化，说明我们的模型稳健性比较好。
3. 根据捕捞量和损失率我们可以得到销售量在前16天为1000公斤，而后5天为1500公斤（见表二），而这两个点刚好是**曲线中的两个间断点的右端点，也就是在同样**下的选取最大允许销售量，从定性角度分析，这样的结果使剩余草鱼种群总量降低，相应地后续天数的草鱼损失量会减少，这样也符合一般的逻辑，比较容易接受。
表二。4. 在最大利润的方案下，每天的捕捞量持续缓慢上升，但第17天出现一个跳跃点。
因为前15天水位高度还能够基本满足鱼群正常生长的需求，适当控制捕捞量以获得较优的售价，而后几天水位较低，鱼群因空间、食物、生长竞争激烈，损失率会急剧升高，因此要加大捕鱼量，适当降低鱼群密度，来减少鱼群基数便于降低鱼群损失。这也可以看出前16天我们捕鱼的决定因素是成本，而后5天的主要影响因素是损失率，且从总体来看，前16的总收益大于后5天的总收益。因此我们认为成本较损失率对收益的影响较大。
六． 优化方向。
1. 上述模型中**相对于**量是分段常函数，这是一种较为简化和粗糙的假设。商品**应围绕价值上下波动，**对**量应比原假设更敏感，如设为连续函数可能更接近市场实际情况，再由于市场本身的特点控制着**的上限我们由计算机模拟出**的定性曲线（如图五所示）。
图五。2. 随着水位下降，损失率增加，模型假设中我们设它们为反比例函数，由于水位的下降，草鱼的密度增加，损失的速率不断加快，在损失率的曲线上表示为斜率的加快，但实际中可能更为复杂，每个生态系统具体情况也不同，。
3. 有结果分析可得，前15天我们捕鱼的决定因素是成本，而后5天的主要影响因素是损失率。因此我们考虑在钱其降低捕鱼成本，可以选择人工加快放水来代替自然放水方式以达到降低成本的途径。
例如，可以在某一天使水位急剧降低到某一个定值l,然后再让水位再自然放水状态下捕捞，抑或者在每天放水量不一的情况下捕捞，通过人工调节水位，都有可能时的最后的收益变得更大。
4. 在第9条假设中，在草鱼**量在1000～1500公斤时，我们将草鱼售价取定位20元/公斤。现在做一下修正，由于草鱼**不会有大幅度降低，大体应该在15～20元/公斤之间分布，于是我们另外分别在15，16，17，18，19元/公斤的售价下进行最优化求解，得到结果列成**，如表三所示：
表三。由表三可以看出，除了单价在16元/公斤时，捕鱼策略稍有变动，第二天捕捞量s2=500公斤，第15天s15=1500公斤之外，其他天数的捕捞量都与单价在20元/公斤时的捕鱼量相同。这从另一个侧面说明我们的第9条假设是合理的。
参考文献：1]谢金星薛毅编著 《优化建模与lindo/lingo软件》清华大学出版社 2023年7月第一版第八章目标规划模型。
2]卢向南，李俊杰寿涌毅编著 《应用运筹学》浙江大学出版社 2023年2月第一版第四章整数规划第三节整数规划的应用。
附录一：程序及结果。
源程序：model:
title 最优捕捞模型；
sets:day/1..21/:r,m,l,q1,q2,q3,y1,y2,y3;
price/1,2,3/:p;
endsets
data:p=30,25,20;
l=15,14.5,14,13.5,13,12.
5,12,11.5,11,10.5,10,9.
5,9,8.5,8,7.5,7,6.
5,6,5.5,5;
m=6,5.85,5.7,5.
55,5.4,5.25,5.
1,4.95,4.8,4.
65,4.5,4.35,4.
2,4.05,3.9,3.
75,3.6 3.45,3.
3,3.15,3;
数学建模案例分析
线性代数建模案例汇编。法正系，思想政治教育13 1 汗克孜 亚森。2015年6月。这里收集了二十个容易理解的案例。和各类数学建模竞赛的题目相比，这些案例确实显得过于简单。但如果学生能通过这些案例加深对线性代数基本概念 理论和方法的理解，培养数学建模的意识，那么我们初步的目的也就达到了。城市道路网中每...
小学数学建模案例
运用，j学数学建模解决此类问题时，要充分发挥学生的自主性，教师需要一步一步地引导学生建立数学模型。解决牛吃草问题的数学模型如下 假定一头牛一天吃草量为 1 草的生长速度 对应的牛头数 吃的较多天数一相应的牛头数x吃的较少天数 原有草量 牛头数x吃的天数一草的生长速度x吃的天数 吃的天数 原有草量 牛...
数学建模案例 上
幻灯片1 经济数学。第7章数学建模案例。7 1 数学建模概述。某饮料生产企业现在需要设计一批容积为v的圆。柱形饮料包装盒，问应怎样设计才能使所用材料最省？讨论 什么是最优设计。易拉罐的形状如何。材料跟易拉罐的什么有关。幻灯片2经济数学。第7章数学建模案例。7 1 数学建模概述。在生活中我们会发现销量...