中国传媒大学《MATLAB基础与应用》课件-第13-14章数值计算方法实际应用案例.pptVIP

将t代入 即可求出则得本问题的一个数学模型： 记 是由 按 构成的向量 example_13_3_1.m 考虑到非线性规划求解对初值依赖性较大，因此我们可在零点处随机生成 若干个初值的计算来获取可能的最优解。 example_13_3_2.m 13.4 最优捕鱼策略 为了求解鱼量稳定情况下的最大捕捞量，我们把捕鱼能力系数k作为一个关键的控制变量，通过建立各相关量与k的关系，以求能最终根据稳定情况下的最大捕捞量条件求得一最佳的k值，并由此得出最优的可持续捕捞的渔业捕获模型。 一般情况下，各龄鱼的数量须经过一段时间才能达到一个稳定的状态，即到平衡时，年末和年初的数量基本不变。 设α为各年龄鱼的自然死亡率(表示单位时间内死亡的鱼的数量与鱼的总量之比)，在 的时间段内，根据死亡率的定义，则有 无捕捞时，鱼群数量s随时间变化的规律 下面考虑有捕捞情况下的情形： 设年初各龄鱼的鱼群数量分别为s10,s20,s30,s40,八月末，经过捕捞及自然死亡后的各龄鱼群数量 ，十二月末，各龄鱼群数量s1,s2,s3,s4，卵的总数量为n，则由前面的式子，t按年计算，有 1-8月为捕捞季节，经过捕捞及自然死亡，八月末，各龄鱼群数量为： 9-12月为产卵季节，此期间不捕捞，则十二月末，各龄鱼群数量为： 再设 分别为3、4龄鱼在产卵期的平均数量，n为3、4龄鱼产卵数量的总和，t按月计算 有捕捞时，鱼群数量s随时间变化的规律 产卵期产卵总量 ，其中a为平均每条4龄鱼产卵个数。设s11,s21,s31,s41表示第2年各龄鱼的数量初值，则有 1、1龄鱼由卵孵化并成活下来的那部分卵子转化而成，即 2、2龄鱼由上一年龄鱼转化而成，即 3、3龄鱼即上一年末2龄鱼 4、4龄鱼即上一年末3龄鱼 对于一种可捕获的鱼来说，设其捕获量为P，则P可表示为 1-8月，捕捞3、4龄鱼的数量 fish_tactics.m 13.5 零件的参数设计 可以假设7个零件参数服从正态分布。根据标定值和容差的定义有： 编写如下程序： mu=[0.1,0.3,0.1,0.1,1.5,16,0.75]; sigma=[0.005/3,0.005,0.005/3,0.005,0.225/3,0.8/3,0.0375/3]; for k=1:7 x(:,k)=normrnd(mu(k),sigma(k),1000,1); end t=(1-2.62*(1-0.36*(x(:,4)./x(:,2)).^(-0.56)).^1.5... .*(x(:,4)./x(:,2)).^1.16)./x(:,6)./x(:,7); y=(x(:,1)./x(:,5)).*(x(:,3)./(x(:,2)-x(:,1))).^0.85; y=174.42*y.*t.^0.5; d=abs(y-1.5); f=sum(9000*(d0.3)+1000*(d=0.3).*(d0.1))/1000 运行结果： 平均损失f =2778 13.6 循环赛排名问题 如何考虑对手的强弱，首先需要回答一个问题 —— 用什么指标刻画对手的强弱？从直观上讲，所谓对手强，就是对手战胜的队多，对手弱就是对手战胜的队少，那么将对手的得分也记录到自己队的得分中，即所谓的二级得分。这样考虑问题要比直接考虑得分更合理，这里计算一下二级得分 按这种计算方法，球队3排名第一，球队一排名第二，…。 这看起来是合理了，但所说的问题并没有根本解决。因为这里只考虑每个队它的对手得分，而没有考虑每个队它的对手的对手得分，也就是说，评判对手强弱的计算公式不合理。为了合理的解决这个问题，需要考虑每个队对手的对手的得分，即三级得分。那么同样的问题，还有四级得分、五级得分……。看一下各队各级的得分。 看来各队名次的排列有点波动，例如球队3和球队1竞争第一名的情况就是这样。如何确定他们的名次呢？再对各队的得分进行分析。将各队比赛的得分情况构造矩阵 同时构造列向量 ，它可以看成各队的实力，比赛前，认为6个队是相同的。那么第一级得分为 s1也可以被看成各队的实力，比赛后就不同了。而二级以后的得分是矩阵与向量的乘积，即 最合理的情况应考虑极限情况，即 式中， 是矩阵A的最大特征值，s是对应于 的特征向量（可以证明：在满足一定的条件下，上述极限是存在的）。 因此，名次排列问题就转化为求得分矩阵的最大特征值和对应的特征向量。对于前面的问题，其最大特征值为 ，相