节点导纳矩阵的建立

2

3

y

y

如上图所示的简单电力系统中，网络各元件参数的标幺值如下：

z12=0.10+j0.40 y120=y210=j0.01528 z13=j0.3,k=1.1 z14=0.12+j0.50 y140=y410=j0.01920 z24=0.08+j0.40 y240=y420=j0.01413

系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点。

节点导纳矩阵的运行程序如下：

clc

Clear

disp('网络各元件参数用标幺值表示');

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

n1=input('请输入支路数:n1=');

l=input('请输入PQ节点的个数=');

for m=1:l

c(m)=input(['请输入第',num2str(m),'个PQ节点的节点号为：']);

end

t=input('请输入PV节点的个数=');

for m=1:t

c(m)=input(['请输入第',num2str(m),'个PV节点的节点号为：']);

end

b=input('请输入平衡节点号:b=');

%%由支路参数形成矩阵B1

disp('各支路连接情况：')

i=1;

for m=1:n1

syms Y N

p=input(['第',num2str(m),'条支路的起始节点']);

q=input(['第',num2str(m),'条支路的终止节点']);

mn=input(['第',num2str(m),'条支路是否有变压器（请输入‘Y’或‘N’）']);

y=0;k=1;

if mn=='Y';

k=input('请输入变压器变比(标幺值):');

z=input(['请输入第',num2str(m),'条支路的线路阻抗']);

else

z=input(['请输入第',num2str(m),'条支路的线路阻抗：']);

y=input(['请输入第',num2str(m),'条支路线路的对地阻抗：']);

end

B1(i,1)=p;B1(i,2)=q;B1(i,3)=z;B1(i,4)=y;B1(i,5)=1/k;

i=i+1;

end

disp('由支路参数形成的矩阵B1')

B1

%求节点导纳矩阵

Y=zeros(N0);

e=zeros(1,N0);

f=zeros(1,N0);

for i=1:n1

p=B1(i,1);q=B1(i,2);

Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));

Y(q,p)=Y(p,q);

Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;

Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;

End

disp('导纳矩阵Y=');

disp(Y)

程序运行结果：