clear;
%各节点参数:节点编号,类型,电压幅值,电压相位,注入有功,注入无功
%类型:1=PQ 节点,2=PV 节点,3=平衡节点
%本程序中将最后一个节点设为平衡节点
R_1=[1 1 1.0 0 0.2 0.2j;
2 1 1.0 0 -0.45 -0.15j;
3 1 1.0 0 -0.45 -0.05j;
4 1 1.0 0 -0.6 -0.1j;
5 3 1.0 0 0 0];
%支路号 首端节点 末端节点 支路导纳
R_2=[1 5 2 1.25-3.75j;
2 2 3 10.00-30.00j;
3 3 4 1.25-3.75j;
4 1 4 2.50-7.50j;
5 1 5 5.00-15.00j;
6 1 2 1.667-5.00j];
n=5;L=6;%需要改变的到此为止
i=0;j=0;a=0;precision=1;k=0;
Y=zeros(n,n);u=zeros(1,n);delt=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);
G=[];B=[];PP=[];uu=[];U=[];dp=[];dq=[];
for a=1:L
i=R_2(a,2);
j=R_2(a,3);
Y(i,j)=-R_2(a,4);
Y(j,i)=Y(i,j);
end
for a=1:n
for b=1:n
if a~=b
Y(a,a)=Y(a,a)+Y(a,b);
end
end
end
for i=1:n
for j=1:n
if i==j
Y(i,j)=-Y(i,j);
end
end
end
Y %形成导纳矩阵
for i=1:n
for j=1:n
G(i,j)=real(Y(i,j));
评论0
最新资源