📄 考虑泄流能力.txt
字号:
// dpView.cpp : implementation of the CDpView class
//
#include "stdafx.h"
#include "dp.h"
#include "dpDoc.h"
#include "dpView.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif
/////////////////////////////////////////////////////////////////////////////
// CDpView
IMPLEMENT_DYNCREATE(CDpView, CView)
BEGIN_MESSAGE_MAP(CDpView, CView)
//{{AFX_MSG_MAP(CDpView)
ON_COMMAND(IDM_DPYH, OnDpyh)
//}}AFX_MSG_MAP
// Standard printing commands
ON_COMMAND(ID_FILE_PRINT, CView::OnFilePrint)
ON_COMMAND(ID_FILE_PRINT_DIRECT, CView::OnFilePrint)
ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView::OnFilePrintPreview)
END_MESSAGE_MAP()
/////////////////////////////////////////////////////////////////////////////
// CDpView construction/destruction
CDpView::CDpView()
{
// TODO: add construction code here
}
CDpView::~CDpView()
{
}
BOOL CDpView::PreCreateWindow(CREATESTRUCT& cs)
{
// TODO: Modify the Window class or styles here by modifying
// the CREATESTRUCT cs
return CView::PreCreateWindow(cs);
}
/////////////////////////////////////////////////////////////////////////////
// CDpView drawing
void CDpView::OnDraw(CDC* pDC)
{
CDpDoc* pDoc = GetDocument();
ASSERT_VALID(pDoc);
// TODO: add draw code for native data here
}
/////////////////////////////////////////////////////////////////////////////
// CDpView printing
BOOL CDpView::OnPreparePrinting(CPrintInfo* pInfo)
{
// default preparation
return DoPreparePrinting(pInfo);
}
void CDpView::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
{
// TODO: add extra initialization before printing
}
void CDpView::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
{
// TODO: add cleanup after printing
}
/////////////////////////////////////////////////////////////////////////////
// CDpView diagnostics
#ifdef _DEBUG
void CDpView::AssertValid() const
{
CView::AssertValid();
}
void CDpView::Dump(CDumpContext& dc) const
{
CView::Dump(dc);
}
CDpDoc* CDpView::GetDocument() // non-debug version is inline
{
ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CDpDoc)));
return (CDpDoc*)m_pDocument;
}
#endif //_DEBUG
/////////////////////////////////////////////////////////////////////////////
// CDpView message handlers
void CDpView::OnDpyh()
{
int i,j,k,l,m,n;//循环变量
int ll,kk;
long int a,b,c;//备用参数
float x,y,z;//备用参数
float x1,x2,x3;
int mark;//标记用
int sign;
float Vu_bottom,Vu_top,hu_bottom=614.0,hu_top=620.0;
float Vd_bottom,Vd_top,hd_bottom=553.5,hd_top=556.5;
int U_st_num,D_st_num;
float U_st[1000],D_st[1000];
readdata_flood_TZY();
////以库容作为状态变量，步长100万m3
long int step=4*10e4;//步长
///////////生成状态序列
///////////上库
Vu_bottom=HtoV_TZY(hu_bottom);
Vu_top=HtoV_TZY(hu_top);
a=Vu_bottom/step;
Vu_bottom=(a+1)*step;//忽略不足步长的库容
a=Vu_top/step;
Vu_top=a*step;//忽略小于步长的库容
U_st_num=(Vu_top-Vu_bottom)/step;
for(i=0;i<U_st_num;i++)
U_st[i]=Vu_bottom+i*step;//上库状态序列
Vd_bottom=HtoV_LD(hd_bottom);
Vd_top=HtoV_LD(hd_top);
a=Vd_bottom/step;
Vd_bottom=(a+1)*step;//忽略不足步长的库容
a=Vd_top/step;
Vd_top=a*step;//忽略小于步长的库容
D_st_num=(Vd_top-Vd_bottom)/step;
for(i=0;i<D_st_num;i++)
D_st[i]=Vd_bottom+i*step;//下库状态序列
///////////////////////////////////////////////////
///////////////////////////////////////////////////
///////////////////////////////////////////////////
FILE *fp=fopen("检验状态序列.txt","w");
for (i=0;i<U_st_num;i++)
{
fprintf(fp,"%4.2f\n",U_st[i]);
}
for (i=0;i<D_st_num;i++)
{
fprintf(fp,"%4.2f\n",D_st[i]);
}
fclose(fp);
////////////////////////////////////////////////////
////////////////////////////////////////////////////
////////////////////////////////////////////////////
////dp优化
float ff[50][50];//阶段优化目标存储数组[U_st][D_st]
float Qo_dp_TZY[50][50];
float tt=3600;
float Qin_dp_LD[50][50],Qo_dp_LD_0[50][25];
float y1;
float Qo_dp_LD_tmp[50][50];
float ff_tmp[50];
float AA=10e10;
////////////第一阶段
for(ll=0;ll<U_st_num;ll++)
for(kk=0;kk<D_st_num;kk++)
{
Qo_dp_TZY[0][ll]=Qin_TZY[0]-(U_st[ll]-U_st[0])/tt;
Qin_dp_LD[0][ll]=Qin_TZY[0]*0.4632+Qo_dp_TZY[0][ll];
Qo_dp_LD_0[ll][kk]=Qin_dp_LD[0][ll]-(D_st[kk]-D_st[0])/tt;
x=VtoQo_TZY(U_st[ll]);
y=VtoQo_LD(D_st[kk]);
x1=VtoH_LD(D_st[kk]);
if((Qo_dp_TZY[0][ll]>=0)&&(Qo_dp_LD_0[ll][kk]>=0))///水量平衡所得流量为正值
if((Qo_dp_TZY[0][ll]<x)&&(Qo_dp_LD_0[ll][kk]<y))///泄流能力约束
if(
((x1<=553.8)&&(Qo_dp_LD_0[ll][kk]<=1800))///（下游）泄流限制，可能要调整
||
(x1>553.8)
)
{
ff[ll][kk]=0.6*Qo_dp_TZY[0][ll]+0.4*Qo_dp_LD_0[ll][kk];
z=10;
}
else
ff[ll][kk]=AA;
z=0;
}
z=0;
//////////中间阶段
float ff_tmp1[50][50][25];
float ff_tmp2[50][25];
int path_tmp[50][50][25],path_u[50][50],path_d[50][25];
for(i=1;i<amount_flood-1;i++)
{
for(l=0;l<U_st_num;l++)
for(ll=0;ll<U_st_num;ll++)
for(k=0;k<D_st_num;k++)
{
for(kk=0;kk<D_st_num;kk++)
{
Qo_dp_TZY[ll][l]=Qin_TZY[i]-(U_st[l]-U_st[ll])/tt;
Qin_dp_LD[ll][l]=Qin_TZY[i]*0.4632+Qo_dp_TZY[ll][l];
Qo_dp_LD_tmp[kk][k]=Qin_dp_LD[ll][l]-(D_st[k]-D_st[kk])/tt;
if(U_st[l]<U_st[ll])
x=VtoQo_TZY(U_st[ll]);
else
x=VtoQo_TZY(U_st[l]);
if(D_st[k]<D_st[kk])
y=VtoQo_LD(D_st[kk]);
else
y=VtoQo_LD(D_st[k]);
x1=VtoH_LD(D_st[k]);
if((Qo_dp_TZY[ll][l]>0)&&(Qo_dp_LD_tmp[kk][k]>0))
if((Qo_dp_TZY[ll][l]<x)&&(Qo_dp_LD_tmp[kk][k]<y))
if(
((x1<553.8)&&(Qo_dp_LD_tmp[kk][k]<1800))
||
(x1>553.8)
)
{
ff_tmp[kk]=0.6*Qo_dp_TZY[ll][l]+0.4*Qo_dp_LD_tmp[kk][k];
}
else
ff_tmp[kk]=AA;
}//kk循环
x=ff_tmp[0];
for(kk=0;kk<D_st_num;kk++)
{
if(x>=ff_tmp[kk])
{
x=ff_tmp[kk];
a=kk;
}
}
z=0;
if(l==2)
{
z=9;
}
if(l==5)
{
z=0;
z=9;
}
ff_tmp1[l][ll][k]=x;
path_tmp[l][ll][k]=a;
}//k,ll,l循环
for(l=0;l<U_st_num;l++)
for(k=0;k<D_st_num;k++)
{
x=ff_tmp1[l][0][k];
for(ll=0;ll<U_st_num;ll++)
{
if(x>=ff_tmp1[l][ll][k])
{
x=ff_tmp1[l][ll][k];
a=ll;
b=path_tmp[l][ll][k];
c=k;
}
}
path_u[i][l]=a;
path_d[i][k]=b;
if(x>ff[a][b])
ff_tmp2[l][k]=x;
else
ff_tmp2[l][k]=ff[a][b];
}
z=0;
for(l=0;l<U_st_num;l++)
for(k=0;k<D_st_num;k++)
{
ff[l][k]=ff_tmp2[l][k];
if(
((ff[l][k]-3649.71)>1)
||
((ff[l][k]-3649.71)<-1)
)
z=1;
}
if(z==0)
{
z=12;
z=11;
}
}//i循环
/////////////////最后阶段
float ff_tmp3[50][25];
float Qo_dp_LD_00[50][25];
for(ll=0;ll<U_st_num;ll++)
for(kk=0;kk<D_st_num;kk++)
{
Qo_dp_TZY[ll][0]=Qin_TZY[i]-(U_st[ll]-U_st[0])/tt;
Qin_dp_LD[ll][0]=0.4623*Qin_TZY[i]+Qo_dp_TZY[ll][0];
Qo_dp_LD_00[ll][kk]=Qin_dp_LD[ll][0]-(D_st[kk]-D_st[0])/tt;
x=VtoQo_TZY(U_st[ll]);
y=VtoQo_LD(D_st[kk]);
if((Qo_dp_TZY[ll][0]>0)&&(Qo_dp_LD_00[ll][kk]>0))
{
if((Qo_dp_TZY[ll][0]<x)&&(Qo_dp_LD_00[ll][kk]<y))
ff_tmp3[ll][kk]=0.6*Qo_dp_TZY[ll][0]+0.4*Qo_dp_LD_00[ll][kk];
else
ff_tmp3[ll][kk]=AA;
}
if(ff_tmp3[ll][kk]>ff[ll][kk])
ff[ll][kk]=ff_tmp3[ll][kk];
}
x=ff[0][0];
for(l=0;l<U_st_num;l++)
for(k=0;k<D_st_num;k++)
if(x>ff[l][k])
{
x=ff[l][k];
path_u[amount_flood-2][0]=l;
path_d[amount_flood-2][0]=k;
}
z=0;
for(l=0;l<U_st_num;l++)
for(k=0;k<D_st_num;k++)
{
if(ff[l][k]==AA)
x3=0;
else
{
z=1;
x2=ff[l][k];
}
}
}
void CDpView::readdata_flood_TZY()
{
int i,j,k;
for(i=0;i<1000;i++)
{
Qin_TZY[i]=-1;
Qin_LD[i]=-1;
Qout_TZY[i]=-1;
Qout_LD[i]=-1;
}///////////////////////////将未赋值的洪水序列初始化为-1；实际洪水序列的长度由amount_flood决定
FILE *fp=fopen("桐子营洪水过程线.txt","r");
fscanf(fp,"%d",&amount_flood);
for (i=0;i<amount_flood;i++)
{
fscanf(fp,"%f",&Qin_TZY[i]);
}
fclose(fp);
fp=fopen("检验桐子营洪水过程线.txt","w");
fprintf(fp,"入库流量\n");
for (i=0;i<amount_flood;i++)
{
fprintf(fp,"%4.2f\n",Qin_TZY[i]);
}
fclose(fp);
}
float CDpView::HtoV_TZY(float h)
{
int i;
float H[15],V[15];
float v;
FILE *fp=fopen("桐子营水位库容曲线.txt","r");
for (i=0;i<15;i++)
{
fscanf(fp,"%f %f",&H[i],&V[i]);
}
fclose(fp);
if((h<590)||(h>625))//最小死水位，最大为所给范围上限
v=-100;//溢出
for(i=0;i<14;i++) //求h水位对应的库容v
if((h>=H[i])&&(h<H[i+1]))
v=V[i]+(V[i+1]-V[i])*(h-H[i])/(H[i+1]-H[i]);
return v;
}
float CDpView::HtoV_LD(float h)
{
int i;
float H[34],V[34];
float v;
FILE *fp=fopen("龙洞水位库容曲线.txt","r");
for (i=0;i<34;i++)
{
fscanf(fp,"%f %f",&H[i],&V[i]);
}
fclose(fp);
if((h<535.5)||(h>557))//最小死水位，最大为所给范围上限
v=-100;//溢出
for(i=0;i<33;i++)
if((h>=H[i])&&(h<H[i+1]))
v=V[i]+(V[i+1]-V[i])*(h-H[i])/(H[i+1]-H[i]);
return v;
}
float CDpView::VtoQo_TZY(float v)
{
int i;
float H[15],V[15];
float Qo;
float h;
FILE *fp=fopen("桐子营水位库容曲线.txt","r");
for (i=0;i<15;i++)
{
fscanf(fp,"%f %f",&H[i],&V[i]);
}
fclose(fp);
if((v<15038000)||(v>101142000))v=-1000;//溢出：最小为死库容，最大为所给范围上限
for(i=0;i<33;i++)
if((v>=V[i])&&(v<V[i+1]))
h=H[i]+(H[i+1]-H[i])*(v-V[i])/(V[i+1]-V[i]);//V插出H
Qo=HtoQo_TZY(h);
return Qo;
}
float CDpView::VtoQo_LD(float v)
{
int i;
float h;
float H[34],V[34];
float Qo;
FILE *fp=fopen("龙洞水位库容曲线.txt","r");
for (i=0;i<34;i++)
{
fscanf(fp,"%f %f",&H[i],&V[i]);
}
fclose(fp);
if((v<27230000)||(v>68480000))v=-1000;//溢出：最小为死库容，最大为所给范围上限
for(i=0;i<33;i++)
if((v>=V[i])&&(v<V[i+1]))
h=H[i]+(H[i+1]-H[i])*(v-V[i])/(V[i+1]-V[i]);//V插出H
Qo=HtoQo_LD(h);
return Qo;
}
float CDpView::HtoQo_TZY(float h)
{
int i;
float H[18],Qout[18];
float Qo;
FILE *fp=fopen("桐子营水位泄流曲线.txt","r");
for (i=0;i<18;i++)
{
fscanf(fp,"%f %f",&H[i],&Qout[i]);
}
fclose(fp);
if((h<608)||(h>620))//所给库泄曲线范围
Qo=-1000;//溢出
for(i=0;i<17;i++) //求h水位对应的泄流能力Qo
if((h>=H[i])&&(h<H[i+1]))
Qo=Qout[i]+(Qout[i+1]-Qout[i])*(h-H[i])/(H[i+1]-H[i]);
return Qo;
}
float CDpView::HtoQo_LD(float h)
{
int i;
float H[14],Qout[14];
float Qo;
FILE *fp=fopen("龙洞水位泄流曲线.txt","r");
for (i=0;i<14;i++)
{
fscanf(fp,"%f %f",&H[i],&Qout[i]);
}
fclose(fp);
if((h<544)||(h>555))//所给库泄曲线范围
Qo=-1000;//溢出
for(i=0;i<13;i++) //求h水位对应的泄流能力Qo
if((h>=H[i])&&(h<H[i+1]))
Qo=Qout[i]+(Qout[i+1]-Qout[i])*(h-H[i])/(H[i+1]-H[i]);
return Qo;
}
float CDpView::VtoH_TZY(float v)
{
int i;
float H[15],V[15];
float h;
FILE *fp=fopen("桐子营水位库容曲线.txt","r");
for (i=0;i<15;i++)
{
fscanf(fp,"%f %f",&H[i],&V[i]);
}
fclose(fp);
if((v<122000)||(v>101142000))//所给范围
h=-100;//溢出
for(i=0;i<14;i++) //求h水位对应的库容v
if((v>=V[i])&&(v<V[i+1]))
h=H[i]+(H[i+1]-H[i])*(v-V[i])/(V[i+1]-V[i]);
return h;
}
float CDpView::VtoH_LD(float v)
{
int i;
float H[34],V[34];
float h;
FILE *fp=fopen("龙洞水位库容曲线.txt","r");
for (i=0;i<34;i++)
{
fscanf(fp,"%f %f",&H[i],&V[i]);
}
fclose(fp);
if((v<27230000)||(v>68480000))//所给范围
h=-100;//溢出
for(i=0;i<33;i++)
if((v>=V[i])&&(v<V[i+1]))
h=H[i]+(H[i+1]-H[i])*(v-V[i])/(V[i+1]-V[i]);
return h;
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -