用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．网络节点的优化

１）静态地按最少出线支路数编号

这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。

3．MATLAB编程应用

Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。

四、设计容

1.设计流程图

2.程序

clear;clc

%重新编号，把原题中的节点1,2,3,4,5重新依次编号为5,1,2,3,4，其中1-4号为PQ节点，5号为平衡节点

y=0;

%输入原始数据，求节点导纳矩阵

y (1,2)=1/(0.06+0.18i); y (1,3)=1/(0.06+0.18i); y (1,4)=1/(0.04+0.12i);

y(1,5)=1/(0.02+0.06i);

y(2,3)=1/(0.01+0.03i);y(2,5)=1/(0.08+0.24i);

y(3,4)=1/(0.08+0.24i);

y(4,5)=0;

for i=1:5

for j=i:5

y(j,i)=y(i,j);

end

end

Y=0;

%求互导纳

for i=1:5

for j=1:5

if i~=j

Y(i,j)=-y(i,j);

end

end

end

%求自导纳

for i=1:5

Y(i,i)=sum(y(i,:));

end

Y %Y 为导纳矩阵

G=real(Y);

B=imag(Y);

%原始节点功率

S(1)=0.2+0.2i;

S(2)=-0.45-0.15i;

S(3)=-0.4-0.05i;

S(4)=-0.6-0.1i;

S(5)=0;

P=real(S);

Q=imag(S);

%赋初值

U=ones(1,5);U(5)=1.06;

e=zeros(1,5);

ox=ones(8,1);fx=ones(8,1);

count=0 %计算迭代次数

while max(fx)>1e-5

for i=1:4

for j=1:4

H(i,j)=0;N(i,j)=0;M(i,j)=0;L(i,j)=0;oP(i)=0;oQ(i)=0;

end

end

for i=1:4

for j=1:5 oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));

oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));

end

oP(i)=oP(i)+P(i); oQ(i)=oQ(i)+Q(i);

end

fx=[oP,oQ]';

%求雅克比矩阵

%当i~=j时候求H,N,M,L 如下：

for i=1:4

for j=1:4

if i~=j H(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));

N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));

L(i,j)=H(i,j);

M(i,j)=-N(i,j);

end

end

end

H,N,M,L

%当i=j 时H,N,M,L如下：

for i=1:4

for j=1:5

if i~=j

H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i, j)*cos (e(i)-e(j))); N(i,i)=N(i,i)-U(i)*U(j)*(G(i, j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));

M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));

L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));

end

end

N(i,i)=N(i,i)-2*(U(i))^2*G(i,i);

L(i,i)=L(i,i)+2*(U(i))^2*B(i,i);

end

J=[H,N;M,L] %J 为雅克比矩阵

ox=-((inv(J))*fx);

for i=1:4

oe(i)=ox(i); oU(i)=ox(i+4)*U(i);

end

for i=1:4

e(i)=e(i)+oe(i); U(i)=U(i)+oU(i);

end

count=count+1;

end

ox,U,e,count

%求节点注入的净功率

i=5;

for j=1:5

P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)))+P(i);

Q(i)=U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)))+Q(i);

end

S(5)=P(5)+Q(5)*sqrt(-1);

S

%求节点注入电流

I=Y*U'

3.运行结果

Y值：

迭代过程：

电压值：

平衡节点注入功率及电流：

五、课程设计心得与体会

在电力系统分析课程中，我们学习了电力网络的数学模型和确定解算方法。在本次课程设计中我们利用前面的基础知识进行制定计算流程和编制计算程序，实践了牛顿-拉夫逊法，体会到其用性。根据潮流计算的基本要求，对潮流计算可归纳为下面几点：（1）计算方法的可靠性或收敛性；（2）对计算速度和存量的要求；（3）计算的方便性和灵活性。基于这几点采用了Matlab进行编程计算。电力系统潮流计算都是基于矩阵的迭代运算，而Matlab 语言正是以处理矩阵见长，实践证明，Matlab 语言在电力系统潮流计算仿真研究中的应用是可行的，而且由于其强大的矩阵处理功能，完全可以应用于电力系统的其它分析计算中，用Matlab语言编程效率高，程序调试十分方便。

电路分析总复习题-分析计算题

