文档库 最新最全的文档下载
当前位置:文档库 › 5节点电力系统牛顿-拉夫逊法潮流计算

5节点电力系统牛顿-拉夫逊法潮流计算

5节点电力系统牛顿-拉夫逊法潮流计算
5节点电力系统牛顿-拉夫逊法潮流计算

二 〇 一 四 年 十 二 月

课 程 论 文 学校代码: 10128 学 号: 20141100304

题 目:五节点系统计算机潮流计算编程 学生姓名:张佳羽

学 院:电力学院

系 别:电力系

专 业:电力系统及其自动化

指导教师:郭力萍

程序设计

% 本程序的功能是用牛顿拉夫逊法进行潮流计算

n=input('请输入节点数:n=');

nl=input('请输入支路数:nl=');

isb=input('请输入平衡母线节点号:isb=');

pr=input('请输入误差精度:pr=');

B1=input('请输入由各支路参数形成的矩阵:B1=');

B2=input('请输入各节点参数形成的矩阵:B2=');

X=input('请输入由节点号及其对地阻抗形成的矩阵:X='); 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);

end

Y(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;

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);

end

for i=1:n

S(i)=B2(i,1)-B2(i,2);

B(i,i)=B(i,i)+B2(i,5);

end

P=real(S);Q=imag(S);

ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;

while IT2~=0

IT2=0;a=a+1;

for i=1:n

if i~=isb

C(i)=0;

D(i)=0;

for j1=1:n

C(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);

end

P1=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)~=3

DP=P(i)-P1;

DQ=Q(i)-Q1;

for j1=1:n

if j1~=isb&j1~=i

X1=-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~=isb

X1=-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;

end

end

else

DP=P(i)-P1;

DV=V(i)^2-V2;

for j1=1:n

if j1~=isb&j1~=i

X1=-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~=isb

X1=-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;

end

end

end

end

end

%求雅可比矩阵

for k=3:N0

k1=k+1;N1=N;

for k2=k1:N1

J(k,k2)=J(k,k2)./J(k,k);

end

J(k,k)=1;

if k~=3;

k4=k-1;

for k3=3:k4

for k2=k1:N1

J(k3,k2)= J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

if k==N0,break;end

for k3=k1:N0

for k2=k1:N1

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

else

for k3=k1:N0

for k2=k1:N1

J(k3,k2)= J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

end

end

for k=3:2:N0-1

L=(k+1)./2;

e(L)=e(L)-J(k,N);

k1=k+1;

f(L)=f(L)-J(k1,N);

end

for k=3:N0

DET=abs(J(k,N));

if DET>=pr

IT2=IT2+1;

end

end

ICT2(a)=IT2;

ICT1=ICT1+1;

for k=1:n

dy(k)=sqrt(e(k)^2+f(k)^2);

end

for i=1:n

Dy(ICT1,i)=dy(i);

end

end

%用高斯消去法解“w=-J*V”disp('迭代次数');

disp(ICT1);

disp('没有达到精度要求的个数'); disp(ICT2);

for k=1:n

V(k)=sqrt(e(k)^2+f(k)^2);

shita(k)=atan(f(k)./e(k))*180/pi;

E(k)=e(k)+f(k)*j;

end

disp('各节点的实际电压标么值E为(节点号从小到大排列):'); disp(E);

disp('各节点的电压大小V为(节点号从小到大排列):');

disp(V);

disp('各节点的电压相角时shita为(节点号从小到大排列):'); disp(shita);

for p=1:n

C(p)=0;

for q=1:n

C(p)=C(p)+conj(Y(p,q))*conj(E(q));

end

S(p)=E(p)*C(p);

end

disp('各节点的功率S为(节点号从小到大排列):');

disp(S);

disp('各条支路的首端功率Si为(顺序同您输入B1时一样):'); 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);

end

Si(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));

end

disp ('各条支路的末端功率Sj为(顺序同您输入B1时一样):');

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);

end

Sj(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));

end

disp('各条支路的功率损耗DS为(顺序同您输入B1时一样):' );

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);

end

DS(i)=Si(p,q)+Sj(q,p);

disp(DS(i));

end

for i=1:ICT1

Cs(i)=i;

end

disp('以下是每次迭代后各节点的电压值(如图所示)');

plot(Cs,Dy),xlabel('迭代次数'),ylabel('电压'),title('电压迭代次数曲线');

运行结果如下所示:

请输入节点数: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 0 0 2;0 2+1i 1 0 0 2;0 1.6+0.8i 1 0 0 2;5 0 1.05 1.05 0 3]

请输入由节点号及其对地阻抗形成的矩阵:X=[1 0;2 0;3 0;4 0;5 0]

导纳矩阵Y

0 -33.3333i 0 +31.7460i 0 0 0

0 +31.7460i 1.5846 -35.7379i -0.8299 + 3.1120i -0.7547 + 2.6415i 0

0 -0.8299 + 3.1120i 1.4539 -66.9808i -0.6240 + 3.9002i 0 +63.4921i

0 -0.7547 + 2.6415i -0.6240 + 3.9002i 1.3787 - 6.2917i 0

