動(dòng)態(tài)規(guī)劃進(jìn)行單庫(kù)優(yōu)化調(diào)度計(jì)算C++程序代碼

1、#ifndef SRDP_H#define SRDP_H class SRDPprivate:/1.階段變量 stage variablelong m_StageNum;double m_Delta_t;double m_DeltaHour;/2.決策變量 decision variable/double *m_OptDeciVar_q;/3.狀態(tài)變量 state variablelong *m_StateNum;/double *m_StateV;/4.指標(biāo)函數(shù) target functiondouble *m_TarFun;/5.目標(biāo)函數(shù)object functiondouble m_Ob

2、jFun;最優(yōu)軌跡記錄數(shù)組long *m_OptArray;原始資料記錄數(shù)組double *m_RunoffArray;double *m_EvapArray;時(shí)間（年月）記錄數(shù)組long *m_YearArray;long *m_MonthArray;計(jì)算結(jié)果double *m_OptVolumnArray;double *m_OptNArray;double *m_OptAbandonQ;階段數(shù)存儲(chǔ)量時(shí)間間隔T=2626560s一個(gè)月的小時(shí)數(shù)：h庫(kù)容離散點(diǎn)數(shù)庫(kù)容離散階段指標(biāo)函數(shù)/目標(biāo)函數(shù)長(zhǎng)系列徑流蒸發(fā)量/最優(yōu)軌跡對(duì)應(yīng)庫(kù)容數(shù)據(jù)數(shù)組序列/最優(yōu)軌跡對(duì)應(yīng)出力數(shù)據(jù)數(shù)組序列/最優(yōu)軌跡對(duì)應(yīng)的棄水量序列d

3、ouble *m_OptUpWaterLevelArray; 最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)上游水位序列 double *m_OptDownWaterLevelArray;/最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)下游水位序列 double *m_OptVarialbe_q;最優(yōu)軌跡對(duì)應(yīng)的發(fā)電流量序列最優(yōu)軌跡對(duì)應(yīng)的預(yù)想出力最優(yōu)軌跡中低于保證出力時(shí)的破壞深度/最優(yōu)軌跡對(duì)應(yīng)的保證率double *m_OptMeanWaterHeadArray; /最優(yōu)軌跡對(duì)應(yīng)的平均發(fā)電水頭 double *m_OptNmaxArray;double *m_OptDamageDepth;double m_OptReliability;double

4、m_OptGenerationProduction; 發(fā)電量double m_OptAnnualMaxPower;/Parameters of Hydropower Station/死水位m/死庫(kù)容.03億m3正常蓄水位835m汛期防洪限制水位m汛期最大庫(kù)容.03億m3/double m_DeadWaterLevel;double m_DeadVolumn;double m_NormalWaterLevel;double m_FloodLimitLevel;double m_FloodVmax;非汛期最大庫(kù)容.37億m3機(jī)組臺(tái)數(shù)裝機(jī)容量*2MW/保證出力.2MWdouble m_NFloodV

5、max;int m_UnitsNum;double m_InstallCapacity;double m_FirmCapacity_N;預(yù)想出力double m_PotentialWaterLevel1;double m_PotentialCapacity_N1;double m_PotentialWaterLevel2;double m_PotentialCapacity_N2;double m_DesignAnnualGeneration;設(shè)計(jì)年發(fā)電量.99xkWe V);上游水位計(jì)算double DownWaterLevelCalculation(double q); 下游水位計(jì)算voi

6、d ReverseSearchCalculation();反向?qū)?yōu)計(jì)算void TerminalOptionalResultCalculation();/最終結(jié)果計(jì)算;#endif類的實(shí)現(xiàn)/ SRDP.cpp: implementation of the SRDP class./#include stdafx.h#include SRDP.h#include #include assert.h#include fstream#include iomanip#include iostream#include using namespace std;/ Construction/Destructi