三、计算分析题 1、图1.5.1所示电路，已知U =3V ，求R 。（2k Ω） 2、图1.5.2所示电路，已知U S ＝3V ，I S ＝2A ，求U AB 和I 。（1V 、5A ） 3、电路如图1.5.5所示，求10V 电压源发出的 功率。 （－35W ） 4、分别计算S 打开与闭合时图1.5.6电路中A 、B 两点的电位。（S 打开：A －10.5V,B －7.5V S 闭合：A 0V ，B 1.6V ） 5、试求图1.5.7所示电路的入端电阻R AB 。（150Ω） U － 图1.5.1 1Ω 图1.5.2 6V 图1.5.5 B －图1.5.6 Ω 图1.5.7

6、试求图2.4.1所示电路的电压U 。 7、已知图2.5.1电路中电压U =4.5V ，试应用已经学过的电路求解法求电阻R 。 （18Ω） 8、求解图2.5.2所示电路的戴维南等效电路。 （U ab =0V ，R 0=8.8Ω） 9、列出图2.5.4所示电路的结点电压方程。 解：画出图2.5.4等效电路图如下： 图2.5.1 9V 图2.5.2 2A Ω U 图2.4.1题电路

对结点A 对结点B 10、应用等效变换求图示电路中的I的值。（10分） 解：等效电路如下： 11、应用等效变换求图示电路中的I的值。

12、应用戴维南定理求解图示电路中的电流I 13、如下图所示，RL等于何值时，能得到最大传输功率P0max？并计算P0max 。

16、图示电路中，开关闭合之前电路已处于稳定状态，已知R1=R2=2Ω 解开关闭合后电感电流iL的全响应表达式。 17、图示电路中，t=0时开关闭合，闭合之前电路已处于稳定状态，请用三要素法求解开关闭合后电容电压uc的全响应表达式。

matlab电力系统潮流计算

华中科技大学 信息工程学院课程设计报告书题目: 电力系统潮流计算 专业：电气工程及其自动化 班级： 学号： 学生姓名： 指导教师： 2015年 11 月 10 日

2015年11月12日

信息工程学院课程设计成绩评定表

摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有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)实例计算

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

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

电力系统分析潮流计算例题

电力系统的潮流计算 西安交通大学自动化学院 2012.10 3．1 电网结构如图3—11所示，其额定电压为10KV 。已知各节点的负荷功率及参数： MVA j S )2.03.0(2 +=， MVA j S )3.05.0(3+=， MVA j S )15.02.0(4+= Ω+=)4.22.1(12j Z ，Ω+=)0.20.1(23j Z ，Ω+=)0.35.1(24j Z 试求电压和功率分布。 解：（1）先假设各节点电压均为额定电压，求线路始端功率。 0068.00034.0)21(103.05.0)(2 2223232232323j j jX R V Q P S N +=++=++=?0019.00009.0)35.1(10 15.02.0)(2 2 224242242424j j jX R V Q P S N +=++=++=?

则： 3068.05034.023323j S S S +=?+= 1519.02009.024424j S S S +=?+= 6587.00043.122423' 12 j S S S S +=++= 又 0346 .00173.0)4.22.1(106587.00043.1)(2 2 212122'12'1212j j jX R V Q P S N +=++=++=? 故： 6933.00216.112'1212 j S S S +=?+= （2） 再用已知的线路始端电压kV V 5.101 =及上述求得的线路始端功率 12 S ，求出线 路 各 点 电 压 。

kV V X Q R P V 2752.05 .104.26933.02.10216.1)(11212121212=?+?=+=? kV V V V 2248.101212=?-≈ kV V V V kV V X Q R P V 1508.100740.0) (24242 2424242424=?-≈?=+=? kV V V V kV V X Q R P V 1156.101092.0) (23232 2323232323=?-≈?=+=? （3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率。 0066.00033.0)21(12.103.05.02 2 223j j S +=++=? 0018.00009.0)35.1(15 .1015.02.02 2 224j j S +=++=? 故 3066.05033.023323j S S S +=?+= 1518.02009.024424j S S S +=?+= 则 6584.00042.122423' 12 j S S S S +=++= 又 0331.00166.0)4.22.1(22 .106584.00042.12 2 212j j S +=++=? 从而可得线路始端功率 6915.00208.112 j S +=

matlab潮流计算

附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % 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 all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=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 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); %对角元K侧 Y(p,p)=Y(p,p)+1、/B1(i,3)+B1(i,4); %对角元1侧 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));

