漫湾水库调度程序

漫湾水库的详细资料及调度程序，可输出，是初学者的不二选择

/*本程序名为“漫湾水电站优化调度”，
本程序有主程序及四个子程序组成。
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]={456789101112123}
Q[T]={64977422502448266621411397744496389355374}
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 jliLT;
int V[100]yue[100];
float Nqfmv0v1hzpzbvzvsvtawxfKNBNX;
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;iscanf（“%f“&Q[i]）;
printf（“请输入对应的月份“）;
for（i=0;iscanf（“%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）/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			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（QT）;
printf（“月初蓄水\t“）;
supt（v0pT）;
printf（“发电流量\t“）;
supt（qpT）;
printf（“弃水量\t\t“）;
supt（wpT）;
printf（“电站出力\t“）;
supt（NpT）;
printf（“电站年出力=%f万kW\n“fm）;
printf（“电站设计中水年年发电量为E=%f亿KW·h\n“（fm*730/10000））;
getch（）;
}
	

float suzp（float v）
{
	float z;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-12-27 11:43  漫湾水电站动态规划调度\
     目录           0  2014-12-17 09:10  漫湾水电站动态规划调度\Debug\
     文件       33792  2014-12-27 09:23  漫湾水电站动态规划调度\Debug\vc60.idb
     文件       45056  2014-12-13 11:32  漫湾水电站动态规划调度\Debug\vc60.pdb
     文件      213088  2014-12-13 11:32  漫湾水电站动态规划调度\Debug\漫湾水电站优化调度.exe
     文件      216764  2014-12-13 11:32  漫湾水电站动态规划调度\Debug\漫湾水电站优化调度.ilk
     文件       14364  2014-12-13 11:32  漫湾水电站动态规划调度\Debug\漫湾水电站优化调度.obj
     文件      209732  2014-12-13 11:32  漫湾水电站动态规划调度\Debug\漫湾水电站优化调度.pch
     文件      484352  2014-12-13 11:32  漫湾水电站动态规划调度\Debug\漫湾水电站优化调度.pdb
     文件        3828  2014-12-13 11:32  漫湾水电站动态规划调度\漫湾水电站优化调度.cpp
     文件        3547  2014-10-10 10:05  漫湾水电站动态规划调度\漫湾水电站优化调度.dsp
     文件         544  2014-10-10 10:06  漫湾水电站动态规划调度\漫湾水电站优化调度.dsw
     文件       41984  2014-12-27 11:43  漫湾水电站动态规划调度\漫湾水电站优化调度.ncb
     文件       48640  2014-12-27 11:43  漫湾水电站动态规划调度\漫湾水电站优化调度.opt
     文件         270  2014-12-27 09:23  漫湾水电站动态规划调度\漫湾水电站优化调度.plg
     文件      196670  2009-06-03 08:53  漫湾水电站动态规划调度\漫湾水电站优化调度结果.exe
     文件     1456640  2010-03-22 20:49  漫湾水电站动态规划调度\漫湾水电站基本资料及数据处理.doc