水资源系统分析作业.docx

1用EXCEL规划求解或Matlab优化工具求解下列随机线性规划问题（10分）目标函数：maxE(z)=E(C1).x1+ E(C2).x2约束条件: P(5 x1+4x2b1)0.975P(2 x1+3x2b2)0.985式中, C1、C2、b1、b2均为正态分布的随机变量C1，N（9，32）；C2，N（8，22）；b1，N（30，82）；b2，N（20，72）（要求附规划求解的屏幕拷贝图，或Matlab程序求解的屏幕拷贝图）解：(1) 目标函数：约束条件：在上述模型中，对于机会约束，查正态分布表得到与和对应的和，于是原约束转化为确定性约束：(2) 在MATLAB中求解，问题如下：Obj: Sb.to: 即目标函数的最大值为25.2514，在x1=3.3886，x2=-0.6557时取得。2. 某水源地可供水量为Q,可以分配给3个用户，分配水量xj给用户j时所产生的效益可近似表示为Ej=ajxj2+bjxj+cj，j=1,2,3。如何分配水量才能使总效益最大？列出数学模型，并用Lagrange乘子法求解。如果Q=19.25，a1=-0.5，a2=-0.4，a3=-0.5，b1=7.65，b2=6.40，b3=6.85，c1=1710，c2=1650，c3=1580，求出具体的水量分配方案（15分）解：(1) 以分配水量获得的总效益最大为目标函数，根据题意建立如下数学模型：目标函数： 约束条件：(2) 构造拉格朗日函数：其驻点满足条件：(3) 解得：考虑到至少有一个为0，则存在以下三种情况。 解得：，不符合约束条件，因而舍去。 此时，约束条件不起作用，解得：，也不符合条件，因而也舍去。 解得：。3一个灌区耕地面积AREA =1500hm2，可用灌溉水量W为600万m3。在安排种植计划时，考虑三种粮食作物A，B，C，其灌溉定额分别为4000m3/hm2、4500 m3/hm2，6000 m3/hm2，净收入分别为4500元/hm2、5000元/hm2、6000元/hm2。问如果希望在保证灌区净收入达到480万元的基础上尽可能多的节约灌溉水量，应如何安排三种作物的种植面积？建立多目标规划模型，并用线性目标规划求解（15分）（要求附MATLAB程序或其他程序求解过程的屏幕拷贝图）解：(1) 依据原问题建立多目标规划模型如下：以作物A、B、C的种植面积为决策变量。目标函数：约束条件：(2) 以作物A、B、C的种植面积为决策变量，以表示灌区净收入与480万元之间的正、负偏差，以表示灌溉水量与600万m3之间的正、负偏差。第一个目标要求净收入达到480万元，即要求尽可能小；第二个目标要求节约灌溉水量最多，即要求尽可能大。原多目标规划模型改为线性目标规划模型为：目标函数： 目标约束： 绝对约束： 非负约束： 利用MATLAB求解上述模型，可得：(3) 求解过程：第一步：求解如下模型：运行结果如下：第二步：求解如下模型运行结果如下：最终得到的结果为：即三种作物的种植面积分别为345.951、363.318、496.878 hm2时能够使净收入达到480万元且节水最大，节水为0 m3。4 为寻求某水库的最优运行策略，将每年划分为3个时段，每个时段的入库水量有两个可能的离散值Qit（i=1，2为离散值编号；t=1，2，3为时段编号），根据历史资料分析，各时段的入库水量相互独立，Qit的取值及其概率Pit见表1。每个时段水库蓄水量St的变化范围为25，有效放水量Rt超过3，St和Rt均间隔1进行离散，各阶段不同放水量Rt下的净效益Bt见表1。如果年初年末水库蓄水量均为2，用随机动态规划方法寻求一个最优运行策略（放水策略）。 (注：时段初水库蓄水量St和时段入库水量Qit为状态变量)。（20分）表1 各时段水库入库水量出现的概率及不同放水量下的净效益时段t入库水量Qit相应概率Pit不同放水量Rt下的净效益Bti=1i=2i=1i=2Rt=0Rt=1Rt=2Rt=31120.20.801015172340.30.701525283230.70.30101213解：(1) 阶段变量：，表示水库年运行期的第t个阶段； (2) 决策变量：第t个阶段水库的有效放水量Rt。 (3) 状态变量：阶段初水库蓄水量St和时段入库水量Qit。 (4) 状态转移方程：水库水量平衡方程 (假设没有蒸发渗漏损失) (5) 指标函数：t阶段的指标函数为该阶段的放水净效益Bt。 (6) 目标函数：调度期内的总净效益最大 (7) 约束条件： (8) 边界约束： 采用顺序法进行递推求解，其基本方程为： 表1 阶段1计算结果S1Qi,1Pi不同R1下的B1EB*1R*1对应的S2弃水WS10123210.20101412020.801015220表2 阶段2计算结果S2Qi,2Pi不同R2下的B2EB*2R*2对应的S3弃水WS20123230.30+1415+1425+1428+144232040.70+1415+1425+1428+14330330.301525282833040.70152528340430.301525282834040.70152528350530.301525282835040.70152528351表3 阶段3计算结果S3Qi,3Pi不同R3下的B3EB*3R*3对应的S4弃水WS30123220.70+4210+4212+4254.322030.30+4210+4212+4213+42320320.70+2810+2812+2813+284132030.30+2810+2812+2813+28321420.70+2810+2812+2813+284132130.30+2810+2812+2813+28322520.70+2810+2812+2813+284132230.30+2810+2812+2813+28323表4 水库最优运行策略时段123根据最优决策确定的净效益入库水量Qt13250放水量Rt132入库水量Qt13351放水量Rt133入库水量Qt14251放水量Rt133入库水量Qt14351放水量Rt133入库水量Qt23255放水量Rt232入库水量Qt23356放水量Rt233入库水量Qt24256放水量Rt233入库水量Qt24356放水量Rt2335投资决策问题。某流域管理局设在今后五年内可用于流域投资的资金总额为900万元，有7个可以考虑的投资项目（表2），假定每个项目只能投资一次，第i个项目所需的投资资金为bi亿元，将会获得的利润为ci亿元，且第4个项目和第5个项目2者只能选其中一个，问如何选择投资项目，才能使获得的总利润最大？试列出该问题的数学模型，并求解。（10分）表2 电站的投资及年利润AiA1A2A3A4A5A6A7Ci/万元2500150030002100270023001800bi/万元220110240140210180130解：引入0-1变量，设第i个项目被选状态为，当时，表示投资该项目；当时，表示不投资该项目。(1) 根据已知条件建立模型目标函数： 约束条件： (2) 采用MATLAB求解，求解结果如下：X=1;1;1;1;0;1;0, Z=1.14亿元，即该管理局未来五年投资项目是第1、2、3、4、6个项目，可得到最大的利润，为1.14亿元。程序编码：6人工神经网络建模：已知14组观测值x1、x2、x3、x4及y（表4），利用BP网络，预测第15组观测值x1、x2、x3、x4取值为122.1、65327、56747、1351.64时，y的值。（10分）（要求附程序，求解过程屏幕拷贝图）表3 试验观测结果变量123456789101112131415x187.1115.6110.877.378.979.5115.5107.7202100.113892.6114.994.4122.1x2423265160652982.55435957552.560746581505644563115651897084466418697747690365327x3239263175632422.53308939847.546606459703613550065526995822456238614946941356747x41357.581357.271356.711356.161355.611355.241354.451353.791353.641353.11352.571352.271351.311351.681351.64y1357.271356.711356.161355.611355.241354.451353.791353.641353.11352.571352.271351.311351.681351.64？解：计算结果为：当时，。程序编码：% 输入X=87.1 115.6 110.8 77.3 78.9 79.5 115.5 107.7 202 100.1 138 92.6 114.9 94.4; 42326 51606 52982.5 54359 57552.5 60746 58150 56445 63115 65189 70844 66418 69774 76903; 23926 31756 32422.5 33089 39847.5 46606 45970 36135 50065 52699 58224 56238 61494 69413 56747; 1357.58 1356.71 1356.16 1355.61 1355.24 1354.45 1353.79 1353.64 1353.1 1352.57 1352.27 1351.31 1351.68;% 期望输出值Y=1357.27 1356.71 1356.16 1355.61 1355.24 1354.45 1353.79 1353.64 1353.1 1352.57 1352.27 1351.31 1351.68 1351.64;%建立BP网络，一层隐含层，隐层神经元数为3，输出为1个单元，训练函数为traingdmnet = newff(minmax(X),3 1,tansig,purelin,traingdm);%设置输入层权值和阈值inputWeights=net.IW1,1;ingputbias=net.b2;%设置训练参数 net.trainParam.lr=0.55; %学习率 net.trainParam.epochs = 6000;%最大训练次数 net.trainParam.goal = 1e-7; %目标误差 net=init(net);%重新初始化%训练网络net = train(net,X,Y);%仿真y = sim(net,X); %将测试数据输入网络进行测试E=Y-y; %计算测试集网络输出和目标的误差mse=MSE(E) %计算均方误差%对得出的网络进行测试X1=122.1;65327;56747;1351.64; y1=sim(net,X1) %用sim仿真7.论述水资源系统分析的一个新理论或新方法（引进时间、方法介绍及应用情况）。（20分）答：(1) 对策论（博弈论）：博弈论，是解决竞争者应该采取何种对策的理论和方法。如果对抗双方可能采取的对策只有有限个，则是有限博弈论；如果可能采取的对策为无限个，则是无限博弈；如果在对抗中获胜的一方和失败的一方得失恰好相等，则是零和对弈。目前博弈论在水资源系统分析中主要应用于水资源配置、解决水资源冲突等方面。(2) 模糊决策方法：以模糊数学为基础发展起来的系统分析方法，是对具有模糊性质的问题提供决策依据的方法，属于不确定数学方法的范畴。模糊决策方法包括隶属度确定方法、模糊聚类分析、模糊数学规划等。模糊决策理论应用于水质模糊综合评价、环境评价、水资源合理配置研究、水资源承载力分析及水资源效益评价体系。(3) 人工神经网络（ANN）人工神经网络是模拟人脑神经网络结构与功能的一种技术系统，用大量的非线性并行处理单元（人工神经元）模拟人脑神经元，用处理器之间错综灵活的连接关系来模拟人工神经元间的突触行为，直接使用样本数据来建立输入与输出间的非线性映射关系。人工神经网络目前在水资源领域也得到了广泛的应用，主要应用于系统建模、预测、模式识别与评价、系统优化等方面，例如，对降水径流等的预测、水文分类预报、水质评价及水资源优化等等。(4) 遗传算法(GA)遗传算法是模拟生物界的遗传和进化过程而建立起来的一种搜索算法，体现着“生存竞争、优胜劣汰、适者生存”的竞争机制，具有高效的随机搜索与全局优化的特点，适于优化问题的求解。目前遗传算法可应用于水文模型参数的优选、水文地质参数的反演、中长期水文预报及水库调度、实时洪水预报、跨流域引水工程优化调度、水资源优化配置、灾情评估等方面。(5) 分形与混沌方法自然界中所有的形状大致可分为两类，一类是具有特征长度的图形，另一类是不具有特征长度的图形，对于后一类图形，尽管其没有特征长度，但是却具有一定的内在分布规律。部分与整体以某种方式相似的形体成为分形，研究分形及其应用的科学成为分形理论。自相似原则和迭代生成原则是分形理论的重要原则。它表征分形在通常的几何变换下具有不变性。目前，分形理论在水文水资源中的应用包括以下几个方面：水系河网结构和流域地形地貌及其演变；河床表面形态及其演变；降雨时空分布；径流过程的分形特征；暴雨时空分布；洪水时空分布；以及产汇流模型中的尺度问题等。(6) 小波分析小波分析是当前应用数学和工程学科中一个迅速