📄 rgwshjts.cpp
字号:
//日光温室热环境模拟程序的函数调试程序
//2007年9月
#include "stdio.h"
#include "math.h"
#include "rgwshjhs.cpp"
#include "rdwsbxhs.cpp"
//-------------------------------------------------------------------------------------------------------
//节点划分函数调试
/*main()
{ int i,i0;
double dw0[10]={0.010,0.011,0.012,0.013,0.014};
// double dw0[10]={0.011,0.012,0.013,0.014,0.015};
// double dw0[10]={0.011,0.001,0.013,0.014,0.015};
// double dw0[10]={0.011,0.001,0.0135,0.014,0.015};
// double dw0[10]={0.011,0.012,0.001,0.014,0.015};
nw=5;
// double dw0[10]={0.010,0.011,0.012,0.013};
// double dw0[10]={0.014,0.011,0.012,0.013};
// double dw0[10]={0.200,0.055,0.032,0.045};
// double dw0[10]={0.200,0.055,0.032,0.600};
// nw=4;
// double dw0[10]={0.004,0.002,0.006};
// double dw0[10]={0.006,0.002,0.004};
// double dw0[10]={0.024,0.006,0.024};
// double dw0[10]={0.240,0.060,0.240};
// double dw0[10]={0.240,0.060,0.120};
// nw=3;
// double dw0[10]={0.002,0.003};
// nw=2;
// double dw0[10]={0.003};
// double dw0[10]={0.004};
// double dw0[10]={0.005};
// nw=1;
n=jdhf(nw,dw0);
printf(" n=%d\n",n);
for(i=0;i<=n;i++)
{
printf(" %2d%7.4f",i,dx[i]);
if(fmod(i,5)==0)printf("\n");
}
scanf("%c",&i0);
}
//-------------------------------------------------------------------------------------------------------
*/
//-------------------------------------------------------------------------------------------------------
//一维导热差分方程组静态参数计算函数调试
/*main()
{ int i,i0;
dtao=60.0;
// double dw0[10]={0.014,0.011,0.012,0.013};
// double dw0[10]={0.200,0.055,0.032,0.600};
nw=4;
// double lmd0[]={0.81,0.032,0.83,0.054};
// double ro0[]={1801.0,22.0,1803.0,24.0};
// double c0[]={881.0,2002.0,883.0,2004.0};
double dw0[10]={0.240,0.060,0.240}; //对于该情况，如条件判断式中不采用e，则存放材料参数时会产生判断错误
// double dw0[10]={0.240,0.060,0.120};
nw=3;
double lmd0[]={0.81,0.032,0.83};
double ro0[]={1801.0,22.0,1803.0};
double c0[]={881.0,2002.0,883.0};
// double dw0[10]={0.080};
// nw=1;
// double lmd0[]={0.81};
// double ro0[]={1801.0};
// double c0[]={881.0};
n=jdhf(nw,dw0);
ywdrcfjc(n,dx,dw0,lmd0,ro0,c0,dtao);
printf(" n=%d\n",n); for(i=0;i<=n;i++)
{printf(" %2d%7.4f",i,dx[i]); if(fmod(i,5)==0)printf("\n");} scanf("%c",&i0);
printf(" lmd=\n"); for(i=0;i<=n;i++)
{printf(" %2d%8.4f",i,lmd[i]); if(fmod(i,5)==0)printf("\n");} scanf("%c",&i0);
printf(" ro=\n"); for(i=0;i<=n;i++)
{printf(" %2d%8.1f",i,ro[i]); if(fmod(i,5)==0)printf("\n");} scanf("%c",&i0);
printf(" c=\n"); for(i=0;i<=n;i++)
{printf(" %2d%8.1f",i,c[i]); if(fmod(i,5)==0)printf("\n");} scanf("%c",&i0);
}
//-------------------------------------------------------------------------------------------------------
*/
/*
//-------------------------------------------------------------------------------------------------------
//日期递进计算函数调试
main()
{ int i,j,i0;
AAI[1]=1; AAI[2]=31; AAD[16]=23.0; dtao=3040.0;
for(i=0; ;i++)
{
j=datecount(AAI[1],AAI[2],AAD[16],dtao);
printf(" %3d月 %3d日 %8.4f时\n",AAI[1],AAI[2],AAD[16]);
if(j==2) scanf("%c",&i0);
}
}
//-------------------------------------------------------------------------------------------------------
*/
//-------------------------------------------------------------------------------------------------------
//计算太阳辐射照度函数调试
/*main()
{ int i,j,i0;
bw=40.0; dj=120.0; bmfw=0.0; bmqx=0.0;
month=12; date=21; timeb=4.0; dtao=1800.0;
for(i=0; ;i++)
{
j=datecount(month,date,timeb,dtao);
month=AAI[1]; date=AAI[2]; timeb=AAD[16];
IDS=TYFZ(bw,dj,month,date,timeb,bmfw,bmqx);
printf("北纬%6.2f° 东经%6.2f°表面方位角%6.2f°表面倾斜角%6.2f° %2d月%2d日\n",bw,dj,bmfw,bmqx,month,date);
printf("%5.2f时 时角%6.2f°太阳赤纬%6.2f°高度%6.2f°方位%6.2f°大气透明度%5.2f\n",timeb,AAD[1],AAD[2],AAD[3]/3.1416*180.0,AAD[4]/3.1416*180.0,AAD[5]);
printf("太阳辐射 I0=%7.2f IDN=%6.2f ID=%6.2f ISH=%6.2f IS=%6.2f I=%7.2f\n\n",AAD[6],AAD[7],AAD[8],AAD[9],AAD[10],IDS);
scanf("%c",&i0);
}
}
*/
//-------------------------------------------------------------------------------------------------------
//网格划分函数与二维静态参数计算的调试
main()
{ int i,j,i0;
// double dw0[10]={0.006,0.002,0.004};
// double dw0[10]={0.024,0.006,0.024};
double dw0[10]={0.240,0.060,0.240};
// double dw0[10]={0.240,0.060,0.120};
nw=3; kd=8.0; dtao=36.0;
wghf(kd,nw,dw0);
nn=31; mm=12;
lmdw[0] =0.800; lmdx[0] =0.030; lmdn[0] =0.600;
lmdw[1] =0.801; lmdx[1] =0.031; lmdn[1] =0.601;
lmdw[2] =0.802; lmdx[2] =0.032; lmdn[2] =0.602;
lmdw[3] =0.803; lmdx[3] =0.033; lmdn[3] =0.603;
lmdw[4] =0.804; lmdx[4] =0.034; lmdn[4] =0.604;
lmdw[5] =0.805; lmdx[5] =0.035; lmdn[5] =0.605;
lmdw[6] =0.806; lmdx[6] =0.036; lmdn[6] =0.606;
lmdw[7] =0.807; lmdx[7] =0.037; lmdn[7] =0.607;
lmdw[8] =0.808; lmdx[8] =0.038; lmdn[8] =0.608;
lmdw[9] =0.809; lmdx[9] =0.039; lmdn[9] =0.609;
lmdw[10]=0.800; lmdx[10]=0.030; lmdn[10]=0.600;
lmdw[11]=0.801; lmdx[11]=0.031; lmdn[11]=0.601;
lmdw[12]=0.802; lmdx[12]=0.032; lmdn[12]=0.602;
row[0]=1800.0; rox[0]=30.0; ron[0]=1600.0;
row[1]=1800.1; rox[1]=30.1; ron[1]=1600.1;
row[2]=1800.2; rox[2]=30.2; ron[2]=1600.2;
row[3]=1800.3; rox[3]=30.3; ron[3]=1600.3;
row[4]=1800.4; rox[4]=30.4; ron[4]=1600.4;
row[5]=1800.5; rox[5]=30.5; ron[5]=1600.5;
row[6]=1800.6; rox[6]=30.6; ron[6]=1600.6;
row[7]=1800.7; rox[7]=30.7; ron[7]=1600.7;
row[8]=1800.8; rox[8]=30.8; ron[8]=1600.8;
row[9]=1800.9; rox[9]=30.9; ron[9]=1600.9;
row[10]=1800.0; rox[10]=30.0; ron[10]=1600.0;
row[11]=1800.1; rox[11]=30.1; ron[11]=1600.1;
row[12]=1800.2; rox[12]=30.2; ron[12]=1600.2;
cw[0]=800.0; cx[0]=3000.0; cn[0]=600.0;
cw[1]=800.1; cx[1]=3000.1; cn[1]=600.1;
cw[2]=800.2; cx[2]=3000.2; cn[2]=600.2;
cw[3]=800.3; cx[3]=3000.3; cn[3]=600.3;
cw[4]=800.4; cx[4]=3000.4; cn[4]=600.4;
cw[5]=800.5; cx[5]=3000.5; cn[5]=600.5;
cw[6]=800.6; cx[6]=3000.6; cn[6]=600.6;
cw[7]=800.7; cx[7]=3000.7; cn[7]=600.7;
cw[8]=800.8; cx[8]=3000.8; cn[8]=600.8;
cw[9]=800.9; cx[9]=3000.9; cn[9]=600.9;
cw[10]=800.0; cx[10]=3000.0; cn[10]=600.0;
cw[11]=800.1; cx[11]=3000.1; cn[11]=600.1;
cw[12]=800.2; cx[12]=3000.2; cn[12]=600.2;
ewdrcfjc(nn,mm,dtx,dty,lmdw,lmdx,lmdn,row,rox,ron,cw,cx,cn,dtao);
/* for(i=0;i<32;i++)
{ printf(" %6.4f",dtx[i]);if(i==5)printf("\n");if(i==9)printf("\n");
if(i==17)printf("\n"); if(i==25)printf("\n");if(i==31)printf("\n\n");
}
for(j=0;j<13;j++)
{ printf(" %6.4f",dty[j]);if(j==6)printf("\n");
}
*/
/* printf("北面外、墙下c=\n");for(j=0;j<13;j++){for(i=0; i<10;i++)printf("%8.2f",cs[i][j]);printf("\n");}
printf("室内北半部c=\n"); for(j=0;j<13;j++){for(i=10;i<18;i++)printf("%8.2f",cs[i][j]);printf("\n");}
printf("室内南半部c=\n"); for(j=0;j<13;j++){for(i=18;i<26;i++)printf("%8.2f",cs[i][j]);printf("\n");}
printf("南面外c=\n"); for(j=0;j<13;j++){for(i=26;i<32;i++)printf("%8.2f",cs[i][j]);printf("\n");}
*/
/* printf("北面外Dtx=\n");for(i=0; i<6; i++)printf("%8.4f",Dtx[i]);printf("\n");
printf("墙下Dtx=\n"); for(i=6; i<11;i++)printf("%8.4f",Dtx[i]);printf("\n");
printf("室内北Dtx=\n");for(i=11;i<19;i++)printf("%8.4f",Dtx[i]);printf("\n");
printf("室内南Dtx=\n");for(i=19;i<27;i++)printf("%8.4f",Dtx[i]);printf("\n");
printf("南面外Dtx=\n");for(i=27;i<33;i++)printf("%8.4f",Dtx[i]);printf("\n");
printf("\n深度方向Dty=\n");for(j=0;j<14;j++)printf("%8.4f",Dty[j]);printf(" \n");
*/
/* printf("北面外lmdim=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.6f",lmdim[i][j]);printf("\n");}
printf("墙下lmdim=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.6f",lmdim[i][j]);printf("\n");}
printf("室内北lmdim=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.6f",lmdim[i][j]);printf("\n");}
printf("室内南lmdim=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.6f",lmdim[i][j]);printf("\n");}
printf("南面外lmdim=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.6f",lmdim[i][j]);printf("\n");}
printf("北面外lmdip=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%8.5f",lmdip[i][j]);printf("\n");}
printf("墙下lmdip=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%8.5f",lmdip[i][j]);printf("\n");}
printf("室内北lmdip=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%8.5f",lmdip[i][j]);printf("\n");}
printf("室内南lmdip=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%8.5f",lmdip[i][j]);printf("\n");}
printf("南面外lmdip=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%8.5f",lmdip[i][j]);printf("\n");}
printf("北面外lmdjm=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.6f",lmdjm[i][j]);printf("\n");}
printf("墙下lmdjm=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.6f",lmdjm[i][j]);printf("\n");}
printf("室内北lmdjm=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.6f",lmdjm[i][j]);printf("\n");}
printf("室内南lmdjm=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.6f",lmdjm[i][j]);printf("\n");}
printf("南面外lmdjm=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.6f",lmdjm[i][j]);printf("\n");}
printf("北面外lmdjp=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%8.5f",lmdjp[i][j]);printf("\n");}
printf("墙下lmdjp=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%8.5f",lmdjp[i][j]);printf("\n");}
printf("室内北lmdjp=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%8.5f",lmdjp[i][j]);printf("\n");}
printf("室内南lmdjp=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%8.5f",lmdjp[i][j]);printf("\n");}
printf("南面外lmdjp=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%8.5f",lmdjp[i][j]);printf("\n");}
printf("北面外roc=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%8.0f",roc[i][j]);printf("\n");}
printf("墙下roc=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%8.0f",roc[i][j]);printf("\n");}
printf("室内北roc=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%8.0f",roc[i][j]);printf("\n");}
printf("室内南roc=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%8.0f",roc[i][j]);printf("\n");}
printf("南面外roc=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%8.0f",roc[i][j]);printf("\n");}
*/
/* printf("北面外AA=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.6f",AA[i][j]);printf("\n");}
printf("墙下AA=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.6f",AA[i][j]);printf("\n");}
printf("室内北AA=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.6f",AA[i][j]);printf("\n");}
printf("室内南AA=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.6f",AA[i][j]);printf("\n");}
printf("南面外AA=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.6f",AA[i][j]);printf("\n");}
printf("北面外BB=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.6f",BB[i][j]);printf("\n");}
printf("墙下BB=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.6f",BB[i][j]);printf("\n");}
printf("室内北BB=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.6f",BB[i][j]);printf("\n");}
printf("室内南BB=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.6f",BB[i][j]);printf("\n");}
printf("南面外BB=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.6f",BB[i][j]);printf("\n");}
printf("北面外KK=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.3f",KK[i][j]);printf("\n");}
printf("墙下KK=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.3f",KK[i][j]);printf("\n");}
printf("室内北KK=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.3f",KK[i][j]);printf("\n");}
printf("室内南KK=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.3f",KK[i][j]);printf("\n");}
printf("南面外KK=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.3f",KK[i][j]);printf("\n");}
printf("北面外PP=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.3f",PP[i][j]);printf("\n");}
printf("墙下PP=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.3f",PP[i][j]);printf("\n");}
printf("室内北PP=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.3f",PP[i][j]);printf("\n");}
printf("室内南PP=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.3f",PP[i][j]);printf("\n");}
printf("南面外PP=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.3f",PP[i][j]);printf("\n");}
*/
/* printf("北面外TS=\n");for(j=0;j<14;j++){for(i=0; i<6; i++)printf("%9.3f",TS[i][j]);printf("\n");}
printf("墙下TS=\n"); for(j=0;j<14;j++){for(i=6; i<11;i++)printf("%9.3f",TS[i][j]);printf("\n");}
printf("室内北TS=\n");for(j=0;j<14;j++){for(i=11;i<19;i++)printf("%9.3f",TS[i][j]);printf("\n");}
printf("室内南TS=\n");for(j=0;j<14;j++){for(i=19;i<27;i++)printf("%9.3f",TS[i][j]);printf("\n");}
printf("南面外TS=\n");for(j=0;j<14;j++){for(i=27;i<33;i++)printf("%9.3f",TS[i][j]);printf("\n");}
*/
/* printf("tfu=\n"); for(i=0;i<=32;i++)printf("%5.1f",tfu[i]);printf("\n");
printf("tfl=\n"); for(j=0;j<=13;j++)printf("%5.1f",tfl[j]);printf("\n");
printf("tfr=\n"); for(j=0;j<=13;j++)printf("%5.1f",tfr[j]);printf("\n");
printf("tfd=\n"); for(i=0;i<=32;i++)printf("%5.1f",tfd[i]);printf("\n");
printf("alfu=\n"); for(i=0;i<=32;i++)printf("%5.0f",alfu[i]);printf("\n");
printf("alfl=\n"); for(j=0;j<=13;j++)printf("%5.0f",alfl[j]);printf("\n");
printf("alfr=\n"); for(j=0;j<=13;j++)printf("%5.0f",alfr[j]);printf("\n");
printf("alfd=\n"); for(i=0;i<=32;i++)printf("%9.0f",alfd[i]);printf("\n");
*/
printf("Pbi0=\n"); for(i=1;i<=nn;i++)printf("%16.2f",Pbi0[i]);printf("\n");
printf("Pbim=\n"); for(i=1;i<=nn;i++)printf("%16.2f",Pbim[i]);printf("\n");
printf("Pbj0=\n"); for(j=1;j<=mm;j++)printf("%16.2f",Pbj0[j]);printf("\n");
printf("Pbjn=\n"); for(j=1;j<=mm;j++)printf("%16.2f",Pbjn[j]);printf("\n");
printf("Pc00=%12.6f Pcn0=%12.6f Pcm0=%12.6f Pcnm=%12.6f",Pc00,Pcn0,Pcm0,Pcnm);
scanf("%c",&i0);
}
//-------------------------------------------------------------------------------------------------------
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -