clear
N1=input('请输入节点数:N1=');
L1=input('请输入支路数:L1=');
B1=input('请输入支路信息:B1=');
%每行是一条支路
%第一列是支路的一个端点I,变压器所在端用负号表示,励磁端用零表示
%第二列是支路的另一个端点J,变压器所在端端用负号表示,励磁端用零表示
%第三列是支路的电阻R
%第四列是支路的电抗X
%第五列是支路的对地电纳或变压器变比k, k在j侧为正
Y=zeros(N1);
G=zeros(N1);
B=zeros(N1);
for m1=1:N1
I=B1(m1,1);
J=B1(m1,2);
R=B1(m1,3);
X=B1(m1,4);
k=B1(m1,5);
if I*J==0 %判断是否是励磁支路
if I==0
G(J,J)=G(J,J)+R;
B(J,J)=B(J,J)+X;
else
G(I,I)=G(I,I)+R;
B(I,I)=B(I,I)+X;
end
elseif I*J>0 %判断是普通支路
B(I,I)=B(I,I)+k;
B(J,J)=B(J,J)+k;
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余1页未读,立即下载