基于MATLAB的电力系统潮流计算

基于MATLAB的电力系统潮流计算 %简单潮流计算的小程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号，将变压器的串联阻抗置于低压侧处理。 %第三列为支路的串列阻抗参数。 %第四列为支路的对地导纳参数。 %第五烈为含变压器支路的变压器的变比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，%“0”为不含有变压器。 %B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中 %“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。 %X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地%参数。 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n); Times=1; %置迭代次数为初始值 %创建节点导纳矩阵 for i=1:n1 if 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)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3);

电力系统分析试题答案(全)

2、停电有可能导致人员伤亡或主要生产设备损坏的用户的用电设备属于（ ）。 A 、一级负荷； B 、二级负荷； C 、三级负荷； D 、特级负荷。 4、衡量电能质量的技术指标是（ ）。 A 、电压偏移、频率偏移、网损率； B 、电压偏移、频率偏移、电压畸变率； C 、厂用电率、燃料消耗率、网损率； D 、厂用电率、网损率、电压畸变率 5、用于电能远距离输送的线路称为（ ）。 A 、配电线路； B 、直配线路； C 、输电线路； D 、输配电线路。 7、衡量电力系统运行经济性的主要指标是（ ）。 A 、燃料消耗率、厂用电率、网损率； B 、燃料消耗率、建设投资、网损率； C 、网损率、建设投资、电压畸变率； D 、网损率、占地面积、建设投资。 8、关于联合电力系统，下述说法中错误的是（ ）。 A 、联合电力系统可以更好地合理利用能源； B 、在满足负荷要求的情况下，联合电力系统的装机容量可以减少； C 、联合电力系统可以提高供电可靠性和电能质量； D 、联合电力系统不利于装设效率较高的大容量机组。 9、我国目前电力系统的最高电压等级是（ ）。 A 、交流500kv ，直流kv 500±； B 、交流750kv ，直流kv 500±； C 、交流500kv ，直流kv 800±；； D 、交流1000kv ，直流kv 800±。 10、用于连接220kv 和110kv 两个电压等级的降压变压器，其两侧绕组的额定电压应为（ ）。 A 、220kv 、110kv ； B 、220kv 、115kv ； C 、242Kv 、121Kv ； D 、220kv 、121kv 。 11、对于一级负荷比例比较大的电力用户，应采用的电力系统接线方式为（ ）。 A 、单电源双回路放射式； B 、双电源供电方式； C 、单回路放射式接线； D 、单回路放射式或单电源双回路放射式。 12、关于单电源环形供电网络，下述说法中正确的是（ ）。 A 、供电可靠性差、正常运行方式下电压质量好； B 、供电可靠性高、正常运行及线路检修（开环运行）情况下都有好的电压质量； C 、供电可靠性高、正常运行情况下具有较好的电压质量，但在线路检修时可能出现电压质量较差的情况； D 、供电可靠性高，但电压质量较差。 13、关于各种电压等级在输配电网络中的应用，下述说法中错误的是（ ）。 A 、交流500kv 通常用于区域电力系统的输电网络； B 、交流220kv 通常用于地方电力系统的输电网络； C 、交流35kv 及以下电压等级通常用于配电网络； D 、除10kv 电压等级用于配电网络外，10kv 以上的电压等级都只能用于输电网络。 14、110kv 及以上电力系统应采用的中性点运行方式为（ ）。 A 、直接接地； B 、不接地； C 、经消弧线圈接地； D 、不接地或经消弧线圈接地。 16、110kv 及以上电力系统中，架空输电线路全线架设避雷线的目的是（ ）。

电力系统分析课后作业题及练习题

第一章 电力系统的基本概念 1-1 什么叫电力系统、电力网及动力系统 1-2 电力线、发电机、变压器和用电设备的额定电压是如何确定的 1-3 我国电网的电压等级有哪些 1-4 标出图1-4电力系统中各元件的额定电压。 1-5 请回答如图1－5所示电力系统中的二个问题： ⑴ 发电机G 、变压器1T 2T 3T 4T 、三相电动机D 、单相电灯L 等各元件的额定电压。 ⑵ 当变压器1T 在+%抽头处工作，2T 在主抽头处工作，3T 在%抽头处工作时，求这些变压器的实际变比。 1-6 图1-6中已标明各级电网的电压等级。试标出图中发电机和电动机的额定电压及变压器的额定变比。 1-7 电力系统结线如图1-7所示，电网各级电压示于图中。试求： ⑴发电机G 和变压器1T 、2T 、3T 高低压侧的额定电压。 习题1-4图

