n=input('请输入节点数:n=');
n1=input('请输入支路数:n1=');
isb=input('请输入平衡节点号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入支路参数:B1=');
B2=input('请输入节点参数:B2=');
X=input('节点号和对地参数:X=');
end
Times=1; %迭代次数
%创建节点导纳矩阵
Y=zeros(n);
for i=1:n1
if B1(i,6)==0 %不含变压器的支路
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-1/B1(i,3);
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);
Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);
else %含有变压器的支路
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-B1(i,5)/B1(i,3);
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3);
Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))/B1(i,3);
end
end
for i=1:n1
Y(i,i)=Y(i,i)+X(i,2); %计及补偿电容电纳
end
disp('导纳矩阵为:');
disp(Y); %显示导纳矩阵
%初始化 OrgS、DetaS
OrgS=zeros(2*n-2,1);
DetaS=zeros(2*n-2,1);
%创建 OrgS,用于存储初始功率参数
h=0;
j=0;
评论1
最新资源