水文分析与计算——频率分析及洪水流量过程

//频率分析.h
void PinLvFenXi()
{
	using namespace std;
	const int NYear = EndYear - 1862 + 1,//调查考证期(1862-2000)
		L = 1,//实测期特大洪水个数
		a = 2,//考证期及实测期特大洪水个数
		YearSuperYearW[a] = {1903, 1975},//特大洪水年份
		SuperYearW[a] = {21000, 12300};//特大洪水值
	int order;//降序的序号
	double AverageMaxYearW[J] = {0},
		HundredYearW[J] = {0},//不同年极值的百年一遇设计值
		K,//模比系数
		Cv = 0;//变差系数近似无偏估计值
		ofstream outfile;
		outfile.open("outfile_YearW_P.txt");


	outfile<>ClassicalMonthQ[i];
	infile.close();
	outfile.open("outfile_ForecastWprocess.txt");

	for(int k = 0; kMaxClassicalYearW[k])
			 {
				MaxClassicalYearW[k] = temp_W;//(m3/s*月)
				date[k] = j;
			 }
		 }
		 cout<<"最大"<>order;
	if(order == 1)
	{
		cout<<"输入由理论频率曲线得到的百年一遇设计值："<>HundredYearW[k];
			if( k == 0) RatioW[k] = HundredYearW[k]/MaxClassicalYearW[k];
			else RatioW[k] = (HundredYearW[k] - HundredYearW[k - 1])/(MaxClassicalYearW[k] - MaxClassicalYearW[k - 1]);
		}
		for(int k = 0; k < J; k++)
		{//求设计洪水过程
			for(int i = date[k]; i < date[k]+NumJ[k]; i++)
			{
				if(k == 0)	ForcastQ[i] = ClassicalMonthQ[i]*RatioW[k];
				else if(i < date[k-1]|| i>=date[k-1]+NumJ[k-1])
					ForcastQ[i] = ClassicalMonthQ[i]*RatioW[k];
			}
		}
		outfile<>Cv;//控制台暂停
	}
}