潮流程序报告 班级：硕721 姓名：胡与非 学号：3107161033 一、概述 .doc

潮流程序报告 班级：硕721 姓名：胡与非 学号：3107161033 一、概述 本程序是PQ分解法中的BX法，本程序中主要用到的三个主要算法概述如下： 1．节点优化编号 本程序使用Tinney-2编号方法，也称为最小度算法，或半动态节点优化编号方法。该方法使用最为广泛，且简单有效，可以大大减少矩阵因子分解过程中产生的注入元的数目。 2．稀疏技术 稀疏技术包括稀疏矩阵技术和稀疏矢量技术。本程序使用了稀疏矩阵技术，利用网络矩阵的稀疏结构，减少了存储量和计算量，提高了求解网络方程的计算速度。 3．潮流计算的快速分解法（FDLF） 快速分解法不是求解高维数的修正方程，而是分别交替求解两个低维数且系数矩阵是常数矩阵的修正方程，因此计算速度快，收敛性好，应用广泛。 二、输入数据格式说明 第一行：支路数 第二行：节点数 第三行：接地支路数 第四行：发电机数 第五行：负荷数 第六行：平衡节点号 第七行：精度 以下是支路数据： 1、支路一端节点号 支路另一端节点号 支路电阻 支路电抗 接地电容(变压器变比) （节点号为负表明变压器支路,负号的节点侧为变压器的非标准变比侧,此时节点电容换为变压器变比；） 2、接地支路号 电阻 电容 以下是发电机数据：（最后一行是平衡机，其余发电机看成PV节点） 节点号 有功功率 无功功率 电压 以下是负荷数据： 节点号 有功功率 无功功率 电压 三、数据结构介绍 //----------------线路数据结构体--------------// struct Line { int I ;//线路一端所连的节点号，节点号绝对值ij,i按绝对值从小到大排列; //节点号为负表示为变压器支路，负号的节点侧为变压器的非标准变比侧;此时节点电容换为变压器变比 int j; //线路另一端所连的节点号 double r;//线路电阻 double x;//线路电抗 double b;//接地电容，若为变压器节点电容换为变压器变比//这里的b为线路导纳的一半b/2 int hh; //标记支路号 用于还原输出 }*line; //----------------接地线路数据结构体--------------// struct Grand_line { int i; //接地支路号 double r;//接地支路电阻 double b;//接地支路电容 int hh; //标记支路号，用于还原输出 }*grand; //----------------发电机数据结构体--------------// struct Generator { int i; //发电机节点号,最后一个发电机为平衡节点 double P;//发电机有功功率 double Q;//发电机无功功率 double V;//发电机机端电压,电压为负表示PV节点 }*generator; struct Pinghengji { int i;//平衡机节点号,最后一个发电机为平衡节点 double P;//平衡机有功功率 double Q;//平衡机无功功率 double V;//平衡机机端电压 }*pinghengji; //----------------负荷数据结构体--------------// struct Load { int i;//负荷节点号 double P;//负荷的有功功率 double Q;//负荷的无功功率 double V;//负荷处的电压(初始为0) }*load; //----------------优化结构体--------------// struct Youhua {int nn;//记录节点的度 int nn1; int *jd;//记录节点号 int *w; }*youhua; //-------------自导纳结构体--------------// struct Yii {double g,b; }*yii; struct Yii1 {double g,b; }*yii1; //------------互导纳结构体----------------// struct Yij {double g,b; int j; }*yij; struct Yij1 {double g,b; int j; }*yij1; //----------------U--------------// struct U_Type {double value; int j; }*U1,*U2,*U; //-------------------节点电压数组-----------------// s