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

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

编程软件：matlab R2010b

程序说明：

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

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

3.标幺值等值电路中，如果变比为1:1,则默认为线路，因此，变压器的非标准变比不可以是1:1。

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

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

程序代码:

a=load('');%从’’中读入数据

[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)<%如果为线路

y(hnode,hnode)=y(hnode,hnode)+b;

y(enode,enode)=y(enode,enode)+b;

end

if (abs(k-1)>%如果为变压器

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 变比

2 3 1

4 2 0 0

5 3 0 0

1 2 1

1 3 0 1

输出数据：

Y=

- + + 0 0

+ + 0 + 0

+ + 0 0 +

0 0 + 0 0 0 0 0 0 + 0 0

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