7、on/SRDP:SRDP() 原始資料長(zhǎng)系列徑流數(shù)據(jù)/蒸發(fā)數(shù)據(jù)this - m_RunoffArray = NULL;this-m_EvapArray = NULL;/動(dòng)態(tài)規(guī)劃模型數(shù)據(jù)階段數(shù)時(shí)間間隔一個(gè)月的小時(shí)數(shù)/狀態(tài)數(shù)階段指標(biāo)/目標(biāo)函數(shù)時(shí)間（年）數(shù)組時(shí)間（月）數(shù)組/最優(yōu)軌跡記錄/最優(yōu)軌跡對(duì)應(yīng)的預(yù)想出力/最優(yōu)軌跡中低于保證出力時(shí)的破壞深度/this-m_StageNum=100;/this-m_Delta_t=2626560;/this-m_DeltaHour=730;this-m_StateNum = NULL;this-m_TarFun = NULL;/this-m_ObjFun=0.00

8、;this-m_YearArray = NULL;this - m_MonthArray = NULL;this-m_OptArray = NULL;this - m_OptNmaxArray = NULL;this-m_OptDamageDepth = NULL;SRDP:SRDP()原始初始資料數(shù)據(jù)析構(gòu)delete this- m_RunoffArray;delete this-m_EvapArray;中間變量析構(gòu)delete this-m_StateNum;delete this-m_YearArray;delete this-m_MonthArray;計(jì)算結(jié)果數(shù)組析構(gòu)delete th

9、is-m_OptVolumnArray;delete this-m_OptNArray;delete this-m_OptAbandonQ;delete this-m_OptUpWaterLevelArray;/最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)上游水位序列delete this-m_OptDownWaterLevelArray; 最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)上游水位序列 delete this-m_OptVarialbe_q;/最優(yōu)軌跡對(duì)應(yīng)的發(fā)電流量序列/最優(yōu)軌跡對(duì)應(yīng)的預(yù)想出力/最優(yōu)軌跡中低于保證出力時(shí)的破壞深度delete this-m_OptMeanWaterHeadArray;delete this-m_Opt

10、NmaxArray;delete this-m_OptDamageDepth;for(int i=0;im_StageNum;i+)/delete this-m_OptDeciVar_qi;/delete this-m_StateVi;delete this-m_OptArrayi;delete this-m_TarFuni;/delete this-m_OptDeciVar_q;/delete this-m_StateV;delete this-m_OptArray;delete this-m_TarFun;/數(shù)據(jù)初始化內(nèi)存分配/void SRDP:InitionalizingSRDP(ch

11、ar *DataFile1, char *DataFile2,char *DataFile3)this-m_StageNum = 516;this-m_BeginMonth=6;this-m_BeginYear=1957;this-m_OptGenerationProduction=0.00;this-m_OptAnnualMaxPower = 0.00;/this-m_DeltaHour=730;/獲取時(shí)間序列this-m_YearArray = new longthis-m_StageNum+1;this- m_MonthArray = new longthis-m_StageNum+1;

12、this-m_StateNum = new longthis-m_StageNum+1;for(int i=0;im_StageNum;i+)this-m_YearArrayi = 0;this-m_MonthArrayi = 0;this-m_StateNumi = 10;this-m_YearArray0=this-m_BeginYear;/水文年起始年份this-m_MonthArray0=this-m_BeginMonth; /水文年起始月份for(int i = 1;im_StageNum;i+ )if(this-m_MonthArrayi-1+1)m_MonthArrayi = t

13、his-m_MonthArrayi-1+1;this-m_YearArrayi = this-m_YearArrayi-1;else if(this-m_MonthArrayi-1+1)12)this-m_MonthArrayi = this-m_MonthArrayi-1-11;this-m_YearArrayi = this-m_YearArrayi-1+1;/coutvvivv m_YearArrayi-m_MonthArrayiendl;/獲取狀態(tài)離散點(diǎn)數(shù)/for(int i=0;im_StageNum;i+)if(this-m_MonthArrayi=6&this-m_MonthAr

14、rayim_StateNumi=1;elsethis-m_StateNumi=51;/coutvvi+1vvtvvthis-m_YearArrayivvVvvthis-m_MonthArrayivvtvvthis-m_StateNumivve ndl;/動(dòng)態(tài)分配內(nèi)存并初始化徑流資料、蒸發(fā)數(shù)據(jù)、狀態(tài)數(shù)據(jù)等/一維數(shù)組this-m_RunoffArray = new doublethis-m_StageNum+1;this-m_EvapArray = new doublethis-m_StageNum+1;/this-m_OptVolumnArray = new doublethis-m_Stage

15、Num+1;this-m_OptNArray = new doublethis-m_StageNum+1;this-m_OptAbandonQ = new doublethis-m_StageNum+1;this-m_OptUpWaterLevelArray = new doublethis-m_StageNum+1;/最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)上游水位序列this-m_OptDownWaterLevelArray = new doublethis-m_StageNum+1; /最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)上游水位 序列this-m_OptVarialbe_q = new doublethis-m_StageN