0 0 0 +63.4921i 0 0 -66.6667i

迭代次数

5

没有达到精度要求的个数

7 8 8 6 0

各节点的实际电压标么值E为(节点号从小到大排列):

1.0500 1.0335 - 0.0774i 1.0260 + 0.3305i 0.8592 - 0.0718i 0.9746 + 0.3907i

各节点的电压大小V为(节点号从小到大排列):

1.0500 1.0364 1.0779 0.8622 1.0500

各节点的电压相角时shita为(节点号从小到大排列):

0 -4.2819 17.8535 -4.7785 21.8433

各节点的功率S为(节点号从小到大排列):

2.5794 + 2.2994i -

3.7000 - 1.3000i -2.0000 - 1.0000i -1.6000 - 0.8000i 5.0000 + 1.8131i

各条支路的首端功率Si为(顺序同您输入B1时一样):

2.5794 + 2.2994i

-1.2774 + 0.2032i

0.1568 + 0.4713i

1.5845 + 0.6726i

5.0000 + 1.8131i

各条支路的末端功率Sj为(顺序同您输入B1时一样): -2.5794 - 1.9745i

1.4155 - 0.2443i

-0.1338 - 0.3909i

-1.4662 - 0.4091i

-5.0000 - 1.4282i

各条支路的功率损耗DS为(顺序同您输入B1时一样): -0.0000 + 0.3249i

0.1381 - 0.0412i

0.0230 + 0.0804i

0.1184 + 0.2635i

0.0000 + 0.3849i

本程序符号说明

(1)clc,clear清屏,清除以前显示结果;

(2)disp 显示待输出的内容,如disp('节点总数为:'),显示:节点总数为:;

(3)for,end;循环语句,用于要循环的地方;

(4)/ 正除;\ 反除;

(5)sum 求和;

(6)conj 复数求共轭;

(7)% 注释

(8):自然数循环如1:5

以下是每次迭代后各节点的电压值(如图1所示):

图1 每次迭代后各节点的电压值

牛顿拉夫逊法潮流计算

摘要 本文,首先简单介绍了基于在MALAB中行潮流计算的原理、意义,然后用具体的实例,简单介绍了如何利用MALAB去进行电力系统中的潮流计算。 众所周知,电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各线的电压、各元件中流过的功率、系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。 此外,在进行电力系统静态及暂态稳定计算时,要利用潮流计算的结果作为其计算的基础;一些故障分析以及优化计算也需要有相应的潮流计算作配合;潮流计算往往成为上述计算程序的一个重要组成部分。以上这些,主要是在系统规划设计及运行方式安排中的应用,属于离线计算范畴。 牛顿-拉夫逊法在电力系统潮流计算的常用算法之一,它收敛性好,迭代次数少。本文介绍了电力系统潮流计算机辅助分析的基本知识及潮流计算牛顿-拉夫逊法,最后介绍了利用MTALAB程序运行的结果。 关键词:电力系统潮流计算,牛顿-拉夫逊法,MATLAB

ABSTRACT This article first introduces the flow calculation based on the principle of MALAB Bank of China, meaning, and then use specific examples, a brief introduction, how to use MALAB to the flow calculation in power systems. As we all know, is the study of power flow calculation of power system steady-state operation of a calculation, which according to the given operating conditions and system wiring the entire power system to determine the operational status of each part: the bus voltage flowing through the components power, system power loss and so on. In power system planning power system design and operation mode of the current study, are required to quantitatively calculated using the trend analysis and comparison of the program or run mode power supply reasonable, reliability and economy. In addition, during the power system static and transient stability calculation, the results of calculation to take advantage of the trend as its basis of calculation; number of fault analysis and optimization also requires a corresponding flow calculation for cooperation; power flow calculation program often become the an important part. These, mainly in the way of system design and operation arrangements in the application areas are off-line calculation. Newton - Raphson power flow calculation in power system is one commonly used method, it is good convergence of the iteration number of small, introduce the trend of computer-aided power system analysis of the basic knowledge and power flow Newton - Raphson method, introduced by the last matlab run results. Keywords:power system flow calculation, Newton – Raphson method, matlab

基于极坐标的牛顿拉夫逊潮流计算

前言 电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行件及系统接线情况确定整个电力系统各部分的运行状态。在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量分析、比较供电方案或运行方式的合理性、可靠性和经济性。本次课程设计任务是闭式网络的潮流计算,用到的方法为牛顿拉夫逊极坐标法潮流计算。 牛顿法是数学中解决非线性方程式的典型方法,有较好的收敛性。解决电力系统潮流计算问题是以导纳距阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数距阵的稀疏性,就可以大大提高牛顿法潮流程序的放率。自从20 世纪60 年代中期利用了最佳顺序消去法以后,牛顿法在收敛性、内存要求、速度方面都超过了阻抗法,成为直到目前仍在广泛采用的优秀方法。

目录 1任务书 (2) 2.模型简介及等值电路 (3) 3.设计原理 (4) 4.修正方程的建立 (7) 5.程序流程图及MATLAB程序编写 (9) 6.结果分析 (15) 7.设计总结 (25) 8.参考文献 (26)

