올려주신 코드를 돌려보았을 때 함수 안에 있는 "nr"이 무엇인지 알 수 없어 말씀하신 에러를 재현할 수 없습니다. 도움을 받기 위해서는 말씀하신 입력 인수 부족에 관한 에러를 재현할 수 있도록 해주셔야 할 것 같습니다.
function[Y]=ybus1(zdata)
nl=zdata(:,1); nt=zdata(:,2); R=zdata(:,3); X=zdata(:,4);
nbr=length(zdata(:,1)); nbus=max(max(nl),max(nr));
Z=R+j*X %branch impedance
y=ones(nbr,1)./Z; %branch admittance
Y=zeros(nbus,nbus); %initialize Y to zero
for k=1:nbr; %formation of th off daiagonal elements
if nl(k)>0&nr(k)>0
Y(nl(k),nr(k))=Y(nl(k),nr(k))-y(k);
Y(nr(k),nl(k))=Y(nl(k),nr(k));
end
end
for n=1:nbus %formation of th off daiagonal elements
for k=1:nbr
if nl(k)==n|nr(k)==n
Y(n,n)=Y(n,n)+y(k);
else, end
end
end
end
% From To R X
z=[ 0 1 0 1
0 2 0 0.8
1 2 0 0.4
1 3 0 0.2
2 3 0 0.2
3 4 0 0.08];
Y=ybus1(z) %bus admittance matrix
Ibus=[-j*1.1; -j*1.25; 0; 0]; %vector of bus currents
Zbus=inv(Y) %bus impedance matrix
Vbus=Zbus*Ibus