16、um+1;/最優(yōu)軌跡對(duì)應(yīng)的發(fā)電流量序列this-m_OptMeanWaterHeadArray= new doublethis-m_StageNum+1;this-m_OptNmaxArray= new doublethis-m_StageNum+1;/最優(yōu)軌跡對(duì)應(yīng)的預(yù)想出力this-m_OptDamageDepth= new doublethis-m_StageNum+1;/最優(yōu)軌跡中低于保證出力時(shí)的破壞深度/二維數(shù)組/this-m_OptDeciVar_q = new double*this-m_StageNum+1;/this-m_StateV = new double*this-m_S

17、tageNum+1;this-m_OptArray = new long*this-m_StageNum+1;this-m_TarFun = new double*this-m_StageNum+1;for(int i=0;im_StageNum;i+)/this-m_OptDeciVar_qi = new doublethis-m_StateNumi;/this-m_StateVi = new doublethis-m_StateNumi;this-m_OptArrayi = new longthis-m_StateNumi;this-m_TarFuni = new doublethis-m

18、_StateNumi;for(int i=0;im_StageNum;i+)this-m_RunoffArrayi = 0.00;this-m_EvapArrayi = 0.00;this-m_OptVolumnArrayi = 0.00;this-m_OptNArray i = 0.00;this-m_OptAbandonQi=0.00;this-m_OptUpWaterLevelArrayi = 0.00;this-m_OptDownWaterLevelArrayi = 0.00;this-m_OptVarialbe_qi = 0.00;this-m_OptMeanWaterHeadArr

19、ayi = -1;this-m_OptNmaxArrayi= -1;最優(yōu)軌跡對(duì)應(yīng)的預(yù)想出力this-m_OptDamageDepthi= 0;/最優(yōu)軌跡中低于保證出力時(shí)的破壞深度f(wàn)or(long j=0;jm_StateNumi;j+)/this-m_OptDeciVar_qij = 0.00;/this-m_StateVij =0.00;this-m_OptArrayij = 0;this-m_TarFunij = -1e+6;/Input System Parameteters from file/ifstream myFileIn1,myFileIn2,myFileIn3;myFileI

20、n1.open(DataFile1);系統(tǒng)特性參數(shù)文件if(!myFileIn1)coutCannot Open System Parameters File this-m_DeadWaterLevelthis-m_DeadVolumnthis-m_NormalWaterLevelthis-m_FloodLimitLevel this-m_FloodVmaxthis-m_NFloodVmax this-m_UnitsNumthis-m_InstallCapacitythis-m_FirmCapacity_Nthis-m_PotentialWaterLevel1this-m_PotentialC

21、apacity_N1this-m_PotentialWaterLevel2this-m_PotentialCapacity_N2this-m_DesignAnnualGenerationthis-m_CCthis-m_LHC this-m_LHmaxthis-m_Delta_t;myFileIn1.close();/*cout”系統(tǒng)特性參數(shù)endl;cout|時(shí)段長(zhǎng)度cout|死水位cout|死庫(kù)容cout|正常蓄水位cout |汛期防洪限制水位cout|汛期最大庫(kù)容cout|非汛期最大庫(kù)容 cout|系統(tǒng)裝機(jī)容量cout|=|endl;| Delta_T = m_Delta_tsendl;|

22、 DeadWaterLevel = m_DeadWaterLevelmendl;| DeadVolumn=m_DeadVolumn億立方米endl;| NormalWaterLevel = m_NormalWaterLevelmendl;| FloodLimitLevel = m_FloodLimitLevelmendl;| FloodPeriodVmax = m_FloodVmax億立方米endl;| NonFloodPeriodVmax = m_NFloodVmax 億立方米endl;| InstallCapacity = m_UnitsNumxm_InstallCapacity/2MWe

23、ndl;cout| 系統(tǒng)保證出力| FirmCapacity = m_FirmCapacity_NMWendl;cout|水電站預(yù)想出力|;cout| 水頭|setw(6)m_PotentialWaterLevel1m|setw(4)m_PotentialWaterLevel2m|endl;cout|出力 |m_PotentialCapacity_N1MW|m_PotentialCapacity_N2MW|endl;cout| 設(shè)計(jì)年發(fā)電量| DesignAnnualGeneration = m_DesignAnnualGenerationkW.hendl;cout|水電站出力系數(shù)cout|水

24、頭損失系數(shù)cout|最大水頭損失cout|徑流系列起始年月| OutputCoefficient=m_CCendl;| HeadLossCoefficient=m_LHCendl;| HeadLossMax=m_LHmaxendl;| BeginTime=m_BeginYear-0m_BeginMonthendl;cout|=|endl; */Input runoff data from file/從文件中讀取長(zhǎng)系列徑流數(shù)據(jù) myFileIn2.open(DataFile2);/徑 流文件 if(!myFileIn2)coutCannot Open Runoff_data File ; exi

