本科课设——电力系统潮流计算

主要内容来自本人大三期间的课程设计，包括最基本的潮流计算，本文内容为五节点的网络结构

// chaoliu.cpp : 定义控制台应用程序的入口点。
//

#include “stdafx.h“
#include “stdio.h“
#include “NEquation.h“
#include “math.h“

#define NODE_TOTAL_NUM 9
#define LINE_TOTAL_NUM 9
#define NODE_DATA_FILENAME	 “D:\\desktop\\chaoliu\\chaoliu\\Node.txt“
#define LINE_DATA_FILENAME	 “D:\\desktop\\chaoliu\\chaoliu\\Line.txt“
#define Y_DATA_FILENAME		 “D:\\desktop\\chaoliu\\chaoliu\\Y.txt“
#define U_DATA_FILENAME		 “D:\\desktop\\chaoliu\\chaoliu\\U.txt“
#define UBL_DATA_FILENAME	 “D:\\desktop\\chaoliu\\chaoliu\\UBL.txt“
#define JCB_DATA_FILENAME		“D:\\desktop\\chaoliu\\chaoliu\\JCB.txt“
#define DELEF_DATA_FILENAME		“D:\\desktop\\chaoliu\\chaoliu\\DELEF.txt“
#define NEWVALUE_DATA_FILENAME “D:\\desktop\\chaoliu\\chaoliu\\NEWVALUE.txt“
#define NETCOM_DATA_FINENAME	“D:\\desktop\\chaoliu\\chaoliu\\NETCOM.txt“
#define eps 0.001
#define pi 3.1415926535

struct  NodeData
{
	unsigned int Index; //node index
	unsigned int Type;  //node type :PQ:1PV:2balance point:0
	float FirstInput;   //PQ or PV:actice powerbalance point:V 
	float SecondInput;  //PQ:reactive powerPV:Vbalance point:angle
};

struct  LineData
{
	unsigned int Index;  //line index
	unsigned int Iftran; //line type: not transformer: 1 transformer: 2
	unsigned int start;
	unsigned int end;
	float R;	 //PQ or PV:actice powerbalance point:V 
	float X;	 //PQ:reactive powerPV:Vbalance point:angle
	float B_2;	 //transformer: k not transformer:B_2
};

struct  NodeData gNodeData[NODE_TOTAL_NUM];
struct  LineData gLineData[LINE_TOTAL_NUM];
double YG[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 } YB[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };   //导纳矩阵
double Gl[LINE_TOTAL_NUM] = { 0 } Bl[LINE_TOTAL_NUM] = { 0 };                                   //支路导纳
double e[NODE_TOTAL_NUM] = { 0 } f[NODE_TOTAL_NUM] = { 0 };						             //电压
double P[NODE_TOTAL_NUM] = { 0 } Q[NODE_TOTAL_NUM] = { 0 };                                     //功率
double DeltaP[NODE_TOTAL_NUM-1] = { 0 } DeltaQU2[NODE_TOTAL_NUM-1] = { 0 };                     //两个小的不平衡列向量
double DeltaPQ[2 * （NODE_TOTAL_NUM - 1）];                                                        //不平衡总列向量
double Jacobi[2 * NODE_TOTAL_NUM][2 * NODE_TOTAL_NUM] = { 0 };                                   //雅可比矩阵
double H[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };                             
double N[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };                             
double J[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };                             
double L[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };                             
double R[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };                             
double S[NODE_TOTAL_NUM][NODE_TOTAL_NUM] = { 0 };                             
double Deltaef[2 * （NODE_TOTAL_NUM -1）] = { 0 };                                                  //总电压偏差量
double Deltaf[NODE_TOTAL_NUM - 1] = { 0 } Deltae[NODE_TOTAL_NUM - 1] = { 0 };                    //横纵电压偏差量
int Count_num=0;                       

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    ..A..H.     35328  2017-06-30 21:14  chaoliu\chaoliu\.vs\chaoliu\v14\.suo

     文件      21015  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\chaoliu.cpp

     文件       8101  2017-06-27 22:24  chaoliu\chaoliu\chaoliu\chaoliu.vcxproj

     文件       1405  2017-06-27 22:24  chaoliu\chaoliu\chaoliu\chaoliu.vcxproj.filters

     文件        837  2017-06-30 21:13  chaoliu\chaoliu\chaoliu\Debug\chaoliu.Build.CppClean.log

     文件       3461  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.log

     文件      68471  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.obj

     文件    3407872  2017-06-30 21:13  chaoliu\chaoliu\chaoliu\Debug\chaoliu.pch

     文件        193  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\chaoliu.lastbuildstate

     文件       1510  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\CL.command.1.tlog

     文件       4714  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\CL.read.1.tlog

     文件        934  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\CL.write.1.tlog

     文件       1222  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\link.command.1.tlog

     文件       2648  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\link.read.1.tlog

     文件        492  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\chaoliu.tlog\link.write.1.tlog

     文件      11740  2017-06-30 21:13  chaoliu\chaoliu\chaoliu\Debug\stdafx.obj

     文件     125952  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\vc140.idb

     文件     151552  2017-06-30 21:14  chaoliu\chaoliu\chaoliu\Debug\vc140.pdb

     文件       3413  2017-06-28 21:26  chaoliu\chaoliu\chaoliu\NEquation.h

     文件       1510  2017-06-26 12:28  chaoliu\chaoliu\chaoliu\ReadMe.txt

     文件        211  2017-06-26 12:28  chaoliu\chaoliu\chaoliu\stdafx.cpp

     文件        234  2017-06-26 12:28  chaoliu\chaoliu\chaoliu\stdafx.h

     文件        240  2017-06-26 12:28  chaoliu\chaoliu\chaoliu\targetver.h

     文件       1303  2017-06-26 12:28  chaoliu\chaoliu\chaoliu.sln

     文件   12189696  2017-06-30 21:14  chaoliu\chaoliu\chaoliu.VC.db

     文件      67072  2017-06-30 21:14  chaoliu\chaoliu\Debug\chaoliu.exe

     文件     374840  2017-06-30 21:14  chaoliu\chaoliu\Debug\chaoliu.ilk

     文件     806912  2017-06-30 21:14  chaoliu\chaoliu\Debug\chaoliu.pdb

     文件       1516  2017-06-30 21:14  chaoliu\chaoliu\DELEF.txt

     文件    3407872  2017-06-28 20:18  chaoliu\chaoliu\ipch\CHAOLIU-bb30bd72\CHAOLIU-bc96da71.ipch

............此处省略22个文件信息