电力系统分析实验指导书.docx

电力系统分析

实验指导书

郑州科技学院电气工程及其自动化 教研室编

目 录

实验一 生成节点导纳矩阵的matlab程序设计 3

实验二 支路追加法生成 节点阻抗矩阵的matlab程序设计 7

实验三 牛顿-拉夫逊法电力系统潮流计算的matlab程序设计 15

实验四 简化模型的静态稳定计算 的matlab程序设计 26

实验五 短路电流计算程序的实现 33

PAGE

PAGE31

实验一 生成节点导纳矩阵的matlab程序设计

一、实验目的：

本实验通过对生成节点导纳矩阵的编程与调试，获得复杂电力系统的节点导纳矩阵，为电力系统的潮流计算程序编制打下基础。通过实验教学加深学生对节点导纳矩阵计算方法的理解，掌握变压器为非标准变比时的修正方法，熟悉MATLAB软件的使用方法，提高编制调试计算机程序的能力。

二、实验器材：

计算机、软件（已安装应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等）。

三、实验内容：

节点导纳矩阵的形成

节点导纳矩阵的对角元Yii(i?1,2,...n)称自导纳，自导纳Yii是节点i以外的所有节点都接地，由节点i向整个网络看到而得到的导纳。在节点i加上单位大小的电压（U?i?1的电压）时，由节点i流向网络的电流就等于节点i的自导纳

?I?Yii?i

?

Ui

U?j?0,j?i

通过计算不难发现，Yii就等于与节点i连接的所有支路导纳的和。

节点导纳矩阵的非对角元Yij(i?1,2,...n；j?1,2,...n；i?

称互导纳，即把节点j以外的节

I?iU?j点全接地，而在节点

I?

i

U?j

Yij?

U?k?0,k?j

同样，Yij其实就是连接节点i和节点j支路的导纳值再加上负号。显然，Yij?Yji。而且，如节点i，j之间没有直接联系，也不计两支路之间，例如两相邻电力线路之间的互感时，

Yij?Yji?0。

互导纳的这些性质决定了节点导纳矩阵是一个对称的稀疏矩阵。而且，由于每个节点所连接的支路数总有一定限度，随着网络中节点数的增加，非零元素数相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素数的比值也就愈来愈高。

导纳矩阵的计算归结如下：

导纳矩阵的阶数等于电力系统网络的节点数。

导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。

导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连支路的导纳之和。

Yii??yij

j?i

式中，yij为节点i与节点j间支路阻抗Zij倒数，符号j?i表示j属于i或与i连接的j，即?内只包括与节点i直接相连的节点j。当节点i有接地支路时，还应包括j?0的情况。

导纳矩阵非对角元素Yij等于节点i与节点j之间的导纳的负数

Yij

??1

zij

??yij

当i、j之间有多条并联支路时，求Yij时应求所有并联支路导纳的代数和的负数

L

Yij

???y(s)

ijs?1

ij

式中，L表示i、j之间并联支路数的条数。

按照以上计算式，对于实际网络均可根据给定的支路参数和连接情况，直观而简单地求出导纳矩阵。可以看出，用以上计算公式求得的导纳矩阵与根据定义得到的导纳矩阵是完全一致的。

变压器为非标准变比时的修正

无论采用有名制或标么制，凡涉及多电压级网络的计算，在精确计算时都必须将网络中所有参数和变量按实际变比归算到同一电压等级。实际上，在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比，也就是不等于标准变比，而且变压器的变比在运行中是可以改变的。这将使每改变一次变比都要重新计算元件参数，很不方便。下面将介绍另一种可等值地体现变压器电压变换功能的模型。

等值变压器模型

?

Z?? I1

ZT 1:K

?

I2 Z?

1

1

?

U1

2

?

U2

图1非标准变比时的修正电路

? ?

I1 KZT I2

1

1

2

?

U1

KZT

K?1

T

1?K

K2Z

?

U2

图2以变压器阻抗表示

增加非零非对角元素为

?Y

??YT

ij

节点i的自导纳，增加一个改变量为

ji K

Y?K?1Y

1Y?Y

ii K T KT T

节点j的自导纳，也增加一个改变量为

Y?1Y

1?KY

?1Y

程序代码及说明

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

jj KT

K2 T

K2T

nl=inp