电力系统分析作业电网节点导纳矩阵的计算机形成_精品文档.doc

电力系统分析作业电网节点导纳矩阵的计算机形成_精品文档.doc

《电力系统分析作业电网节点导纳矩阵的计算机形成_精品文档.doc》由会员分享，可在线阅读，更多相关《电力系统分析作业电网节点导纳矩阵的计算机形成_精品文档.doc（2页珍藏版）》请在冰豆网上搜索。

电力系统分析作业——电网节点导纳矩阵的计算机形成

编程软件：

matlabR2010b

程序说明：

1.如果已经输入i-j支路的信息，则不可再输入j-i支路的信息。

2.变压器支路的第一个节点编号默认为变压器一次侧，即变压器的等值电路中的阻抗归算侧，亦即变压器非标准变比的1:

k中的‘1’。

3.标幺值等值电路中，如果变比为1:

1,则默认为线路，因此，变压器的非标准变比不可以是1:

1。

5.如果变压器支路也有导纳B不为零，则说明此导纳就是励磁导纳，与线路的导纳B/2不同含义，只算作变压器原边的自导纳。

4.由于程序执行的是复数运算，所以即使实部为零时，也会输出实部‘0’。

程序代码:

a=load（'data.txt'）;%从’data.txt’中读入数据

[m,n]=size（a）;

w=1i;

u=1;

while（u<=m）

hnode=a（u,1）;

enode=a（u,2）;

z=a（u,3）+a（u,4）*w;

b=a（u,5）*w;

k=a（u,6）;

y（hnode,enode）=-1/（k*z）;

y（enode,hnode）=-1/（k*z）;

y（hnode,hnode）=y（hnode,hnode）+1/（k*z）+（k-1）/（k*z）;

y（enode,enode）=y（enode,enode）+1/（k*z）+（1-k）/（k*k*z）;

if（abs（k-1）<0.0001）%如果为线路

y（hnode,hnode）=y（hnode,hnode）+b;

y（enode,enode）=y（enode,enode）+b;

end

if（abs（k-1）>0.0001）%如果为变压器

y（hnode,hnode）=y（hnode,hnode）-b;

end

u=u+1;

end

[m,n]=size（y）;

disp（‘Y=’）;

disp（y（1:

m,1:

n））;

clear;

算例

输入数据：

首端编号末端编号电阻电抗电纳/2变比

230.080.300.251

4200.01501.05

5300.0301.05

120.040.250.251

130.10.3501

输出数据：

Y=

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

-0.6240+3.9002i1.4539-66.9808i-0.8299+3.1120i0+63.4921i0

-0.7547+2.6415i-0.8299+3.1120i1.5846-35.7379i00+31.7460i

00+63.4921i00-66.6667i0

000+31.7460i00-33.3333i

经手算校验，程序结果准确。