⑵设变压器1T 工作于+%抽头， 2T 工作于主抽头，3T 工作于-5%抽头，求这些变压器的实际变比。 1-8 比较两种接地方式的优缺点，分析其适用范围。 1-9 什么叫三相系统中性点位移它在什么情况下发生中性点不接地系统发生单相接地时，非故障相电压为什么增加3倍 1-10 若在变压器中性点经消弧线圈接地，消弧线圈的作用是什么 第二章 电力系统各元件的参数及等值网络 2-1 一条110kV 、80km 的单回输电线路，导线型号为LGJ —150，水平排列，其线间距离为4m ，求此输电线路在40℃时的参数，并画出等值电路。 2-2 三相双绕组变压器的型号为SSPL —63000/220，额定容量为63000kVA ，额定电压为242/，短路损耗404=k P kW ，短路电压45.14%=k U ，空载损耗93=o P kW ，空载电流 41.2%=o I 。求该变压器归算到高压侧的参数，并作出等值电路。 2-3 已知电力网如图2-3所示： 各元件参数如下： 变压器：1T ：S =400MVA ，12%=k U ， 242/ kV 2T ：S =400MVA ，12%=k U ， 220/121 kV 线路：2001=l km, /4.01Ω=x km (每回路) 602=l km, /4.01Ω=x km 115kV T 1 T 2 l 1 l 2 习题2-3图

基于MATLAB的潮流计算源程序代码(优.选)

%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;

电力系统分析练习题

电力系统分析练习题 一、单项选择题 1.对电力系统的基本要求就是(A) A.保证供电可靠性、保证良好的电能质量、保证系统运行的经济性 B.保证供电可靠性、保证良好的电能质量、保证继电保护动作的选择性 C.保证供电可靠性与良好的电能质量 D.保证供电可靠性与系统运行的经济性 2.下图所示的电力系统中,变压器T1的额定变比应为(B) A.242/121/11kV B.220/121/11kV C.220/110/10kV D.242/121/10kV 3.连接220kV 电力系统与110kV 电力系统的降压变压器,其额定变比应为(D) A.220/110kV B.220/115、5kV C.242/121kV D.220/121kV 4.我国110kV 以上电力系统中性点通常采用的运行方式就是(B) A 、不接地 B 、直接接地 C 、经消弧线圈接地 D 、经电容器接地 5.电力系统中性点经消弧线圈接地时,应采用的补偿方式为(A) A.过补偿 B.欠补偿 C.全补偿 D.全补偿或欠补偿 6.采用同一型号导线的三相架空输电线路,相间距离增大时,其电容(B) A.增大 B.减小 C.不变 D.无法确定 7.架空输电线路采用分裂导线的目的就是(A) A 、减小线路电抗 B 、减小线路电阻 C 、减小线路电容 D 、增大线路电抗 8.三相架空输电线路导线全换位的目的就是(A) A.减小三相参数的不平衡 B.减小线路电抗 C.减小线路电容 D.提高输电线路的输送容量 9.变压器参数B T 由试验数据 确定。(2章)(B) A.%K U B.%0I C.0P ? D.K P ? 10.变压器的电抗参数 T X ,由实验数据 确定。(A) A 、%K U B 、%0I C 、0P ? D 、K P ? 11.变压器的电纳参数T G ,由实验数据 确定。(C) A.%K U B.%0I C.0P ? D.K P ? 12.变压器的电阻参数T R ,由实验数据 确定。(D) A.%K U B.%0I C.0P ? D.K P ? 13.如果三相功率基准值为b S 、线电压基准值为b U ,则阻抗基准值为(D) A.b b U S / B. b b U S 3/ C.2/b b U S D.b b S U /2 14.输电线路运行时,其始端电压电压与末端电压的大小关系就是(C) A.始端电压一定高于末端电压 B. 始端电压始终等于末端电压 C.在不计电压降落横分量影响情况下,PR+QX>0时,始端电压高于末端电压 D.在不计电压降落横分量影响情况下,PR+QX>0时,始端电压低于末端电压 15.输电线路的电压降落就是指(A) A 、线路始端电压与末端电压的相量差 B 、线路始端电压与末端电压的数值差 C 、线路始端电压与额定电压的数值差 D 、线路末端电压与额定电压的数值差

