|
24.10.2012., 19:28 | #1 |
Registered User
Datum registracije: Nov 2009
Lokacija: Novi Sad
Postovi: 5
|
Matlab
Pomoc oko koda iz matlaba: Code:
for n = 1:nbus busprt = 0; for L = 1:nbr; if busprt == 0 fprintf(' \n'), fprintf('%6g', n), fprintf(' %9.3f', P(n)*basemva) fprintf('%9.3f', Q(n)*basemva), fprintf('%9.3f\n', abs(S(n)*basemva)) busprt = 1; else, end if nl(L)==n k = nr(L); %In - Nominalna struja In = (V(n) - a(L)*V(k))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(n); %y= ones(nbr,1)./Z; Ik = (V(k) - V(n)/a(L))*y(L) + Bc(L)*V(k); %Admitansa grane L Snk = V(n)*conj(In)*basemva; Skn = V(k)*conj(Ik)*basemva; SL = Snk + Skn; SLT = SLT + SL; elseif nr(L)==n k = nl(L); In = (V(n) - V(k)/a(L))*y(L) + Bc(L)*V(n); Ik = (V(k) - a(L)*V(n))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(k); Snk = V(n)*conj(In)*basemva; Skn = V(k)*conj(Ik)*basemva; SL = Snk + Skn; SLT = SLT + SL; else, end if nl(L)==n | nr(L)==n fprintf('%12g', k), fprintf('%9.3f', real(Snk)), fprintf('%9.3f', imag(Snk)) fprintf('%9.3f', abs(Snk)), fprintf('%9.3f', real(SL)), if nl(L) ==n & a(L) ~= 1 fprintf('%9.3f', imag(SL)), fprintf('%9.3f\n', a(L)) else, fprintf('%9.3f\n', imag(SL)) end else, end end end SLT = SLT/2; fprintf(' \n'), fprintf(' Total loss ') fprintf('%9.3f', real(SLT)), fprintf('%9.3f\n', imag(SLT)) clear Ik In SL SLT Skn Snk Code:
if nl(L)==n k = nr(L); %In - Nominalna struja In = (V(n) - a(L)*V(k))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(n); %y= ones(nbr,1)./Z; Ik = (V(k) - V(n)/a(L))*y(L) + Bc(L)*V(k); %Admitansa grane L Snk = V(n)*conj(In)*basemva; Skn = V(k)*conj(Ik)*basemva; SL = Snk + Skn; SLT = SLT + SL; elseif nr(L)==n k = nl(L); In = (V(n) - V(k)/a(L))*y(L) + Bc(L)*V(n); Ik = (V(k) - a(L)*V(n))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(k); Snk = V(n)*conj(In)*basemva; Skn = V(k)*conj(Ik)*basemva; SL = Snk + Skn; SLT = SLT + SL; Zadnje izmijenjeno od: Neshone. 24.10.2012. u 19:53. |
|
|
Oglas
|
Oglasni prostor
|
Uređivanje | |
|
|