25、t(0);elsewhile(myFileIn2.good()for(int i=0;im_StageNum;i+)myFileIn2m_RunoffArrayi;/*coutsetiosflags(ios:left)setw(6)m_RunoffArrayi;if(i+1)%12=0) coutendl; */ break;myFileIn2.close();/Input evaporation data from file/從文件中讀取蒸發(fā)數(shù)據(jù)myFileIn3.open(DataFile3);/ 蒸發(fā)數(shù)據(jù)文件if(!myFileIn3)coutCannot Open Evaporatio

26、n_data File ;exit(0);elsewhile(myFileIn3.good()for(int i=0;im_StageNum;i+)myFileIn3this-m_EvapArrayi;if(0=(i+1)%12)myFileIn3.close();myFileIn3.open(DataFile3);break;myFileIn3.close();/運(yùn)行計(jì)算/void SRDP:RunningSRDP()反向逆推計(jì)算尋優(yōu)this-ReverseSearchCalculation();正向順推計(jì)算this-TerminalOptionalResultCalculation();/

27、保存計(jì)算結(jié)果/void SRDP:SaveCaluResults(char *FileName)ofstream myFileOut;myFileOut.open(FileName);/FileNam e 是保存數(shù)據(jù)的文件名myFileOutsetiosflags(ios:fixed)setprecision(3);myFileOut=endl;myFileOutsetw(35)”系統(tǒng)特性參數(shù)endl;myFileOut=endl;myFileOut時(shí)段長(zhǎng)度m_Delta_tsendl;myFileOut死水位m_DeadWaterLevelmendl;myFileOut死庫(kù)容m_DeadVo

28、lumn億立方米endl;myFileOut正常蓄水位m_NormalWaterLevelmendl;myFileOut汛期防洪限制水位m_FloodLimitLevelmendl;myFileOut汛期最大庫(kù)容”m_FloodVmax”億立方米endl;myFileOut非汛期最大庫(kù)容m_NFloodVmax億立方米endl;myFileOut系統(tǒng)裝機(jī)容量m_UnitsNumx m_InstallCapacity/2MW endl;myFileOut系統(tǒng)保證出力m_FirmCapacity_NMWendl;myFileOut水電站預(yù)想出力水頭(m)setw(6)m_PotentialWate

29、rLevel1setw(8)m_PotentialWaterLevel2endl;myFileOut出力(MW)setw(7)m_PotentialCapacity_N1setw(6)m_PotentialCapacity_N2endl;myFileOut設(shè)計(jì)年發(fā)電量m_DesignAnnualGenerationkW.hendl;myFileOut水電站出力系數(shù)m_CCendl;myFileOut水頭損失系數(shù)m_LHCendl;myFileOut最大水頭損失m_LHmaxendl;myFileOut徑流系列起始年月m_BeginYear/0m_BeginMonthendl;myFileOut

30、=endl;myFileOutendl;myFileOut= =endl;myFileOutsetw(35)最優(yōu)化計(jì)算結(jié)果endl;myFileOut=endl;myFileOutendl;myFileOutsetiosflags(ios:fixed)setprecision(2);myFileOut總發(fā)電量：m_OptGenerationProduction/100000000億KW.hendl;myFileOut多年平均發(fā)電量：m_OptGenerationProduction/43.0/100000000”億KW.hendl;myFileOut可靠保證率：m_OptReliability

31、*100%endl;myFileOut年最大發(fā)電量：m_OptAnnualMaxPower/100000000億KW.hendl;myFileOutendl;myFileOutsetiosflags(ios:fixed)setiosflags(ios:right)setw(10)時(shí)間(年/月)setw(11)來(lái)流量”setw(8)蒸發(fā)量”setw(10)初庫(kù)容”setw(10)末庫(kù)容”setw(10)初水位”setw(10)末水位”setw(10)”下游水位”setw(10)發(fā)電水頭”setw(10)發(fā)電流量setw(10)棄水量”setw(10)出力值”setw(13)預(yù)想出力值”setw(

32、11)”破壞深度”endl;for(int i=0;im_StageNum-1;i+)時(shí)間(年/月)輸出myFileOutsetw(5)m_YearArrayim_MonthArrayi10)myFileOut0m_MonthArrayit;elsemyFileOutm_MonthArrayit;數(shù)據(jù)輸出myFileOutsetiosflags(ios:fixed)setprecision(2)setw(6)m_RunoffArrayisetw(7)m_EvapArrayi;myFileOutsetiosflags(ios:fixed)setprecision(3)setiosflags(io

33、s:right)setw(11)m_OptVolumnArrayisetw(10)m_OptVolumnArrayi+1setw(10)m_OptUpWaterLevelArrayisetw(10)m_OptUpWaterLevelArrayi+1setw(10)m_OptDownWaterLevelArrayisetw(10)m_OptMeanWaterHeadArrayisetw(10)m_OptVarialbe_qisetw(10)m_OptAbandonQisetw(10)m_OptNArrayisetw(12)m_OptNmaxArrayisetw(10)m_OptDamageDep

34、thi;myFileOutsetw(5)m_YearArrayim_MonthArrayi10)myFileOut0 m_MonthArrayi;elsemyFileOutm_MonthArrayi;myFileOutendl;myFileOut.close();ofstream myfileout;myfileout.open(WaterLevel.txt);for(int i=0;im_StageNum-1;i+)myfileoutsetw(5)m_YearArrayim_MonthArrayi10)myfileout0m_MonthArrayit;elsemyfileoutm_Month