基于matlab--psat软件的电力系统潮流计算课程设计

东北电力大学课程设计改革试用任务书: 电力系统潮流计算课程设计任务书 设计名称：电力系统潮流计算课程设计 设计性质：理论计算，计算机仿真与验证 计划学时：两周 一、设计目的 1.培养学生独立分析问题、解决问题的能力； 2.培养学生的工程意识，灵活运用所学知识分析工程问题的能力 3.编制程序或利用电力系统分析计算软件进行电力系统潮流分析。 二、原始资料 1、系统图：IEEE14节点。 2、原始资料：见IEEE14节点标准数据库 三、课程设计基本内容： 1.采用PSAT仿真工具中的潮流计算软件计算系统潮流； 1)熟悉PSAT仿真工具的功能； 2)掌握IEEE标准数据格式内容； 3)将IEEE标准数据转化为PSAT计算数据； 2.分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况， 分析两种方法计算潮流的优缺点； 3.分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。

4.选择以下内容之一进行分析： 1)找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施 的特点，并仿真验证； 2)找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证； 3)找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较 各种措施的特点，并仿真验证； 5.任选以下内容之一作为深入研究：（不做要求） 1)找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有 功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗； 2)找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功 功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗； 3)找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压， 并有效调整发电机端电压，提高该节点电压水平； 四、课程设计成品基本要求： 1.绘制系统潮流图，潮流图应包括： 1)系统网络参数 2)节点电压幅值及相角 3)线路和变压器的首末端有功功率和无功功率 2.撰写设计报告，报告内容应包括以下几点： 1)本次设计的目的和设计的任务； 2)电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点； 3)对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性； 4)找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施； 5)分析各种调整措施的特点并比较它们之间的差异； 6)结论部分以及设计心得； 五、考核形式 1.纪律考核：学生组织出勤情况和工作态度等； 2.书面考核：设计成品的完成质量、撰写水平等； 3.答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩； 4.采用五级评分制：优、良、中、及格、不及格五个等级。

第三章简单电力系统的潮流计算汇总

第一章 简单电力系统的分析和计算 一、 基本要求 掌握电力线路中的电压降落和功率损耗的计算、变压器中的电压降落和功率损耗的计 算；掌握辐射形网络的潮流分布计算；掌握简单环形网络的潮流分布计算；了解电力网络的简化。 二、 重点内容 1、电力线路中的电压降落和功率损耗 图3-1中，设线路末端电压为2U 、末端功率为222~jQ P S +=，则 （1）计算电力线路中的功率损耗 ① 线路末端导纳支路的功率损耗： 222 2* 222~U B j U Y S Y -=?? ? ??=? ……………（3-1） 则阻抗支路末端的功率为： 222~~~Y S S S ?+=' ② 线路阻抗支路中的功率损耗： ()jX R U Q P Z I S Z +'+'==?2 2 22222 ~ ……（3-2） 则阻抗支路始端的功率为： Z S S S ~ ~~21?+'=' ③ 线路始端导纳支路的功率损耗： 2121* 122~U B j U Y S Y -=?? ? ??=? …………（3-3） 则线路始端的功率为： 111~ ~~Y S S S ?+'= ~~~图3-3 变压器的电压和功率 ~2 ? U （2）计算电力线路中的电压降落 选取2U 为参考向量，如图3-2。线路始端电压 U j U U U δ+?+=2 1 其中 2 2 2U X Q R P U '+'= ? ； 222U R Q X P U '-'=δ ……………（3-4） 则线路始端电压的大小： ()()2 221U U U U δ+?+= ………………（3-5） 一般可采用近似计算： 2 2 2221U X Q R P U U U U '+'+ =?+≈ ………………（3-6）

电力系统分析计算题.doc

