讲解R语言、R开发解析、辅导R编程、R语言编程讲解

B1=b1+(b2.*b0)./(b2+b0);
c=a./B1;
B2=b0./(b2+b0);
B0=b2./(b2+b0);
i1=sum(c);
i2=-i1*B2(11);
i0=-i1*B0(11);
u1=n-i1*b1;
u2=-i2*b2;
u0=-i0*b0;
U1 = sparse( LineI, LineJ, u1(LineI)-u1(LineJ),22,22);I1=U1./Z1;
U2 = sparse( LineI, LineJ, u2(LineI)-u2(LineJ),22,22);I2=U2./Z2;
U0 = sparse( LineI, LineJ, u0(LineI)-u0(LineJ),22,22);
k=1;
for j=1:BranchNum
if LineX0(k)==0
U0( LineI(k), LineJ(k))=0;
end
k=k+1;
end
Z=zeros(22,22);
k=1;
for j=1:BusNum
if ismember(LineI(k), Linei)==1 || ismember(LineJ(k), Linei)==1
Z( LineI(k),LineJ(k))=LineX0(k)*1i ;
end
k=k+1;
end
Z=sparse(Z);
Z=Z+Z0;
m=1;
for i=1:BusNum
x=1;
for j=1:BusNum
if Z(x,m)==0
Z(x,m)=1;
end
x=x+1;
end
m=m+1;
end
I0=U0./Z;
Ia=I1+I2+I0;
Ib=p*I1+q*I2+I0;
Ic=q*I1+p*I2+I0;
Ua=u1+u2+u0;
Ub=p*u1+q*u2+u0;
Uc=q*u1+p*u2+u0;
If=3*i0;