【正文】 double z。 if(vV0[0]) z=Zv[0]+(Zv[1]Zv[0])/(V0[1]V0[0])*(vV0[0])。 else { if(vV0[m1]) z=Zv[m2]+(Zv[m1]Zv[m2])/(V0[m1]V0[m2])*(vV0[m2])。 else { for(int i=0。im1。i++) { if(vV0[i]) index++。 } z=Zv[index1]+(Zv[index]Zv[index1])/(V0[index]V0[index1])*(vV0[index1])。 } } return z。 } double ZtoV(double z) { int index=0。 double v。 if(zZv[0]) v=V0[0]+(V0[1]V0[0])/(Zv[1]Zv[0])*(zZv[0])。 else { if(zZv[m1]) v=V0[m2]+(V0[m1]V0[m2])/(Zv[m1]Zv[m2])*(zZv[m2])。 else { for(int i=0。im1。i++) { if(zZv[i]) index++。 } v=V0[index1]+(V0[index]V0[index1])/(Zv[index]Zv[index1])*(zZv[index1])。 } } return v。 } double ZtoQ(double z) { int index=0。 double q。 if(zZq[0]) q=Q0[0]+(Q0[1]Q0[0])/(Zq[1]Zq[0])*(zZq[0])。 else { if(zZq[n1]) q=Q0[n2]+(Q0[n1]Q0[n2])/(Zq[n1]Zq[n2])*(zZq[n2])。 else { for(int i=0。in1。i++) { if(zZq[i]) index++。 } q=Q0[index1]+(Q0[index]Q0[index1])/(Zq[index]Zq[index1])*(zZq[index1])。 } } return q。 } //演算防洪高水位 void calculateZ1() { fstream out1。 (5%,ios::out)。 double Q[period],q[period],V[period]。 double max=0。 double Z1[period]。 fstream input。 (5%.txt,ios::in)。 for(int i=0。iperiod。i++) inputQ[i]。 V[0]=ZtoV(Z0)。 max=V[0]。 for(i=0。iperiod。i++) { q[i]=Q[i]。 if(q[i]Qsafe) q[i]=Qsafe。 } for(int j=1。jperiod。j++) { V[j]=V[j1]+*(Q[j]+Q[j1])*dt/*(q[j]+q[j1])*dt/1000000。 if(maxV[j]) max=V[j]。 } Z1max=VtoZ(max)。 out1Q q V Zendl。 coutQq V Zendl。 for(i=0,j=0。iperiod。i++,j++) { Z1[i]=VtoZ(V[j])。 out1Q[i] q[i] V[i] Z1[i]endl。 } for(i=0。iperiod。i++) coutQ[i] q[i] V[i] Z1[i]endl。 out1防洪高水位為 Z1maxendl。 cout防洪高水位為 Z1maxendl。 } //演算設計洪水位 void calculateZ2() { fstream out2。 (1%,ios::out)。 double t=0。 double Q[period],q[period],V[period]。 double Z2[period]。 double max=0,temp=0。 fstream input。 (1%.txt,ios::in)。 for(int i=0。iperiod。i++) inputQ[i]。 V[0]=ZtoV(Z0)。 max=V[0]。 for(i=0。iperiod。i++) { q[i]=Q[i]。 if(q[i]Qsafe) q[i]=Qsafe。 } for(int j=1。jperiod。j++) { if(VtoZ(V[j1])Z1max) t=VtoZ(V[j1])。 if(tZ1max) { while(true) { V[j]=V[j1]+*(Q[j]+Q[j1])*dt/*(q[j]+q[j1])*dt/1000000。 temp=ZtoQ(VtoZ(V[j]))。 if(fabs(q[j]temp)e) break。 else q[j]=*(q[j]+temp)。 } } else V[j]=V[j1]+*(Q[j]+Q[j1])*dt/*(q[j]+q[j1])*dt/1000000。 if(VtoZ(V[j])Z0) q[j]=Q[j]。 if(maxV[j]) max=V[j]。 } Z2max=VtoZ(max)。 out2Q q V Zendl。 coutQq V Zendl。 for(i=0,j=0。iperiod。i++,j++) { Z2[i]=VtoZ(V[j])。 out2Q[i] q[i] V[i] Z2[i]endl。 } for(i=0。iperiod。i++) coutQ[i] q[i] V[i] Z2[i]endl。 out2設計洪水位為 Z2maxendl。 cout設計洪水位為 Z2maxendl。 } //演算校核洪水位 void calculateZ3() { fstream out3。 (%,ios::out)。 double a=0,b=0。 double Q[period],q[period],V[period]。 double Z3[period]。 double max=0,temp=0。 fstream input。 (%.txt,ios::in)。 for(int i=0。iperiod。i++) inputQ[i]。 V[0]=ZtoV(Z0)。 max=V[0]。 for(i=0。iperiod。i++) { q[i]=Q[i]。 if(q[i]Qsafe) q[i]=Qsafe。 } for(int j=1。jperiod。j++) { b=VtoZ(V[j1])。 if(bZ2max) { while(true) { V[j]=V[j1]+*(Q[j]+Q[j1])*dt/*(q[j]+q[j1])*dt/1000000。 temp=ZtoQ(VtoZ(V[j]))+100。 if(fabs(q[j]temp)e) break。 else q[j]=*(q[j]+temp)。 } } else { if(VtoZ(V[j1])Z1max) a=VtoZ(V[j1])。 if(aZ1max) { while(true) { V[j]=V[j1]+*(Q[j]+Q[j1])*dt/*(q[j]+q[j1])*dt/1000000。 temp=ZtoQ(VtoZ(V[j]))。 if(fabs(q[j]temp)e) break。 else q[j]=*(q[j]+temp)。 } } else V[j]=V[j1]+*(Q[j]+Q[j1])*dt/*(q[j]+q[j1])*dt/1000000。 } if(VtoZ(V[j])Z0) q[j]=Q[j]。 if(maxV[j]) max=V[j]。 } Z3max=VtoZ(max)。 out3Q q V Zendl。 coutQq V Zendl。 for(i=0,j=0。iperiod。i++,j++) { Z3[i]=VtoZ(V[j])。 out3Q[i] q[i] V[i] Z3[i]endl。 } for(i=0。iperiod。i++) coutQ[i] q[i] V[i] Z3[i]endl。 out3校核洪水位為 Z3maxendl。 cout校核洪水位為 Z3maxendl。 } void main() { inputData()。 calculateZ1()。 calculateZ2()。 calculateZ3()。 }