电力系统分析计算题 1.假设电压U 1固定不变，试就图示系统分析为什么投入电容器C 后可以降低线损和提高电压U 2。 2. 试就图示系统分析为什么对中枢点U 1进行逆调压可以对负荷点电压U 2进行控制的原理。 3. 试就图示系统分析（a ）、（b ）两种情况下线路的电能损耗ΔA ，你的结论是什么？ 4.某负荷由发电厂经电压为110kV 的输电线路供电。 Ω, X=41Ω，C Q 2 1 =1.7Mvar ， 发 线路的参数为：R=17 电厂高压母线电压U 1=116kV, 线路末端负荷为20+j10 MVA,求输电线路的功率损耗和末端电压U 2（计及电压降落的横分量)。(12分) 5、某降压变电所装有两台并联工作的有载调压变压器，电压为110±5×2.5/11 kV ，容量为31.5MVA, 已知最大负荷时:高压母线电压为103 kV,两台变压器并列运行时的电压损耗为5.849 kV,最小负荷时: 高压母线电压为108.5 kV, 两台变压器并列运行时的电压损耗为2.631 kV 。变电所低压母线要求逆调压，试选择有载调压变压器分接头。(12分) 6.一双电源电力系统如下图所示，如在f 点发生b 、 P+jQ P+jQ S ~ （h ） （a ） （b ） 20+j10 MVA U 2 2 T 1 1 MVA X X j E d G 3013.013.01"2., ,. 1. === U K %=10.5 60MV A U K %=10.5 31.5MV A 80km X 1=0.4Ω/km X 0=3.5 X 1 MVA X X j E d G 60125.0125.01".2, ,. 2. ===

MATLAB下的潮流计算实现-稀疏技术毕业设计

毕业设计（论文）MATLAB下的潮流计算实现-稀疏技术

毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作者签名：日期： 指导教师签名：日期： 使用授权说明 本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：日期：

摘要 电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。因此潮流计算是研究电力系统的一种很重要和很基础的计算。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流问题的方程式阶数越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。 本文旨在于研究潮流计算的牛顿—拉夫逊法的基本原理，在Matlab环境中实现牛顿—拉夫逊法潮流计算的数学模型，程序流程以及编制相应程序，并在程序中融合了节点优化编号和稀疏技术，以提高计算效率。最后用IEEE-3O节点标准测试系统验证所编程序。 关键词：潮流计算Newtom-Raphson法节点优化稀疏技术Matlab ABSTRACT Power flow calculation is fundanmental of analysis. Network reconfiguration,fault management,state estimator etc also need the data of electrial system power flow.There is important significance to develop power flow calculation in allusion to traits of distribution network. This paper introduces the principle of Newtom-Raphson algorithm, which is developed for calculation of power flow calculation ,where zero sequence network is open.With this algorithm,the three-phase load is resolved into positive/negative sequence power and coupling power,thus,decoupling three phase power flow into sequencet component power flow.The power flow can be obtained by just finding the positive sequence power flow and then finding the negative sequent component from the coupling https://www.wendangku.net/doc/726214587.html,pared with the existing methods,the jacobian matrix with the proposed algorithm is of much lower order,thus substantially reducing the computation burden.The proposed algorithm,together with a reference algorithm,has been simulated on an actual IEEE-30 system using statistic load date.And then it will

电力系统潮流计算课程设计(终极版)

目录 摘要................................................. - 1 - 1.设计意义与要求..................................... - 2 - 1.1设计意义 ...................................... - 2 - 1.2设计要求（具体题目）........................... - 2 - 2.题目解析........................................... - 3 - 2.1设计思路 ...................................... - 3 - 2.2详细设计 ...................................... - 4 - 2.2.1节点类型.................................. - 4 - 2.2.2待求量 ................................... - 4 - 2.2.3导纳矩阵.................................. - 4 - 2.2.4潮流方程.................................. - 5 - 2.2.5牛顿—拉夫逊算法.......................... - 6 - 2.2.5.1牛顿算法数学原理：................... - 6 - 2.2.5.2修正方程............................. - 7 - 2.2.5.3收敛条件............................. - 9 - 3.结果分析.......................................... - 10 - 4.小结.............................................. - 11 - 参考文献............................................ - 12 -

Matlab牛拉法潮流计算程序

%本程序的功能是用牛顿——拉夫逊法进行潮流计算 % 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=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);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 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; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧 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节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量 end %=================================================================== P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2