MATLAB在水资源动态规划中运用与改进.docVIP

MATLAB在水资源动态规划中运用与改进 　　【摘 要】水资源是人类赖以生存的宝贵资源，搞好水资源优化配置对实现水资源的可持续利用尤为重要。本文结合一个实例，详细介绍了用MATLAB解决水资源动态规划的思路与方法，同时与传统的逆序法对比，突出了MATLAB的优越性。并且，参照了之前学者所编写的程序，给出了自己的程序优化方案，使程序运行速度、可读性等方面得到了提高。 　　【关键词】动态规划；水资源优化配置；MATLAB 　　0.引言 　　动态规划是解决多阶段决策过程最优化问题的一种方法，该方法是由美国数学家贝尔曼（R.Bellman）等人在20世纪50年代提出。他们针对多阶段决策问题的特点，提出了解决这类问题的最优化原理，并成功解决了生产管理、资源分配等方面的许多实际问题，从而建立了运筹学的分支——动态规划[1]。MATLAB是一个功能强大的用于矩阵运算的数值计算软件，是决策系统优化计算和设计的有力工具，将MATLAB运用到动态规划中可以达到计算简便的目的[2]。现在利用MATLAB解决动态规划问题的程序不少，但是，如何充分利用MATLAB的特点来优化程序确实一个值得深思的问题。 　　1.MATLAB程序设计 　　MATLAB有2种工作方式：一是交互式的命令工作方式，直接在Command Window中输入指令求解。另一种是M脚本文件的工作方式，这种工作方式用来解决指令较多和所用指令结构较为复杂的问题[3]。 　　由于动态规划问题的特殊性，其涉及的实际问题均需要反复的计算求解，所以在使用MATLAB语言编程中采用最多的是循环结构。MATLAB提供了2种实现循环结构的语句——for语句和while语句[4]。 　　2.应用实例 　　有一引水渠，其设计最大流量为6m3/s，为甲、乙、丙、丁4个地区供水。据统计，在给四个地区供水时，不同的供水量产生的效益见表1（效益单位为万元）。那么，如何分配供水量可以使产生的总效益最大。（本实例取自文献[5]） 　　表1 供水量与增产效益的关系 　　2.1逆序法计算 　　采用逆序法表格计算，可以得到最优分配方案有3个，分别是（1，2，1，2），（1，2，2，1），（2，2，1，1），最大效益为19万元。 　　2.2 MATLAB编程计算 　　2.2.1思路分析 　　由于本例较???简单，所以不采用M脚本文件的工作方式，直接在Command Window中输入指令即可。 　　Step1.输入参数 　　建立A、B、C、D 4个数组，分别来储存供给甲、乙、丙、丁4个地区不同水量时产生的效益。考虑到在动态规划时，某个城市的供水量可以为0，对应的效益也为0，所以A、B、C、D四个数组都应该加入一个元素0。这里采用单下标表示法，分别用i、j、k、 　　l表示数组A、B、C、D中元素的下标。A（3）=2，表示当i=3时，A中元素为2，即向A供水2m3/s。Step2.找出所有可行解，储存其分配方法和对应的效益，要找出所有的可行解，最简便的的方法是采用循环结构枚举。 　　①确定循环条件： 　　循环的限制条件应该是最大设计流量。由于最大设计流量是6 m3/s，现在任意假设一例分配方法以确定循环条件。假设所有的水量全部供给了丁，即l=7，那么此时可有i=j=k=1，得i+j+k+l=10。由于无论是那种分配方法都要满足最大设计流量的限制，所以i+j+k+l=10即为循环条件。由于循环次数是已知的，所以采用for语句循环。 　　②储存分配方法： 　　建立一个二维数组E，采用全下标的方法将分配方法存入其中，每一行为一种分配方式。 　　③储存效益： 　　建立一个数组d，常用单下标的方式将所有分配方法对应的效益存入其中。 　　Step3.确定最大效益的值 　　所有方案的效益都存储在数组d中，所以传统的方法就是采用循环结构遍历数组d中的所有元素以找到最大值。但是在编程中应该尽量避免使用循环结构。因为在循环结构中，要反复进行存储变量间的“存放”操作和算符调用操作，消耗计算时间。MATLAB中为用户提供了大量的现成函数，尽可能多的采用现成函数编程可以使所编程序更可靠、更快速、可读性更好[6]。 　　①采用循环结构确定最大效益 　　程序如下： 　　MAX=d（1）； 　　for n=1：length（d） 　　if d（n）MAX 　　MAX=d（n）； 　　end 　　end 　　②采用max函数确定最大效益 　　由于max（X）求的是数组X中各列的最大值，所以应该先对数组d转置。 　　程序如下： 　　MAX=max（d） 　　可以清楚的看到，采用matlab中现成的max函数比用循环结构更好，不仅使程序的输入得到了简化，同时也使程序的可读性更好。