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

//频率分析.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<<setw(10)<<"年份"<<setw(10)<<"最大1日"
			<<setw(10)<<"降序序号"<<setw(10)<<"经验频率"<<endl;
	outfile<<"调查考证期："<<endl;
	for(int M = 0; M < a; M++)//最大1日特大洪水经验频率
		outfile<<setw(10)<<YearSuperYearW[M]<<setw(10)<<SuperYearW[M]
			<<setw(10)<<M+1<<setw(10)<<(double)(M+1)/(NYear+1)<<endl;
	outfile<<"实