%% Nombre del Estudiante %%
clear all
clc
%% Zprim= matriz primitiva de impedancias 
Zprim=[0.4013+1.4133*i 0.0953+0.8515*i 0.0953+0.7266*i 0.0953+0.7524*i;...
    0.0953+0.8515*i 0.4013+1.4133*i 0.0953+0.7802*i 0.0953+0.7865*i;...
    0.0953+0.7266*i 0.0953+0.7802*i 0.4013+1.4133*i 0.0953+0.7674*i;...
    0.0953+0.7524*i 0.0953+0.7865*i 0.0953+0.7674*i 0.6873+1.5465*i];
%% Definición de sub-matrices de la matriz primitiva de impedancias
Zij= Zprim(1:3,1:3);
Zin= Zprim(1:3,4);
Znj= Zprim(4,1:3);
Znn= Zprim(4,4);
%% Determinar la Matriz de impedancias de fase Zabc en Ohm/milla

%% Matrices generalizadas a,b,c y d (Se desprecia la admitancia Shunt)
a=eye(3);
b=ZabcT;
c=zeros(3,3);
d=a;
%% Datos de las tensiones y corrientes en la subestación (Nodo n)
t=(cosd(120)+i*sind(120));
Vabcn=7200*[1;conj(t);t];
Iabcn=[258*(cosd(-20)+i*sind(-20));288*(cosd(-147)+i*sind(-147));324*(cosd(86)+i*sind(86))];
%% Determinar las tensiones y corrientes la carga (Nodo m)
Vabcm=
Iabcm=
%% Determinar las matrices generalizadas para el regulador
ar=
dr=
%% Determinar las tensiones y corrientes en el secundario del regulador
Vrabc=
Irabc=
%% Tensiones en la carga
Vload_abc=