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

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

编程软件：matlab R2010b

程序说明：

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

2 3 0.08 0.30 0.25 1

4 2 0 0.01

5 0 1.05

5 3 0 0.03 0 1.05

1 2 0.04 0.25 0.25 1

1 3 0.1 0.35 0 1

输出数据：

Y=

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

2.6415i 0 0 -0.6240 +

3.9002i 1.4539 -66.9808i -0.8299 + 3.1120i 0 +63.4921i 0 -0.7547 + 2.6415i -0.8299 + 3.1120i 1.5846 -35.7379i 0 0 +31.7460i 0 0 +63.4921i 0 0 -66.6667i 0 0 0 0 +31.7460i 0 0 -33.3333i 经手算校验，程序结果准确。