《电力系统分析》 课程设计任务书 题目极坐标表示的牛顿拉夫逊法潮流计算程序设计学生姓名学号专业班级

设计内容与要求1. 设计要求 掌握MATLAB语言编程方法;理解和掌握运用计算机进行潮流计算的基本算法原理;针对某一具体电网,进行潮流计算程序设计。 其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解,培养学生运用所学知识分析和解决问题的能力。 2. 内容 1)学习并掌握MATLAB语言。 2)掌握变压器非标准变比概念及非标准变比变压器的等值电路。掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。 3)掌握电力系统功率方程、变量和节点分类。 4)掌握利用极坐标表示的牛-拉法进行潮流计算的方法和步骤。 5)选择一个某一具体电网,编制程序流程框图。 6)利用MATLAB语言编写该模型的潮流计算程序,并上机调试程序,对计算结果进行分析。 7)整理课程设计论文。 起止时间2013 年7 月 4 日至2013 年7月10日指导教师签名年月日 系(教研室)主任签 名 年月日学生签名年月日 2 模型简介及等值电路 2.1课程设计模型:模型3

电力系统分析课程设计-潮流计算

目录 摘要 (1) 1.任务及题目要求 (2) 2.计算原理 (3) 2.1牛顿—拉夫逊法简介 (3) 2.2牛顿—拉夫逊法的几何意义 (7) 3计算步骤 (7) 4.结果分析 (9) 小结 (11) 参考文献 (12) 附录:源程序 (13) 本科生课程设计成绩评定表 (32)

摘要 电力系统的出现,使高效,无污染,使用方便,易于调控的电能得到广泛应用,推动了社会生产各个领域的变化,开创了电力时代,发生率第二次技术革命。电力系统的规模和技术水准已经成为一个国家经济发展水平的标志之一。 电力系统稳态分析包括潮流计算和静态安全分析。电力系统潮流计算是电力系统最基本的计算,也是最重要的计算。所谓潮流计算,就是已知电网的接线方式与参数及运行条件,计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统,通过潮流计算可以判断电网母线电压、支路电流和功率是否越限,如果有越限,就应采取措施,调整运行方式。对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 在电力系统运行方式和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中,则采用在线潮流计算。 关键词:电力系统潮流计算牛顿-拉夫逊法

用matlab电力系统潮流计算

题目:潮流计算与matlab 教学单位电气信息学院姓名 学号 年级 专业电气工程及其自动化指导教师 职称副教授

摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算,潮流计算是电力网络设计与运行中最基本的运算,对电力网络的各种设计方案及各种运行方式进行潮流计算,可以得到各种电网各节点的电压,并求得网络的潮流及网络中的各元件的电力损耗,进而求得电能损耗。本位就是运用潮流计算具体分析,并有MATLAB仿真。 关键词:电力系统潮流计算 MATLAB Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation

基于MATLAB牛顿拉夫逊法进行潮流计算

>> %本程序的功能是用牛顿拉夫逊法进行潮流计算n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); isb=input('请输入平衡母线节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); 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); end Y(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; 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); end for i=1:n S(i)=B2(i,1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5); end P=real(S);Q=imag(S); ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; while IT2~=0 IT2=0;a=a+1; for i=1:n if i~=isb C(i)=0; D(i)=0; for j1=1:n C(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); end

电力系统课程设计-牛顿拉夫逊法潮流计算

课程设计说明书 题目电力系统分析系(部) 专业(班级) 学号 指导教师 起止日期

电力系统分析课程设计任务书系(部):专业:指导教师:

目录 一、潮流计算基本原理 1.1潮流方程的基本模型 1.2潮流方程的讨论和节点类型的划分 1.3、潮流计算的意义 二、牛顿-拉夫逊法 2.1牛顿-拉夫逊法基本原理 2.2节点功率方程 2.3修正方程 2.4牛顿法潮流计算主要流程 三、收敛性分析 四、算例分析 总结 参考文献

电力系统分析潮流计算 一、潮流计算基本原理 1.1潮流方程的基本模型 电力系统是由发电机、变压器、输电线路及负荷等组成,其中发电机及负荷是非线性元件,但在进行潮流计算时,一般可以用接在相应节点上的一个电流注入量来代表。因此潮流计算所用的电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成,并用集中参数表示的串联或并联等值支路来模拟。结合电力系统的特点,对这样的线性网络进行分析,普通采用的是节点法,节点电压与节点电流之间的关系 V Y I = (1-1) 其展开式为 j n j ij i V Y I ∑==1 ),,3,2,1 (n i = (1-2) 在工程实际中,已经的节点注入量往往不是节点电流而是节点功率,为此必须应用联 系节点电流和节点功率的关系式 i i i i V jQ P I * -= ),,3,2,1(n i = (1-3) 将式(1-3)代入式(1-2)得到 j n j ij i i i V Y V jQ P ∑=* =-1 ),,3,2,1(n i = (1-4) 交流电力系统中的复数电压变量可以用两种极坐标来表示 i j i i e V V θ= (1-5) 或 i i i jf e V += (1-6) 而复数导纳为

电力系统潮流计算

电力系统潮流计算 Company number:【0089WT-8898YT-W8CCB-BUUT-202108】

电力系统 课程设计题目: 电力系统潮流计算 院系名称:电气工程学院 专业班级:电气F1206班 学生姓名: 学号: 指导教师:张孝远 1 2 节点的分类 (5) 3 计算方法简介 (6) 牛顿—拉夫逊法原理 (6) 牛顿—拉夫逊法概要 (6) 牛顿法的框图及求解过程 (8) MATLAB简介 (9) 4 潮流分布计算 (10)

系统的一次接线图 (10) 参数计算 (10) 丰大及枯大下地潮流分布情况 (14) 该地区变压器的有功潮流分布数据 (15) 重、过载负荷元件统计表 (17) 5 设计心得 (17) 参考文献 (18) 附录:程序 (19) 原始资料 一、系统接线图见附件1。 二、系统中包含发电厂、变电站、及其间的联络线路。500kV变电站以外的系统以一个等值发电机代替。各元件的参数见附件2。 设计任务 1、手动画出该系统的电气一次接线图,建立实际网络和模拟网络之间的联系。 2、根据已有资料,先手算出各元件的参数,后再用Matlab表格核算出各元件的参数。 3、潮流计算 1)对两种不同运行方式进行潮流计算,注意110kV电网开环运行。 2)注意将电压调整到合理的范围 110kV母线电压控制在106kV~117kV之间; 220kV母线电压控制在220 kV~242kV之间。 附件一:

72 水电站2 水电站1 30 3x40 C 20+8 B 2x8 A 2x31.5 D 4x7.5 水电站5 E 2x10 90+120 H 12.5+31.5 F G 1x31.5 水电站3 24 L 2x150 火电厂 1x50 M 110kV线路220kV线路课程设计地理接线示意图 110kV变电站220kV变电站牵引站火电厂水电站500kV变电站

电力系统分析潮流计算

电力系统分析潮流计算报告

目录 一.配电网概述 (3) 1.1 配电网的分类 (3) 1.2 配电网运行的特点及要求 (3) 1.3 配电网潮流计算的意义 (4) 二.计算原理及计算流程 (4) 2.1 前推回代法计算原理 (4) 2.2 前推回代法计算流程 (7) 2.3主程序清单: (9) 2.4 输入文件清单: (11) 2.5计算结果清单: (12) 三.前推回代法计算流程图 (13) 参考文献 (14)

一.配电网概述 1.1 配电网的分类 在电力网中重要起分配电能作用的网络就称为配电网; 配电网按电压等级来分类,可分为高压配电网(35—110KV),中压配电网(6—10KV,苏州有20KV的),低压配电网(220/380V); 在负载率较大的特大型城市,220KV电网也有配电功能。 按供电区的功能来分类,可分为城市配电网,农村配电网和工厂配电网等。 在城市电网系统中,主网是指110KV及其以上电压等级的电网,主要起连接区域高压(220KV及以上)电网的作用。 配电网是指35KV及其以下电压等级的电网,作用是给城市里各个配电站和各类用电负荷供给电源。 从投资角度看,我国与国外先进国家的发电、输电、配电投资比率差异很大,国外基本上是电网投资大于电厂投资,输电投资小于配电投资。我国刚从重发电轻供电状态中转变过来,而在供电投资中,输电投资大于配电投资。从我国城网改造之后,将逐渐从输电投资转入配电建设为主。 本文是基于前推回代法的配电网潮流分析计算的研究,研究是是以根节点为10kV的电压等级的配电网。 1.2 配电网运行的特点及要求 配电系统相对于输电系统来说,由于电压等级低、供电范围小,但与用户直接相连,是供电部门对用户服务的窗口,因而决定了配电网运行有如下特点和基本要求:

电力系统潮流计算方法分析

电力系统潮流分析 —基于牛拉法和保留非线性的随机潮流 , 姓名:*** 学号:***

1 潮流算法简介 常规潮流计算 常规的潮流计算是在确定的状态下。即:通过已知运行条件(比如节点功率或网络结构等)得到系统的运行状态(比如所有节点的电压值与相角、所有支路上的功率分布和损耗等)。 常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法。当初始值和方程的精确解足够接近时,该方法可以在很短时间内收敛。下面简要介绍该方法。 牛顿拉夫逊方法原理 对于非线性代数方程组式(1-1),在待求量x 初次的估计值(0)x 附近,用泰勒级数(忽略二阶和以上的高阶项)表示它,可获得如式(1-2)的线性化变换后的方程组,该方程组被称为修正方程组。'()f x 是()f x 对于x 的一阶偏导数矩阵,这个矩阵便是重要的雅可比矩阵J 。 12(,,,)01,2, ,i n f x x x i n == (1-1) (0)'(0)(0)()()0f x f x x +?= (1-2) ' 由修正方程式可求出经过第一次迭代之后的修正量(0)x ?,并用修正量(0)x ?与估计值(0) x 之和,表示修正后的估计值(1)x ,表示如下(1-4)。 (0)'(0)1(0)[()]()x f x f x -?=- (1-3) (1)(0)(0)x x x =+? (1-4) 重复上述步骤。第k 次的迭代公式为: '()()()()()k k k f x x f x ?=- (1-5) (1)()()k k k x x x +=+? (1-6) 当采用直角坐标系解决潮流方程,此时待解电压和导纳如下式: i i i ij ij ij V e jf Y G jB =+=+ (1-7) 假设系统的网络中一共设有n 个节点,平衡节点的电压是已知的,平衡节点表示如下。 n n n V e jf =+ (1-8) }

