下载文档原格式
电气工程及其自动化专业（二学位）人才培养方案学科门类工学代码08
类别电气信息类代码0806
专业名称电气工程及其自动化代码080601 一、学制与学位
学制：二年
学位：第二学士学位
二、培养目标
培养基础扎实，知识面宽，具有较强的综合素质和一定的创新精神，能够从事电气工程及其自动化领域相关的工程设计、生产制造、系统运行、系统分析、技术开发、教育科研、经济管理等方面工作的特色鲜明的复合型高级工程技术人才。
三、专业培养基本要求
1.了解本专业领域的理论前沿和发展动态；
2.获得较好的工程实践训练，具有综合解决工程实际问题的能力；
3.具有较强的知识获取与运用能力，具备创新意识和从事科学研究、科技开发的能力；
4.具有较强的工作适应性、人际交往能力和团队协作精神，具有一定的组织管理才能。
四、学时与学分
五、专业主干课程
1.学科技术基础骨干课：工程制图、电路理论、模拟电子技术基础、数字电子技术基础、网络与通信技术
2.学科专业基础骨干课：电机学、工程电磁场、电力电子技术
3.学科专业骨干课：电力系统分析基础、电力系统暂态分析、电力系统继电保护原理、发电厂电气部分、高电压技术
Table of Teaching Schedule for Required Course
Table of Teaching Schedule for Specialty Elective Course
独立开设的实践教学环节进程
Table of Teaching Schedule for Major Practical Training。
《电力系统潮流上机》课程设计报告院系班级：学号：学生姓名：指导教师：设计周数成绩：日期：年月日一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细内容见附录）1.手算： 要求应用牛顿-拉夫逊法或P-Q 分解法手算求解，要求精度为0.001MW 。
节点1为平衡节点，电压︒∠=00.11U ，节点2为PQ 节点，负荷功率6.08.0~2j S +=，节点3是PV 节点，1.1,4.033==U P ，两条支路分别为04.001.013j Z +=，2.005.012j Z +=，对地支路33.030j y =。
2.计算机计算：编写潮流计算程序，要求如下：2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器模型，电压等级的归算，标幺值的计算；2.2基础数据的计算机存储：节点数据，支路数据（包括变压器）； 2.3用牛顿-拉夫逊法计算；2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV 节点无功功率和平衡节点功率，数据输出模块； 2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图简单系统如下图所示，支路数据如下：41.01.012j z +=，3.013j z =，5.012.014j z +=，40.008.024j z +=01528.01,202,10j y y ==，0192.01,404,10j y y ==，01413.02,404,20j y y == 1.1=k节点数据如下：18.030.01j S --= ，13.055.02j S --= ， 5.03=S ，10.13=U ，o U 005.14∠=1)节点导纳阵#include <stdio.h> #include <math.h> #include <fstream.h> #include "LF.h"//form node conductance matrixintMakeY( intnB, intnL, Line* sL, double** YG, double** YB ) { inti,j,l;double r,x,d1,g,b,t; for(i=0;i<nB;i++) for(j=0;j<nB;j++) { YG[i][j]=0.0;YB[i][j]=0.0;} for(i=0;i<nL;i++) {r=sL[i].R; x=sL[i].X; g=r/(r*r+x*x); b=-x/(r*r+x*x); switch(sL[i].Type){case 1://Linebreak;case 2://Transformerg*=1/sL[i].K;b*=1/sL[i].K;break;}YG[sL[i].NumI][sL[i].NumI]+=g;YG[sL[i].NumJ][sL[i].NumJ]+=g;YG[sL[i].NumI][sL[i].NumJ]-=g;YG[sL[i].NumJ][sL[i].NumI]-=g;YB[sL[i].NumI][sL[i].NumI]+=b+sL[i].B;YB[sL[i].NumJ][sL[i].NumJ]+=b+sL[i].B;YB[sL[i].NumI][sL[i].NumJ]-=b;YB[sL[i].NumJ][sL[i].NumI]-=b; }printf("实部：\n");for(i=0;i<nB;i++){for(j=0;j<nB;j++)printf("%lf\t",YG[i][j]);printf("\n");}printf("虚部：\n");for(i=0;i<nB;i++){for(j=0;j<nB;j++)printf("%lf\t",YB[i][j]);printf("\n");}/* Check the Y matrix */ofstreamfout("out.txt");fout<< "--------------Y Matrix--------------------" <<endl;for(i=0;i<nB;i++){for(j=0;j<nB;j++)fout<< YG[i][j] << "+j" << YB[i][j] << "\t";fout<<endl;}fout.close();return 0;}2）计算功率不平衡量#include <stdio.h>#include <math.h>#include <fstream.h>#include "LF.h"//form delta p and delta qintCalDeltaPQ( intnpv, intnpq, Bus* bus, double** YG, double** YB, int* p_Jtobus, double* deltaf ){intk,i,j;for(k=0;k<npv+npq*2;k++){ i=p_Jtobus[k];if(k<npv){ deltaf[k]=bus[i].GenP-bus[i].LoadP;for(j=0;j<npv+npq+1;j++){deltaf[k]-=bus[i].Volt*bus[j].Volt*(YG[i][j]*cos(bus[i].Phase-bus[j].Phase)+YB[i][j ]*sin(bus[i].Phase-bus[j].Phase));}printf("PV节点%d的有功功率是%lf\n",i,deltaf[k]);}if(k<npq+npv&&k>=npv){ deltaf[k]=bus[i].GenP-bus[i].LoadP;for(j=0;j<npv+npq+1;j++){deltaf[k]-=bus[i].Volt*bus[j].Volt*(YG[i][j]*cos(bus[i].Phase-bus[j].Phase)+YB[i][j ]*sin(bus[i].Phase-bus[j].Phase));}printf("PQ节点%d的有功功率是%lf\n",i,deltaf[k]);}if(k<npq*2+npv&&k>=npv+npq){deltaf[k]=bus[i].GenQ-bus[i].LoadQ;for(j=0;j<npv+npq+1;j++){deltaf[k]-=bus[i].Volt*bus[j].Volt*(YG[i][j]*sin(bus[i].Phase-bus[j].Phase)-YB[i][j ]*cos(bus[i].Phase-bus[j].Phase));}printf("PQ节点%d的无功功率是 %lf\n",i,deltaf[k]);}}return 0;}3）雅各比矩阵的计算/*Purpose: for undergraduate courseTask: Load FlowCopyright @ NCEPU, Liu Chongru*/#include <stdio.h>#include <math.h>#include <fstream.h>#include "LF.h"//form Jacobian matrixintFormJacobian( intnpv, intnpq, Bus* bus, double** YG, double** YB, int* p_Jtobus, double** Jac ){ intnp = npv+npq,j,k,i,m;//TODOdouble a[14],q[14];for(k=0;k<npv+npq*2;k++){i=p_Jtobus[k];a[i]=0;q[i]=0;if(k<np)//H N{for(j=0;j<np+1;j++)if(j!=i){a[i]+=bus[j].Volt*(YG[i][j]*sin(bus[i].Phase-bus[j].Phase)-YB[i][j]*cos(bus[i].Phase-bu s[j].Phase));q[i]+=bus[j].Volt*(YG[i][j]*cos(bus[i].Phase-bus[j].Phase)+YB[i][j]*sin(bus[i].Phase-bu s[j].Phase));}for(m=0;m<npv+npq*2;m++){j=p_Jtobus[m];if(j!=i){if(m<np)Jac[k][m]=bus[i].Volt*bus[j].Volt*(YG[i][j]*sin(bus[i].Phase-bus[j].Phase)-YB[i][j]*cos (bus[i].Phase-bus[j].Phase));//Form HelseJac[k][m]=bus[i].Volt*bus[j].Volt*(YG[i][j]*cos(bus[i].Phase-bus[j].Phase)+YB[i][j]*sin (bus[i].Phase-bus[j].Phase));//Form N}else if(j==i){if(m<np)Jac[k][m]=-bus[i].Volt*a[i];//Form HelseJac[k][m]=bus[i].Volt*q[i]+2*bus[i].Volt*bus[i].Volt*YG[i][j];//Form N}}}else{for(j=0;j<np+1;j++)if(j!=i){a[i]+=bus[j].Volt*(YG[i][j]*sin(bus[i].Phase-bus[j].Phase)-YB[i][j]*cos(bus[i].Phase-bu s[j].Phase));q[i]+=bus[j].Volt*(YG[i][j]*cos(bus[i].Phase-bus[j].Phase)+YB[i][j]*sin(bus[i].Phase-bu s[j].Phase));}for(m=0;m<npv+npq*2;m++){j=p_Jtobus[m];if(j!=i){if(m<np)Jac[k][m]=-bus[i].Volt*bus[j].Volt*(YG[i][j]*cos(bus[i].Phase-bus[j].Phase)+YB[i][j]*si n(bus[i].Phase-bus[j].Phase)); //Form JelseJac[k][m]=bus[i].Volt*bus[j].Volt*(YG[i][j]*sin(bus[i].Phase-bus[j].Phase)-YB[i][j] *cos(bus[i].Phase-bus[j].Phase)); //Form L }else if(j==i){ if(m<np)Jac[k][m]=bus[i].Volt*q[i];elseJac[k][m]=bus[i].Volt*a[i]-2*bus[i].Volt*bus[i].Volt*YB[i][j];} }}}for(i=0;i<np+npq;i++){for(int j=0;j<np+npq;j++){printf("%d %d %f ",i,j,Jac[i][j]);}printf("\n");}//Output the matrix to check the Jacobian matrixofstreamfout("out.txt",ios::app);fout<< "--------------Jacobian Matrix--------------------" <<endl;for(i=0; i<np+npq;i++ ){for(j=0; j<np+npq; j++ )fout<<Jac[i][j] << "\t";fout<<endl;}fout.close();return 0;}4）线路损耗//8.calculate the power flowdouble* p_Pij, *p_Qij, *p_Pji, *p_Qji;p_Pij = new double[nL];p_Qij = new double[nL];p_Pji = new double[nL];p_Qji = new double[nL];int x1,x2;for( i=0; i<nL; i++ ){ x1=line[i].NumI;x2=line[i].NumJ;if(line[i].Type==1){p_Pij[i]=bus[x1].Volt*bus[x1].Volt*(-YG[x1][x2])-bus[x1].Volt*bus[x2].Volt*((-YG[x1][x2 ])*cos(bus[x1].Phase-bus[x2].Phase)+(-YB[x1][x2])*sin(bus[x1].Phase-bus[x2].Phase));p_Qij[i]=-bus[x1].Volt*bus[x1].Volt*(line[i].B+(-YB[x1][x2]))-bus[x1].Volt*bus[x2]. Volt*((-YG[x1][x2])*sin(bus[x1].Phase-bus[x2].Phase)-(-YB[x1][x2])*cos(bus[x1].Phase-bu s[x2].Phase));p_Pji[i]=bus[x2].Volt*bus[x2].Volt*(-YG[x2][x1])-bus[x2].Volt*bus[x1].Volt*((-YG[x2][x1 ])*cos(bus[x2].Phase-bus[x1].Phase)+(-YB[x2][x1])*sin(bus[x2].Phase-bus[x1].Phase));p_Qji[i]=-bus[x2].Volt*bus[x2].Volt*(line[i].B+(-YB[x2][x1]))-bus[x2].Volt*bus[x1]. Volt*((-YG[x2][x1])*sin(bus[x2].Phase-bus[x1].Phase)-(-YB[x2][x1])*cos(bus[x2].Phase-bu s[x1].Phase));}else{p_Pij[i]=bus[x1].Volt*bus[x1].Volt*(-YG[x1][x2])/line[i].K-bus[x1].Volt*bus[x2].Vol t*((-YG[x1][x2])*cos(bus[x1].Phase-bus[x2].Phase)+(-YB[x1][x2])*sin(bus[x1].Phase-bus[x 2].Phase));p_Qij[i]=-bus[x1].Volt*bus[x1].Volt*((-YB[x1][x2])/line[i].K+line[i].B)-bus[x1].Volt*bus[x2].Volt*((-YG[x1][x2])*sin(bus[x1].Phase-bus[x2].Phase)-(-YB[x1][x2])*cos(bus[x1].Ph ase-bus[x2].Phase));p_Pji[i]=bus[x2].Volt*bus[x2].Volt*(-YG[x2][x1]*line[i].K)-bus[x2].Volt*bus[x1].Volt*((-YG[x2][x1])*cos(bus[x2].Phase-bus[x1].Phase)+(-YB[x2][x1])*sin(bus[x2].Phase-bus[x1].P hase));p_Qji[i]=-bus[x2].Volt*bus[x2].Volt*((-YB[x2][x1])*line[i].K+line[i].B)-bus[x2].Volt*bu s[x1].Volt*((-YG[x2][x1])*sin(bus[x2].Phase-bus[x1].Phase)-(-YB[x2][x1])*cos(bus[x2].Ph ase-bus[x1].Phase)); }}//p and q of PH bus and PV bus int s=0;double p[9],q[9],Ps[9],Qs[9],PS=0,QS=0; for( i=0; i<nB; i++ ){p[i]=0;q[i]=0; for(int j=0; j<nB; j++ ){p[i]+=(bus[j].Volt*(YG[i][j])*cos(bus[j].Phase)-bus[j].Volt*(YB[i][j])*sin(bus[j].Phas e));q[i]-=(bus[j].Volt*(YG[i][j])*sin(bus[j].Phase)+bus[j].Volt*(YB[i][j])*cos(bus[j].Phase ));}Ps[i]=bus[i].Volt*cos(bus[i].Phase)*p[i]-bus[i].Volt*sin(bus[i].Phase)*q[i]; Qs[i]=bus[i].Volt*cos(bus[i].Phase)*q[i]+bus[i].Volt*sin(bus[i].Phase)*p[i];} for(i=0;i<nB;i++) { PS+=Ps[i]; QS+=Qs[i];}printf("PS=%7.7f,QS=%7.7f\n",PS,QS);}//lossdoublePloss=0, Qloss=0; for( i=0; i<nB; i++ ){Ploss+=p_Pij[i]+p_Pji[i]; Qloss+=p_Qij[i]+p_Qji[i];} 5）程序流图如下6--------------Y Matrix--------------------0+j-17.3611 0+j0 0+j0 0+j17.3611 0+j0 0+j0 0+j0 0+j0 0+j00+j0 0+j-16 0+j0 0+j0 0+j0 0+j0 0+j16 0+j0 0+j00+j0 0+j0 0+j-17.0648 0+j0 0+j0 0+j0 0+j0 0+j0 0+j17.06480+j17.3611 0+j0 0+j0 3.30738+j-39.3089 -1.36519+j11.6041 -1.94219+j10.5107 0+j0 0+j0 0+j0 0+j0 0+j0 0+j0 -1.36519+j11.6041 2.55279+j-17.3382 0+j0 -1.1876+j5.97513 0+j0 0+j00+j0 0+j0 0+j0 -1.94219+j10.5107 0+j0 3.2242+j-15.8409 0+j0 0+j0 -1.28201+j5.588240+j0 0+j16 0+j0 0+j0 -1.1876+j5.97513 0+j0 2.80473+j-35.4456 -1.61712+j13.698 0+j00+j0 0+j0 0+j0 0+j0 0+j0 0+j0 -1.61712+j13.698 2.77221+j-23.3032 -1.15509+j9.784270+j0 0+j0 0+j17.0648 0+j0 0+j0 -1.28201+j5.58824 0+j0 -1.15509+j9.78427 2.4371+j-32.1539 --------------Jacobian Matrix--------------------16.4 0 0 0 0 -16.4 0 0 0 0 0 0 0 00 17.4915 0 0 0 0 0 -17.4915 0 0 0 0 0 00 0 40.1703 -11.6041 -10.5107 0 0 0 3.30738 -1.36519 -1.94219 0 0 00 0 -11.6041 17.5792 0 -5.97513 0 0 -1.36519 2.55279 0 -1.1876 0 00 0 -10.5107 0 16.0989 0 0 -5.58824 -1.94219 0 3.2242 0 0 -1.28201-16.4 0 0 -5.97513 0 36.0731 -13.698 0 0 -1.1876 0 2.80473 -1.61712 00 0 0 0 0 -13.698 23.4822 -9.78427 0 0 0 -1.61712 2.77221 -1.155090 -17.4915 0 0 -5.58824 0 -9.78427 32.864 0 0 -1.28201 0 -1.15509 2.43710 0 -3.30738 1.36519 1.94219 0 0 0 38.4474 -11.6041 -10.5107 0 0 00 0 1.36519 -2.55279 0 1.1876 0 0 -11.6041 17.0972 0 -5.97513 0 00 0 1.94219 0 -3.2242 0 0 1.28201 -10.5107 0 15.5829 0 0 -5.588240 0 0 1.1876 0 -2.80473 1.61712 0 0 -5.97513 0 34.8181 -13.698 00 0 0 0 0 1.61712 -2.77221 1.15509 0 0 0 -13.698 23.1242 -9.784270 0 0 0 1.28201 0 1.15509 -2.4371 0 0 -5.58824 0 -9.78427 31.4437--------------Jacobian Matrix--------------------16.9269 0 0 0 0 -16.9269 0 0 0 0 0 1.68793 0 00 18.1691 0 0 0 0 0 -18.1691 0 0 0 0 0 0.8836270 0 41.9297 -12.1301 -11.1536 0 0 0 3.54272 -1.0628 -1.76646 0 0 00 0 -12.0455 18.0609 0 -6.01539 0 0 -1.78138 1.30819 0 -2.10262 0 00 0 -11.0484 0 16.8144 0 0 -5.76607 -2.33608 0 2.42598 0 0 -1.97778-16.9269 0 0 -6.36224 0 37.9476 -14.6585 0 0 -0.357534 0 3.05959 -0.930027 00 0 0 0 0 -14.4721 24.8873 -10.4152 0 0 0 -2.509 1.86088 -1.473890 -18.1691 0 0 -6.05157 0 -10.4721 34.6928 0 0 -0.733327 0 -0.9919732.66270 0 -3.52149 1.0628 1.76646 0 0 0 42.0299 -12.1301 -11.1536 0 0 00 0 1.78138 -3.884 0 2.10262 0 0 -12.0455 17.2037 0 -6.01539 0 00 0 2.33608 0 -4.31386 0 0 1.97778 -11.0484 0 16.2993 0 0 -5.766071.68793 0 0 0.357534 0 -2.97549 0.930027 0 0 -6.36224 0 38.3226 -14.6585 00 0 0 0 0 2.509 -3.98289 1.47389 0 0 0 -14.4721 24.2355 -10.41520 0.883627 0 0 0.733327 0 0.991973 -2.60893 0 0 -6.05157 0 -10.4721 34.8585 --------------Jacobian Matrix--------------------16.7457 0 0 0 0 -16.7457 0 0 0 0 0 1.63043 0 00 18.0388 0 0 0 0 0 -18.0388 0 0 0 0 0 0.8501960 0 41.3695 -11.8919 -10.9686 0 0 0 3.48069 -1.02775 -1.73712 0 0 00 0 -11.8057 17.6918 0 -5.8861 0 0 -1.7602 1.28091 0 -2.0217 0 00 0 -10.8651 0 16.5476 0 0 -5.68251 -2.29737 0 2.40655 0 0 -1.91027-16.7457 0 0 -6.21183 0 37.3041 -14.3465 0 0 -0.382862 0 2.95313 -0.937485 00 0 0 0 0 -14.1704 24.4052 -10.2348 0 0 0 -2.42909 1.86079 -1.433530 -18.0388 0 0 -5.94693 0 -10.2872 34.273 0 0 -0.757656 0 -0.9892082.598470 0 -3.48089 1.02775 1.73712 0 0 0 41.3703 -11.8919 -10.9686 0 0 00 0 1.7602 -3.78189 0 2.0217 0 0 -11.8057 16.6941 0 -5.8861 0 00 0 2.29737 0 -4.20764 0 0 1.91027 -10.8651 0 15.9488 0 0 -5.682511.63043 0 0 0.382862 0 -2.95077 0.937485 0 0 -6.21183 0 37.3083 -14.3465 00 0 0 0 0 2.42909 -3.86262 1.43353 0 0 0 -14.1704 23.7059 -10.23480 0.850196 0 0 0.757656 0 0.989208 -2.59706 0 0 -5.94693 0 -10.2872 34.2743 --------------Jacobian Matrix--------------------16.7435 0 0 0 0 -16.7435 0 0 0 0 0 1.63 0 00 18.0374 0 0 0 0 0 -18.0374 0 0 0 0 0 0.850 0 41.3625 -11.8888 -10.9664 0 0 0 3.48016 -1.02713 -1.73662 0 0 00 0 -11.8026 17.6871 0 -5.8845 0 0 -1.76008 1.28053 0 -2.02045 0 00 0 -10.8628 0 16.5444 0 0 -5.68158 -2.29703 0 2.40632 0 0 -1.90929-16.7435 0 0 -6.20987 0 37.296 -14.3426 0 0 -0.383399 0 2.95114 -0.937742 00 0 0 0 0 -14.1667 24.3994 -10.2326 0 0 0 -2.42794 1.86097 -1.433020 -18.0374 0 0 -5.94567 0 -10.285 34.2681 0 0 -0.758139 0 -0.9892032.597340 0 -3.48016 1.02713 1.73662 0 0 0 41.3625 -11.8888 -10.9664 0 0 00 0 1.76008 -3.78053 0 2.02045 0 0 -11.8026 16.6871 0 -5.8845 0 00 0 2.29703 0 -4.20632 0 0 1.90929 -10.8628 0 15.9444 0 0 -5.681581.63 0 0 0.383399 0 -2.95114 0.937742 0 0 -6.20987 0 37.296 -14.3426 00 0 0 0 0 2.42794 -3.86097 1.43302 0 0 0 -14.1667 23.6994 -10.23260 0.85 0 0 0.758139 0 0.989203 -2.59734 0 0 -5.94567 0 -10.285 34.2681--------------iteration---------------iteration = 4--------------voltage magnitude and angle--------------------1.04 0 01.025 0.161967 9.280011.025 0.0814153 4.664761.02579 -0.0386902 -2.216790.995631 -0.0696178 -3.988811.01265 -0.0643572 -3.68741.02577 0.064921 3.71971.01588 0.0126979 0.7275371.03235 0.0343257 1.96672-------------bus P and Q------------------1 0.71641 0.2704592 1.63 0.06653663 0.85 -0.1085974 0 05 -1.25 -0.56 -0.9 -0.37 0 08 -1 -0.359 0 0----------------line flow--------------------NUM------i-------j-----------begin------------end--------------1 4 1 -0.71641+j-0.239232 0.71641+j0.270462 7 2 -1.63+j0.0917816 1.63+j0.06653653 9 3 -0.85+j0.149553 0.85+j-0.1085974 7 8 0.763799+j-0.00797398 -0.759046+j-0.1070415 9 8 0.241834+j0.0311946 -0.240954+j-0.2429586 7 5 0.866201+j-0.0838079 -0.843202+j-0.1131287 9 6 0.608166+j-0.180748 -0.594627+j-0.1345668 5 4 -0.406798+j-0.386872 0.409373+j0.2289319 6 4 -0.305372+j-0.165433 0.307037+j0.0102993--------------------Ploss and Qloss-----------------Ploss = 0.0471901 Qloss = -0.9574833.思考题3.1潮流计算的方法有哪些？各有何特点？答：潮流计算分为手算和机算两大类，其中机算又有高斯-赛德尔迭代法、牛顿-拉夫逊迭代法、P-Q 分解法等算法。
潮流能发电及潮流能发电装置戴庆忠摘要潮流能发电是利用潮汐动能的一种发电方式。
由于潮流能发电不需要筑坝拦水，具有对环境影响小等许多优点。
因此，近年来潮流能发电引起许多国家重视，潮流能发电技术发展很快。
本文从分析潮流能的特点入手，介绍了国内外潮流能发电的近况，重点介绍目前出现的各种潮流能发电装置，包括水平轴潮流能水轮机、竖井潮流能水轮机、振荡水翼式潮流能装置等。
关键词潮汐潮流能潮流能水轮机潮流能发电1前言1.1潮流能的特点潮流主要是指伴随潮汐现象而产生的有规律的海水流，潮流每天两次改变其大小和方向。
而潮流能发电则是直接利用涨落潮水的水流冲击叶轮等机械装置进行发电。
众所周知，潮汐是海水在月球、太阳等引力作用下形成的周期性海水涨落现象。
潮汐现象伴随两种运动形态：一是涨潮和落潮引起的海水垂直升降，即通常所指的潮汐；二是海水的水平运动，即潮流。
前者(海水垂直升降)所携带的能量(潮汐能)为势能；而后者所携带的能量(潮流能)为动能。
可以说，两者是与潮汐涨落相伴共生的孪生兄弟。
对前者，可以采用类似河川水力发电的方式，筑坝蓄水发电；而对本文所介绍的潮流能，可以采用类似于海流发电方式，利用潮流的动能发电。
与常规能源比较，潮流能有以下特点：(1) 潮流能是一种可再生的清洁能源。
(2) 潮流能的能量密度较低(但远大于风能和太阳能)，但总储量较大。
(3) 与海流能不同，潮流能是一种随时间、空间而变化的能源，但其变化有规律可循，并可提前预测预报。
(4) 潮流能发电不拦海建坝，且发电机组通常浸没在海中，对海洋生物影响较小，也不会对环境产生三废污染，不存在常规水电建设中头疼的占用农田、移民安置等诸多问题。
(5) 与陆地电力建设相比，潮流能开发环境恶劣，一次性投资大，设备费用高，安装维护和电力输送等都存在一系列关键技术问题。
1.2 潮流能水轮机输出功率的计算潮流能机组输出功率的计算公式为：a :AV 3p=一2式中P——功率，W3P ——海水密度，1025kg/mA ――潮流水轮机转子扫掠面积，卅V ――潮流速度，m/sn ----- 效率从上述可以看出，潮流能机组的输出功率很大程度决定于潮流速度。
潮流能发电原理潮流能发电是一种利用海洋潮汐能转化为电能的可再生能源技术。
潮汐是由于地球受到月球和太阳的引力作用而产生的规律性海洋水位变化。
潮汐是一种相对固定的自然现象，能够稳定地提供能量。
因此，利用潮汐能发电具有持续稳定的特点，被视为一种可靠的清洁能源。
潮流能发电的原理是利用海洋潮汐能转化为动力能，再将动力能转化为电能。
具体来说，主要包括潮汐动力装置和发电装置两部分。
潮汐动力装置是将潮汐能转化为机械能的装置。
一种常见的潮汐动力装置是潮流涡轮机。
潮流涡轮机通过将水流引入涡轮机，利用水流的动能带动涡轮旋转，进而驱动涡轮机内部的发电机发电。
涡轮机通常设在海底，以最大限度地利用潮流能。
发电装置是将机械能转化为电能的装置。
一般情况下，潮流涡轮机会通过传动系统将动力传递给内部的发电机。
发电机利用电磁感应原理，将旋转的机械能转化为电能。
由于潮流涡轮机的运行受到潮汐的影响，因此发电装置需要具备一定的调节和控制系统，以保证发电的稳定性和可控性。
潮流能发电具有一些独特的优势。
首先，潮汐能是一种可再生的能源，不会像化石能源一样耗尽。
其次，潮汐能是一种稳定的能源，潮汐的周期性和规律性使得发电系统能够提前计划和调度发电。
再次，潮流能发电是一种清洁能源，不会产生二氧化碳等温室气体和污染物，对环境没有直接的负面影响。
此外，潮流能发电的设备通常位于海洋中，不会占用大量陆地资源。
然而，潮流能发电也存在一些挑战和问题。
首先，由于潮汐能的利用需要特定的地理环境和条件，因此潮流能发电的可行性和适用性受到一定限制。
其次，潮流能发电设备的建设和维护成本较高，需要大量的投资。
此外，潮流能发电还面临着对海洋生态环境的影响和冲击，需要进行充分的环境评估和监测。
潮流能发电技术还处于发展初期，尚未广泛商业化应用。
但随着清洁能源的需求和发展，潮流能发电作为一种可持续的能源选择，具有重要的发展前景。
未来，随着技术的进一步突破和经济的进一步成熟，潮流能发电有望成为海洋能源领域的重要组成部分。
智能电网信息工程专业人才培养方案Undergraduate Program for Smart Grid InformationEngineering Major学科门类工学代码08Discipline Type: Engineering Code:08专业类电气类代码0806 Type: Electrical Information Code: 0806 专业名称智能电网信息工程代码080602T Title of the Major: Smart Grid Information Engineering Code: 080602T一、学制与学位Length of Schooling and Degree学制：四年Duration：Four years学位：工学学士Degree：Bachelor of Engineering二、培养目标Educational Objectives培养具有扎实的专业理论和专业技能，具备较强的综合素质和一定的创新精神，掌握信息采集和处理的基本理论和电力系统通信技术，掌握电力系统生产、运行的规律和特点，并对智能电网体系结构和关键技术有一定认识，可以在信息化、自动化、互动化的电力系统领域从事研究、开发、设计、制造、运行维护与管理等工作的复合型高级工程技术人才。
To cultivate advanced and vertisale enginerring talents with solid theoretical knowledge and professional skills, having strong comprehensive qualities and a certain spirit of innovation, mastersing basic theory of information collection and processing, power system communication techniques and the characteristics of power production and operation. At the same time, the architecture and key techniques of smart grid are also learned. Upon the graduation, the students are expected to be comprehensive senior engineers with dominant speciality characteristic, and be capable for all information, automation, interactive areas of electrical engineering such as engineering design, electrical production, technology development, education and research etc.三、专业培养基本要求Skills Profile本专业学生主要学习电工电子技术, 计算机应用技术, 信息通信技术，同时掌握电气工程的基本知识和理论，使学生具备跨学科的综合素质和专业技能。
潮流能发电技术（一）Tidal Current Power Generation(1) 在浅海、海峡、海湾或河口一带，涨潮与退潮会引起较强的潮流，水流速度较高，可直接利用潮流前进的动能来推动水轮机发电，方法类似于风力发电机，称之为潮流发电机。
风力机可用的形式基本上都可用于潮汐动能发电，比较典型的的是采用与水平轴风力机相似的结构与原理，下面举两个例子介绍一下。
敞水型水平轴桨叶潮流发电装置图1是一种潮汐动能发电机示意图，与顺风式风力机工作原理相同，也叫桨叶式潮流发电机。
转轮的叶片与风力机类似，由于水的密度约为空气的800多倍，相同功率下的叶片面积与长度可大大缩小，不过潮汐水流速一般没有风速快。
叶轮通过增速齿轮箱与发电机连接，一同安装在机舱内，机舱安装在支柱上，可绕支柱旋转，支柱固定在海底,属坐底式安装。
水流推动叶轮旋转，并保持叶轮面与水流方向垂直。
由于叶轮在敞开的水中旋转，属敞水式结构。
图1 桨叶式潮流发电机图2 水下的桨叶式潮流发电机桨叶式潮流发电机可安装多个组成潮流发电站，图2是六个桨叶式潮流发电机的水下模拟图。
桨叶式潮流发电机在国外已有成功的产品，比较典型的是英国研制的一种桨叶式潮流发电机，名为“SeaGen”，形似倒置的风车，由一对涡轮组成，每个涡轮直径20米，容量1200千瓦。
涡轮安装在柱桩上，柱桩固定在海底。
图3是该潮流发电机的图片。
图3 SeaGen桨叶式潮流发电机（图片来自网络）扩张型涡轮潮流发电装置还有一种是把叶轮安装在水平扩张导流管内，扩张导流管在风力机中称为扩散器，利用导流管对水流进行加速，可明显提高叶轮的输出功率。
叶轮由多个叶片组成，这是一种导流管增速涡轮机结构，或称为导流管增速水轮机，见图4。
图4 导管涡轮式潮流发电机图中浅蓝色箭头线表示水流走向，水流推动发电机发电，涡轮机采用双向对称结构，反方向水流同样推动发电机发电，适用于潮流发电。
导流管增速涡轮发电机可安装多个组成潮流发电站，图5是五个导管涡轮式潮流发电机水下模拟图。
大学生的个人SWOT分析。
第一篇：大学生的个人SWOT分析。
大学生的个人SWOT分析SWOT分析方法是一种企业内部分析方法，即根据企业自身的既定内在条件进行分析，找出企业的优势，劣势及核心竞争力之所在。
其中，S代表strength(优势)，W代表weakness(弱势)，O代表opportunity(机会)，T代表threat(威胁)，其中，S、W是内部因素，O、T是外部因素。
近来，SWOT分析已广被应用在许多领域上，如学校的自我分析，个人的能力自我分析等方面，本文将对自身进行SWOT分析，从而正确认识自己，了解自己，完善自己。
一、个人优劣势分析1、个人优点(Strengths)：性格上，做事认真，待人真诚，乐于助人，有强烈的责任感，言行举止得体，做人安分守己；能力上，有一定的分析能力，并有寻根究底的兴趣，一定要将事情想清楚；生活上，勤俭节约，不与别人攀比，生活态度比较乐观积极向上，自控能力强，意志坚定；学习上，努力勤奋，不易受外界环境的影响，成绩良好，上课专心致志，已经通过计算机二级。
2、个人缺点(Weakness)：性格大体倾向于内向，做事不果断，经常犹豫不决；在大众面前显得不够热情大方，说话不稳重，也没幽默感；有时不够积极参加学校和学院的各项活动。
在大学习没有结合长远目标，创新能力有待提高；写作能力差，语言表达不够流利和简洁；竞争意识不强，不够主动地利用校园资源，比如学校图书馆和电脑室等，很少借阅书和到电脑室利用英语资源。
二、环境外界分析1、环境优势(Opportunities)：现在我处于大三的第二学期期末，还有一年的时间来改进自身的缺点和锻炼各方面的能力。
正所谓二十一世纪是效率至上的社会，这说明作为管理学院的一员，我在二十一世纪将有一个大好的发展前途。
另一方面，国家振兴物流产业政策的实施，这又为我的就业提供了好机会。
国家出台了相关的政策来解决大学生就业的问题，我同样也受益在其中。
我住在一个和谐的宿舍里，舍友们个个勤奋好学，有他们的帮助与支持。
大学本科各大专业及学习课程详细介绍大全1.【专业名称】中药资源与开发本专业主要学习中药的基本理论和基本知识，学习获得调查分析中药资源和中药材生产加工、中药资源的综合开发和利用以及保护更新等方面的知识和技能。
学生将受到系统的中药资源与开发的专业基本训练，具有中药资源调查，中药原料的生产、加工，中药资源的综合开发和中药新药开发的基本能力。
毕业后能从事中药资源调查、中药材培养、中药材鉴定、中药原料采购、中药资源的综合开发和合理利用、中药新药研究开发等方面工作。
主要课程：计算机应用基础、程序设计语言、大学英语,高等数学(ⅠⅡ)、线性代数、概率统计、物理学、无机化学、分析化学、有机化学、物理化学、法律基础、中医药基础理论、基础生物学、药用植物学、药用动物学、植物生理和生态学、生药学（中药鉴定学）、药用植物栽培学、中药资源学、天然药物资源化学、植物化学分类学、中药分析化学、中药生物技术、中药药理学、中药材加工和炮制学、中药药剂学、中药新药和保健食品开发概论、药事法规。
2.【专业名称】包装工程业务培养目标：本专业培养具备包装系统设计与管理等方面的能力，能在商品生产与流通部门、包装企业、科研机构、外贸、商检等部门从事包装系统设计、质量检测、技术管理和科学研究的高级工程技术人才。
业务培养要求：本专业学生要学习保护产品、方便流通、促进销售的包装基础理论．包装设计原理和方法、包装材料、包装印刷、包装测试、包装艺术设计、包装设计等基本知识。
毕业生应获得以下几方面的知识和能力：1．掌握工程力学、材料学、生物学、设计美学等包装工程的基础理论；2．掌握包装工艺、包装结构设计方法和包装测试、包装管理技术；3．具有制定包装工艺、合理选择包装材料和包装设备的初步能力；4．熟悉国家有关包装的方针、政策和法规；5．了解包装学前沿和发展趋势；6．掌握文献检索、资料查询的基本方法，具有分析解决包装工程技术问题，研究、开发包装新材料、新结构、新工艺、新设备和技术管理的初步能力。
课程设计报告( 2011-- 2012年度第1学期)名称：电力系统潮流上机院系：电气与电子工程学院班级：电网1101 学号：学生姓名：指导教师：王莉丽设计周数：两周成绩：日期：2014年1月《电力系统综合仿真》课程设计一、 目的与要求培养学生编程实现电力系统仿真计算的能力，掌握潮流计算、短路计算、暂态稳定计算的相关知识二、 主要内容1． 编写潮流计算程序，要求如下：2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器模型，电压等级的归算，标幺值的计算；2.2基础数据的计算机存储：节点数据，支路数据（包括变压器）； 2.3用牛顿-拉夫逊法计算；2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV 节点无功功率和平衡节点功率，数据输出模块； 2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图2． 利用潮流计算结果编写故障计算程序，要求如下：2.1 发电机参数一律为P N =200MW ，cos φ=0.85, 19.0=''dx ,21x x = 2.2 变压器一律为Yd11接线，其中三角形侧接发电机，星形侧都直接接地 2.3 线路零序参数一律为103x x =2.4 负荷都采用Y/Y 型变压器接入，中性点与系统相连侧部接地。
负荷按照对称负荷考虑，纯阻抗形式。
2.5 计算5节点发生三相对称短路、单相接地短路、两相短路和两相接地短路时的故障电流，以及各支路电流和节点电压。
3．思考题3.1潮流计算和短路电流计算的元件模型有哪些不同点？节点导纳矩阵有哪些不同？计算过程有哪些不同点？3.2 短路电流计算中进行了哪些假设？3.3如果交给你一个任务，请你用已有的潮流计算软件计算北京城市电网的潮流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析）3.4比较过去的计算机语言上机学习和本课程上机有哪些不同？设计中遇到的问题和解决的办法。
三、进度计划四、设计成果要求1．提交任务书2．计算机计算潮流程序流程图；3．计算机计算短路电流程序流程图；4．完整的仿真程序(包括变量的定义，子程序的调用等，计算结果)，并有程序注释。
5．回答思考题五、考核方式总成绩=考勤+设计报告+面试成绩学生姓名：指导教师：2014年1月附录：1.手算：2.计算机计算潮流程序流程图3.程序#include"stdafx.h"#include"fstream"#include"iostream"#include"DataForm.h"#include"NEquation.h"#include"math.h"using namespace System;using namespace std;#define PI 3.1415926void DataRead();void FormY();void V oltageInitial();void SetUnbalance();void FormJacobi();void GetRevised();void GetNewV oltageValue();void BusBalance();void BranchPQ();void BranceLoss();int main(array<System::String ^> ^args){DataRead();FormY();V oltageInitial();for(int i=0;i<50;i++){SetUnbalance();FormJacobi();GetRevised();GetNewV oltageValue();//break()//if(fabs(H[i][j])>=1e-5||(N[i][j])>=1e-5||(H[i][j])>=1e-5||(J[i][j])>=1e-5||(L[i][j])>=1e-5) }BusBalance();BranchPQ();BranceLoss();Console::WriteLine(L"Hello World");return 0;}void DataRead() //节点、支路数据读入子程序{ ifstream infile1;infile1.open("BranchData2.txt");for(int i=0;i<BranchNum;i++)infile1>>Branch[i].Num>>Branch[i].BusFirst>>Branch[i].BusEnd>>Branch[i].R>>Branch[i].X>>Branch[i].B>>Branch[i].K;infile1.close();for(int i=0;i<BranchNum;i++)cout<<Branch[i].Num<<" "<<Branch[i].BusFirst<<" "<<Branch[i].BusEnd<<" "<<Branch[i].R<<" "<<Branch[i].X<<" "<<Branch[i].B<<" "<<Branch[i].K<<endl;ifstream infile2;infile2.open("BusData2.txt");for(int i=0;i<BusNum;i++)infile2>>Bus[i].Num>>Bus[i].BusType>>Bus[i].V oltage>>Bus[i].Phase>>Bus[i].Pg>>Bus[i].Qg>>Bus[i].Pl>>Bus[i].Ql;infile2.close();for(int i=0;i<BusNum;i++)cout<<Bus[i].Num<<Bus[i].BusType<<Bus[i].V oltage<<Bus[i].Phase<<Bus[i].Pg<<Bus[i].Qg<<Bus[i].Pl<<Bus[i].Ql<<endl;}void FormY() //形成节点导纳矩阵子程序{for(int i=0;i<BusNum;i++){ for(int j=0;j<BusNum;j++){ G[i][j]=0;B[i][j]=0;}}int m,n;//设置两个变量分别存储之路手摸节点的编号for(int i=0;i<BranchNum;i++)//以支路循环{float GG,BB,B0,G1,G2,B1,B2;m=Branch[i].BusFirst-1;//n=Branch[i].BusEnd-1;//以下为变压器变比不唯Branch[i].R=Branch[i].K*Branch[i].R;Branch[i].X=Branch[i].K*Branch[i].X;//将变压器变比引起的对地之路转换成电哪里？G1=(1-Branch[i].K)/(Branch[i].K*Branch[i].K)/(Branch[i].R*Branch[i].R+Branch[i].X*Branch[i].X)*Branch[i].R;B1=-(1-Branch[i].K)/(Branch[i].K*Branch[i].K)/(Branch[i].R*Branch[i].R+Branch[i].X*Branch[i].X)*Branch[i].X;//将变压器变比引起的对地之路换成导纳G2=(Branch[i].K-1)/Branch[i].K/(Branch[i].R*Branch[i].R+Branch[i].X*Branch[i].X)*Branch[i].R;B2=-(Branch[i].K-1)/Branch[i].K/(Branch[i].R*Branch[i].R+Branch[i].X*Branch[i].X)*Branch[i].X;//变压器非标准变比处理完毕GG=Branch[i].R/(Branch[i].R*Branch[i].R+Branch[i].X*Branch[i].X);BB=-Branch[i].X/(Branch[i].R*Branch[i].R+Branch[i].X*Branch[i].X);B0=Branch[i].B;//接地数据G[m][m]+=GG+G1;G[m][n]-=GG;G[n][m]-=GG;G[n][n]+=GG+G2;B[m][m]+=BB+B0+B1;B[m][n]-=BB;B[n][m]-=BB;B[n][n]+=BB+B0+B2;//B[][]+=}ofstream outfile1("result.txt");for(int i=0;i<BusNum;i++){ for(int j=0;j<BusNum;j++){outfile1<<G[i][j]<<"\t"<<B[i][j]<<"\t";}outfile1<<endl;}outfile1.close();}void V oltageInitial() //设置节点电压子程序{for(int i=0;i<BusNum;i++){if(Bus[i].BusType==3){Ue[i]=1.04;Uf[i]=0.0;}else{Ue[i]=1.0;Uf[i]=0.0;}cout<<Ue[i]<<" "<<Uf[i]<<endl;}}void SetUnbalance() //计算功率不平衡量子程序{for(int i=1;i<BusNum;i++){ P[i]=0;Q[i]=0;for(int j=0;j<BusNum;j++){ P[i]=P[i]+Ue[i]*(G[i][j]*Ue[j]-B[i][j]*Uf[j])+Uf[i]*(G[i][j]*Uf[j]+B[i][j]*Ue[j]);Q[i]=Q[i]+Uf[i]*(G[i][j]*Ue[j]-B[i][j]*Uf[j])-Ue[i]*(G[i][j]*Uf[j]+B[i][j]*Ue[j]);}//cout<<P[i]<<endl;//cout<<Q[i]<<endl;if(Bus[i].BusType==1) //PQ{ DeltaP[i]=Bus[i].Pg-Bus[i].Pl-P[i];DeltaQU2[i]=Bus[i].Qg-Bus[i].Ql-Q[i];}else if(Bus[i].BusType==2) //PV{ DeltaP[i]=Bus[i].Pg-Bus[i].Pl-P[i];DeltaQU2[i]=Bus[i].V oltage*Bus[i].V oltage-(Ue[i]*Ue[i]+Uf[i]*Uf[i]);}//cout<<DeltaP[i]<<endl;//cout<<DeltaQU2[i]<<endl;}for(int i=0;i<BusNum-1;i++){ Delta[2*i]=DeltaP[i+1];Delta[2*i+1]=DeltaQU2[i+1];} //形成delta值ofstream outfile2("Delta2.txt");for(int j=0;j<2*(BusNum-1);j++){ outfile2<<Delta[j]<<endl;}outfile2.close();}void FormJacobi() //形成雅各比矩阵子程序{float H[BusNum-1][BusNum-1],N[BusNum-1][BusNum-1],J[BusNum-1][BusNum-1],L[BusNum-1][BusNum-1];float a[BusNum],b[BusNum];for(int i=1;i<BusNum;i++){a[i]=G[i][i]*Ue[i]-B[i][i]*Uf[i];b[i]=G[i][i]*Uf[i]+B[i][i]*Ue[i];for(int j=0;j<BusNum;j++){ if(j!=i){a[i]=a[i]+(G[i][j]*Ue[j]-B[i][j]*Uf[j]);b[i]=b[i]+(G[i][j]*Uf[j]+B[i][j]*Ue[j]);}}}for(int i=1;i<BusNum;i++){for(int j=1;j<BusNum;j++){ if(i==j){ H[i][i]=-B[i][i]*Ue[i]+G[i][i]*Uf[i]+b[i];N[i][i]=G[i][i]*Ue[i]+B[i][i]*Uf[i]+a[i];if(Bus[i].BusType==1){ J[i][i]=-G[i][i]*Ue[i]-B[i][i]*Uf[i]+a[i];L[i][i]=-B[i][i]*Ue[i]+G[i][i]*Uf[i]-b[i];}else if(Bus[i].BusType==2){ J[i][i]=2*Uf[i];}}else if(i!=j){ H[i][j]=-B[i][j]*Ue[i]+G[i][j]*Uf[i];N[i][j]=G[i][j]*Ue[i]+B[i][j]*Uf[i];if(Bus[i].BusType==1){ J[i][j]=-G[i][j]*Ue[i]-B[i][j]*Uf[i];L[i][j]=-B[i][j]*Ue[i]+G[i][j]*Uf[i];}else if(Bus[i].BusType==2){ J[i][j]=0;L[i][j]=0;}}}}for(int i=0;i<BusNum-1;i++){ for(int j=0;j<BusNum-1;j++){ Jacobi[2*i][2*j]=H[i+1][j+1];Jacobi[2*i][2*j+1]=N[i+1][j+1];Jacobi[2*i+1][2*j]=J[i+1][j+1];Jacobi[2*i+1][2*j+1]=L[i+1][j+1];}}ofstream outfile4;outfile4.open("Jacobi2.txt");for(int i=0;i<2*(BusNum-1);i++){for(int j=0;j<2*(BusNum-1);j++){outfile4<<Jacobi[i][j]<<"\t ";}outfile4.close();}void GetRevised() //求解修正方程得到电压增量子程序{NEquation Object;//建立一个对象Object.SetSize(2*(BusNum-1));//设置矩阵的维数ifstream infile1;infile1.open("Jacobi2.txt");for(int i=0;i<2*(BusNum-1);i++){ for(int j=0;j<2*(BusNum-1);j++){ Object.Data(i,j)=Jacobi[i][j];infile1>>Object.Data(i,j);}}infile1.close();ifstream infile2;infile2.open("Delta2.txt");for(int i=0;i<2*(BusNum-1);i++)infile2>>Object.Value(i);infile2.close();Object.Run();for(int i=0;i<2*(BusNum-1);i++){ DeltaU[i]=Object.Value(i);}//结果输出ofstream outfile5("revised equation2.txt");for(int j=0;j<2*(BusNum-1);j++){ outfile5<<DeltaU[j]<<endl;}outfile5.close();}for(int i=0;i<BusNum;i++)if(Bus[i].BusType!=3){Ue[i]+=DeltaU[2*i-1];Uf[i]+=DeltaU[2*i-2];}ofstream outfile6("new voltage value2.txt");for(int i=0;i<BusNum;i++){ //if(Bus[i].BusType==1||Bus[i].BusType==2)outfile6<<Ue[i]<<" "<<Uf[i]<<endl;}outfile6.close();}void BusBalance() //计算平衡节点功率子程序{for(int i=0;i<BusNum;i++){if(Bus[i].BusType==3){P[i]=0.0;Q[i]=0.0;for(int j=0;j<BusNum;j++){P[i]+=Ue[i]*(G[i][j]*Ue[j]-B[i][j]*Uf[j])+Uf[i]*(G[i][j]*Uf[j]+B[i][j]*Ue[j]);//有功Q[i]+=Uf[i]*(G[i][j]*Ue[j]-B[i][j]*Uf[j])-Ue[i]*(G[i][j]*Uf[j]+B[i][j]*Ue[j]);//无功}}}ofstream outfile7("BusBalance2.txt");for(int i=0;i<BusNum;i++)if(Bus[i].BusType==3){ outfile7<<P[i]<<" \t"<<Q[i]<<endl;}outfile7.close();//for(int i=0;i<BusNum;i++)//if(Bus[i].BusType==3)//{ cout<<P[i]<<" "<<Q[i]<<endl;}}void BranchPQ() //计算支路功率子程序{ float B0;for(int i=0;i<BusNum;i++){for(int j=0;j<BusNum;j++){if(i!=j&(G[i][j])!=0&(B[i][j])!=0){B0=Branch[i].B;BranchP[i][j]=Ue[i]*(-Uf[i]*B0+Ue[i]*0+(Ue[i]-Ue[j])*G[i][j]-(Uf[i]-Uf[j])*B[i][j])+Uf[i]*(Ue[i]*B0+Uf[i]*0+(Ue[i]-Ue [j])*B[i][j]+(Uf[i]-Uf[j])*G[i][j]);BranchQ[i][j]=Uf[i]*(-Uf[i]*B0+Ue[i]*0+(Ue[i]-Ue[j])*G[i][j]-(Uf[i]-Uf[j])*B[i][j])-Ue[i]*(Ue[i]*B0+Uf[i]*0+(Ue[i]-Ue[j])*B [i][j]+(Uf[i]-Uf[j])*G[i][j]);BranchP[j][i]=Ue[j]*(-Uf[j]*B0+Ue[j]*0+(Ue[j]-Ue[i])*G[j][i]-(Uf[j]-Uf[i])*B[j][i])+Uf[j]*(Ue[j]*B0+Uf[j]*0+(Ue[j]-Ue[i])*B [j][i]+(Uf[j]-Uf[i])*G[j][i]);BranchQ[j][i]=Uf[j]*(-Uf[j]*B0+Ue[j]*0+(Ue[j]-Ue[i])*G[j][i]-(Uf[j]-Uf[i])*B[j][i])-Ue[j]*(Ue[j]*B0+Uf[j]*0+(Ue[j]-Ue[i])*B [j][i]+(Uf[j]-Uf[i])*G[j][i]);}}}ofstream outfile8("BranchPQ2.txt");for(int i=0;i<BusNum;i++){ for(int j=0;j<BusNum;j++)//if(i!=j&(G[i][j])!=0&(B[i][j])!=0){ outfile8<<BranchP[i][j]<<" \t"<<BranchQ[i][j]<<" \t";}outfile8<<endl;}outfile8.close();}void BranceLoss() //计算线路上损耗子程序{ float LossP[BusNum][BusNum],LossQ[BusNum][BusNum],LossPP,LossQQ;LossPP=0;LossQQ=0;for(int i=0;i<BusNum;i++){for(int j=0;j<BusNum;j++){if(i!=j&(G[i][j])!=0&(B[i][j])!=0){LossP[i][j]=BranchP[i][j]+BranchP[j][i];LossQ[i][j]=BranchQ[i][j]+BranchQ[j][i];//cout<<LossP[i]<<" "<<LossQ[i]<<endl;}}}for(int i=0;i<BusNum;i++)for(int j=0;j<BusNum;j++){if(i!=j&(G[i][j])!=0&(B[i][j])!=0){LossPP+=LossP[i][j];LossQQ+=LossQ[i][j];}}//cout<<0.5*LossPP<<" "<<0.5*LossQQ;ofstream outfile9("BranchLoss2.txt");//for(int i=0;i<BusNum;i++){ outfile9<<0.5*LossPP<<"\t "<<0.5*LossQQ<<endl;} outfile9.close();4.思考题1潮流计算的方法有哪些？各有何特点？潮流计算目前主要的计算方法有：高斯-赛德尔迭代法、牛顿拉夫逊法和P-Q分解法。
