/*本程序名為“漫灣水電站優(yōu)化調(diào)度”，
本程序有主程序及四個子程序組成。
1）其中四個子程序的功能說明如下：
SUZP--水位庫容關(guān)系子程序；
SUZB--下游水位流量關(guān)系子程序；
SUPT--輸出計算結(jié)果子程序
PT--輸出資料子程序。
2）簡單變量標(biāo)識符：
T--總時段數(shù)；
L--總格點數(shù)；
j--階段變量數(shù)（時段數(shù)）；
l--當(dāng)前時段格點數(shù)
i--（有時為臨時變量）下一時段格點數(shù)；
a--單位格點庫容增量，億m3；
vz--水庫正常庫容，億m3；
vs--水庫死庫容，億m3；
NX--水庫總裝機(jī)（限制出力），萬kW；
NB--水庫保證出力，萬kW；
K--出力系數(shù)；
N--水電站出力，萬kW；
h--水電站發(fā)電水頭，m；
q--發(fā)電流量，m3/s；
w--棄水流量，m3/s；
v0--末庫容，億m3；
v1--初庫容，億m3；
zp--水庫上游平均水位；
zb--下游水位；
fm--全時期總出力，萬kW；
x--余留時期出力，萬kW；
f--已確定最優(yōu)余留時期出力，萬kW。
3）數(shù)組標(biāo)識符：
V[T]--最優(yōu)化庫容格點；
yue[T]--時段標(biāo)號；
Q[T]--時段來水量，m?/s；
N1[L+1][?L+1]--本時段各組合庫容出力，萬kW；
N2[L+1]--余留期效益，萬kW；
v0p[T]--最優(yōu)化后各時段末水位，億m?；
qp[T]--最優(yōu)化后各時段發(fā)電流量，m?/s；
wp[T]--最優(yōu)化后各時段棄水量，m?/s；
Np[T]--最優(yōu)化后各時段出力，萬kW；
4）漫灣水電站設(shè)計中水年的資料為:
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漫灣水電站興利調(diào)度程序\n“）;
printf（“輸入階段數(shù)，狀態(tài)數(shù)（格點數(shù)），末狀態(tài)及其他原始數(shù)據(jù)\n“）;
printf（“請輸入階段數(shù)T=“）;
scanf（“%d“&T）;
printf（“請輸入格點數(shù)L=“）;
scanf（“%d“&L）;
printf（“請輸入出力系數(shù)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（“請輸入對應(yīng)的月份“）;
for（i=0;iscanf（“%d“&yue[i]）;
printf（“請輸入末狀態(tài)“億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漫灣水電站優(yōu)化調(diào)度軌跡線\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（“發(fā)電流量\t“）;
supt（qpT）;
printf（“棄水量\t\t“）;
supt（wpT）;
printf（“電站出力\t“）;
supt（NpT）;
printf（“電站年出力=%f萬kW\n“fm）;
printf（“電站設(shè)計中水年年發(fā)電量為E=%f億KW·h\n“（fm*730/10000））;
getch（）;
}
	

float?suzp（float?v）
{
	float?z;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-12-27?11:43??漫灣水電站動態(tài)規(guī)劃調(diào)度\
?????目錄???????????0??2014-12-17?09:10??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\
?????文件???????33792??2014-12-27?09:23??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\vc60.idb
?????文件???????45056??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\vc60.pdb
?????文件??????213088??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\漫灣水電站優(yōu)化調(diào)度.exe
?????文件??????216764??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\漫灣水電站優(yōu)化調(diào)度.ilk
?????文件???????14364??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\漫灣水電站優(yōu)化調(diào)度.obj
?????文件??????209732??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\漫灣水電站優(yōu)化調(diào)度.pch
?????文件??????484352??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\Debug\漫灣水電站優(yōu)化調(diào)度.pdb
?????文件????????3828??2014-12-13?11:32??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度.cpp
?????文件????????3547??2014-10-10?10:05??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度.dsp
?????文件?????????544??2014-10-10?10:06??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度.dsw
?????文件???????41984??2014-12-27?11:43??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度.ncb
?????文件???????48640??2014-12-27?11:43??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度.opt
?????文件?????????270??2014-12-27?09:23??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度.plg
?????文件??????196670??2009-06-03?08:53??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站優(yōu)化調(diào)度結(jié)果.exe
?????文件?????1456640??2010-03-22?20:49??漫灣水電站動態(tài)規(guī)劃調(diào)度\漫灣水電站基本資料及數(shù)據(jù)處理.doc