disp('VALORES BASE')
Sb=100e6
Vb=220e3
disp('IMPEDANCIA BASE')
Zb=((Vb^2)/Sb) % PARA HALLAR LA IMPEDANCIA BASE
disp('VALORES POR UNIDAD DE LAS IMPEDANCIA SERIE')
Z12=3.2+15.5j;
Z13=2.5+12.4j;
Z14=2.7+8.3j;
Z25=2.9+14.5j;
Z35=4.2+20.7j;
Z45=1.7+8.3j;
Zpu12=(Z12/Zb)
Zpu13=(Z13/Zb)
Zpu14=(Z14/Zb)
Zpu25=(Z25/Zb)
Zpu35=(Z35/Zb)
Zpu45=(Z45/Zb)
% UNA VEZ OBTENIDA LA MATRIZ ADMITANCIA "Y" (no pongo la mia porque cambia dependiendo % del unifilar)
Y % Matriz
% ---- CALCULO DE P y Q----
% def de parametros:
disp('CALCULAR P y Q')
V=[1,1,1,1,1]; % PERFIL PLANO DE TENSIONES
d=[0,0,0,0,0]; % PERFIL PLANO DE TENSIONES son los desfases
N=5; % porque mi matriz es 5X5
P=zeros(1,N); % para inicializar en cero
disp('P:')
for ii=1:N
for kk=1:N
P(1,ii)=P(1,ii)+ (abs(V(ii))*abs(V(kk))*(real(Y(ii,kk))*cos(d(ii)-d(kk))+imag(Y(ii,kk))*sin(d(ii)-d(kk))) );
end
end
P % para imprimir P fuera del bucle.
disp('Q:')
Q=zeros(1,N);
for ii=1:n
for kk=1:n
Q(1,ii)=Q(1,ii)+ (abs(V(ii))*abs(V(kk))*(real(Y(ii,kk))*sin(d(ii)-d(kk))-imag(Y(ii,kk))*cos(d(ii)-d(kk))) );
end
end
Q