电力系统的潮流分布计算计划书
电力系统分析课程设计——牛顿拉夫逊潮流计算
电力系统编程潮流计算1 设计任务及初步分析1.1 设计任务条件：节点数：3 支路数：3 计算精度：0.00010支路1： 0.0300+j0.09001┠—————□—————┨2支路2： 0.0200+j0.09002┠—————□—————┨3支路3： 0.0300+j0.09003┠—————□—————┨1节点1：PQ 节点，S （1）=-0.5000-j0.2000 节点2：PQ 节点，S （2）=-0.6000-j0.2500 节点3：平衡节点，U （3）=1.0000∠0.0000 要求：编写程序计算潮流1.2 初步分析潮流计算在数学上可归结为求解非线性方程组，其数学模型简写如下：⎪⎪⎩⎪⎪⎨⎧=⋅⋅⋅⋅⋅⋅=⋅⋅⋅=⋅⋅⋅0)(0)(0)(21212211n n n n x x x f x x x f x x x f ，，，，，，，，，2 牛顿-拉夫逊法简介2.1概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。
2．2 一般概念对于非线性代数方程组()0=x f即 ()0,,,21=n i x x x f ()n i ,2,1= (2－1)在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组()()()()()0000=∆'+x x f x f (2－2)上式称之为牛顿法的修正方程式。
电力系统分析课程设计——电力系统潮流计算
信息工程学院课程设计报告书题目: 电力系统潮流计算专业：电气工程及其自动化班级：0310406学号：031040635学生姓名：陈代才指导教师：钟建伟2013年 4 月15 日信息工程学院课程设计任务书2013年4月15日目录1 任务提出与方案论证 (2)2 总体设计 (3)2.1潮流计算等值电路 (3)2.2建立电力系统模型 (3)2.3模型的调试与运行 (3)3 详细设计 (4)3.1 计算前提 (4)3.2手工计算 (7)4设计图及源程序 (11)4.1MA TLAB仿真 (11)4.2潮流计算源程序 (11)5 总结 (19)参考文献 (20)1 任务提出与方案论证潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法,是电力系统规划和运营中不可缺少的一个重要组成部分。
可以说,它是电力系统分析中最基本、最重要的计算,是系统安全、经济分析和实时控制与调度的基础。
常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。
潮流计算的结果是电力系统稳定计算和故障分析的基础。
在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。
同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。
因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。
在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。
是电力系统研究人员长期研究的一个课题。
它既是对电力系统规划设计和运行方式的合理性、可靠性及经济性进行定量分析的依据，又是电力系统静态和暂态稳定计算的基础。
潮流计算经历了一个由手工到应用数字电子计算机的发展过程，现在的潮流算法都以计算机的应用为前提用计算机进行潮流计算主要步骤在于编制计算机程序，这是一项非常复杂的工作。
电力系统分布式潮流计算
目录1.课题背景 (1)1.1 理论背景 (1)1.1.1 电力系统静态等值理论 (1)1.1.2 电力系统分块计算理论 (1)1.1.3 分布式潮流计算 (1)1.1.4 计算节点（导纳阵）拼接 (2)1.1.5 思路的特点比较 (2)1.2工程背景 (3)2.研究内容 (3)2.1 平等地位的分布式潮流计算研究内容 (4)2.2计算节点拼接思路研究内容 (4)3.研究思路 (4)3.1 平等地位的分布式潮流计算研究思路 (4)3.1.1 使用导纳矩阵来简化协调计算的思路 (4)3.1.1.1 思路来源 (5)3.1.1.2 思路介绍 (5)3.1.1.3实验结果 (5)3.1.2使用最小二乘估计理论来简化协调计算的思路 (5)3.1.2.1 思路来源 (5)3.1.2.2思路介绍 (6)3.1.2.3 初步实验及结果 (6)3.1.2.4 思路优缺点 (7)3.1.3使用灵敏度分析减少外层迭代次数和计算结果误差的思路 (8)3.1.3.1 思路来源 (8)3.1.3.2 算法介绍 (8)3.1.3.3 算法难点 (10)3.2 计算节点（导纳阵）拼接研究思路 (10)4.课题工作进度安排 (11)附：大扰动情况下的数值实验 (12)互联电网分布式动态潮流计算的算法研究1. 课题背景1.1 理论背景随着电力系统互联不断加强，为提高系统稳态分析水平，研究现有技术条件下，如何在分层调度中心之间交换少量信息，就能使子网分别独立潮流计算并达到或接近和全网一体化潮流计算同样的效果，是一个有实际意义的课题。
如何使子网分别独立计算达到或接近和全网一体化计算同样的效果，目前看来有如下几种思路。
1.1.1 电力系统静态等值理论使用电力系统静态等值的思路。
常见的有W ARD等值法和REI等值法。
目前广泛使用的W ARD等值法又主要包括W ARD节点注入法、解藕ＷＡＲＤ等值、扩展ＷＡＲＤ等值、缓冲网等值等。
这些方法或者是对部分重要的ＰＶ点进行了保留，或者是通过等值支路和虚拟ＰＶ点，或者通过保留一个或若干个母线层来模拟基态和内网发生扰动后外网的功率响应，在一些场合下取得了令人满意的效果。
潮流计算课程设计-电力系统潮流计算
课程设计(论文)题目名称电力系统潮流计算课程名称电力系统稳态分析学生姓名学号系、专业指导教师2014 年 1 月 5 日课程设计（论文）任务书注：1．此表由指导教师填写，经系、教研室审批，指导教师、学生签字后生效；2．此表1式3份，学生、指导教师、教研室各1份。
指导教师（签字）：学生（签字）：邵阳学院课程设计（论文）评阅表学生姓名学号系专业班级题目名称潮流计算课程设计课程名称二、指导教师评定注：1、本表是学生课程设计（论文）成绩评定的依据，装订在设计说明书（或论文）的“任务书”页后面；2、表中的“评分项目”及“权重”根据各系的考核细则和评分标准确定。
目录摘要 (I)1课题内容要求及目的 (1)1.1课题背景 (1)1.2课题意义 (2)1.3课题要求及内容 (3)1.4课题目的 (4)2 潮流计算步骤与原理 (5)2.1潮流计算流程图 (5)2.2潮流计算步骤 (6)3 方案设计 (7)3.1系统框图设计 (7)4 仿真调试及实验分析 (8)总结 (11)参考文献 (12)致谢 (13)摘要电力系统的主体结构有电源、电力网络和负荷中心。
电源指各类发电厂，它将一次能源转换成电能；电力网络由电源的升压变电所、输电线路、负荷中心变电所、配电线路等构成。
它的功能是将电源发出的电能升压到一定等级后输送到负荷中心变电所，再降压至一定等级后，经配电线路与用户相联。
电力系统中网络结点千百个交织密布，有功潮流、无功潮流、高次谐波、负序电流等以光速在全系统范围传播。
它既能输送大量电能，创造巨大财富，也能在瞬间造成重大的灾难性事故。
实际电力系统的潮流计算主要采用牛顿-拉夫逊法。
按电压的不同表示方法，牛顿-拉夫逊潮流计算分为直角坐标形式和极坐标形式两种。
本次计算采用极坐标形式下的牛顿-拉夫逊法，牛顿-拉夫逊法有很好的收敛性，但要求有合适的初值。
本设计采用电力系统仿真软件PSCAD，可以直观地看出电力系统运行时的潮流分布，从而完成课程设计的要求。
第三章简单电力系统的潮流分布计算(李庚银新教材)
❖ 有发出,有吸收,构成无功功率平衡.
❖ 同步电机发出感性无功,异步电机吸收感性无 功.因此电力系统中,发电机是同步的,电动机 是异步的.
❖ 如果发电机、电动机都做出异步的，或都做 出同步的会怎么样？
《辞海》对无功功率的解释
❖ 具有电感和电容的交流电路中，电感的磁场 或电容的电场在一个周期内的一部分时间内 从电源吸收能量，另一部分时间内将能量返 回电源。在整个周期内平均功率为零，也就 是没有能量消耗。但能量是在电源和电感或 电容之间来回交换的。能量交换率的最大值 叫做无功功率。
电容上的无功叫消耗容性无功，或发出感性无功。
❖ 无功功率是电感或电容吞吐有功功率的最大 值,是标量,相当于容器的容量;
❖ 但是,由于电感与电容吞吐有功功率步调相反, 因此将无功功率赋予方向:
电容元件吞吐有功功率的过程,称为发出感性无 功功率,是感性无功功率电源;
电感元件吞吐有功功率的过程,称为吸收感性无 功功率,是感性无功功率负荷.
S3 U3 式①
S23
式② y30 U3
S4 U 4 U34 式③ δU34
QyT
I0 % 100
SN
S2 SN
PZT
Pk 1000
QZT
Uk % 100
SN
PyT
P0 1000
QyT
I0 % 100
SN
S1 SN
PZT
Pk 1000
QZT
Uk % 100
SN
PyT
P0 1000
QyT
I0 % 100
SN
实际功率是S时
ST
PyT
jQyT
PZT
S SN
3Z
S 2
3U
电力系统潮流分析与计算设计
电力系统潮流分析与计算设计电力系统潮流分析是电力系统运行和规划的重要工具之一，用于计算和分析电力系统中各个节点的电压幅值和相角。
潮流分析可以帮助我们了解电力系统中电流、功率等重要变量的分布和流动情况，进而指导电力系统的运行和调度，提高电力系统的稳定性和经济性。
在潮流分析中，我们常常使用P-Q分解法来计算电力系统中各个节点的电压幅值和相角。
P-Q分解法是一种迭代的方法，可以逐步计算出节点的电压幅值和相角，并满足节点的功率平衡条件。
首先，我们需要列出发电机节点和负荷节点的功率平衡方程。
对于发电机节点，功率平衡方程可以表示为：P_g-P_d-P_l=0其中，P_g表示发电机的有功出力，P_d表示发电机的有功损耗，P_l 表示节点的有功负荷。
对于负荷节点，功率平衡方程可以表示为：-P_l+P_g-2BVL=0其中，B表示节点的导纳，V表示节点的电压幅值，L表示节点的电流幅值。
然后，我们需要列出节点的电压-相角方程。
对于每个节点，电压-相角方程可以表示为：V_i ∠θi = Vj ∠θj - (Rij + Xij)Iij其中，V_i和θ_i表示节点i的电压幅值和相角，V_j和θ_j表示节点j的电压幅值和相角，R_ij和X_ij表示节点i和节点j之间的电阻和电抗，I_i_j表示节点i到节点j的电流。
在实际计算中，我们通常从平衡节点开始，逐步计算各个节点的电压幅值和相角，直到所有节点的电压幅值和相角满足要求。
在进行计算前，我们需要先给定节点的电压幅值和相角的初值。
对于平衡节点，我们可以直接给定其电压幅值和相角的固定值，而对于其他节点，我们可以先给定其电压幅值的初值，然后通过P-Q分解法来计算其相角。
在每一步计算中，我们首先根据电压-相角方程计算出节点之间的电流，然后再根据功率平衡方程计算出节点的有功出力和有功负荷，最后根据节点的有功出力和有功负荷来更新节点的电压幅值和相角。
通过多次迭代计算，直到节点的电压幅值和相角满足要求，我们就可以得到电力系统的潮流分析结果。
电力系统分析潮流计算最终完整版
电力系统分析潮流计算实验报告姓名：XXXXXX 学号：XXXXXXXXXX 班级：XXXXXXXX一、实验目的掌握潮流计算计算机算法的方法，熟悉MATLAB的程序调试方法。
二、实验准备根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。
三、实验要求每人一组，在实验课时内，调试和修改运行程序，用算例计算输出潮流结果。
四、程序流程五、实验程序%本程序的功能是用牛拉法进行潮流计算%原理介绍详见鞠平著《电气工程》%默认数据为鞠平著《电气工程》例8.4所示数据%B1是支路参数矩阵%第一列和第二列是节点编号。
节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号%第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗%第四列为支路的对地导纳参数，含变压器支路此值不代入计算%第五烈为含变压器支路的变压器的变比，变压器非标准电压比%第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器%B2为节点参数矩阵%第一列为节点注入发电功率参数%第二列为节点负荷功率参数%第三列为节点电压参数%第四列%第五列%第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数%X为节点号和对地参数矩阵%第一列为节点编号%第二列为节点对地参数%默认算例% n=4;% n1=4;% isb=4;% pr=0.00001;% B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];% B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];% X=[1 0;2 0.05i;3 0;4 0];clear;clc;num=input('是否采用默认数据？(1-默认数据;2-手动输入)');if num==1n=4;n1=4;isb=4;pr=0.00001;B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];X=[1 0;2 0.05i;3 0;4 0];elsen=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X=');endTimes=1; %迭代次数%创建节点导纳矩阵Y=zeros(n);for i=1:n1if B1(i,6)==0 %不含变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);else %含有变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-B1(i,5)/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3);Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))/B1(i,3);endendfor i=1:n1Y(i,i)=Y(i,i)+X(i,2); %计及补偿电容电纳enddisp('导纳矩阵为：');disp(Y); %显示导纳矩阵%初始化OrgS、DetaSOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%创建OrgS，用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2 %不是平衡点&是PQ点h=h+1;for j=1:n%公式8-74%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j ,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3 %不是平衡点&是PV点h=h+1;for j=1:n%公式8-75-a%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend%创建PVU 用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;for i=1:nif B2(i,6)==3t=t+1;PVU(t,1)=B2(i,3);endend%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1); %delPiDetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1); %delQiendendt=0;for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,1))-OrgS(2*h-1,1); %delPiDetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2; %delUi endend% DetaS%创建I，用于存储节点电流参数i=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));%conj求共轭endend%创建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %对PQ节点的处理if B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendk=0;for i=1:n %对PV节点的处理if B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendenddisp('初始雅可比矩阵为：');disp(Jacbi);%求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS; %inv矩阵求逆% DetaU%修正节点电压j=0;for i=1:n %对PQ节点处理if B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:n %对PV节点的处理if B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endend% B2%开始循环**********************************************************************while abs(max(DetaU))>prOrgS=zeros(2*n-2,1);h=0;j=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:nif i~=isb&B2(i,6)==3h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend% OrgS%创建DetaSh=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);endendt=0;for i=1:nif i~=isb&B2(i,6)==3h=h+1;t=t+1;% DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h-1,1)=real(B2(i,1))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;endend% DetaS%创建Ii=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));endend% I%创建JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:nif B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendend% JacbiDetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;% DetaU%修正节点电压j=0;for i=1:nif B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:nif B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endend% B2Times=Times+1; %迭代次数加1enddisp('迭代次数为：');disp(Times);disp('收敛时电压修正量为：：');disp(DetaU);for k=1:nE(k)=B2(k,3);e(k)=real(E(k));f(k)=imag(E(k));V(k)=sqrt(e(k)^2+f(k)^2);sida(k)=atan(f(k)./e(k))*180./pi;end%=============== 计算各输出量=========================== disp('各节点的实际电压标幺值E为(节点号从小到大排列)：'); disp(E); %显示各节点的实际电压标幺值E用复数表示disp('-----------------------------------------------------')disp('各节点的电压大小V为(节点号从小到大排列)：');disp(V); %显示各节点的电压大小V的模值disp('-----------------------------------------------------');disp('各节点的电压相角sida为(节点号从小到大排列)：');disp(sida); %显示各节点的电压相for p=1:nfor q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点的注入电流的共轭值endS(p)=E(p)*C(p); %计算各节点的功率S = 电压X 注入电流的共轭值enddisp('各节点的功率S为(节点号从小到大排列)：');disp(S); %显示各节点的注入功率Sline=zeros(n1,5);disp('-----------------------------------------------------');disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');for i=1:n1p=B1(i,1);q=B1(i,2);Sline(i,1)=B1(i,1);Sline(i,2)=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*((1-B1(i,5))/B1(i,3))+(conj(E(p))-conj(E(q)))*(B1(i,5)/B1(i,3)));Siz(i)=Si(p,q);endSSi(p,q)=Si(p,q);Sline(i,3)=Siz(i);ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');for i=1:n1p=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*((B1(i,5)*(B1(i,5)-1))/B1(i,3))+(conj(E(q))-conj(E(p)))*(B1(i,5)/B1(i,3)));Sjy(i)=Sj(q,p);endSSj(q,p)=Sj(q,p);Sline(i,4)=Sjy(i);ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');for i=1:n1p=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);DDS(i)=DS(i);Sline(i,5)=DS(i);ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];disp(ZF);enddisp('-----------------------------------------------------');disp('各支路首端编号末端编号首端功率末端功率线路损耗');disp(Sline);六、运行结果及其分析是否采用默认数据？(1-默认数据;2-手动输入)1导纳矩阵为：2.9056 -11.5015i 0.0000 + 5.3173i -1.6606 +3.1617i -1.2450 + 2.3710i0.0000 + 5.3173i 0.0000 - 4.6633i 0.0000 + 0.0000i 0.0000 + 0.0000i-1.6606 + 3.1617i 0.0000 + 0.0000i 2.4904 - 4.7039i -0.8298 + 1.5809i-1.2450 + 2.3710i 0.0000 + 0.0000i -0.8298 + 1.5809i 2.0749 - 3.9089i初始雅可比矩阵为：11.1267 2.7603 -5.3173 0 -3.1617 -1.6606-3.0509 11.8762 0 -5.3173 1.6606 -3.1617-5.3173 0 5.3173 0 0 00 -5.3173 0 4.0092 0 0-3.3198 -1.7436 0 0 4.8217 2.69800 0 0 0 0 2.1000迭代次数为：4收敛时电压修正量为：：1.0e-05 *0.0349-0.2445-0.0101-0.5713-0.0931-0.0073各节点的实际电压标幺值E为(节点号从小到大排列)：0.9673 - 0.0655i 1.0252 - 0.1666i 1.0495 - 0.0337i 1.0500 + 0.0000i -----------------------------------------------------各节点的电压大小V为(节点号从小到大排列)：0.9695 1.0387 1.0500 1.0500-----------------------------------------------------各节点的电压相角sida为(节点号从小到大排列)：-3.8734 -9.2315 -1.8419 0各节点的功率S为(节点号从小到大排列)：-0.0000 + 0.0000i -0.5000 - 0.3000i 0.2000 + 0.1969i 0.3277 + 0.0443i -----------------------------------------------------各条支路的首端功率Si为(顺序同您输入B1时一致)：S(1,2)=-0.5-0.30713iS(1,3)=-0.24266-0.197iS(1,4)=-0.25734-0.11013iS(3,4)=-0.055551+0.0017528i-----------------------------------------------------各条支路的末端功率Sj为(顺序同您输入B1时一致)：S(2,1)=0.5+0.24606iS(3,1)=0.25555+0.1952iS(4,1)=0.2712+0.1014iS(4,3)=0.056496-0.057061i-----------------------------------------------------各条支路的功率损耗DS为(顺序同您输入B1时一致)：DS(1,2)=0-0.06107iDS(1,3)=0.012892-0.0018014iDS(1,4)=0.013863-0.0087295iDS(3,4)=0.00094545-0.055308i-----------------------------------------------------各支路首端编号末端编号首端功率末端功率线路损耗1.0000 + 0.0000i2.0000 + 0.0000i -0.5000 - 0.3071i 0.5000 + 0.2461i 0.0000 - 0.0611i 1.0000 + 0.0000i3.0000 + 0.0000i -0.2427 - 0.1970i 0.2556 + 0.1952i 0.0129 - 0.0018i 1.0000 + 0.0000i4.0000 + 0.0000i -0.2573 - 0.1101i 0.2712 + 0.1014i 0.0139 - 0.0087i3.0000 + 0.0000i4.0000 + 0.0000i -0.0556 + 0.0018i 0.0565 - 0.0571i 0.0009 - 0.0553i七、实验体会及感悟通过这次实验，首先让我对matlab软件有了初步的了解，对它强大的矩阵运算能力有了更深的体会，同时掌握了设置断点和断点调试的一般方法，结合课本上的程序流程图和参考资料上的例子单步跟踪调试，再一次的熟悉了牛顿拉夫逊法潮流计算的一般方法和步骤，对计算机计算潮流计算有了更进一步的认识，在学习潮流计算时，虽然依次学习了节点导纳矩阵，功率方程、雅可比矩阵，但不能将它们联系起来，更不知道其中的原委，通过程序的编写，知道了其中的联系，也知道了每个方程、矩阵在计算中的作用。
电力系统分析课程设计——电力系统潮流计算
信息工程学院课程设计报告书题目: 电力系统潮流计算专业：电气工程及其自动化班级： 0310406学号： *********学生姓名：***指导教师：***2013年 4 月 15 日信息工程学院课程设计任务书2013年4月15日目录1 任务提出与方案论证 (4)2 总体设计 (5)2.1潮流计算等值电路 (5)2.2建立电力系统模型 (5)2.3模型的调试与运行 (6)3 详细设计 (6)3.1 计算前提 (6)3.2手工计算 (9)4设计图及源程序 (13)4.1MATLAB仿真 (13)4.2潮流计算源程序 (13)5 总结 (18)参考文献 (19)1 任务提出与方案论证潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法,是电力系统规划和运营中不可缺少的一个重要组成部分。
可以说,它是电力系统分析中最基本、最重要的计算,是系统安全、经济分析和实时控制与调度的基础。
常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。
潮流计算的结果是电力系统稳定计算和故障分析的基础。
在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。
同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。
因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。
在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。
是电力系统研究人员长期研究的一个课题。
它既是对电力系统规划设计和运行方式的合理性、可靠性及经济性进行定量分析的依据，又是电力系统静态和暂态稳定计算的基础。
潮流计算经历了一个由手工到应用数字电子计算机的发展过程，现在的潮流算法都以计算机的应用为前提用计算机进行潮流计算主要步骤在于编制计算机程序，这是一项非常复杂的工作。
电力系统潮流计算课设论文
变电所1变电所2变电所3变电所435kV 母线35kV 母线35kV 母线10kV 母线一次侧电压220kV 一次侧电压220kV母线1母线3母线2线路长为100km线路长为90km线路长为80km线路长为100km线路长为80km2*QFQ-50-2 3*QFs-50-2TQN-100-22*TQN-100-2。
线路长为70km。
线路长为70km电力系统潮流计算课程设计任务书一、题目原始资料1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。
2、发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为350MW ，母线3为机压母线，机压母线上装机容量为100MW ，最大负荷和最小负荷分别为40MW 和25MW ；发电厂二总装机容量为200MW 。
3、变电所资料：(一) 变电所1、2、3、4低压母线的电压等级分别为：35KV 35KV 35KV 10KV (二) 变电所的负荷分别为：（1）50MW 50MW 70MW 60MW （2）60MW 70MW 50MW 45MW （3）50MW 40MW 50MW 70MW电厂一 电厂二（4）70MW 50MW 60MW 60MW （5）70MW 60MW 50MW 50MW （6）45MW 55MW 60MW 70MW （7）70MW 45MW 50MW 60MW （8）50MW 40MW 50MW 60MW （9）40MW 40MW 50MW 60MW （10）50MW 40MW 40MW 60MW (三) 每个变电所的功率因数均为cos φ=0.85；(四) 变电所1和变电所2分别配有两台容量为75MV A 的变压器，短路损耗414KW ，短路电压（%）=16.7；变电所3和变电所4分别配有两台容量为63MV A 的变压器，短路损耗为245KW ，短路电压（%）=10.5；4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。
第3章 简单电力系统的潮流分布计算
3.3.2 对多端网络的处理
1.变电所的运算负荷功率 2.发电厂的运算电源功率
3.3.2 对多端网络的处理
图3-9 多端网络
1.变电所的运算负荷功率
如图3-10a所示，变电所两侧连接着线路，此变 电器低压侧负荷使用的功率S称为负荷功率。ΔS 为变压器的功率损耗，这里将变电器低压侧的 负荷功率加上变压器的功率损耗，称为变电所 的等效负荷功率，用S′表示，S′=S+ΔS。可见 等效负荷功率即是变电所从网络中吸取的功率。
图3-15 例3-1
(示意图)
3.4 环形网中的潮流分布计算
3.4.1 3.4.2 环形网中的初步功率分布 环形网的分解及潮流分布
3.4 环形网中的潮流分布计算
图3-16 简单环形网 a)环形网接线图 b)简化等效电路
3.4.1 环形网中的初步功率分布
1.两端电压相等时的功率分布 2.两端电压不等时的功率分布
3.2 变压器运行状况的分析与计算
3.2.1 3.2.2 变压器的功率损耗和电压降落 变压器的电能损耗
3.2.1 变压器的功率损耗和电压降落
1.变压器的功率损耗 2.变压器的电压降落
3.2.1 变压器的功率损耗和电压降落
图3-7 变压器的等效电路
1.变压器的功率损耗
阻抗支路中的功率损耗ΔSZT为 ΔSZT=S′2U22ZT=P′22+Q′22U22(RT+jXT) =P′22+Q′22U22RT+jP′22+Q′22U22XT =ΔPZT+jΔQZT(3-14) 励磁支路中的功率损耗ΔSyT为 ΔSyT=U21Y∧T=U21(GT+jBT)=U21GT+jU21BT=ΔPy T+jΔQyT(3-15)
电力系统的潮流计算讲解
课程设计（ 2012级本科）题目：给定电力系统的潮流分析学院：物理与机电工程学院专业：电气工程及其自动化作者姓名：李凤朋指导教师：刘永科职称：副教授完成日期： 2015 年 6 月 25 日二〇一五年六月二十五日表5. 负荷数据表节点号有功负荷（MW）无功负荷（MVA）5 125 506 90 308 100 35三、要求：1、完成调试潮流计算与短路计算程序调试；计算程序以《电力系统计算机辅助分析》所用教材中的例程为基础，要求掌握算法原理，建议进行适当改写以利于本次设计应用。
2、根据设计资料准备计算数据；整理原始数据，使其符合计算程序的输入条件。
要求比较采用变压器Π型模型前后两种情况的差异。
3、为指定电网安排合理的稳态运行方式；合理设置目标电网中的控制变量，以1号机为平衡机，将潮流分布调整为可行的运行方式（要求所有节点的电压不得低于1.0 p.u.，也不能高于1.05p.u.）绘制潮流分布图。
分析该运行方式的特点，提出进一步降低网损的改进的措施，并利用计算结果简单证明之。
具体内容参考设计资料所附思考题。
4、列表分析节点7发生各种短路故障时的短路电流，以及运行方式对故障电流的影响；5、完成设计说明书。
要求结构清晰，即从电网模型建立，到初态潮流分布计算，到问题分析，再到调整后的方案计算。
故障电流计算也类似。
三、论文(设计)进度安排阶段论文（设计）各阶段名称起止日期1 熟悉设计任务书、设计题目及设计背景资料 5.20～5.252 查阅有关资料 5.26～5.273 阅读设计要求必读的参考资料 5.28～5.294 书写设计说明书 5.30～6.15摘要潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。
电力系统稳态分析包括潮流计算和静态安全分析。
本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。
电分潮流计算课程设计
电分潮流计算课程设计一、课程目标知识目标：1. 掌握电力系统潮流计算的基本概念、原理和方法；2. 了解电力系统中各元件参数的物理意义及其对潮流计算的影响；3. 学会使用标准算法进行电力系统的潮流计算。
技能目标：1. 能够正确建立电力系统的数学模型，并进行潮流计算；2. 能够运用相关软件工具（如PSS/E、DIgSILENT PowerFactory等）进行潮流计算，并分析计算结果；3. 能够识别并解决潮流计算过程中出现的问题，优化计算结果。
情感态度价值观目标：1. 培养学生严谨的科学态度，注重理论与实践相结合；2. 激发学生对电力系统分析的兴趣，提高学生主动学习的积极性；3. 培养学生团队协作意识，增强沟通与交流能力。
课程性质：本课程为电力系统分析领域的基础课程，旨在帮助学生掌握电力系统潮流计算的基本理论和方法，提高解决实际问题的能力。
学生特点：学生已具备一定的电力系统基础知识，具备一定的数学和物理背景，但可能对实际应用中的复杂问题分析能力较弱。
教学要求：结合学生特点，注重理论与实践相结合，以案例教学为主，培养学生实际操作和解决问题的能力。
通过本课程的学习，使学生能够达到上述课程目标，并为后续相关课程的学习打下坚实基础。
二、教学内容1. 引言：电力系统潮流计算的基本概念、意义及其应用范围；教材章节：第一章2. 潮流计算基本原理：- 交流电力系统的基本方程；- 节点功率方程；- 输电线路的潮流方程；教材章节：第二章3. 潮流计算方法：- 牛顿-拉夫逊法；- 高斯-赛德尔法；- 快速分解法；教材章节：第三章4. 潮流计算中的特殊情况处理：- 不平衡功率处理；- 无解或多解现象；- 非线性方程的求解；教材章节：第四章5. 潮流计算软件应用：- PSS/E、DIgSILENT PowerFactory等软件介绍；- 潮流计算案例分析；- 软件操作与结果分析；教材章节：第五章6. 潮流计算在实际工程中的应用：- 电力系统稳定性分析；- 电压无功优化；- 短路计算；教材章节：第六章本教学内容安排和进度根据课程目标和教学要求制定，确保学生能够系统掌握电力系统潮流计算的理论和方法，注重理论与实践相结合，提高学生的实际操作能力。
电力系统的潮流分布计算项目计划书
电力系统的潮流分布计算计划书第一章电力系统潮流计算简述电力系统的潮流分布计算（牛顿——拉夫逊直角坐标法）第一节潮流计算的意义(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。
(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。
(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。
(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。
总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。
同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。
因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。
在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。
第二节潮流计算的发展史利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。
此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。
对潮流计算的要求可以归纳为下面几点：（1）算法的可靠性或收敛性（2）计算速度和内存占用量（3）计算的方便性和灵活性电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。
因此其数学模型不包含微分方程，是一组高阶非线性方程。
非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。
随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已达到几千阶甚至上万阶，对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。
电力系统课程设计-潮流计算
摘要潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。
对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。
潮流计算是电力系统分析最基本的计算。
除它自身的重要作用之外，潮流计算还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础。
实际电力系统的潮流计算主要采用牛顿-拉夫逊法。
按电压的不同表示方法，牛顿-拉夫逊潮流计算分为直角坐标形式和极坐标形式两种。
本次计算采用直角坐标形式下的牛顿-拉夫逊法，牛顿-拉夫逊法有很好的收敛性，但要求有合适的初值。
传统的潮流计算程序缺乏图形用户界面，结果显示不直接难与其他分析功能集成。
网络原始数据输入工作大量且易于出错。
本文采用MATLAB语言运行WINDOWS操作系统的潮流计算软件。
目前MATLAB已成为国际控制界最流行、使用最广泛的语言了。
它的强大的矩阵处理功能给电力系统的分析、计算带来很多方便，而且采用MATLAB界面直观，运行稳定，计算准确。
所以本次课程设计程序设计采用MATLAB计算。
关键词：电力系统潮流计算牛顿—拉夫逊法潮流计算 MATLAB目录一、概述1.1设计目的与要求 (3)1.1.1 设计目的 (3)1.1.2 设计要求 (3)1.2 设计题目 (3)1.3 设计内容 (3)二电力系统潮流计算概述 (4)2.1 电力系统简介 (4)2.2 潮流计算简介 (4)2.3 潮流计算的意义及其发展..................... . (5)三潮流计算设计题目 (6)3.1 潮流计算题目 (6)3.2 对课题的分析及求解思路 (7)四潮流计算算法及手工计算 (7)4.1 极坐标下P-Q法的算法 (7)4.2 节点电压方程 (8)4.3节点导纳矩阵 (9)4.4 导纳矩阵在潮流计算 (10)4.5 潮流计算的手工计算 (12)五 Matlab概述 (13)5.1 Matlab简介 (14)5.2 Matlab的应用 (14)5.3 矩阵的运算 (14)5.3.1 与常数的运算 (14)5.3.2 基本数学运算 (14)5.3.3 逻辑关系运算 (14)5.4 Matlab中的一些命令 (15)六潮流计算流程图及源程序 (18)6.1 潮流计算流程图 (18)6.2 潮流计算源程序图 (19)6.3 运行计算结果 (27)七总结 (29)八参考文献 (29)第一章系统概述1.1 设计目的与要求1.1.1设计目的1.掌握电力系统潮流计算的基本原理；2.掌握并能熟练运用一门计算机语言（MATLAB语言或C语言或C++语言）；3.采用计算机语言对潮流计算进行计算机编程。
电力系统潮流分布计算
第二章 电力系统潮流分布计算例1：某110kV 输电线，长100km ，r=0.21Ω/km,x=0.409Ω/km,b=62.74/10S km -⨯,线路末端功率10M ，cos 0.95φ=滞后，已知末端电压为110kV ，试计算始端电压大小和相角，始端功率，并作相量图。
例2、某110kV 输电线，长100km ，r=0.21Ω/km ，x=0.409Ω/km ，b=62.74/10S km -⨯,线路末端功率10MW ，cos 0.95φ=超前，已知始端电压为112kV ，试计算末端电压大小和相角、始端功率，并作相量图。
解：线路的等值电路及参数见例1解。
由于cos 0.95φ=超前，故18.195φ=-︒，末端功率为：思考题与习题一、思考题：1、输电线路和变压器的功率损耗如何计算?它们在各导纳支路上损耗的无功功率有什么不同?2、输电线路和变压器阻抗元件上的电压降落如何计算?电压降落的大小主要由什么决定?电压降落的相位主要由什么决定?什么情况下会出现线路末端电压高于首端电压的情况?3、电压降落、电压损耗、电压偏移各是如何定义的?4、运算功率是什么?运算负荷是什么?如何计算升压变电站的运算功率和降压变电站的运算负荷?5、辐射形网络潮流计算可以分为哪两种类型?分别怎样计算?6、简单闭式网络主要有哪两种类型?其潮流计算的主要步骤是什么?7、节点阻抗矩阵如何形成?是否具有稀疏性?其各元素的物理意义是什么?二、习题⨯s/km.1、某110kv输电线路、长80km, r=0.2l/kmΩ, b=2.746Ω，x=0.409/km10-线路未端功率10MW. cos φ端电压为110kV ,试计算始端电压大小和相角、始端功率.并作相量图。
2、一双绕组变压器，型号SFL1-10000,电压35±5%/11kV ,KP= 58.29kW,P=11.75kW,%k u =7.5，0%I =1.5，低压侧负荷10MW. cos φ=0.85,低压侧电压10kV ,变压器抽头电压+5%,求： （1）功率分布; （2）高压侧电压。