两机五节点网络潮流计算—牛拉法
电力系统稳态分析课程设计说明书
电力系统稳态分析课程设计说明书题目：电力系统的潮流分布计算（牛顿——拉夫逊直角坐标法）学生姓名:高峰专业：电气工程及其自动化班级：电气2009——1班内蒙古科技大学课程设计任务书12系统接线图其中节点1为平衡节点，节点2、3、4、5为PQ节点。
前言一内蒙古科技大学课程设计任务书 (2)第一章电力系统潮流分布简述 (6)第一节潮流计算的发展趋势 (6)第二节潮流计算的发展史 (6)第三节潮流计算的意义 (8)第二章设计程序 (8)第一节潮流计算的基本步骤 (8)第二节潮流计算设计资料及参数 (10)第三节P—Q法潮流计算流程图 (10)第四节潮流计算的matlab程序 (11)第五节潮流计算程序运行结果显示 (16)附录感想 (39)参考 (39)第一章电力系统潮流计算简述电力系统的潮流分布计算（牛顿——拉夫逊直角坐标法）第一节潮流计算的意义(1）在电网规划阶段，通过潮流计算,合理规划电源容量及接入点,合理规划网架，选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。
（2）在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。
（3）正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制，指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求.(4）预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。
总结为在电力系统运行方式和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。
同时,为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。
因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。
在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。
关于直角坐标牛拉法系统潮流分布计算的答辩
关于直角坐标牛拉法系统潮流分布计算的答辩直角坐标牛拉法系统潮流分布计算是电力系统分析中的一种重要方法。
在电力系统中，潮流计算是评估系统的稳定性和可靠性的关键步骤之一、本次答辩将介绍直角坐标牛拉法系统潮流分布计算的原理、优势和应用。
直角坐标牛拉法是一种解决复杂电力系统潮流计算问题的有效方法。
在进行潮流计算时，我们需要根据功率平衡方程和节点电压相位差方程来建立数学模型。
直角坐标牛拉法通过将复数形式的电压和注入功率转换为直角坐标形式，大大简化了计算过程。
通过迭代求解非线性方程组，可以得到系统中各节点的电压和相角。
直角坐标牛拉法系统潮流分布计算有以下几个优势：1.精度高：直角坐标牛拉法可以得到较为准确的潮流计算结果。
通过改变迭代次数和误差限制条件，可以实现较高的计算精度。
2.适用范围广：直角坐标牛拉法适用于各种不同类型的电力系统，包括输电网、配电网和混合电力系统等。
无论是小型还是大型系统，均可使用该方法进行潮流计算。
3.可拓展性强：直角坐标牛拉法可以根据需要进行拓展，例如考虑输电线路的参数变化、不同节点的发电机和负荷模型等。
这使得在复杂电力系统的分析中，直角坐标牛拉法能够应对更多的预测和实际运行问题。
直角坐标牛拉法系统潮流分布计算在电力系统的多个方面应用广泛：1.安全评估：通过潮流计算，我们可以对电力系统进行稳定性和可靠性评估。
这对于电力系统的规划、运行和维护非常重要。
2.输电线路调度：潮流计算可以提供输电线路的实时状态和导纳情况，从而帮助进行输电线路的调度安排。
这有助于优化电力系统的运行效率和功率传输能力。
3.发电机出力控制：通过潮流计算，我们可以获得发电机的出力和负载之间的关系。
这对于发电机的运行和控制至关重要，以确保电力系统的平衡和稳定。
4.变电站优化：潮流计算可以帮助确定变电站的最佳参数配置和运行方案。
通过优化变电站的调度计划，可以降低电力损耗和提高系统的供电质量。
总之，直角坐标牛拉法系统潮流分布计算是电力系统分析的重要方法之一、它通过简化计算过程和提高计算精度，为电力系统的规划和运行提供了有力的支持。
简单的潮流计算-本科毕业设计
摘要本文分别对4节点电力网络，5节点电力网络和30节点的电力网络进行简单的电力网络的潮流计算，采用C++软件平台编制潮流程序，分别采用牛拉法和PQ快速解耦法编制了两个程序，通过对比两种算法的计算结果、迭代精度和迭代次数，进而分析得出两种潮流计算机方法的特点。
电力系统的潮流计算是用来监测和研究有关于电力系统稳态运行状况的基本计算，是电力系统规划设计、运行与分析的基本工具。
为提高电力网络系统自动化操作的快速性及准确性，而且让系统可以进行安全、系统地进行实时分析，电力网络的潮流变化信息是一项重要的资源，所以对于电力网络进行潮流计算是现在潮流计算的发展的一个方面。
本文的主要内容就是基于简单电力网络地潮流计算设计，结合牛顿-拉夫逊法潮流计算，采用PQ快速解耦发潮流计算程序，完成基于C语言的电力系统潮流计算程序的设计。
在电力系统在正常运行的情况下，通过对电力系统的潮流（功率、电压、相角等）来对电力系统网络的某一个特定时间点的运行状态进行描述。
本文基于以上的基础知识对于电力系统网络潮流计算进行了设计、仿真，并且分析了在不同场景两种计算方法的结果精确度。
关键词：C++编程；潮流计算；牛顿-拉夫逊算法；PQ快速解耦法；迭代比较ABSTRACTElectric power network, this article respectively to four node 5 node power network with 30 nodes of power network to carry on the simple electric network power flow calculation, using c + + software platform trend program, with rafa and PQ fast decoupling method created two programs, by comparing the calculation results of two algorithms, iterative precision and the number of iterations, then analysis the characteristics of the two methods of trend of the computer。
牛拉法
一、牛顿拉夫逊法进行潮流计算的MATLAB程序%本程序的功能是用牛顿拉夫逊法进行潮流计算%n=input('请输入节点数:n=');%nl=input('请输入支路数:nl=');%isb=input('请输入平衡母线节点号:isb=');%pr=input('请输入误差精度:pr=');%B1=input('请输入由各支路参数形成的矩阵:B1=');%B2=input('请输入各节点参数形成的矩阵:B2=');n=5; %输入节点数nl=5; %输入支路数isb=1; %输入平衡母线节点号pr=0.00001; %输入误差精度B1=[1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1];%输入由各支路参数形成的矩阵B2=[0 0 1.05 1.05 0 1;0 3.7+1.3i 1.05 0 0 2;0 2+1i 1.05 0 0 2;0 1.6+0.8i 1.05 0 0 2;5 0 1.05 1.05 0 3];%输入各节点参数形成的矩阵Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);for i=1:nl %形成节点导纳矩阵if B1(i,6)==0 %判断变压器变比是否在低压侧p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %节点互导纳Y(q,p)=Y(p,q); %节点互导纳Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%节点导纳=自导纳+X修正量Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;end%求导纳矩阵disp('导纳矩阵Y=');disp(Y);G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部for i=1:n %给定各节点初始电压的实部和虚部e(i)=real(B2(i,3)); %取出实部f(i)=imag(B2(i,3)); %取出虚部V(i)=B2(i,4); %PV节点电压给定模值endfor i=1:n %给定各节点注入功率S(i)=B2(i,1)-B2(i,2);%i节点注入功率SG-SLB(i,i)=B(i,i)+B2(i,5);%i节点无功补偿量endP=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2~=0IT2=0;a=a+1;for i=1:nif i~=isbC(i)=0;D(i)=0;for j1=1:nC(i)= C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);D(i)= D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); endP1=C(i)*e(i)+f(i)*D(i);Q1=f(i)*C(i)-D(i)*e(i);V2=e(i)^2+f(i)^2;if B2(i,6)~=3DP=P(i)-P1;DQ=Q(i)-Q1;for j1=1:nif j1~=isb&j1~=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);X2=B(i,j1)*e(i)-G(i,j1)*f(i);X3=X2;X4=-X1;p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; elseif j1==i&j1~=isbX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; endendelseDP=P(i)-P1;DV=V(i)^2-V2;for j1=1:nif j1~=isb&j1~=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);X2=B(i,j1)*e(i)-G(i,j1)*f(i);X5=0;X6=0;p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2; elseif j1==i&j1~=isbX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X5=-2*e(i);X6=-2*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;endendendendend%求雅可比矩阵for k=3:N0k1=k+1;N1=N;for k2=k1:N1J(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=1;if k~=3;k4=k-1;for k3=3:k4for k2=k1:N1J(k3,k2)= J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endif k==N0,break;endfor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefor k3=k1:N0for k2=k1:N1J(k3,k2)= J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendend %119页zhengtong打%for k=3:2:N0-1L=(k+1)./2;e(L)=e(L)-J(k,N);k1=k+1;f(L)=f(L)-J(k1,N);endfor k=3:N0DET=abs(J(k,N));if DET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;for k=1:ndy(k)=sqrt(e(k)^2+f(k)^2);endfor i=1:nDy(ICT1,i)=dy(i);endend%用高斯消去法解“w=-J*V”disp('迭代次数');disp(ICT1);disp('没有达到精度要求的个数');disp(ICT2);for k=1:nE(k)=sqrt(e(k)^2+f(k)^2);O(k)=atan(f(k)./e(k))*180./pi;endE=e+f*j;disp('各节点的实际电压标么值E为(节点号从小到大排列):'); disp(E);disp('各节点的电压大小V为(节点号从小到大排列):');disp(V);disp('各节点的电压相角时O为(节点号从小到大排列):'); disp(O);for p=1:nC(p)=0;for q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);end %120页打disp('各节点的功率S为(节点号从小到大排列):');disp(S);disp('各条支路的首端功率Si为(顺序同您输入B1时一样):');for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSi(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))));disp(Si(p,q));enddisp ('各条支路的末端功率Sj为(顺序同您输入B1时一样):');for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSj(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))));disp(Sj(q,p));enddisp('各条支路的功率损耗DS为(顺序同您输入B1时一样):' );for i=1:nlif B1(i,6)==0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i));endfor i=1:ICT1Cs(i)=i;enddisp('以下是每次迭代后各节点的电压值（如图所示）');plot(Cs,Dy),xlabel('迭代次数'),ylabel('电压'),title('电压迭代次数曲线');二、程序输出结果导纳矩阵Y=Columns 1 through 40 -33.3333i 0 +31.7460i 0 00 +31.7460i 1.5846 -35.7379i -0.8299 + 3.1120i -0.7547 + 2.6415i0 -0.8299 + 3.1120i 1.4539 -66.9808i -0.6240 + 3.9002i0 -0.7547 + 2.6415i -0.6240 + 3.9002i 1.3787 - 6.2917i0 0 0 +63.4921i 0 Column 50 +63.4921i0 -66.6667i迭代次数5没有达到精度要求的个数7 8 8 7 0各节点的实际电压标么值E为(节点号从小到大排列):Columns 1 through 41.0500 1.0335 - 0.0774i 1.0260 + 0.3305i 0.8592 - 0.0718i Column 50.9746 + 0.3907i各节点的电压大小V为(节点号从小到大排列):1.0500 0 0 0 1.0500各节点的电压相角时O为(节点号从小到大排列):0 -4.2819 17.8535 -4.7785 21.8433各节点的功率S为(节点号从小到大排列):Columns 1 through 42.5794 + 2.2994i -3.7000 - 1.3000i -2.0000 - 1.0000i -1.6000 - 0.8000i Column 55.0000 + 1.8131i各条支路的首端功率Si为(顺序同您输入B1时一样):2.5794 + 2.2994i-1.2774 + 0.2032i0.1568 + 0.4713i1.5845 + 0.6726i5.0000 + 1.8131i各条支路的末端功率Sj为(顺序同您输入B1时一样): -2.5794 - 1.9745i1.4155 - 0.2443i-0.1338 - 0.3909i-1.4662 - 0.4091i-5.0000 - 1.4282i各条支路的功率损耗DS为(顺序同您输入B1时一样):0 + 0.3249i0.1381 - 0.0412i0.0230 + 0.0804i0.1184 + 0.2635i-0.0000 + 0.3849i以下是每次迭代后各节点的电压值（如图所示）三、程序说明1、牛拉法进行潮流计算的计算机方法电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。
第7章潮流计算的数学模型及基本解法
代入(7-10)式，经整理后有
Vn D I n Ys Vs L Vn U Vn 1
2013-1-10
（7-11）
12
考虑到电流和功率的关系式，上式写成迭代格式为
（ k 1）
Vi
i -1 n （ k） （ k） 1 Si Ys Vs Yij V j Yij V j i 1,2, ,n Yii j1 ji 1 k Vi
2013-1-10
7
综上所述，若选第N个节点为平衡节点，剩下n个节点（n=N-1） 中有r个节点是PV节点，则有n-r个节点是PQ节点。因此除了平衡 节点外，有n个节点注入有功功率，n-r个节点注入无功功率以及r 个节点的电压幅值是已知量。 在直角坐标系，待求的状态变量共2n个，用
x e
T
f
T T
2013-1-10 13
0
（7-12）
考所以，用 V j 代替 V j 可出得到更好的收敛果。 这就是高斯，赛德尔(Gauss-Seidel)选代的基本思 想，即一旦求出电压新值，在随后的迭代中立即使 用。这种方法的选代格式是
（ k 1）
（ k 1）
（ k）
Vi
（ k） （ k） i -1 n 1 Si Ys Vs Yij V j Yij V j i 1,2, ,n Yii j1 ji 1 k Vi
（7-13）
考高斯一赛德尔法比高斯迭代法收敛性要好。 考在导纳矩阵法的迭代公式中，导纳矩阵高度稀 疏，每行只有少数几个是非零元素，非对角非零元 素个数与和节点j相联的支路数相等。所以，上一次
电力系统稳态分析第四章1
节点导纳矩阵的修改（5）
原网络节点i、j之间的导纳由
• 相当于先切除导纳为yij的支路，再增加一条导纳为 • 修改i、j节点的自导纳的增量 • 修改i、j节点间互导纳的增量 的支路
节点导纳矩阵的修改（6）
原网络节点i、j之间变压器的变比由k 原网络 k'
新网络 切除变比为k的变压器 增加变比为k'的变压器
n+1
节点导纳矩阵的修改（3）
在原网络节点i、j之间增加一条支路
• 节点导纳矩阵阶数不变 • 修改i、j节点的自导纳 • 修改i、j节点间互导纳
节点导纳矩阵的修改（4）
在原网络节点i、j之间切除一条导纳为yij的支路
• 相当于增加一条(-yij)的支路 • 节点导纳矩阵阶数不变 • 修改i、j节点的自导纳 • 修改i、j节点间互导纳
节点电压方程
设图中各量均为 标么值，以下不 再区分单相/三相 电路。
节点电压方程
节点注入电流 相量列向量
节点导 纳矩阵
节点电压相 量列向量
N节点电压方程
N节点网络的节点电压方程
n×n阶节点 导纳矩阵
节点注入电流 列向量
节点电压 列向量
• n是网络中的独立（非地）节点数，大地节点编号为0 • 自导纳Yii (YB的对角元) 与节点i直接相连的各支路导纳之和 • 互导纳Yij (YB的非对角元, i≠j) 直接连接于i、j节点之间的导纳的相反数
PV节点
P、 U
Q 、δ
平衡节点
U、δ
P、Q
设置平衡节点的必要性
系统的功率损耗在潮流计算完成之前是未知的， 即功率损 耗是状态变量的函数。必须设置至少一个节点来平衡全网 的功率。 功率方程中节点电压相位是以相对相位（相位差） 的形式 出现的，要求节点电压的绝对相位，必须有一个相位参考 节点。
牛顿-拉夫逊迭代法电力网潮流计算方法与程序
牛顿-拉夫逊迭代法电力网潮流计算方法与程序编写 佘名寰牛顿-拉夫逊迭代法在电力网潮流计算中因其收敛性较好获得广泛运用，该算法的难点是需反复计算功率方程中雅可比矩阵各个元素表达式。
本文简叙了牛顿-拉夫逊迭代法潮流计算的基本公式，通过例题介绍了用牛顿-拉夫逊法计算电力网潮流电压的MATLAB 程序。
程序采用MATLAB 语言的符号矩阵简化了雅可比矩阵系数的计算。
本文可供电力系统电气技术人员和大专院校电力类专业师生参考。
2．牛顿-拉夫逊迭代法潮流计算2.1 极坐标型式节点功率方程式由电源注入或从负载流出节点的电流统称节点电流，由节点电压和节点电流所求得的节点注入功率为：S ti =V i I ̂iI i =∑Y ij V j (I,j=1, 2, ….n)P ti =V i ∑V j n 1(G ij cos θij +B ij sin θij )Q ti =V i ∑V j n 1(G ij sin θij −B ij cos θij ) (i,j=1, 2…..n) （2-01）式中 P ti, Q tI ——节点注入的有功和无功功率Vi ，V j ——节点 i, j 电压幅值θij =θi -θj 节点 i, j 电压的相角差G ij , B ij 节点导纳矩阵的元素，Y IJ =G IJ +jB IJ节点功率平衡关系为：P gi-P lI= P tiQ gi-Q li= Q tiΔP i=P gi-P lI-P ti=0ΔQ i=Q gi-Q li-Q ti=0 （2-02）P gi, Q gi——节点i发电机输入有功和无功功率P li , Q li——节点i负荷有功和无功功率ΔP i，ΔQ i--节点i不平衡功率不平衡功率的微分d(ΔP i), d(ΔQ i)d(∆pi )=−(∂p ti∂v1∆v1+ ∂p ti∂v2∆v2…+∂p ti∂ϑ1∆ϑ1+∂p ti∂ϑ2∆ϑ2….)（i=1,2,…,n）d(∆qi )=−(∂q ti∂v1∆v1+ ∂q∂v2∆v2…+∂q ti∂ϑ1∆ϑ1+∂q ti∂ϑ2∆ϑ2….)（2-03）对于n个节点系统可得如下矩阵形式修正方程式；[∆P1∆P2:∆P n∆Q1∆Q2:∆Q n]=[∂∆P1∂ϑ1∂∆P2∂ϑ1:∂∆P n∂ϑ1∂∆Q1∂ϑ1∂∆Q2∂ϑ1:∂∆Q n∂ϑ1∂∆P1∂ϑ2∂∆P2∂ϑ2:∂∆P n∂ϑ2∂∆Q1∂ϑ2∂∆Q2∂ϑ2:∂∆Q n∂ϑ2⋯⋯⋯⋯⋯⋯⋯⋯∂∆P1∂ϑn∂∆P2∂ϑn:∂∆P n∂ϑn∂∆Q1∂ϑn∂∆Q2∂ϑn:∂∆Q n∂ϑn∂∆P1∂v1∂∆P2∂v1:∂∆P n∂v1∂∆Q1∂v1∂∆Q2∂v1:∂∆Q n∂v1∂∆P1∂v2∂∆P2∂v2:∂∆P n∂v2∂∆Q1∂v2∂∆Q2∂v2:∂∆Q n∂v2⋯⋯⋯⋯⋯⋯⋯⋯∂∆P1∂v n∂∆P2∂v n:∂∆P n∂v n∂∆Q1∂v n∂∆Q2∂v n:∂∆Q n∂v n][Δϑ1∆ϑ2:∆ϑn∆v1∆v2:∆v n]（2-04）式中偏微分矩阵为雅可比矩阵。
电力系统牛拉法潮流计算概要PPT文档43页
6、最大的骄傲于最大的自卑都表示心灵的最软弱无力。——斯宾诺莎 7、自知之明是最难得的知识。——西班牙 8、勇气通往天堂，怯懦通往地狱。——塞内加 9、有时候读书是一种巧妙地避开思考的方法。——赫尔普斯 10、阅读一切好书如同和过去最杰出的人谈话。——统牛拉法潮流计算概要
51、没有哪个社会可以制订一部永远 适用的 宪法， 甚至一 条永远 适用的 法律。 ——杰 斐逊 52、法律源于人的自卫本能。——英 格索尔
53、人们通常会发现，法律就是这样 一种的 网，触 犯法律 的人， 小的可 以穿网 而过， 大的可 以破网 而出， 只有中 等的才 会坠入 网中。 ——申 斯通 54、法律就是法律它是一座雄伟的大 夏，庇 护着我 们大家 ；它的 每一块 砖石都 垒在另 一块砖 石上。 ——高 尔斯华 绥 55、今天的法律未必明天仍是法律。 ——罗·伯顿
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
两机五节点网络潮流计算—牛拉法
基于Matpower2机5节点系统的潮流计算方法
引言
Matpower是基于MatlabM文件的组建包，主要用来解决电力潮流和优化潮流的问题。
Matpower的设计理念是尽可能简单易懂，它可以执行电力常规潮流运算，如牛顿拉夫逊法，P-Q分解法等，也可以执行最优潮流程序。
本文主要对执行常规的潮流计算进行分析。
1Matpower的简介
内蒙古科技大学课程设计说明书
是导线则为0，如果支路元件为变压器，则该变比为fbu侧母线的基准电压与tbu侧母线的基准电压之比；angle用来设置支路的相位角度，如果支路元件为变压器，其值就是变压器的转角，如果支路元件不是变压器，则相位角度为0°。
可以处理高价机组停机的直流OPF；runopf.m是运行最优潮流计算程序，可运行一个潮流计算程序；runuopf.m运行可以处理高价机组停机的OPF。
本文应用的是runpf.m，即将编写好的程序存为caei(i是自己设定的序号)，并存在Matpower的文件夹下，接着在Matlab的命令窗口输入runpf(‘caei’)，然后回车即可得到结果。
2潮流计算
本文以两机五节点系统为例来对其电网潮流情况进行分析。
系统接线图如下：
2
内蒙古科技大学课程设计说明书
G0.45+j0.150.4+j0.0510.08+j0.2430.01+j0.0340.02+j0.0620.04+j0.1 2G-(0.2+j0.2)50.6+j0.10.08+j0.24
j0.18+60.0.180.06+j0系统接线图
（其中节点1为平衡节点，节点2、3、4、5为PQ节点。
）
该系统的基准容量是100MVA，有5个节点数，每个节点就是一条母线，母线下有各供电线路，支路数为3条。
电源点的基准电压设为
1.05(p.u.)，电源点有功容量为500MW，电源点无功容量为0Mvar。
其已知的节点参数、发电机参数和支路参数如图1所示。
将已知参数输入Matpower软件系统，通过运行分析，所得到的各节点和支路的计算参数如图2所示。
图一
3
内蒙古科技大学课程设计说明书
图二（1）
4
内蒙古科技大学课程设计说明书
5。