用Matlab计算潮流计算电力系统分析

《电力系统潮流上机》课程设计报告 院系:电气工程学院 班级:电088班 学号: 0812002221 学生姓名:刘东昇 指导教师:张新松 设计周数:两周 日期:2010年 12 月 25 日

一、课程设计的目的与要求 目的:培养学生的电力系统潮流计算机编程能力,掌握计算机潮流计算的相关知识 要求:基本要求: 1.编写潮流计算程序; 2.在计算机上调试通过; 3.运行程序并计算出正确结果; 4.写出课程设计报告 二、设计步骤: 1.根据给定的参数或工程具体要求(如图),收集和查阅资料;学习相关软件(软件自选:本设计选择Matlab进行设计)。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 三、设计原理 1.牛顿-拉夫逊原理 牛顿迭代法是取x0 之后,在这个基础上,找到比x0 更接近的方程的跟,一步一步迭代,从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0 的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根。电力系统潮流计算,一般来说,各个母线所供负荷的功率是已知的,各个节点电压是未知的(平衡节点外)可以根据网络结构形成节点导纳矩阵,然后由节点导纳矩阵列写功率方程,由于功率方程里功率是已知的,电压的幅值和相角是未知的,这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组,需要将上述功率方程改写成功率平衡方程,并对功率平衡方程求偏导,得出对应的雅可比矩阵,给未知节点赋电压初值,一般为

额定电压,将初值带入功率平衡方程,得到功率不平衡量,这样由功率不平衡量、雅可比矩阵、节点电压不平衡量(未知的)构成了误差方程,解误差方程,得到节点电压不平衡量,节点电压加上节点电压不平衡量构成新的节点电压初值,将新的初值带入原来的功率平衡方程,并重新形成雅可比矩阵,然后计算新的电压不平衡量,这样不断迭代,不断修正,一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤: (1)形成各节点导纳矩阵Y。 (2)设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 (3)计算各个节点的功率不平衡量。 (4)根据收敛条件判断是否满足,若不满足则向下进行。 (5)计算雅可比矩阵中的各元素。 (6)修正方程式个节点电压 (7)利用新值自第(3)步开始进入下一次迭代,直至达到精度退出循环。 (8)计算平衡节点输出功率和各线路功率 2.网络节点的优化 1)静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数,然后,按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时,则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中,出线支路数最少的节点所对应的行中非零元素也2)动态地按增加出线支路数最少编号在上述的方法中,各节点的出线支路数是按原始网络统计出来的,在编号过程中认为固定不变的,事实上,在节点消去过程中,每消去一个节点以后,与该节点相连的各节点的出线支路数将发生变化(增加,减少或保持不变)。因此,如果每消去一个节点后,立即修正尚未编号节点的出线支路数,然后选其中支路数最少的一个节点进行编号,就可以预期得到更好的效果,动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3.MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写,主要包括:一般数值分析,矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致,所以不像学习高级语言那样难于掌握,而且编程效率和计算效率极高,还可在计算机上直接输出结果和精美的图形拷贝,所以它的确为一高效的科研助手。 四、设计内容

牛顿法潮流计算综述

潮流例题:根据给定的参数或工程具体要求(如图),收集和查阅资料;学习相关软件(软件自选:本设计选择Matlab进行设计)。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1.牛顿-拉夫逊原理 牛顿迭代法是取x0 之后,在这个基础上,找到比x0 更接近的方程的跟,一步一步迭代,从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0 的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根。电力系统潮流计算,一般来说,各个母线所供负荷的功率是已知的,各个节点电压是未知的(平衡节点外)可以根据网络结构形成节点导纳矩阵,然后由节点导纳矩阵列写功率方程,由于功率方程里功率是已知的,电压的幅值和相角是未知的,这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组,需要将上述功率方程改写成功率平衡方程,并对功率平衡方程求偏导,得出对应的雅可比矩阵,给未知节点赋电压初值,一般为额定电压,将初值带入功率平衡方程,得到功率不平衡量,这样由功率不平衡量、雅可比矩阵、节点电压不