35、Arrayit; myfileoutm_OptUpWaterLevelArrayim_TarFunthis-m_StageNum0=0.00; for(int i=this-m_StageNum-1;i=0;i-) 對(duì)庫(kù)容進(jìn)行離散 if(1=this-m_StateNumi) Delta_V1 = this-m_FloodVmax - this-m_DeadVolumn; elseDelta_V1 = (this-m_NFloodVmax - this-m_DeadVolumn)/(this-m_StateNumi-1); if(1=this-m_StateNumi+1)Delta_V2 =

36、this-m_FloodVmax - this-m_DeadVolumn; elseDelta_V2 = (this-m_NFloodVmax - this-m_DeadVolumn )/(this-m_StateNumi+1-1);/*coutvvivvsetw(4)vvthis-m_MonthArrayivvsetw(10)vvDelta_V1vvsetw(4)vvthis-m_StateNumivv setw(10)vvthis-m_MonthArrayi+1vvsetw(10)vvDelta_V2vvsetw(4)vvthis-m_StateNumi+1vvendl;*/ 進(jìn)行逆推迭代

37、計(jì)算for(int j=0;jm_StateNumi;j+) V1 = this-m_DeadVolumn+j*Delta_V1;單位：億立方米for(int k=0;km_StateNumi+1;k+)V2 = this-m_DeadVolumn+k*Delta_V2;單位：億立方米水量平衡驗(yàn)算求出發(fā)電流量q=this-m_RunoffArrayi-(V2 -V1)*100000000/this-m_Delta_t-this-m_EvapArrayi;/ 單 位：立方米/秒if(qUpWaterLevelCalculation(V1)+this-UpWaterLevelCalculation

38、(V2)/2.0; 水位流量關(guān)系曲線求下游水位值Z_down = this-DownWaterLevelCalculation(q); 水頭損失計(jì)算和發(fā)電水頭計(jì)算 HeadLoss= this-m_LHC*q*q; if(HeadLossthis-m_LHmax)HeadLoss = this-m_LHmax; GWH = Z_up-Z_down-HeadLoss; if(GWH=0)N = 0.00; else 最大出力值計(jì)算if(GWH=57) Nmax = 120.0;elseNmax=2.94*GWH-47.58;/coutsetw(10)Nmaxm_CC*q*GWH/1000;if(

39、Nm_FirmCapacity_N)/低于保證出力 N=N-2000*(this-m_FirmCapacity_N-N); else if(N=Nmax)/大于最大出力限制時(shí)按最大出力計(jì)算N = Nmax; /couti”t”j”t”k”t”Nm_TarFuni+1k;if(tempTarFun=this-m_TarFunij)this-m_TarFunij=tempTarFun;this-m_OptArrayij = k;this-m_ObjFun = this-m_TarFun00;/coutm_ObjFunm_OptVolumnArray0=this-m_DeadVolumn;for(i

40、nt i=0;im_StageNum;i+)this-m_OptUpWaterLevelArrayi=this-UpWaterLevelCalculation(this-m_OptVolumnArrayi);對(duì)庫(kù)容進(jìn)行離散if(1=this-m_StateNumi+1)Delta_V = this-m_FloodVmax - this-m_DeadVolumn;elseDelta_V = (this-m_NFloodVmax - this-m_DeadVolumn)/(this-m_StateNumi+1-1);m=this-m_OptArrayim;this-m_OptVolumnArray

41、i+1=this-m_DeadVolumn+m*Delta_V;時(shí)段初末水位值計(jì)算Z1=this-UpWaterLevelCalculation(this-m_OptVolumnArrayi);Z2=this-UpWaterLevelCalculation(this-m_OptVolumnArrayi+1);水量平衡計(jì)算單位：立方米/秒q=this-m_RunoffArrayi-this-m_EvapArrayi-(this-m_OptVolumnArrayi+1-this-m_OptVolumnArrayi)*100000000/this-m_Delta_t;if(qm_OptNArrayi

42、=0.00;this-m_OptVarialbe_qi=0.00;this-m_OptAbandonQi=0.00;this-m_OptDamageDepthi=-10000;count+; else /平均發(fā)電水頭計(jì)算Z_up = (Z1+Z2)/2;Z_down = this-DownWaterLevelCalculation(q);Delta_H = this-m_LHC*q*q;if(Delta_H5.59)Delta_H = 5.59;GWH = Z_up-Z_down-Delta_H;預(yù)想出力值計(jì)算if(GWH=57)Nmax = 120.0;elseNmax=2.94*GWH-4

43、7.58;N=this-m_CC*q*GWH/1000; 出力計(jì)算if(NNmax) 出力大于最大出力double q0,q1,q2;q0=q; q1=0.00; q2=q;while(abs(N-Nmax)0.0001)q=(q1+q2)/2;Z_down = this-DownWaterLevelCalculation(q);Delta_H = this - m_LHC*q*q;if(Delta_H5.59)Delta_H = 5.59;GWH = Z_up-Z_down-Delta_H;/預(yù)想出力值計(jì)算if(GWH=57)Nmax = 120.0;elseNmax=2.94*GWH-47

44、.58;N=this-m_CC*q*GWH/1000; 出力計(jì)算if(NNmax)q2=q;elseq1=q;計(jì)算結(jié)果this-m_OptNArrayi=N;this-m_OptAbandonQi=q0-q;this-m_OptVarialbe_qi=q;/最優(yōu)軌跡對(duì)應(yīng)出力數(shù)據(jù)數(shù)組序列/最優(yōu)軌跡對(duì)應(yīng)的棄水量序列/最優(yōu)軌跡對(duì)應(yīng)的發(fā)電流量序列this-m_OptMeanWaterHeadArrayi=GWH; /最優(yōu)軌跡對(duì)應(yīng)的平均發(fā)電水頭 this-m_OptNmaxArrayi=Nmax;/最優(yōu)軌跡對(duì)應(yīng)的預(yù)想出力this-m_OptDownWaterLevelArrayi=Z_down最優(yōu)軌跡對(duì)應(yīng)的水庫(kù)下游水位序列 else if(Nm_FirmCapa