水电站优化调度C语言

漫湾水电站优化调度C++程序!

/*本程序名为“漫湾水电站优化调度”， 本程序有主程序及四个子程序组成。 1）其中四个子程序的功能说明如下： SUZP--水位库容关系子程序； SUZB--下游水位流量关系子程序； SUPT--输出计算结果子程序 PT--输出资料子程序。 2）简单变量标识符： T--总时段数； L--总格点数； j--阶段变量数（时段数）； l--当前时段格点数 i--(有时为临时变量)下一时段格点数； a--单位格点库容增量，亿m3； vz--水库正常库容，亿m3； vs--水库死库容，亿m3； NX--水库总装机（限制出力），万kW； NB--水库保证出力，万kW； K--出力系数； N--水电站出力，万kW； h--水电站发电水头，m； q--发电流量，m3/s； w--弃水流量，m3/s； v0--末库容，亿m3； v1--初库容，亿m3； zp--水库上游平均水位； zb--下游水位； fm--全时期总出力，万kW； x--余留时期出力，万kW； f--已确定最优余留时期出力，万kW。 3）数组标识符： V[T]--最优化库容格点； yue[T]--时段标号； Q[T]--时段来水量，m?/s； N1[L+1][ L+1]--本时段各组合库容出力，万kW； N2[L+1]--余留期效益，万kW； v0p[T]--最优化后各时段末水位，亿m?； qp[T]--最优化后各时段发电流量，m?/s； wp[T]--最优化后各时段弃水量，m?/s； Np[T]--最优化后各时段出力，万kW； 4）漫湾水电站设计中水年的资料为: NB=26.584 NX=125 K=8.4366 T=12 yue[12]={4,5,6,7,8,9,10,11,12,1,2,3} Q[T]={649,774,2250,2448,2666,2141,1397,744,496,389,355,374} vz=9.2 vs=6.63 vt=6.63*/ #include "conio.h" #include "stdio.h" float suzp(float); float suzb(float); void supt(float *,int); void main() { int j,l,i,L,T; int V[100],yue[100]; float N,q,fm,v0,v1,h,zp,zb,vz,vs,vt,a,w,x,f,K,NB,NX; float N1[100][100],N2[100]; float v0p[100],qp[100],wp[100],Np[100]; float Q[100]; printf("\t\t漫湾水电站兴利调度程序\n"); printf("输入阶段数，状态数（格点数），末状态及其他原始数据\n"); printf("请输入阶段数T="); scanf("%d",&T); printf("请输入格点数L="); scanf("%d",&L); printf("请输入出力系数K="); scanf("%f",&K); printf("请输入正常库容“亿m3”vz="); scanf("%f",&vz); printf("请输入死库容“亿m3”vs="); scanf("%f",&vs); printf("请输入保证出力“万KW·h”NB="); scanf("%f",&NB); printf("请输入限制出力“万KW·h”NX="); scanf("%f",&NX); printf("请输入月来水量“m3/s”(按顺序输入，末库容确定的为最后一个月)"); for(i=0;i<T;i++) scanf("%f",&Q[i]); printf("请输入对应的月份"); for(i=0;i<T;i++) scanf("%d",&yue[i]); printf("请输入末状态“亿m3”vt="); scanf("%f",&vt); V[T-1]=0; fm=0; f=0; x=0; a=(vz-vs)/L; for(j=T;j>=1;j--) { for(l=0;l<=L;l++) { for(i=0;i<=L;i++) { v0=vs+a*l; if(j==T) v1=vt; else v1=vs+a*i; zp=suzp((v0+v1)/

漫湾水电站优化调度C++程序!

2); q=(v0-v1)/0.0263+Q[j-1]; zb=suzb(q); h=zp-zb; N=K*h*q/10000; N1[l][i]=N; if(N<NB) break; if(N>NX) N1[l][i]=NX; if(j==T) break; if(fm<=(N1[l][i]+N2[i]+f)) { fm=N1[l][i]+N2[i]+f; x=N2[i]+f; V[j-1]=i; } } if(j==1) break; } for(l=0;l<=L;l++) N2[l]=N1[l][V[j-1]]; f=x; } for(j=1;j<=T;j++) { v0=vs+a*V[j-2]; if(j==1) v0=vs; v1=vs+a*V[j-1]; zp=suzp((v0+v1)/2); q=(v0-v1)/0.0263+Q[j-1]; zb=suzb(q); h=zp-zb; N=K*h*q/10000; w=0; if(N>NX) { N=NX; w=q-N/K/h*10000; q=N/K/h*10000; } v0p[j-1]=v0; qp[j-1]=q; wp[j-1]=w; Np[j-1]=N; } printf("\t\t\t\t\t漫湾水电站优化调度轨迹线\n"); printf("月份\t\t"); for(j=1;j<=T;j++) printf("%d\t",yue[j-1]); printf("\n"); printf("入库流量\t"); supt(Q,T); printf("月初蓄水\t"); supt(v0p,T); printf("发电流量\t"); supt(qp,T); printf("弃水量\t\t"); supt(wp,T); printf("电站出力\t"); supt(Np,T); printf("电站年出力=%f万kW\n",fm); printf("电站设计中水年年发电量为E=%f亿KW·h\n",(fm*730/10000)); getch(); } float suzp(float v) { float z; z=934.4+9.1185*v-v*v*0.2857; return z; } float suzb(float q) { float z; z=(89077+0.32*q-q*q*0.00001)*0.01; return z; } void supt(float *a,int T) { int j; for(j=1;j<=T;j++) printf("%3.2f\t",a[j-1]); printf("\n"); }