平衡量(未知的)构成了误差方程,解误差方程,得到节点电压不平衡量,节点电压加上节点电压不平衡量构成新的节点电压初值,将新的初值带入原来的功率平衡方程,并重新形成雅可比矩阵,然后计算新的电压不平衡量,这样不断迭代,不断修正,一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤: (1)形成各节点导纳矩阵Y。 (2)设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 (3)计算各个节点的功率不平衡量。 (4)根据收敛条件判断是否满足,若不满足则向下进行。 (5)计算雅可比矩阵中的各元素。 (6)修正方程式个节点电压 (7)利用新值自第(3)步开始进入下一次迭代,直至达到精度退出循环。 (8)计算平衡节点输出功率和各线路功率 2.网络节点的优化 1)静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数,然后,按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时,则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中,出线支路数最少的节点所对应的行中非零元素也2)动态地按增加出线支路数最少编号在上述的方法中,各节点的出线支路数是按原始网络统计出来的,在编号过程中认为固定不变的,事实上,在节点消去过程中,每消去一个节点以后,与该节点相连的各节点的出线支路数将发生变化(增加,减少或保持不变)。因此,如果每消去一个节点后,立即修正尚未编号节点的出线支路数,然后选其中支路数最少的一个节点进行编号,就可以预期得到更好的效果,动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3.MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写,主要包括:一般数值分析,矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致,所以不像学习高级语言那样难于掌握,而且编程效率和计算效率极高,还可在计算机上直接输出结果和精美的图形拷贝,所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

电力系统稳态分析牛顿拉夫逊法

0引言 潮流是配电网络分析的基础,用于电网调度、运行分析、操作模拟和设计规划,同时也是电压优化和网络接线变化所要参考的内容。潮流计算通过数值仿真的方法把电力系统的详细运行情况呈现给工作人员,从而便于研究系统在给定条件下的稳态运行特点。随着市场经济的发展,经济利益是企业十分看重的,而线损却是现阶段阻碍企业提高效益的一大因素。及时、准确的潮流计算结果,可以给出配电网的潮流分布、理论线损及其在网络中的分布,从而为配电网的安全经济运行提供参考。从数学的角度来看,牛顿-拉夫逊法能有效进行非线性代数方程组的计算且具有二次收敛的特点,具有收敛快、精度高的特点,在输电网中得到广泛应用。随着现代计算机技术的发展,利用编程和相关软件,可以更好、更快地实现配电网功能,本文就是结合牛顿-拉夫逊法的基本原理,利用C++程序进行潮流计算,计算结果表明该方法具有良好的收敛性、可靠性及正确性。 1牛顿-拉夫逊法基本介绍 1.1潮流方程 对于N 个节点的电力网络(地作为参考节点不包括在内),如果网络结构和元件参数已知,则网络方程可表示为: =&& YV I (1-1) 式中,Y 为N*N 阶节点导纳矩阵;&V 为N*1维节点电压列向量;&I 为N*1维节点注入电流列向量。如果不计网络元件的非线性,也不考虑移相变压器,则Y 为对称矩阵。 电力系统计算中,给定的运行变量是节点注入功率,而不是节点注入电流,这两者之间有如下关系: ??=&&&EI S (1-2) 式中,&S 为节点的注入复功率,是N*1维列矢量;? &S 为&S 的共轭;

??i diag ??=???? &&E V 是由节点电压的共轭组成的N*N 阶对角线矩阵。 由(1-1)和(1-2),可得: ??=&&&S EYV 上式就是潮流方程的复数形式,是N 维的非线性复数代数方程组。将其展开,有: ?i i i ij j j i P jQ V Y V ∈-=∑&& j=1,2,….,N (1-3) 式中,j i ∈表示所有和i 相连的节点j ,包括j i =。 将节点电压用极坐标表示,即令i i i V V θ=∠&,代入式(1-3)中则有: ()i i i i ij ij j j j i P jQ V G jB V θθ∈-=∠-+∠∑ ()()cos sin i j ij ij ij ij j i V V G jB j θθ∈=+-∑ 故有: () ()cos sin sin cos i i j ij ij ij ij j i i i j ij ij ij ij j i P V V G B Q V V G B θθθθ∈∈?=+?? =-?? ∑∑ i=1,2,…,N (1-4) 式(1-4)是用极坐标表示的潮流方程。 而节点功率误差: (cos sin )θθ∈?=-+∑SP i i i j ij ij ij ij j i P P V V G B (1-5) (cos sin )θθ∈?=--∑SP i i i j ij ij ij ij j i Q Q V V G B (1-6) 式中:SP i P ,SP i Q 为节点i 给定的有功功率及无功功率。 1.2牛顿-拉夫逊法基本原理 1.2.1牛拉法的一般描述 牛拉法是把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程,即非线性问题通过线性化逐步近似,这就是牛拉法的核心。下面以非线性方程式的求解过程来进行说明。 设电力网络的节点功率方程一般形式如下:

(完整word版)9节点电力系统潮流计算

电力系统分析课程设计 设计题目9节点电力网络潮流计算 指导教师 院(系、部)电气与控制工程学院 专业班级 学号 姓名 日期

电气工程系课程设计标准评分模板

