clear R=[0.715,1.235,1.33,1.242,0.966,0.552,0.69,0.828]; S=[0.05,0.030,0.100,0.030,0.063,0.1]; S1=35.460+26.410j; P=[0.0015,0.0008,0.0023,0.0008,0.0014,0.002]; P_0=[0.00019,0.00015,0.00038,0.00015,0.00022,0.00032]; t=555; K=1.08; U=10; N=6; S_H=0; P_0H=0; P_H=0; Irms=K*sqrt(real(S1)^2+imag(S1)^2) /(t*U*sqrt(3)); for m=1:N; S_H=S_H+S(m); end for m=1:N; P_H=P_H+P(m); end R_T=[(U/S_H)^2]*P_H; R_L=[R(4)*S(6)^2+R(8)*S(3)^2+R(3)*[S(3)+S(6)]^2+R(7)*[S(4)+S(5)]^2+R(2)*[S(3)+S(4)+S(5)+S(6)] ^2+R(5)*S(2)^2+R(6)*S(1)^2+R(1)*S_H^2]/S_H^2; for m=1:N; P_0H=P_0H+P_0(m); end R_T R_L W=[3*Irms^2*(R_T+R_L)+P_0H]*t W0=W*100/real(S1)