电力系统分析实验指导书.docx
- 文档编号:707702
- 上传时间:2023-04-29
- 格式:DOCX
- 页数:41
- 大小:383.70KB
电力系统分析实验指导书.docx
《电力系统分析实验指导书.docx》由会员分享,可在线阅读,更多相关《电力系统分析实验指导书.docx(41页珍藏版)》请在冰点文库上搜索。
电力系统分析
实验指导书
郑州科技学院电气工程及其自动化 教研室编
目 录
实验一 生成节点导纳矩阵的matlab程序设计 3
实验二 支路追加法生成 节点阻抗矩阵的matlab程序设计 7
实验三 牛顿-拉夫逊法电力系统潮流计算的matlab程序设计 15
实验四 简化模型的静态稳定计算 的matlab程序设计 26
实验五 短路电流计算程序的实现 33
实验一 生成节点导纳矩阵的matlab程序设计
一、实验目的:
本实验通过对生成节点导纳矩阵的编程与调试,获得复杂电力系统的节点导纳矩阵,为电力系统的潮流计算程序编制打下基础。
通过实验教学加深学生对节点导纳矩阵计算方法的理解,掌握变压器为非标准变比时的修正方法,熟悉MATLAB软件的使用方法,提高编制调试计算机程序的能力。
二、实验器材:
计算机、软件(已安装应用软件MATLAB等)、移动存储设备(学生自备,软盘、U盘等)。
三、实验内容:
1.节点导纳矩阵的形成
节点导纳矩阵的对角元Yii(i=1,2,...n)称自导纳,自导纳Yii是节点i以外的所有节点都接地,由节点i向整个网络看到而得到的导纳。
在节点i加上单位大小的电压(U&i=1的电压)时,由节点i流向网络的电流就等于节点i的自导纳
39
&
I&Yii=i
Ui
U&j=0,j¹i
通过计算不难发现,Yii就等于与节点i连接的所有支路导纳的和。
节点导纳矩阵的非对角元Yij(i=1,2,...n;j=1,2,...n;i¹
j)称互导纳,即把节点j以外的节
I&
i
U&j
点全接地,而在节点j加以单位电压时,由节点i流向j的电流加上负号的值。
Yij=
U&k=0,k¹j
同样,Yij其实就是连接节点i和节点j支路的导纳值再加上负号。
显然,Yij=Yji。
而且,如节点i,j之间没有直接联系,也不计两支路之间,例如两相邻电力线路之间的互感时,
Yij=Yji=0。
互导纳的这些性质决定了节点导纳矩阵是一个对称的稀疏矩阵。
而且,由于每个节点所连接的支路数总有一定限度,随着网络中节点数的增加,非零元素数相对愈来愈少,节点导纳矩阵的稀疏度,即零元素数与总元素数的比值也就愈来愈高。
导纳矩阵的计算归结如下:
(1)导纳矩阵的阶数等于电力系统网络的节点数。
(2)导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。
(3)导纳矩阵的对角元素,即各节点的自导纳等于相应节点所连支路的导纳之和。
Yii=åyij
jÎi
式中,yij为节点i与节点j间支路阻抗Zij倒数,符号jÎi表示j属于i或与i连接的j,即S内只包括与节点i直接相连的节点j。
当节点i有接地支路时,还应包括j=0的情况。
(4)导纳矩阵非对角元素Yij等于节点i与节点j之间的导纳的负数
Yij
=-1
zij
=-yij
当i、j之间有多条并联支路时,求Yij时应求所有并联支路导纳的代数和的负数
L
Yij
=-åy(s)
ij
s=1
式中,L表示i、j之间并联支路数的条数。
按照以上计算式,对于实际网络均可根据给定的支路参数和连接情况,直观而简单地求出导纳矩阵。
可以看出,用以上计算公式求得的导纳矩阵与根据定义得到的导纳矩阵是完全一致的。
2.变压器为非标准变比时的修正
无论采用有名制或标么制,凡涉及多电压级网络的计算,在精确计算时都必须将网络中所有参数和变量按实际变比归算到同一电压等级。
实际上,在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比,也就是不等于标准变比,而且变压器的变比在运行中是可以改变的。
这将使每改变一次变比都要重新计算元件参数,很不方便。
下面将介绍另一种可等值地体现变压器电压变换功能的模型。
等值变压器模型
·
ZI¢ I1
ZT 1:
K
·¢
I2 ZC
1
·
U1
2
·
U2
图1非标准变比时的修正电路
· ·
I1 KZT I2
1
2
·
U1
KZT
K-1
T
1-K
K2Z
·
U2
图2以变压器阻抗表示
增加非零非对角元素为
Y=Y
=-YT
ij
节点i的自导纳,增加一个改变量为
ji K
Y=K-1Y
+1Y=Y
ii K T KT T
节点j的自导纳,也增加一个改变量为
Y=1Y
+1-KY
=1Y
3.程序代码及说明
n=input('请输入节点数:
n=');
jj KT
K2 T
K2T
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);
%zeros(m,n)产生m×n的零矩阵,zeros(n)产生n×n的全0方阵。
fori=1:
n %节点数
ifX(i,2)~=0;%matlab关系符号:
>(大于),>=(大于等于),<(小于),<=(小于等于),
==(等于)~=(不等于)p=X(i,1);
Y(p,p)=1./X(i,2); %接地支路,“./”点除代表矩阵对应元素相除
end
end
fori=1:
nl %支路数
ifB1(i,6)==0 %折算到哪一侧的标志,0非标准变比在q侧,1非标准变比在在p侧p=B1(i,1);q=B1(i,2); %B1(i,1),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));
%计算节点p和q间的互导纳,B1(i,3)为支路阻抗,B1(i,5)为变压器支路变比k,B1(i,5)=1时不是变压器支路
Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
%对角元,节点q的自导纳,B1(i,4)为支路对地容抗Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;%对角元,节点p的自导纳
end
%输出导纳矩阵
disp('导纳矩阵Y=');disp(Y);
四、实验数据:
电力系统接线图如图3所示,根据图中数据可得B1和X矩阵:
P5=5
j0.015
1:
1.05 ③
0.8+j0.30
②1.05:
1
j0.03
⑤
U5=1.05
2+j1
j0.25 j0.25
0.04+j0.25
3.7+j1.3
①
U1=1.05
δ1=0°
j0.25 0.1+j0.35
j0.25 ④
1.6+j0.8
图3电力系统接线图
B1=[120.03i01.050;230.08+0.3i0.5i10;240.1+0.35i010;34
0.04+0.25i0.5i10;350.015i01.051];
X=[10;20;30;40;50];
五、实验过程
1.手算图3所示网络的导纳矩阵。
2.认真阅读源程序,建立matlab的M程序,输入节点导纳矩阵计算程序。
3.输入实验数据,用B1和X矩阵代替B1和X输入语句并调试程序。
B1矩阵
B1(i,1)为支路始端编号p;
B1(i,2)为支路末端编号q,且p B1(i,3)为支路阻抗; B1(i,4)为支路对地容抗,无填0; B1(i,5)为变压器支路变比,非变压器支路填1; B1(i,6)为变压器变比折算到哪一侧的标志,0非标准变比在q侧,1非标准变比在在p侧,非变压器支路填0; X矩阵 X(i,1)为节点编号; X(i,2)为该节点的接地支路阻抗,如果节点没有接地支路填”0”; 4.正确输入节点数和支路数进行计算,将计算结果与手算结果进行比较。 六、要求 认真书写实验报告并上交实验数据 实验二 支路追加法生成 节点阻抗矩阵的matlab程序设计 一、实验目的: 本实验通过对生成节点阻抗矩阵的编程与调试,获得复杂电力系统的节点阻抗矩阵,为采用阻抗法的电力系统的短路电流计算程序编制打下基础。 通过实验教学加深学生对节点阻抗矩阵计算方法的理解,掌握支路追加法生成节点阻抗矩阵的方法,熟悉MATLAB软件的使用方法,提高编制调试计算机程序的能力。 二、实验器材: 计算机、软件(已安装应用软件MATLAB等)、移动存储设备(学生自备,软盘、U盘等)。 三、实验内容: 1.节点阻抗矩阵的形成 目前求节点阻抗矩阵的方法,归纳起来主要有两种: 一种是导纳矩阵求逆,间接求出阻抗矩阵;另一种是用支路追加法,直接形成节点阻抗矩阵。 支路追加法形成节点阻抗矩阵可以从网络的接线图中直接形成,它从某一个与地相连的支路开始,以后每次追加一条支路,直至形成整个网络的节点阻抗矩阵。 该方法矩阵形成的规律性很强,易于理解和记忆。 用支路追加法形成节点阻抗矩阵时,所追加的支路可以分成4类,即追加接地树支,追加树支,追加接地连支,追加连支。 现以图1所示的含有3个独立节点的网络为例,别就上述4种情况阐述节点阻抗的形成。 先将节点的电压、电流关系表示如下: éU&1ù éZ11 Z12 Z13ùéI&ù ê ú ê úê1ú êU&2ú=êZ21 Z22 Z23úêI&2ú U ê&ú êZ Z Zúê&ú 3 ë31 32 33ûI ë û ë3û 图1原始网络 (1)追加接地支路(0,4) 如图2所示,由于网络增加了一个节点4,所以原网络矩阵增加一阶,但是由于新增节点4与原网络矩阵增加一阶,但由于新增节点4与原网络3个节点1、2、3无电气直接联系,所以原来 3各方程不收影响,只是新增一个方程U&4=zI&4,其中z是新增支路的阻抗。 这样节点方程修改为 式 (1),即节点阻抗矩阵原有的各元素均不变,新增的行、列元素均为零,只有新增的对角元素为z。 图2追加接地树支 éU&1ù éZ11 Z12 Z13 0ùéI&ù ê ú ê úê1ú êU&2ú=êZ21 Z22 Z23 0úêI&ú êU&ú êZ31 Z32 Z33 2 0úêI&ú (1) ê3ú ê úê3ú êU&4ú ë0 0 0 zûêI&ú ë û (2)追加树支(2,4) ë4û 图3追加树支 4 2 4 如图3所示。 由于网络增加了一个节点4,所以矩阵增加一阶,设节点4的注入电流为I&,从原网络看进去,节点2的注入电流变为I&+I&,而其他节点注入电流不变,而新增了一个方程: 4 U&4=U&2+zI&,所以节点方程改为 I I I I 111 U&1=Z I&+Z (I& +I&)+Z I&=Z I&+Z I&+Z I&+ZI& I I I I 12 2 4 13 3 112 12 2 13 3 12 4 U&2=Z & I+Z 211 22(I& +I&)+Z23& & I =Z 212 +Z22& +Z23& +Z22& 2 4 3 2 3 4 U&3=Z & I+Z 311 32(I& +I&)+Z33& & I I =Z 312 +Z32& +Z33& +Z32& 2 4 4 3 2 3 4 U&4=U&2+zI& & I+Z =Z 211 & I 223 +Z23& +(Z22+z)I& 4 3 写成矩阵形式 (2)。 éU&1ù éZ11 Z12 Z13 Z12 ùéI&ù ê ú ê úê1ú êU&2ú=êZ21 Z22 Z23 Z22 úêI&ú 2 êU&ú êZ Z Z Z úêI&ú (2) ê3ú ê31 32 33 32 úê3ú êU&ú ëZ21 Z22 Z23 Z22+zûêI&ú ë4û ë4û 由上式可以看出,追加树支时,矩阵增加一阶,新增的行、列元素分别等于树支所接的原网络节点2所对应的行、列元素,新对角元素等于树支所接的节点2的对角元素加上新增支路的阻抗值。 电力网络中包含有许多变压器。 在追加变压器支路时,一般用一个等值阻抗同一个理想变压器相串联的支路来表示,具体可分为追加树支和追加连支两种情况。 下面先讨论追加树支的情况,如图4所示,追加变压器树支(2,4),网络增加了一个节点4, 图4追加变压器树支 4 2 4 2 4 I 3 设其注入电流为I&,从原网络看,节点2的注入电流为(I&+KI&),所以有 I 1 U&1=Z11& +Z12 (I& +KI& )+Z13& 另外还有 U&2=Z21& I I 1 1 U&3=Z31& +Z22 +Z32 (I& 2 2 (I& +KI& 4 4 +KI& )+Z23& I I 3 3 )+Z33& 4 U&4=K(U&2+KzI&) U&=K{[ZI&+Z (I& +KI&)+ZI&]+KzI&} U&4整理后,有 4 211 22 2 4 233 4 éU&1ù éZ11 Z12 Z13 KZ12 ùéI&ù ê ú ê úê1ú êU&2ú=êZ21 Z22 Z23 KZ22 úêI&ú 2 êU&ú êZ31 Z32 Z33 KZ32 úêI&ú ê3ú ê úê3ú ê&ú KZ KZ KZ K2(Z +z)ê&ú ëU4û ë 21 22 23 22 ûëI4û 由此可见,追加变压器树支和追加普通树支支路相似,只是在新增行、列的元素分别乘以变比K,新对角元素乘以变比K2。 (3)追加接地连支(0,2) 图5追加接地连支 如图5所示,由于追加接地连支,网络的节点数没有变化,故矩阵的阶次不变。 对原网络来 说,节点2的注入电流变为I&¢=I&-I&,其他节点注入电流不变,则各节点电压方程变为 2 2 I I I I I 1 U&1=Z11& +Z12 (I& -I&)+Z13& =Z11& +Z12& +Z13& -Z12I& 2 I I I I I I I I I 3 1 2 3 1 3 1 2 3 U&2=Z21& +Z22 (I& -I&)+Z23& =Z21& +Z22& +Z23& -Z22I& 2 I I I 3 1 2 3 1 2 3 U&3=Z31& +Z32 (I& -I&)+Z33& =Z31& +Z32& +Z33& -Z32I& 2 I 1 0=-U&2+zI&=-Z21& 写成矩阵形式为 -Z22& -Z33& +(Z22 +z)I& éU&1ù éZ Z Z-Z ùéI&ù 1 ê ú ê 11 12 13 12úê ú êU&2ú=êZ21 Z22 Z23 -Z22 úêI&2ú êU&ú êZ Z Z -Z úêI&ú (3) ê3ú ê 31 32 33 32úê3ú êë0úû ë-Z21 -Z22 -Z23 Z22+zûêëI&úû 可见,矩阵可暂时增加一阶,原矩阵元素不变,新行、列元素分别等于该追加连支的非零节点所对应的行列元素的负值;新对角元素等于该点的自阻抗加上连支阻抗z。 形成了暂时增加一阶 的节点阻抗矩阵以后,用高斯消去法消去矩阵的暂增行与列,原矩阵的元素Zij为 Zi¢j =Zij -Zi2Z2j(i,j=1,2,3) Z (4) 44 (4)追加连支(2,3) 图6追加连支 如图6所示,由于追加连支,所以矩阵阶次不变,与采用追加接地连支同样的处理方法,矩阵可暂时增加一阶。 设连支电流为I&,由节点3流向节点2,从原网络看进去,节点2的注入电流 为(I&+I&),节点3的注入电流变为(I&-I&)。 则节点电压方程为 2 3 éU&1ù é Z Z Z (Z -Z) ùéI&ù ê ú ê 11 12 13 12 13 úê1ú êU&2ú=ê Z21 Z22 Z23 (Z22-Z32) úêI&2ú êU&ú ê Z31 Z32 Z33 (Z32-Z33) úêI&ú ê3ú ê úê3ú êë0úû ë(Z21-Z31) (Z22-Z32) (Z23-Z33) (Z22+Z33-Z23-Z32+z)ûêëI&úû (5) 可见,追加连支后,原矩阵的元素不变,暂时增加的行、列元素分别等于该追加连支的两个节点所对应的行元素之差和列元素之差;新增对角元素为这两个节点自阻抗之和减去相互间的互阻抗之和加上该连支阻抗。 形成了暂时增加一阶的节点阻抗矩阵之后,用高斯消去法消去暂增行、列,即得追加连支的节点阻抗矩阵,消元公式同式(4)。 图7追加变压器连支 如图7,追加变压器连支(2,3)。 由节点3流入理想变压器的电流为I& 则由变压器支路流 2 入节点2的电流为KI&,从原网络看节点2的注入电流变为(I&+KI&),节点3的注入电流变为 3 (I&+I&),则: 111 U&1=Z I&+Z (I& +KI&)+Z(I& -I&) 12 2 13 3 U&2=Z & I+Z 211 22(I& +KI&)+Z23 (I& -I&) 2 3 U&3=Z & I+Z 311 32(I& +KI&)+Z33 (I& -I&) 2 3 另有: U&3=K(U&2+zKI&) 2 3 即,KU&-U&+zK2I&=0 把U&2,U&3代入整理后得 (KZ21 -Z31 )I&+(KZ22 -Z32 )I& +(KZ23 -Z33 )I& +(K2Z -KZ32 -KZ23 +Z33 +K2z)I&=0 1 2 3 22 所以得 éU&1ù é ê ú ê Z11 Z12 Z13 KZ12-Z13 ùéI&ù 2 úê1ú êU&2ú=ê Z21 Z22 Z23 KZ22-Z23 úêI&ú U ê&ú ê ê3ú ê Z31 Z32 Z33 KZ32-Z33 úêI&ú úê3ú ê0ú (KZ -Z) (KZ -Z) (KZ -Z) (K2Z +KZ -KZ-KZ +K2z) ê&ú ë û ë 21 31 22 32 23 33 22 33 23 32 ûëIû 前面已经讨论的都是变压器的漏抗归算至低压侧。 如果变压器的漏抗归算至高压侧,则只需将变比变为1/K即可。 2.程序代码及说明 n=input('请输入节点数: n='); nl=input('请输入支路数: nl='); B=input('请输入由各支路参数形成的矩阵: B='); B=[012i010;024i010;122i010;0320i010;238i010;1 35i010]; m=0;Z=zeros(n);%zeros(n)产生n×n的全0方阵。 fork1=1: nl %支路数 p=B(k1,1);q=B(k1,2); %B(k1,1),B(k1,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 分析 实验 指导书
![提示](https://static.bingdoc.com/images/bang_tan.gif)
文档标签
- 电力系统分析实验
- 电力系统分析综合实验
- 电力系统仿真实验指导书
- 信号系统分析实验指导书
- 电力系统保护实验指导书
- 电力系统分析实验指导书
- 多元统计分析实验指导书
- 系统分析设计实验指导书
- 电力系统分析计算实验
- 电力系统实验指导
- 电渗析实验指导书
- 电力系统自动化实验指导书
- 电力系统过电压实验指导
- 电力系统实验指导全面
- mcgs电力系统实验指导书
- 电力系统分析实验讲义
- 电路分析实验指导书
- 分布式系统实验指导书
- 0301电力系统分析实验
- 电力系统分析春北交电力系统分析
- 电力系统分析西工大电力系统分析
- 电力系统分析Ⅰ电力系统分析同步
- 电力系统分析大工19电力系统
- 电力系统分析Ⅱ20电力系统分析
- 电力系统分析实验报告金科电力系统分析
- 电力系统分析电力系统分析阶段