目录 1 PSASP软件简介 (1) 1.1 PSASP平台的主要功能和特点 (1) 1.2 PSASP的平台组成 (2) 2 牛顿拉夫逊潮流计算简介 (3) 2.1 牛顿—拉夫逊法概要 (3) 2.2 直角坐标下的牛顿—拉夫逊潮流计算 (5) 2.3 牛顿—拉夫逊潮流计算的方法 (6) 3 九节点系统单线图及元件数据 (7) 3.1 九节点系统单线图 (7) 3.2 系统各项元件的数据 (8) 4 潮流计算的结果 (10) 4.1 潮流计算后的单线图 (10) 4.2 潮流计算结果输出表格 (10) 5 结论 (14)

电力系统分析课程设计任务书9节点系统单线图如下: 基本数据如下:

表3 两绕组变压器数据 负荷数据

1 PSASP软件简介 “电力系统分析综合程序”(Power System Analysis Software Package,PSASP)是一套历史悠久、功能强大、使用方便的电力系统分析程序,是高度集成和开发具有我国自主知识产权的大型软件包。 基于电网基础数据库、固定模型库以及用户自定义模型库的支持,PSASP可进行电力系统(输电、供电和配电系统)的各种计算分析,目前包括十多个计算机模块,PSASP的计算功能还在不断发展、完善和扩充。 为了便于用户使用以及程序功能扩充,在PSASP7.0中设计和开发了图模一体化支持平台,应用该平台可以方便地建立电网分析的各种数据,绘制所需要的各种电网图形(单线图、地理位置接线图、厂站主接线图等);该平台服务于PSASP 的各种计算,在此之外可以进行各种分析计算,并输出各种计算结果。 1.1PSASP平台的主要功能和特点 PSASP图模一体化支持平台的主要功能和特点可概括为: 1. 图模支持平台具备MDI多文档操作界面,是一个单线图图形绘制、元件数据录入编辑、各种计算功能、结果显示、报表和曲线输出的集成环境。用户可以方便地建立电网数据、绘制电网图形、惊醒各种分析计算。人机交互界面全部汉化,界面良好,操作方便。 2. 真正的实现了图模一体化。可边绘图边建数据,也可以在数据已知的情况下进行图形自动快速绘制;图形、数据自动对应,所见即所得。 3. 应用该平台可以绘制各种电网图形,包括单线图、地理位置接线图、厂站主接线图等。 ●所有图形独立于各种分析计算,并为各计算模块所共享; ●可在图形上进行各种计算操作,并在图上显示各种计算结果; ●同一系统可对应多套单线图,多层子图嵌套; ●单线图上可细化到厂站主接线结构;

c语言编写的牛顿拉夫逊法解潮流程序

c语言编写的牛顿拉夫逊法解潮流程序 闲来无事,最近把牛拉法用c语言重写一遍,和matlab相比,c语言编写潮流程序最大的难点在于矩阵求逆,我使用的求逆方法是初等行变换法,程序段如下: #include #define N 3 void main() { int i,j,k; float t; float Jacob[N][N]={{1,2,2},{1,3,4},{2,3,4}};//欲进行求逆的矩阵 float inv_J[N][N];//逆矩阵存储于此 //初始化inv_J[N][N] for(i=0;i

} //输出逆矩阵 for(i=0;i #include #define N 4 //节点数 #define n_PQ 2 //PQ节点数 #define n_PV 1 //PV节点数 #define n_br 5 //串联支路数 void main() { void disp_matrix(float *disp_p,int disp_m,int disp_n); //矩阵显示函数 float Us[2*N]={1.0,0,1.0,0,1.05,0,1.05,0}; //电压初值 float Ps[N]={0,-0.5,0.2}; //有功初值 float Qs[N]={0,-0.3}; //无功初值float G[N][N],B[N][N]; //各几点电导电纳 struct //阻抗参数 { int nl; //左节点 int nr; //右节点 float R; //串联电阻值 float X; //串联电抗值 float Bl; //左节点并联电导 float Br; //右节点并联电纳 }ydata[n_br]={ {1,2,0,0.1880,-0.6815,0.6040}, {1,3,0.1302,0.2479,0.0129,0.0129}, {1,4,0.1736,0.3306,0.0172,0.0172}, {3,4,0.2603,0.4959,0.0259,0.0259}, {2,2,0,0.05,0,0} };

第3章作业答案电力系统潮流计算(已修订)

