谁知道一个简单电力系统的牛顿拉夫逊法的分析!要一个简单的,有具体过程和编程!
kuaidi.ping-jia.net 作者：佚名 更新日期：2024-06-15
考虑负荷静态特性时,牛顿法潮流计算的每次迭代过程中,都需要根据当前电压值对负荷的有功和无功值进行修正
1.3MATLAB概述
目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单。
MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。
MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。
另外，MATLAB提供了一种特殊的工具：工具箱（TOOLBOXES）.这些工具箱主要包括：信号处理（SIGNAL PROCESSING）、控制系统（CONTROL SYSTEMS）、神经网络（NEURAL NETWORKS）、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。
MATLAB设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。
原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。
2．1 电力系统的基本概念
2．1．1电力系统
（1）电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。
（2）电力网：电力系统中除发电机和用电设备外的部分。
（3）动力系统：电力系统和“动力部分”的总和。
“动力部分”：包括火力发电厂的锅炉、汽轮机、热力网和用电设备，水力发电厂的水库和水轮机，核电厂的反应堆等。
2．1．2电力系统的负荷和负荷曲线
（1）电力系统的负荷：系统中千万个用电设备消费功率的总和，包括异步电动机、同步电动机、电热炉、整流设备、照明设备等若干类。
（2）电力系统的供电负荷：综合用电负荷加上电力网中损耗的功率。
（3）电力系统的发电负荷：供电负荷加上发电厂本身的消耗功率。
（4）各用电设备的有功功率和无功功率随受电电压和系统频率的变化而变化，其变化规律不尽相同，综合用电负荷随电压和频率的变化规律是各用电负荷变化规律的合成。
（5）负荷曲线：某一时间段内负荷随时间而变化的规律。
（6）按负荷种类可分有功功率负荷和无功功率负荷；按时间长短可分为日负荷和年负荷曲线；按计量地点可分为个别用户、电力线路、变电所、发电厂以至整个系统的负荷曲线。将上述三种分类相结合，就确定了某一种特定的负荷曲线。不同行业的有功功率日负荷曲线差别很大。负荷曲线对电力系统的运行又很重要的意义，它是安排日发电计划，确定各发电厂发电任务以及确定系统运行方式等的重要依据。
2．2 电力系统的基本元件
2．2．1 发电机
现代电力工业中，无论是火力发电、水力发电或核能发电，几乎全部采用同步交流发电机。电机的电枢布置在定子上，励磁绕组布置在转子上，作为旋转式磁极。同步发电机的转速（转/MIN）和系统频率f（HZ）之间有着严格的关系，即n=60f/p式中p为电机的极对数。
根据转子结构型式的不同，分为隐极式和凸极式发电机，前者转子没有显露出来的磁极，后者则有。
转子的励磁型式有直流励磁系统和可控硅励磁系统，后者利用同轴交流励磁机或由同步发电机本身发出的交流电，经整流后供给转子。直流励磁机有换向问题，故其制造容量受到限制，所以，在大容量发电机中均可采用可控硅励磁系统。
2．2．2 电力变压器
电力变压器是电力系统中广泛使用的升压和降压设备。据统计，电力系统中变压器的安装总容量约为发电机安装容量的6－8倍。按用途，电力变压器可分为升压变压器、降亚变压器、配电变压器和联络变压器。按相数分，变压器可分为单相式和三相式。按每相线圈分，又有双绕组和三绕组之分。按线圈耦合的方式，可分为普通变压器和自耦变压器。
2．2．3 电力线路
（1）架空线路：由导线、避雷针、杆塔、绝缘子和金具等构成。
（2）电缆线路：由导线、绝缘层、包护层等构成。
2．2．4 无功功率补偿设备
主要的无功功率补偿设备有同步调相机、电力电容器和静止补偿器。
2．3 电力系统元件的数学模型
2．3．1 电力线路的等值电路
在电力系统分析中，一般只考虑电力线路两侧端口的电压和电流，把电力线路作为无源双口网络处理。
线路的双口网络方程：
Z=B=*L*
2．3．2 变压器的等值电路
（1）双绕组变压器等值电路
（2）三绕组变压器等值电路
2．3．3 同步发电机的数学模型
2．3．4 电力系统负荷
2．3．5 多级电压电力系统的等值电路
2．4 电力系统稳态运行分析
2．4．1 电力线路的电压损耗与功率损耗
2．4．2 变压器中的功率损耗与电压损耗
2．4．3 辐射形网络的分析计算
辐射形电力网的特点是各条线路有明确的始端与末端。辐射形电力网的分析计算就是利用已知的负荷、节点电压来求取未知的节点电压、线路功率分布、功率损耗及始端输出功率。
辐射形电力网的分析计算，根据已知条件的不同分两种
1 已知末端功率与电压：即 从末端逐级往上推算，直至求得各要求的量。
2 已知末端功率、始端电压：末端可理解成一负荷点，始端为电源点或电压中枢点。采用迭代法。
（1）假设末端电压为线路额定电压，利用第一种方法求得始端功率及全网功率分布。
（2）用求得的线路始端功率和已知的线路始端电压，计算线路末端电压和全网功率分布。
（3）用第（2）步求得的线路末端电压计算线路始端功率和全网功率分布，如求得的各线路功率与前一次相同计算的结果相差小于允许值，就可以认为本步求得的线路电压和全网功率分布为最终计算结果。否则，返回第二步重新进行计算。
2．4．4 复杂电力系统潮流计算
电力系统潮流计算始对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标始求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。因此，潮流计算是电力系统计算分析中的一种最基本的计算。
潮流计算结果的用途，例如用于电力系统稳定研究、安全估计或最优潮流等也对潮流计算的模型和方法有直接影响。
2．5 电力系统潮流计算机算法
2．5．1电力系统潮流计算机算法概述
2．5．1．1 导纳矩阵的形成
2．5．1．2 节点类型
（1）PV节点：柱入有功功率P为给定值，电压也保持在给定数值。
（2）PQ节点：诸如有功功率和无功功率是给定的。
（3）平衡节点：用来平衡全电网的功率。选一容量足够大的发电机担任平衡全电网功率的职责。
平衡节点的电压大小与相位是给定的，通常以它的相角为参考量，即取其电压相角为0。一个独立的电力网中只设一个平衡点。
2．5．1．3 高斯迭代法
2．5．2 牛顿－拉夫逊法
2．5．2．1 原理
2．5．2．2 基本步骤
基本步骤：
（1）形成节点导纳矩阵
（2）将各节点电压设初值U，
（3）将节点初值代入相关求式，求出修正方程式的常数项向量
（4）将节点电压初值代入求式，求出雅可比矩阵元素
（5）求解修正方程，求修正向量
（6）求取节点电压的新值
（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步
（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。
2．5．2．3 注意事项
2．5．2．4 程序流程框图
2．6 软件设计
2．6．1 方案选择及说明
2．6．2 方案求解
2．6．3 MATLAB编程说明及元件描述
2．6．4 程序
#include<stdio.h>
struct powernode
{
float pi;
float qi;
int i;
float vi;
};
struct powernode wg[20];
struct powernode wl[20];
struct linedata
{
int i;
int j;
float r;
float x;
float y; /*包括变压器变比*/
float k; /*只用作标析变压器，变压器变比仍在y中*/
};
struct linedata zl[20];
struct linedata t3; /*临时数组*/
static double y[][3]; /*在matrixform中应用*/
int t=0;
int t2,ti,tj; /*临时记数单元*/
float temp;
float tx,tr,YK; /*中间工作单元(在matrixform中应用)*/
double GIJ,BIJ; /*中间工作单元(在matrixform中应用)*/
int N; /*总节点数*/
int zls;
int Q,V,PVS,PVD;
int GS;
int LS;
float vo;
float Eps;
static double GII[]={0},BII[]={0},YDS[]={0},YDZ[]={0},B[]={0};/*添加数组*/
/*因子表形成时定义的数据*/
struct pvdata
{
float vis;
int i;
};
static struct pvdata pv[]={0};
datain()
{
clrscr();
printf("program runningn" );
printf("n");
printf("please input the aggregate to the system note");/*总节点数*/
scanf("%d",&N);
printf("n");
printf(" PQ note IN ALL?");/*总节点数*/
scanf("%d",&Q);
PVS=(N-Q)-1;
printf("n");
printf("them input the aggregate to the system power line");
scanf("%d",&zls);/*输电线路数和变压器的总数*/
printf("n");
printf("electromotor node in all :?");/*发电机节点总数*/
scanf("%d",&GS);
printf("n");
printf("load node in all : ?");/*负荷节点总数*/
scanf("%d",&LS);
printf("n");
printf("average electric voltage");/*平均电压*/
scanf("%f",vo);
printf("n");
printf("n");
printf("please input the date messagen");
printf("follow the format like it: i,j,r,x,y,kn");
do{
t++;
scanf("%d,%d,%f,%f,%f",&zl[t].i,&zl[t].j,&zl[t].r,&zl[t].x,&zl[t].y,&zl[t].k);
printf("processing....n");
if(zl[t].i>zl[t].j)
{
temp=zl[t].i;
zl[t].i=zl[t].j;
zl[t].j=temp;
/* if(zl[t].k!=1) */ /*要考虑归算问题不？？？？*/
}
printf("data you input is:n " );
printf("%d,%d,%f,%f,%f",zl[t].i,zl[t].j,zl[t].r,zl[t].x,zl[t].y,zl[t].k);
}while(zl[t].i!=0&&zl[t].j==0);
for(t2=t;t>0;t--) /*冒泡法排序*/
{
for(;t2>0;t2--)
{
if(zl[t2].i<zl[t2-1].i)
{
t3.i=zl[t2].i;t3.j=zl[t2].j;t3.r=zl[t2].r;t3.x=zl[t2].x;t3.y=zl[t2].y;t3.k=zl[t2].k;
zl[t2].i=zl[t2-1].i;zl[t2].j=zl[t2-1].j;zl[t2].r=zl[t2-1].r;zl[t2].x=zl[t2-1].x;zl[t2].y=zl[t2-1].y;zl[t2].k=zl[t2-1].k;
zl[t2-1].i=t3.i;zl[t2-1].j=t3.j;zl[t2-1].r=t3.r;zl[t2-1].x=t3.x;zl[t2-1].y=t3.y;zl[t2-1].k=t3.k;
}
else if(zl[t2].i==zl[t2-1].i)
{if(zl[t2].j<zl[t2-1].j)
{
t3.i=zl[t2].i;t3.j=zl[t2].j;t3.r=zl[t2].r;t3.x=zl[t2].x;t3.y=zl[t2].y;t3.k=zl[t2].k;
zl[t2].i=zl[t2-1].i;zl[t2].j=zl[t2-1].j;zl[t2].r=zl[t2-1].r;zl[t2].x=zl[t2-1].x;zl[t2].y=zl[t2-1].y;zl[t2].k=zl[t2-1].k;
zl[t2-1].i=t3.i;zl[t2-1].j=t3.j;zl[t2-1].r=t3.r;zl[t2-1].x=t3.x;zl[t2-1].y=t3.y;zl[t2-1].k=t3.k;
}
}
}
}
printf("n");
t=0;
printf("please input wg~!n");
do
{
scanf("%f,%f,%d,%f",&wg[t].pi,&wg[t].qi,&wg[t].i,&wg[t].vi);
t++;
}while(t!=GS);ti=0;
for(t=0;t<GS;t++){if(wg[t].vi<o){pv[ti].vis=labs(wg[t].vi);pv[ti].i=wg[t].i;ti++;}}
t2=0;
printf("please input WL~!n");
do
{
scanf("%f,%f,%d,%f",&wl[t2].pi,&wl[t2].qi,&wl[t2].i,&wl[t2].vi);
t2++;
}while(t2!=LS);
for(t=0;t<LS;t++){if(wl[t].vi<o){pv[ti].vis=labs(wl[t].vi);pv[ti].i=wl[t].i;ti++;}}
}
matrixform()
{
for(t=1;t<N;t++)
{
GII[t]=0;
BII[t]=0;
YDS[t]=0;
}
for(t2=1;t<zls;t2++)
{
ti=labs(zl[t2].i);
tj=labs(zl[t2].j);
tr=zl[t2].r;
tx=zl[t2].x;
temp=ldexp(tr,1)+ldexp(tx,1);
GIJ=tr/temp; BIJ=tx/temp;
y[t2][1]=-GIJ;
y[t2][2]=-BIJ;
y[t2][3]=tj;
GII[ti]=GII[ti]+GIJ; BII[ti]=BII[ti]+BIJ;
GII[tj]=GII[tj]+GIJ; BII[tj]=BII[tj]+BIJ;
YDS[ti]=YDS[ti]+1;
}
YDZ[1]=1;
for(t=1;t<N-1;t++)
{
YDZ[t+1]=YDZ[t]+YDS[t];
} /*矩阵型成第一部完成*/
/*矩阵型成第二部开始*/
for(t2=1;t<zls;t2++)
{ /*.k只用作变压器的标析，变压器变比仍在y中*/
ti=zl[t2].i;tj=zl[t2].j;YK=zl[t2].y;
if(ti<0||tj<0)
{ if(ti<0)
ti=labs(ti);
else
ti=labs(tj);
GIJ=y[t2][1];BIJ=y[t2][2];
GII[t2]=GII[t2]+(1-1/YK/YK)*GIJ;
BII[t2]=BII[t2]+(1-1/YK/YK)*BIJ;
y[t2][1]=GIJ/YK;
y[t2][2]=BIJ/YK;
}
else
GIJ=0;
BIJ=YK/2;
SY(tr); /*这个东东要调用，实现节点累计自导纳*/
SY(tj); /*SY的过程是完成向一个节点累计相应自导纳的实部和虚部*/
}
}
int sign,ld,k2,x,im,ai; /*k2控制台开关，负荷静态特性开关*/
static float fd[]={0};
unsigned AF[1];
static int u[]={0}; /*???????????怎么实现？来自那里???????*/
divisorform()
{
/*暂时不知道LD PVD 等的作用……待善*/
PVD=pv[0].i;
ld=wl[0].i;
t=0;
do{
t2++;
if(sign==1&&t2==PVD)
{t=t+1;pvd=pv[t].i;fd[t2]=0;di[t2]=0;
if(k2==0&&t2==ld)
{t2=t2+1;ld=wl[t2].i;}
}continue;
else
B[t2]=BII[t2];
if(k2==0&&sing==1&&t2==ld)
{
B[t2]=B[t2]+AF[1]*wl[t2].qi/wl[t2].vi/wl[t2].vi;t2=t2+1;ld=wl[t2].i;
}
for(temp=YDZ[t2];temp<YDZ[t2+1]-1;temp++)
{
tj=Y[temp][3];B[tj]=Y[temp][2];
}
if(sign=1)
{for(temp=1;temp<PVS;temp++)
tj=pv[temp][2];
B[tj]=0;
}
x=2;im=1;
do{im++;
if(im>t2-1)
break;
else
temp=1;
for(;temp!>fd[im];){if(u[x+1]!=1){temp=temp+1;x=x+2;}else ai=u[x]/} /*u[]未完成*/
continue;
}
}while(t2!=N-1);
}
dataout()
{
clrscr();
printf("note 1 voltagen");
printf("(.639696730300784) + j (1.832939) = 1.94136001255537 ∠ 70.7609880529659°n");
printf("87u& 婾[1]??u?孢???�u
--------------------------------------------------------------------------------
??虍鉧C&8u謤蛝髻??n");
}
main()
{
datain(); /*数据输入及处理*/
matrixform(); /*矩阵的形成*/
/* divisorform(); */ /*因子表的形成*/
matrixsolve(); /*矩阵线形方程的求解*/
/* nodepower(); */ /*迭代过程中节点功率的计算*/
/* iterate(); */ /*迭代*/
dataout(); /*数据输出及支路功率计算*/
}
牛顿拉夫逊法和PQ分解法的区别与联系是什么?求高人指点
答：1、60年代中期，基于导纳矩阵的牛顿—拉夫逊法。牛顿一拉夫逊法(简称牛顿法)是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中...
电力系统里的潮流如何计算? 有何方便的软件?需要输入什么参数?
答：电力系统中的潮流计算归结为求解一个用功率作为注入的非线性方程组。牛顿拉夫逊法和PQ分解法都是求解潮流方程的典型算法。这方面的软件有很多，国外、国内都有。在国内应用较广的如电科院的PSASP，中国版BPA都有潮流计算的功能。不同的软件在输入方式、数据文件格式和输入内容上都有差别。但直接与功率...
潮流计算论文开题报告
答：4. 直角坐标和极坐标的牛顿拉夫逊法在电力系统潮流计算中的应用. 5. 潮流计算的在MATLAB实现方法.6. 利用MATLAB软件编写计算机程序进行牛顿拉夫逊法潮流计算的分析. 7. 对编制好的程序选用算例进行检验，对检验的结果进行详细的分析总结 四、研究方案 1.编制形成电力网节点导纳矩阵的框图和程序 在掌握...
有谁给我讲下潮流计算?谢谢
答：从数学上说，潮流计算是求解一组由潮流方程（ 2）描述的非线性代数方程组。牛顿-拉夫逊方法是解非线性代数方程组的一种基本方法，在潮流计算中也得到应用。当采用了稀疏矩阵技术和节点优化编号技术后，牛顿-拉夫逊潮流算法成为电力系统潮流计算中的优秀算法，至今仍是各种潮流算法的基础。此外，还有各种快速...
交流电力系统牛额-拉夫逊法潮流计算的雅可比矩阵,一般具有下述哪些性 ...
答：【答案】：B、D 主要考查雅可比矩阵的特点:△f=J△x,称为雅可比矩阵。雅可比矩阵也是稀疏矩阵:但不是对称矩阵。由于节点电压有极坐标和直角坐标两种表示方式，因此牛顿-方程式也有两种形式，从而潮流计算的具体方法也有两种. C 选项是错的。雅克比矩阵是非奇异矩阵。
如何用matlab分析潮流计算
答：收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 三、设计原理 1． 牛顿-拉夫逊原理 牛顿迭代法是取x0&...
电力系统分析目录
答：2. 电力系统潮流计算方法 2.1 功率方程和网络方程 2.2 牛顿-拉夫逊潮流算法 2.3 快速解耦潮流算法和直流潮流算法 思考题：潮流算法在实际问题中的应用 3. 经济运行与优化 3.1 经济运行概念 3.4 最优潮流分析与计算 思考题：电力系统成本效益分析 4. 同步电机数学模型 ...
电力系统分析,利用 P—Q 分解法和牛顿—拉夫逊法进行潮流计算,二者的收 ...
答：收敛速度是指迭代次数，牛拉法的迭代次数比PQ法少，所以收敛速度快。不同情况两种方法收敛速度不同。牛顿—拉夫逊法比较通用，但是收敛速度不高，但基本所有问题都通用；P—Q 分解法适用于有P-Q能分解开的情况，适用面没有牛顿—拉夫逊法广，但是一旦可以适用，则收敛速度比较快。
潮流计算是什么?
答：潮流是电力系统中电压（各节点）、功率（有功、无功）（各支路）的稳态分布。潮流的分布是运行调度单位和维修部门所必须知道的事项。而潮流计算，是指给定电网中一些参数、已知值和未知值中假设的初始值，通过重复迭代，最终求出潮流分布的精确值，常用方法有牛顿-拉夫逊法和PQ分解法。在发电机母线上功率...
牛顿拉夫逊如何求功率判断收敛
答：牛顿拉夫逊使用牛顿迭代法求功率判断收敛。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x)=0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
对的吧。考虑负荷特性时，一般系统的频率保持不变，因此只考虑负荷的电压特性，即负荷吸收的有功和无功功率随其端电压的波动而变化。所以考虑负荷静特性时，采用把负荷功率当做该点电压的线性函数和非线性函数。潮流计算中，计算功率误差时的P0和修正方程式系数矩阵中有关元素都应是根据负荷节点电压而改变的
牛顿-拉夫逊法早在50年代末就已应用于求解电力系统潮流问题，但作为一种实用的，有竞争力的电力系统潮流计算方法，则是在应用了稀疏矩阵技巧和高斯消去法求修正方程后。牛顿-拉夫逊法是求解非线性代数方程有效的迭代计算。 1.3MATLAB概述 目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单。 MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。 MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。 另外，MATLAB提供了一种特殊的工具：工具箱（TOOLBOXES）.这些工具箱主要包括：信号处理（SIGNAL PROCESSING）、控制系统（CONTROL SYSTEMS）、神经网络（NEURAL NETWORKS）、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。 MATLAB设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。 原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。 2．1 电力系统的基本概念 2．1．1电力系统 （1）电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。 （2）电力网：电力系统中除发电机和用电设备外的部分。 （3）动力系统：电力系统和“动力部分”的总和。 “动力部分”：包括火力发电厂的锅炉、汽轮机、热力网和用电设备，水力发电厂的水库和水轮机，核电厂的反应堆等。 2．1．2电力系统的负荷和负荷曲线 （1）电力系统的负荷：系统中千万个用电设备消费功率的总和，包括异步电动机、同步电动机、电热炉、整流设备、照明设备等若干类。 （2）电力系统的供电负荷：综合用电负荷加上电力网中损耗的功率。 （3）电力系统的发电负荷：供电负荷加上发电厂本身的消耗功率。 （4）各用电设备的有功功率和无功功率随受电电压和系统频率的变化而变化，其变化规律不尽相同，综合用电负荷随电压和频率的变化规律是各用电负荷变化规律的合成。 （5）负荷曲线：某一时间段内负荷随时间而变化的规律。 （6）按负荷种类可分有功功率负荷和无功功率负荷；按时间长短可分为日负荷和年负荷曲线；按计量地点可分为个别用户、电力线路、变电所、发电厂以至整个系统的负荷曲线。将上述三种分类相结合，就确定了某一种特定的负荷曲线。不同行业的有功功率日负荷曲线差别很大。负荷曲线对电力系统的运行又很重要的意义，它是安排日发电计划，确定各发电厂发电任务以及确定系统运行方式等的重要依据。 2．2 电力系统的基本元件 2．2．1 发电机 现代电力工业中，无论是火力发电、水力发电或核能发电，几乎全部采用同步交流发电机。电机的电枢布置在定子上，励磁绕组布置在转子上，作为旋转式磁极。同步发电机的转速（转/MIN）和系统频率f（HZ）之间有着严格的关系，即n=60f/p式中p为电机的极对数。 根据转子结构型式的不同，分为隐极式和凸极式发电机，前者转子没有显露出来的磁极，后者则有。 转子的励磁型式有直流励磁系统和可控硅励磁系统，后者利用同轴交流励磁机或由同步
牛顿-拉夫逊法早在50年代末就已应用于求解电力系统潮流问题，但作为一种实用的，有竞争力的电力系统潮流计算方法，则是在应用了稀疏矩阵技巧和高斯消去法求修正方程后。牛顿-拉夫逊法是求解非线性代数方程有效的迭代计算。1.3MATLAB概述
目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单。
MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。
MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。
另外，MATLAB提供了一种特殊的工具：工具箱（TOOLBOXES）.这些工具箱主要包括：信号处理（SIGNAL PROCESSING）、控制系统（CONTROL SYSTEMS）、神经网络（NEURAL NETWORKS）、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。
MATLAB设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。
原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。
2．1 电力系统的基本概念
2．1．1电力系统
（1）电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线路及各种用电设备等联系在一起组成的统一整体称为电力系统。
（2）电力网：电力系统中除发电机和用电设备外的部分。
（3）动力系统：电力系统和“动力部分”的总和。
“动力部分”：包括火力发电厂的锅炉、汽轮机、热力网和用电设备，水力发电厂的水库和水轮机，核电厂的反应堆等。
2．1．2电力系统的负荷和负荷曲线
（1）电力系统的负荷：系统中千万个用电设备消费功率的总和，包括异步电动机、同步电动机、电热炉、整流设备、照明设备等若干类。
（2）电力系统的供电负荷：综合用电负荷加上电力网中损耗的功率。
（3）电力系统的发电负荷：供电负荷加上发电厂本身的消耗功率。
（4）各用电设备的有功功率和无功功率随受电电压和系统频率的变化而变化，其变化规律不尽相同，综合用电负荷随电压和频率的变化规律是各用电负荷变化规律的合成。
（5）负荷曲线：某一时间段内负荷随时间而变化的规律。
（6）按负荷种类可分有功功率负荷和无功功率负荷；按时间长短可分为日负荷和年负荷曲线；按计量地点可分为个别用户、电力线路、变电所、发电厂以至整个系统的负荷曲线。将上述三种分类相结合，就确定了某一种特定的负荷曲线。不同行业的有功功率日负荷曲线差别很大。负荷曲线对电力系统的运行又很重要的意义，它是安排日发电计划，确定各发电厂发电任务以及确定系统运行方式等的重要依据。
2．2 电力系统的基本元件
2．2．1 发电机
现代电力工业中，无论是火力发电、水力发电或核能发电，几乎全部采用同步交流发电机。电机的电枢布置在定子上，励磁绕组布置在转子上，作为旋转式磁极。同步发电机的转速（转/MIN）和系统频率f（HZ）之间有着严格的关系，即n=60f/p式中p为电机的极对数。
根据转子结构型式的不同，分为隐极式和凸极式发电机，前者转子没有显露出来的磁极，后者则有。
转子的励磁型式有直流励磁系统和可控硅励磁系统，后者利用同轴交流励磁机或由同步发电机本身发出的交流电，经整流后供给转子。直流励磁机有换向问题，故其制造容量受到限制，所以，在大容量发电机中均可采用可控硅励磁系统。
2．2．2 电力变压器
电力变压器是电力系统中广泛使用的升压和降压设备。据统计，电力系统中变压器的安装总容量约为发电机安装容量的6－8倍。按用途，电力变压器可分为升压变压器、降亚变压器、配电变压器和联络变压器。按相数分，变压器可分为单相式和三相式。按每相线圈分，又有双绕组和三绕组之分。按线圈耦合的方式，可分为普通变压器和自耦变压器。
2．2．3 电力线路
（1）架空线路：由导线、避雷针、杆塔、绝缘子和金具等构成。
（2）电缆线路：由导线、绝缘层、包护层等构成。
2．2．4 无功功率补偿设备
主要的无功功率补偿设备有同步调相机、电力电容器和静止补偿器。
2．3 电力系统元件的数学模型
2．3．1 电力线路的等值电路
在电力系统分析中，一般只考虑电力线路两侧端口的电压和电流，把电力线路作为无源双口网络处理。
线路的双口网络方程：
Z=B=*L*
2．3．2 变压器的等值电路
（1）双绕组变压器等值电路
（2）三绕组变压器等值电路
2．3．3 同步发电机的数学模型
2．3．4 电力系统负荷
2．3．5 多级电压电力系统的等值电路
2．4 电力系统稳态运行分析
2．4．1 电力线路的电压损耗与功率损耗
2．4．2 变压器中的功率损耗与电压损耗
2．4．3 辐射形网络的分析计算
辐射形电力网的特点是各条线路有明确的始端与末端。辐射形电力网的分析计算就是利用已知的负荷、节点电压来求取未知的节点电压、线路功率分布、功率损耗及始端输出功率。
辐射形电力网的分析计算，根据已知条件的不同分两种
1 已知末端功率与电压：即 从末端逐级往上推算，直至求得各要求的量。
2 已知末端功率、始端电压：末端可理解成一负荷点，始端为电源点或电压中枢点。采用迭代法。
（1）假设末端电压为线路额定电压，利用第一种方法求得始端功率及全网功率分布。
（2）用求得的线路始端功率和已知的线路始端电压，计算线路末端电压和全网功率分布。
（3）用第（2）步求得的线路末端电压计算线路始端功率和全网功率分布，如求得的各线路功率与前一次相同计算的结果相差小于允许值，就可以认为本步求得的线路电压和全网功率分布为最终计算结果。否则，返回第二步重新进行计算。
2．4．4 复杂电力系统潮流计算
电力系统潮流计算始对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标始求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。因此，潮流计算是电力系统计算分析中的一种最基本的计算。
潮流计算结果的用途，例如用于电力系统稳定研究、安全估计或最优潮流等也对潮流计算的模型和方法有直接影响。
2．5 电力系统潮流计算机算法
2．5．1电力系统潮流计算机算法概述
2．5．1．1 导纳矩阵的形成
2．5．1．2 节点类型
（1）PV节点：柱入有功功率P为给定值，电压也保持在给定数值。
（2）PQ节点：诸如有功功率和无功功率是给定的。
（3）平衡节点：用来平衡全电网的功率。选一容量足够大的发电机担任平衡全电网功率的职责。
平衡节点的电压大小与相位是给定的，通常以它的相角为参考量，即取其电压相角为0。一个独立的电力网中只设一个平衡点。
2．5．1．3 高斯迭代法
2．5．2 牛顿－拉夫逊法
2．5．2．1 原理
2．5．2．2 基本步骤
基本步骤：
（1）形成节点导纳矩阵
（2）将各节点电压设初值U，
（3）将节点初值代入相关求式，求出修正方程式的常数项向量
（4）将节点电压初值代入求式，求出雅可比矩阵元素
（5）求解修正方程，求修正向量
（6）求取节点电压的新值
（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步
（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。
2．5．2．3 注意事项
2．5．2．4 程序流程框图
2．6 软件设计
2．6．1 方案选择及说明
2．6．2 方案求解
2．6．3 MATLAB编程说明及元件描述
2．6．4 程序
#include<stdio.h>
struct powernode
{
float pi;
float qi;
int i;
float vi;
};
struct powernode wg[20];
struct powernode wl[20];
struct linedata
{
int i;
int j;
float r;
float x;
float y; /*包括变压器变比*/
float k; /*只用作标析变压器，变压器变比仍在y中*/
};
struct linedata zl[20];
struct linedata t3; /*临时数组*/
static double y[][3]; /*在matrixform中应用*/
int t=0;
int t2,ti,tj; /*临时记数单元*/
float temp;
float tx,tr,YK; /*中间工作单元(在matrixform中应用)*/
double GIJ,BIJ; /*中间工作单元(在matrixform中应用)*/
int N; /*总节点数*/
int zls;
int Q,V,PVS,PVD;
int GS;
int LS;
float vo;
float Eps;
static double GII[]={0},BII[]={0},YDS[]={0},YDZ[]={0},B[]={0};/*添加数组*/
/*因子表形成时定义的数据*/
struct pvdata
{
float vis;
int i;
};
static struct pvdata pv[]={0};
datain()
{
clrscr();
printf("program runningn" );
printf("n");
printf("please input the aggregate to the system note");/*总节点数*/
scanf("%d",&N);
printf("n");
printf(" PQ note IN ALL?");/*总节点数*/
scanf("%d",&Q);
PVS=(N-Q)-1;
printf("n");
printf("them input the aggregate to the system power line");
scanf("%d",&zls);/*输电线路数和变压器的总数*/
printf("n");
printf("electromotor node in all :?");/*发电机节点总数*/
scanf("%d",&GS);
printf("n");
printf("load node in all : ?");/*负荷节点总数*/
scanf("%d",&LS);
printf("n");
printf("average electric voltage");/*平均电压*/
scanf("%f",vo);
printf("n");
printf("n");
printf("please input the date messagen");
printf("follow the format like it: i,j,r,x,y,kn");
do{
t++;
scanf("%d,%d,%f,%f,%f",&zl[t].i,&zl[t].j,&zl[t].r,&zl[t].x,&zl[t].y,&zl[t].k);
printf("processing....n");
if(zl[t].i>zl[t].j)
{
temp=zl[t].i;
zl[t].i=zl[t].j;
zl[t].j=temp;
/* if(zl[t].k!=1) */ /*要考虑归算问题不？？？？*/
}
printf("data you input is:n " );
printf("%d,%d,%f,%f,%f",zl[t].i,zl[t].j,zl[t].r,zl[t].x,zl[t].y,zl[t].k);
}while(zl[t].i!=0&&zl[t].j==0);
for(t2=t;t>0;t--) /*冒泡法排序*/
{
for(;t2>0;t2--)
{
if(zl[t2].i<zl[t2-1].i)
{
t3.i=zl[t2].i;t3.j=zl[t2].j;t3.r=zl[t2].r;t3.x=zl[t2].x;t3.y=zl[t2].y;t3.k=zl[t2].k;
zl[t2].i=zl[t2-1].i;zl[t2].j=zl[t2-1].j;zl[t2].r=zl[t2-1].r;zl[t2].x=zl[t2-1].x;zl[t2].y=zl[t2-1].y;zl[t2].k=zl[t2-1].k;
zl[t2-1].i=t3.i;zl[t2-1].j=t3.j;zl[t2-1].r=t3.r;zl[t2-1].x=t3.x;zl[t2-1].y=t3.y;zl[t2-1].k=t3.k;
}
else if(zl[t2].i==zl[t2-1].i)
{if(zl[t2].j<zl[t2-1].j)
{
t3.i=zl[t2].i;t3.j=zl[t2].j;t3.r=zl[t2].r;t3.x=zl[t2].x;t3.y=zl[t2].y;t3.k=zl[t2].k;
zl[t2].i=zl[t2-1].i;zl[t2].j=zl[t2-1].j;zl[t2].r=zl[t2-1].r;zl[t2].x=zl[t2-1].x;zl[t2].y=zl[t2-1].y;zl[t2].k=zl[t2-1].k;
zl[t2-1].i=t3.i;zl[t2-1].j=t3.j;zl[t2-1].r=t3.r;zl[t2-1].x=t3.x;zl[t2-1].y=t3.y;zl[t2-1].k=t3.k;
}
}
}
}
printf("n");
t=0;
printf("please input wg~!n");
do
{
scanf("%f,%f,%d,%f",&wg[t].pi,&wg[t].qi,&wg[t].i,&wg[t].vi);
t++;
}while(t!=GS);ti=0;
for(t=0;t<GS;t++){if(wg[t].vi<o){pv[ti].vis=labs(wg[t].vi);pv[ti].i=wg[t].i;ti++;}}
t2=0;
printf("please input WL~!n");
do
{
scanf("%f,%f,%d,%f",&wl[t2].pi,&wl[t2].qi,&wl[t2].i,&wl[t2].vi);
t2++;
}while(t2!=LS);
for(t=0;t<LS;t++){if(wl[t].vi<o){pv[ti].vis=labs(wl[t].vi);pv[ti].i=wl[t].i;ti++;}}
}
matrixform()
{
for(t=1;t<N;t++)
{
GII[t]=0;
BII[t]=0;
YDS[t]=0;
}
for(t2=1;t<zls;t2++)
{
ti=labs(zl[t2].i);
tj=labs(zl[t2].j);
tr=zl[t2].r;
tx=zl[t2].x;
temp=ldexp(tr,1)+ldexp(tx,1);
GIJ=tr/temp; BIJ=tx/temp;
y[t2][1]=-GIJ;
y[t2][2]=-BIJ;
y[t2][3]=tj;
GII[ti]=GII[ti]+GIJ; BII[ti]=BII[ti]+BIJ;
GII[tj]=GII[tj]+GIJ; BII[tj]=BII[tj]+BIJ;
YDS[ti]=YDS[ti]+1;
}
YDZ[1]=1;
for(t=1;t<N-1;t++)
{
YDZ[t+1]=YDZ[t]+YDS[t];
} /*矩阵型成第一部完成*/
/*矩阵型成第二部开始*/
for(t2=1;t<zls;t2++)
{ /*.k只用作变压器的标析，变压器变比仍在y中*/
ti=zl[t2].i;tj=zl[t2].j;YK=zl[t2].y;
if(ti<0||tj<0)
{ if(ti<0)
ti=labs(ti);
else
ti=labs(tj);
GIJ=y[t2][1];BIJ=y[t2][2];
GII[t2]=GII[t2]+(1-1/YK/YK)*GIJ;
BII[t2]=BII[t2]+(1-1/YK/YK)*BIJ;
y[t2][1]=GIJ/YK;
y[t2][2]=BIJ/YK;
}
else
GIJ=0;
BIJ=YK/2;
SY(tr); /*这个东东要调用，实现节点累计自导纳*/
SY(tj); /*SY的过程是完成向一个节点累计相应自导纳的实部和虚部*/
}
}
int sign,ld,k2,x,im,ai; /*k2控制台开关，负荷静态特性开关*/
static float fd[]={0};
unsigned AF[1];
static int u[]={0}; /*???????????怎么实现？来自那里???????*/
divisorform()
{
/*暂时不知道LD PVD 等的作用……待善*/
PVD=pv[0].i;
ld=wl[0].i;
t=0;
do{
t2++;
if(sign==1&&t2==PVD)
{t=t+1;pvd=pv[t].i;fd[t2]=0;di[t2]=0;
if(k2==0&&t2==ld)
{t2=t2+1;ld=wl[t2].i;}
}continue;
else
B[t2]=BII[t2];
if(k2==0&&sing==1&&t2==ld)
{
B[t2]=B[t2]+AF[1]*wl[t2].qi/wl[t2].vi/wl[t2].vi;t2=t2+1;ld=wl[t2].i;
}
for(temp=YDZ[t2];temp<YDZ[t2+1]-1;temp++)
{
tj=Y[temp][3];B[tj]=Y[temp][2];
}
if(sign=1)
{for(temp=1;temp<PVS;temp++)
tj=pv[temp][2];
B[tj]=0;
}
x=2;im=1;
do{im++;
if(im>t2-1)
break;
else
temp=1;
for(;temp!>fd[im];){if(u[x+1]!=1){temp=temp+1;x=x+2;}else ai=u[x]/} /*u[]未完成*/
continue;
}
}while(t2!=N-1);
}
dataout()
{
clrscr();
printf("note 1 voltagen");
printf("(.639696730300784) + j (1.832939) = 1.94136001255537 ∠ 70.7609880529659°n");
printf("87u& 婾[1]??u?孢???�u
--------------------------------------------------------------------------------
??虍鉧C&8u謤蛝髻??n");
}
main()
{
datain(); /*数据输入及处理*/
matrixform(); /*矩阵的形成*/
/* divisorform(); */ /*因子表的形成*/
matrixsolve(); /*矩阵线形方程的求解*/
/* nodepower(); */ /*迭代过程中节点功率的计算*/
/* iterate(); */ /*迭代*/
dataout(); /*数据输出及支路功率计算*/
}
答：1、60年代中期，基于导纳矩阵的牛顿—拉夫逊法。牛顿一拉夫逊法(简称牛顿法)是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中...
答：电力系统中的潮流计算归结为求解一个用功率作为注入的非线性方程组。牛顿拉夫逊法和PQ分解法都是求解潮流方程的典型算法。这方面的软件有很多，国外、国内都有。在国内应用较广的如电科院的PSASP，中国版BPA都有潮流计算的功能。不同的软件在输入方式、数据文件格式和输入内容上都有差别。但直接与功率...
答：4. 直角坐标和极坐标的牛顿拉夫逊法在电力系统潮流计算中的应用. 5. 潮流计算的在MATLAB实现方法.6. 利用MATLAB软件编写计算机程序进行牛顿拉夫逊法潮流计算的分析. 7. 对编制好的程序选用算例进行检验，对检验的结果进行详细的分析总结 四、研究方案 1.编制形成电力网节点导纳矩阵的框图和程序 在掌握...
答：从数学上说，潮流计算是求解一组由潮流方程（ 2）描述的非线性代数方程组。牛顿-拉夫逊方法是解非线性代数方程组的一种基本方法，在潮流计算中也得到应用。当采用了稀疏矩阵技术和节点优化编号技术后，牛顿-拉夫逊潮流算法成为电力系统潮流计算中的优秀算法，至今仍是各种潮流算法的基础。此外，还有各种快速...
答：【答案】：B、D 主要考查雅可比矩阵的特点:△f=J△x,称为雅可比矩阵。雅可比矩阵也是稀疏矩阵:但不是对称矩阵。由于节点电压有极坐标和直角坐标两种表示方式，因此牛顿-方程式也有两种形式，从而潮流计算的具体方法也有两种. C 选项是错的。雅克比矩阵是非奇异矩阵。
答：收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 三、设计原理 1． 牛顿-拉夫逊原理 牛顿迭代法是取x0&...
答：2. 电力系统潮流计算方法 2.1 功率方程和网络方程 2.2 牛顿-拉夫逊潮流算法 2.3 快速解耦潮流算法和直流潮流算法 思考题：潮流算法在实际问题中的应用 3. 经济运行与优化 3.1 经济运行概念 3.4 最优潮流分析与计算 思考题：电力系统成本效益分析 4. 同步电机数学模型 ...
答：收敛速度是指迭代次数，牛拉法的迭代次数比PQ法少，所以收敛速度快。不同情况两种方法收敛速度不同。牛顿—拉夫逊法比较通用，但是收敛速度不高，但基本所有问题都通用；P—Q 分解法适用于有P-Q能分解开的情况，适用面没有牛顿—拉夫逊法广，但是一旦可以适用，则收敛速度比较快。
答：潮流是电力系统中电压（各节点）、功率（有功、无功）（各支路）的稳态分布。潮流的分布是运行调度单位和维修部门所必须知道的事项。而潮流计算，是指给定电网中一些参数、已知值和未知值中假设的初始值，通过重复迭代，最终求出潮流分布的精确值，常用方法有牛顿-拉夫逊法和PQ分解法。在发电机母线上功率...
答：牛顿拉夫逊使用牛顿迭代法求功率判断收敛。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x)=0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
