【优化调度】基于协方差矩阵自适应贪心搜索实现配水系统优化附matlab代码+仿真结果和运行方法.zip

# Covariance Matrix Adaptation Greedy Search Applied to Water Distribution System Optimization # 01/08/2019 ================================================================ The Matlab hybrid framework consists of three componenets including 1. Covariance Matrix Adaptation Evolution Strategy (CMA-ES): Applied for minimizing the cost of the pipes. The recommended version of CMA-ES : https://www.lri.fr/~hansen/cmaes.m 2. UpWard Greedy Search: Applied for handling the violations of constraints 3. DownWard Greedy Search: Applied for reducing the extra cost of pipes. ----------------------------------------------------------------------- Please pay attention to some essential points : a. main.m is the main file for running this framework, but before running, please check the below code line in your Matlab prompt mex - setup if the below error message is shown: please follow the recommended steps: ====================================================================== # Functions # |Function|Description| |---------|---------| |set_Network_name|initializing the network name | |set_dim|initializing the diameters of the applied pipes| |set_fitfun|setting the cost function name | |Def_Iter | used for an iterative network like DNYTP and 50NYTP| |Def_Constraint | determining the Nodal pressure head constraints| |set_Bestfmin | estimating of the near-optimum network cost design | |initialize_xstart | initializing the first population| |Def_TolFun | defining the tolerance of the fitness function as a stopping criterian| |set_Max_Eval | initializing the maximum evaluation number| |cmaes | CMAES function| |Finding_best | recording the best continuous, discrete and possible pipe configurations | |Downward_Greedy_Search| applying the strategy of downward greedy search | |Upward_Greedy_Search | applying the strategy of upward greedy search| |rep | counting the same pipe configurations| |Dim_Correction | converting the continuous designs to the nearest possible ones| # The structure of recorded results # |Variables |Description| |---------|---------| |Pipe| The main structure| |UGSflag| a flag (0/1) shows the upward greedy search is run or not | |DGSflag| a flag (0/1) shows the downward greedy search is run or not| |Continuous_Sols| a matrix of all continuous pipe designs (design number * network pipe number)| |Continuous_Cost| an array of the continuous pipe costs | |Continuous_Violation| an array of total violations of all continuous pipes| |Rounded_Sols| a matrix of all possible pipe designs (design number * network pipe number)| |Rounded_Cost| an array of the possible pipe costs | |Rounded_Violation| an array of total violations of all possible pipes| |UGS_Sols| a matrix of all possible pipe designs after applying the upward greedy search (design number * network pipe number)| |UGS_Cost| an array of the possible pipe costs after applying the upward greedy search| |XminU| a variable shows the best pipe design of upward greedy search| |FminU | a variable shows the best pipe cost of upward greedy search| |XminC | a variable shows the best continuous pipe design | |FminC | a variable shows the best continuous pipe cost | |XminD | a variable shows the best discrete pipe design | |FminD | a variable shows the best discrete pipe cost| |UGS_Sols| a matrix of all possible pipe designs after applying the downward greedy search (design number * network pipe number)| |UGS_Cost | an array of the possible pipe costs after applying the downward greedy search | In the following there are the details of EPANET Matlab version. EPANET-MATLAB-Toolkit ================================== The `EPANET-Matlab Toolkit` is an open-source software, originally developed by the [KIOS Research Center for Intelligent Systems and Networks of the University of Cyprus](http://www.kios.ucy.ac.cy/) which operates within the Matlab environment, for providing a programming interface for the latest version of [EPANET](https://github.com/OpenWaterAnalytics/epanet), a hydraulic and quality modeling software created by the US EPA, with Matlab, a high-level technical computing software. The goal of the EPANET Matlab Toolkit is to serve as a common programming framework for research and development in the growing field of smart water networks. The `EPANET-Matlab Toolkit` features easy to use commands/wrappers for viewing, modifying, simulating and plotting results produced by the EPANET libraries. For support, please use the OWA community forum : http://community.wateranalytics.org/ # How to use the Toolkit # Download the folder `Download ZIP`, set the run path in Matlab within the saved folder, and run `RunTests.m`. This will execute all the commands which have been implemented in the Class. Example: d=epanet('Net1.inp') d.getNodeCount d.getNodeElevations # How to fix/report bugs # To fix a bug `Fork` the `EPANET-Matlab Toolkit`, `Edit` the code and make the appropriate change, and then `Pull` it so that we evaluate it. Keep in mind that some bugs may exist in the `EPANET` libraries, in case you are not receiving the expected results. # List of Matlab Class Functions # |Function|Description| |---------|---------| |epanet| Load Input file and open the EPANET Toolkit system| |unload|Unload library and close the EPANET Toolkit system| |loadEPANETFile| Open the EPANET Toolkit system| |getError| Returns the description of an error code| |getComputedHydraulicTimeSeries|Computed Hydraulic Time Series| |getComputedQualityTimeSeries|Computed Quality Time Series| |getComputedTimeSeries|Computed Hydraulic & Quality Time Series via ENepanet binary file| |getConnectivityMatrix|Return connectivity matrix of the network| |getControlRulesCount|Retrieves the number of control rules| |getControls|Retrieves the controls| |getCurveCount|Retrieves the number of curves| |getCurveIndex|Retrieves index of curve with specific ID| |getCurveLengths|Retrieves number of points in a curve| |getCurveNameID|Retrieves curve id| |getCurveValue|Retrieves (x,y) values of specific curve index| |getCurveXY|Retrieves (x,y) values of all curves| |getENfunctionsImpemented|Retrieves the epanet functions that have been developed| |getFlowUnits|Retrieves the units used to express all flow rates| |getLibFunctions|Retrieves the functions of DLL| |getLinkBulkReactionCoeff|Retrieves the value of all link bulk reaction coefficients| |getLinkCount|Retrieves the number of links| |getLinkDiameter|Retrieves the value of all link diameters| |getLinkEnergy|Retrieves the value of all computed energy in kwatts| |getLinkEfficiency|Retrieves the value of all computed efficiency| |getLinkFlows|Retrieves the value of all computed link flow rates| |getLinkHeadloss|Retrieves the value of all computed link headloss| |getLinkIndex|Retrieves the indices of all links, or the indices of an ID set of links| |getLinkInitialSetting|Retrieves the value of all link roughness for pipes or initial speed for pumps or initial setting for valves| |getLinkInitialStatus|Retrieves the value of all link initial status| |getLinkLength|Retrieves the value of all link lengths| |getLinkMinorLossCoeff|Retrieves the value of all link minor loss coefficients| |getLinkNameID|Retrieves the ID label(s) of all links, or the IDs of an index set of links| |getLinkNodesIndex|Retrieves the indexes of the from/to nodes of all links| |getLinkPipeCount|Retrieves the number of pipes| |getLinkPipeIndex|Retrieves the indices of pipes| |getLinkPipeNameID|Retrieves the pipe IDs| |getLinkPumpCount|Retrieves the number of pumps| |getLinkPumpHeadCurveIndex|Retrieves index of a head curve for specif