第三章 电力系统的潮流计算 3-1 电力系统潮流计算就是对给定的系统运行条件确定系统的运行状态。系 统运行条件是指发电机组发出的有功功率和无功功率(或极端电压),负荷的有 功功率和无功功率等。运行状态是指系统中所有母线(或称节点)电压的幅值和 相位,所有线路的功率分布和功率损耗等。 3-2 电压降落是指元件首末端两点电压的相量差。 电压损耗是两点间电压绝对值之差。当两点电压之间的相角差不大时, 可以近似地认为电压损耗等于电压降落的纵分量。 电压偏移是指网络中某点的实际电压同网络该处的额定电压之差。电压 偏移可以用kV 表示,也可以用额定电压的百分数表示。 电压偏移= %100?-N N V V V 功率损耗包括电流通过元件的电阻和等值电抗时产生的功率损耗和电压 施加于元件的对地等值导纳时产生的损耗。 输电效率是是线路末端输出的有功功率2P 与线路首端输入的有功功率 1P 之比。 输电效率= %1001 2 ?P P 3-3 网络元件的电压降落可以表示为 ()? ? ? ? ? +=+=-2221V V I jX R V V δ? 式中,?2V ?和? 2V δ分别称为电压降落的纵分量和横分量。 从电压降落的公式可见,不论从元件的哪一端计算,电压降落的纵、横分量计算公式的结构都是一样的,元件两端的电压幅值差主要有电压降落的纵分量决定,电压的相角差则由横分量决定。在高压输电线路中,电抗要远远大于电阻,即R X ??,作为极端的情况,令0=R ,便得 V QX V /=?,V PX V /=δ 上式说明,在纯电抗元件中,电压降落的纵分量是因传送无功功率而产生的,而电压降落的横分量则是因为传送有功功率产生的。换句话说,元件两端存在电压幅值差是传送无功功率的条件,存在电压相角差则是传送有功功率的条件。 3-4 求解已知首端电压和末端功率潮流计算问题的思路是,将该问题转化成 已知同侧电压和功率的潮流计算问题。

《电力系统分析》习题第3-6章(1)

3 简单电力系统潮流计算 3.1 思考题、习题 1)电力线路阻抗中的功率损耗表达式是什么?电力线路始、末端的电容功率表达式是什 么? 2)电力线路阻抗中电压降落的纵分量和横分量的表达式是什么? 3)什么叫电压降落、电压损耗、电压偏移、电压调整及输电效率? 5)对简单开式网络、变电所较多的开式网络和环形网络潮流计算的内容及步骤是什么? 6)变压器在额定状况下,其功率损耗的简单表达式是什么? 9)为什么要对电力网络的潮流进行调整控制?调整控制潮流的手段主要有哪些? 10)欲改变电力网络的有功功率和无功功率分布,分别需要调整网络的什么参数? 16)110kV 双回架空线路,长度为150kM ,导线型号为LGJ-120,导线计算外径为15.2mm , 三相导线几何平均距离为5m 。已知电力线路末端负荷为30+j15MVA ,末端电压为106kV ,求 始端电压、功率,并作出电压向量图。 17)220kV 单回架空线路,长度为200kM ,导线型号为LGJ-300,导线计算外径为24.2mm , 三相导线几何平均距离为7.5m 。已知电力线路始端输入功率为120+j50MVA ,始端电压为 240kV ,求末端电压、功率,并作出电压向量图。 18)110kV 单回架空线路,长度为80kM ,导线型号为LGJ-95,导线计算外径为13.7mm ,三 相导线几何平均距离为5m 。已知电力线路末端负荷为15+j10MVA ,始端电压为116kV ,求末 端电压和始端功率。 19)220kV 单回架空线路,长度为220kM ,电力线路每公里的参数分别为: kM S b kM x kM r /1066.2,/42.0,/108.06111-?=Ω=Ω=、 线路空载运行,当线路末端电压为205kV ,求线路始端的电压。 20)有一台三绕组变压器,其归算至高压侧的等值电路如图3-1所示,其中 ,68~,45~,8.3747.2,5.147.2,6547.232321MVA j S MVA j S j Z j Z j Z T T T +=+=Ω+=Ω-=Ω+=当变压器变比 为110/38.5(1+5%)/6.6kV ,U 3=6kV 时,试计算高压、中压侧的实际电压。 图3- 1 图3-2

牛顿——拉夫逊法进行潮流计算

%本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳 % 5、支路的变比;6、支路首端处于K侧为1,1侧为0 % B2矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值% 4、PV节点电压V的给定值;5、节点所接的无功补偿设备的容量% 6、节点分类标号:1为平衡节点(应为1号节点);2为PQ节点;% 3为PV节点; clear; n=10;%input('请输入节点数:n='); nl=10;%input('请输入支路数:nl='); isb=1;%input('请输入平衡母线节点号:isb='); pr=0.00001;%input('请输入误差精度:pr='); B1=[1 2 0.03512+0.08306i 0.13455i 1 0; 2 3 0.0068+0.18375i 0 1.02381 1; 1 4 0.05620+0.13289i 0.05382i 1 0; 4 5 0.00811+0.24549i 0 1.02381 1; 1 6 0.05620+0.13289i 0.05382i 1 0; 4 6 0.04215+0.09967i 0.04037i 1 0; 6 7 0.0068+0.18375i 0 1.02381 1; 6 8 0.02810+0.06645i 0.10764i 1 0; 8 10 0.00811+0.24549i 0 1 1; 8 9 0.03512+0.08306i 0.13455i 1 0] B2=[0 0 1.1 1.1 0 1; 0 0 1 0 0 2; 0 0.343+0.21256i 1 0 0 2; 0 0 1 0 0 2; 0 0.204+0.12638i 1 0 0 2; 0 0 1 0 0 2; 0 0.306+0.18962i 1 0 0 2; 0 0 1 0 0 2; 0.5 0 1.1 1.1 0 3; 0 0.343+0.21256i 1 0 0 2] ;%input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); % % %--------------------------------------------------- for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end

相关文档
相关文档 最新文档