开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 行业应用 > 能源行业(电力石油煤炭) > 查看源码
PCL.cpp
资源名称：PCL.rar　[点击查看]
上传用户：hyjxc386
上传日期：2021-09-01
资源大小：2k
文件大小：7k
源码类别：

能源行业(电力石油煤炭)

开发平台：

Visual C++

PCL.cpp：源码内容
  1. #include<iostream.h>
  2. #include<iomanip.h>
  3. #include<math.h>
  4. #include<fstream.h>
  5. #include <stdio.h> 
  6. #include <stdlib.h>
  7. #define N 50
  8. void guss(double **a,int n)
  9. { 
  10. int i,j,k;double c; 
  11. for(k=0;k<n-1;k++)
  12. { 
  13. c=a[k][k]; 
  14. for(j=k;j<=n;j++) a[k][j]/=c; 
  15. for(i=k+1;i<n;i++) 
  16. { 
  17. c=a[i][k]; 
  18. for(j=k;j<=n;j++) a[i][j]-=c*a[k][j]; 
  19. } 
  20. } 
  21. a[n-1][n]/=a[n-1][n-1]; 
  22. for(k=n-2;k>=0;k--) 
  23. for(j=k+1;j<n;j++) a[k][n]-=a[k][j]*a[j][n]; 
  24. } 
  25. void main()
  26. {
  27. ifstream pfile;
  28. pfile.open("test.txt");
  29. if (!pfile)
  30. {cout<<"can't open!n"<<endl;}
  31. cout<<"%%%%%%潮流计算%%%%%%"<<endl;
  32. int n=1,m,sum,sumPQ=0,sumPV=0,sumdetaV=0;
  33. double k;
  34. int type[N];
  35. type[0]=0;
  36. double V[N],P[N],Q[N],P1[50];
  37. double x[50][50],r[N][N],bb[N][N],K[N][N],B1[N][N],mo[N][N],G1[N][N],KK[N][N];
  38. for(int i=0;i<=N;i++)
  39. {P[i]=0;Q[i]=0;P1[i]=0;V[i]=1;}
  40. for(i=0;i<N;i++)
  41. for(int j=0;j<N;j++)
  42. {x[i][j]=0;r[i][j]=0;bb[i][j]=0;K[i][j]=1;KK[i][j]=1;B1[i][j]=0;G1[i][j]=0;mo[i][j]=1;}
  43. i=1;
  44. for(n=1;n!=0;)
  45. {
  46.  pfile>>n;
  47.  if(n) 
  48.  {
  49.  pfile>>type[n];
  50.  if(type[n]==1) 
  51.  {pfile>>P[n]>>Q[n];sumPQ++;}
  52.  if(type[n]==2) 
  53.  {pfile>>P[n]>>k;sumPV++;}
  54.  if(type[n]==3) 
  55.  {pfile>>k>>k;sumdetaV++;}
  56.  sum=i++;
  57.  }
  58.  }
  59. for(n=1;n!=0;)
  60.  {
  61.  pfile>>n;
  62.  if(type[n]==2) pfile>>P1[n]>>V[n];
  63.  if(type[n]==3) pfile>>k>>V[n];
  64. }
  65. for(n=1;n!=0;)
  66. {
  67. pfile>>n;
  68. if(n)
  69. {
  70. pfile>>m;
  71. pfile>>r[n][m]>>x[n][m]>>bb[n][m]>>K[m][n];
  72. r[m][n]=r[n][m];x[m][n]=x[n][m];KK[n][m]=K[m][n];
  73. bb[m][n]=bb[n][m];
  74. }
  75. }
  76. for(n=1;n<=sum;n++)
  77. for(m=1;m<=sum;m++) 
  78. {
  79. mo[n][m]=x[n][m]*x[n][m]+r[n][m]*r[n][m];
  80. if(mo[n][m]==0) mo[n][m]=1;
  81. }
  82. for(n=1;n<=sum;n++)
  83. for(m=1;m<=sum;m++)
  84. {
  85. B1[n][n]+=bb[n][m]-x[n][m]*K[n][m]*K[n][m]/mo[n][m];
  86.     G1[n][n]+=r[n][m]*K[n][m]*K[n][m]/mo[n][m];
  87. }
  88. for(n=1;n<=sum;n++)
  89. for(m=1;m<=sum;m++)
  90. if(K[n][m]==1) {K[n][m]=K[m][n];KK[m][n]=-1;}
  91. for(n=1;n<=sum;n++)
  92. for(m=1;m<=sum;m++)
  93. {
  94. if(m!=n)
  95. {
  96. B1[n][m]+=x[n][m]*K[n][m]/mo[n][m];
  97. G1[n][m]+=-r[n][m]*K[n][m]/mo[n][m];
  98. }
  99. }
  100. double deta[N],detaP[N],detaV[N],PS[N],Pi[N],Ddeta[N],detaQ[N],Qi[N];
  101. double aa[N][N],cc[N][N];
  104. for(n=1;n<=sum;n++)
  105. {deta[n]=0;detaP[n]=0;PS[n]=0;Pi[n]=0;detaQ[n]=0;Qi[n]=0;}
  106. for(n=1;n<=sum-1;n++)
  107. PS[n]=-P[n]+P1[n];
  108. int Kp=1,Kq=1,maxP=0,maxQ=0;
  109. k=0;
  110. cout<<"**********************************************************************"<<endl;
  111. P:for(n=1;n<=sum-1;n++)
  112. {
  113. for(m=1;m<=sum;m++)
  114. Pi[n]+=V[m]*(G1[n][m]*cos(deta[n]-deta[m])+B1[n][m]*sin(deta[n]-deta[m]));
  115. Pi[n]*=V[n];
  116. }
  117. double maxdetaP=0,maxdetaQ=0;
  118. for(n=1;n<=sum-1;n++)
  119. {
  120. detaP[n]=PS[n]-Pi[n];
  121. Pi[n]=0;
  122. }
  123.     for(n=0;n<sum-1;n++)
  124. for(m=0;m<sum-1;m++) aa[n][m]=-B1[n+1][m+1];
  125. for(n=0;n<sum;n++)
  126. aa[n][sum-1]=detaP[n+1]/V[n+1];
  127. double *a[N];
  128. for (n=0;n<sum-1;n++)
  129. a[n]=aa[n];
  130. for(n=1;n<=sum-1;n++)
  131. if(fabs(maxdetaP)<fabs(detaP[n])) {maxdetaP=detaP[n];maxP=n;}
  132. if(fabs(maxdetaP)<1e-5) 
  133. {
  134. Kp=0;
  135. if(Kq==0) goto R;
  136. if(Kq!=0) goto Q;
  137. }
  138. guss(a,sum-1);
  139. for(n=0;n<sum-1;n++)
  140. {Ddeta[n+1]=aa[n][sum-1]/V[n+1];deta[n+1]=deta[n+1]+Ddeta[n+1];}
  141. Kq=1;
  142. Q:for(n=1;n<=sumPQ;n++)
  143. {
  144. for(m=1;m<=sum;m++)
  145. Qi[n]+=V[m]*(G1[n][m]*sin(deta[n]-deta[m])-B1[n][m]*cos(deta[n]-deta[m]));
  146. Qi[n]*=V[n];
  147. }
  148.  
  149. for(n=1;n<=sumPQ;n++)
  150. {
  151. detaQ[n]=-Q[n]-Qi[n];
  152. Qi[n]=0;
  153. }
  154. for(n=0;n<sumPQ;n++)
  155. for(m=0;m<sumPQ;m++) cc[n][m]=-B1[n+1][m+1];
  156. for(n=0;n<sumPQ;n++)
  157. cc[n][sumPQ]=detaQ[n+1]/V[n+1];
  158.  for(n=1;n<=sumPQ;n++)
  159.   ;
  160. double *c[N];
  161. for (n=0;n<sumPQ;n++)
  162. c[n]=cc[n];
  163. for(n=1;n<=sumPQ;n++)
  164. if(fabs(maxdetaQ)<fabs(detaQ[n])) {maxdetaQ=detaQ[n];maxQ=n;}
  165. if(fabs(maxdetaQ)<1e-5) 
  166. {
  167. Kq=0;
  168. if(Kp==0) goto R;
  169.     if(Kq!=0) goto T;
  170. }
  171. guss(c,sumPQ);
  172. for(n=0;n<sumPQ;n++)
  173. {detaV[n+1]=cc[n][sumPQ];V[n+1]=V[n+1]+detaV[n+1];}
  174. T: if(fabs(maxdetaQ)<fabs(maxdetaP)) {maxdetaQ=maxdetaP;maxQ=maxP;
  175. cout<<"第"<<k<<"次迭代后误差:"<<setw(15)<<maxdetaQ<<" 最大的功率误差在节点:"<<maxQ<<" 最大误差是detaP"<<endl;}
  176. else
  177. cout<<"第"<<k<<"次迭代后误差:"<<setw(15)<<maxdetaQ<<" 最大的功率误差在节点:"<<maxQ<<" 最大误差是detaQ"<<endl;
  178. k=k+1;
  179.     goto P;
  180. R:cout<<"*******************************************************************************"<<endl;
  181.   cout<<"节点数据:"<<endl;
  182.   for(n=1;n<=sum;n++)
  183. {
  184. cout<<"节点:"<<n<<"电压:"<<setw(10)<<V[n]<<"  相角:"<<setw(12)<<deta[n]*180/3.1415926;//转换为角度制
  185. if(type[n]==1) cout<<"发电有功:"<<setw(10)<<"0.0"<<"  发电无功:"<<setw(10)<<"0.0"<<"负荷有功:"<<setw(10)<<-PS[n]<<"  负荷无功:"<<setw(10)<<Q[n]<<endl;//PQ节点
  186. else if(type[n]==2)
  187. {
  188. for(m=1;m<=sum;m++)
  189. Q[n]+=V[m]*(G1[n][m]*sin(deta[n]-deta[m])-B1[n][m]*cos(deta[n]-deta[m]));
  190. Q[n]*=V[n];
  191. cout<<"发电有功:"<<setw(10)<<PS[n]<<"  发电无功:"<<setw(10)<<Q[n]<<"负荷有功:"<<setw(10)<<"0.0"<<"  负荷无功:"<<setw(10)<<"0.0"<<endl;
  192. }
  193. else
  194. {
  195. for(m=1;m<=sum;m++)
  196. PS[n]+=V[m]*(G1[n][m]*cos(deta[n]-deta[m])+B1[n][m]*sin(deta[n]-deta[m]));
  197. PS[n]*=V[n];
  198. for(m=1;m<=sum;m++)
  199. Q[n]+=V[m]*(G1[n][m]*sin(deta[n]-deta[m])-B1[n][m]*cos(deta[n]-deta[m]));
  200. Q[n]*=V[n];
  201. cout<<"发电有功:"<<setw(10)<<PS[n]<<"  发电无功:"<<setw(10)<<Q[n]<<"负荷有功:"<<setw(10)<<"0.0"<<"  负荷无功:"<<setw(10)<<"0.0"<<endl;//输出浮点型的0.0
  202. }
  203. }
  204.   cout<<"-------------------------------------------------------------------------------"<<endl;
  205.   cout<<"支路数据:"<<endl;
  206.   double LP[N][N],LQ[N][N];
  207.   for(n=1;n<N;n++)
  208.   for(m=n;m<N;m++)
  209.   LP[n][m]=LP[m][n]=LQ[n][m]=LQ[m][n]=0;
  210.   cout<<"支路功率流动"<<endl;
  211.   for(n=1;n<=sum;n++)
  212. for(m=1;m<=sum;m++)
  213. {
  214. if(r[n][m]||x[n][m]||bb[n][m])
  215. {
  216. LP[n][m]=-V[n]*V[n]*((1-K[n][m])*G1[n][m]/KK[n][m])-V[n]*cos(deta[n])*(G1[n][m]*(V[n]*cos(deta[n])-V[m]*cos(deta[m]))-B1[n][m]*(V[n]*sin(deta[n])-V[m]*sin(deta[m])))
  217. -V[n]*sin(deta[n])*(G1[n][m]*(V[n]*sin(deta[n])-V[m]*sin(deta[m]))+B1[n][m]*(V[n]*cos(deta[n])-V[m]*cos(deta[m])));
  218. LQ[n][m]=V[n]*V[n]*(-bb[n][m]+(1-K[n][m])*B1[n][m]/KK[n][m])-V[n]*sin(deta[n])*(G1[n][m]*(V[n]*cos(deta[n])-V[m]*cos(deta[m]))-B1[n][m]*(V[n]*sin(deta[n])-V[m]*sin(deta[m])))
  219. -V[n]*cos(deta[n])*(-G1[n][m]*(V[n]*sin(deta[n])-V[m]*sin(deta[m]))-B1[n][m]*(V[n]*cos(deta[n])-V[m]*cos(deta[m])));
  220. cout<<n<<"到"<<m<<"的有功:"<<setw(10)<<LP[n][m]<<" ";
  221. cout<<n<<"到"<<m<<"的无功:"<<setw(10)<<LQ[n][m]<<endl;
  222. }
  223. }
  224.   cout<<"支路功率损耗"<<endl;
  225.   for(n=1;n<=sum;n++)
  226. for(m=n+1;m<=sum;m++)
  227. {
  228. if(LP[n][m]&&LP[m][n]) cout<<"节点"<<n<<"到"<<"节点"<<m<<"有功损耗"<<setw(15)<<fabs(LP[n][m]+LP[m][n])<<" ";
  229. if(LQ[n][m]&&LQ[m][n]) cout<<"节点"<<n<<"到"<<"节点"<<m<<"无功损耗"<<setw(15)<<fabs(LQ[n][m]+LQ[m][n])<<endl;
  230. }
  232.   cout<<"-------------------------------------------------------------------------------"<<endl;
  233.   cout<<"全网数据:"<<endl;
  234.   double GenerateP=0,LoadP=0,LossP,MaxV=0,MinV=100000;
  235.   int NumbermaxV=1,NumberminV=1;
  236.   for(n=1;n<=sum;n++)
  237.   {
  238. if(type[n]!=1) GenerateP+=PS[n];
  239. else LoadP+=-PS[n];
  240. if(V[n]>MaxV) {MaxV=V[n];NumbermaxV=n;}
  241. else if(V[n]<MinV) {MinV=V[n];NumberminV=n;}
  242.   }
  243.   LossP=GenerateP-LoadP;
  244.   cout<<"发电有功:"<<GenerateP<<" 负荷有功:"<<LoadP<<" 有功损耗:"<<LossP<<" 最高电压:"<<MaxV<<" 最高电压节点号:"<<NumbermaxV<<" 最低电压:"<<MinV<<" 最低电压节点号:"<<NumberminV<<endl;
  245. }