% Programa para Construccion de la Matriz Admitancia de Barra 
% Realizado por : Francisco M. Gonzalez-Longatt
% fglongatt@cantv.net, fglongatt@ieee.org
% Curso de Sistemas de Potencia II
% Caracas, Universidad Central de Venezuela
% Abril de 2003
% Propiedad de Francisco M. Gonzalez-Longatt
% Derechos de autor reservados, Copyright 2003
% Solo para ser empleados con fines de investigacion o academicos (C) 2008
clear all;
clc; op = 8;
while (op ~= 7)  % BEGIN WHILE MAIN
    % PRESENTACION Y MENU DE INICIO DEL PROGRAMA
    clc
    disp('     Programa para Construccion de la Matriz Admitancia de Barra')
    disp('           Realizado por : Francisco M. Gonzalez-Longatt')
    disp('                  Curso de Sistemas de Potencia II')
    disp('     Caracas, Universidad Central de Venezuela, Abril de 2003   ')
    disp(' --------------------------------------------------------------')
    disp('        ')
    disp('     MENU PRINCIPAL')
    disp(' --------------------------------------------------------------')
    disp('     [1] Ingresas Datos del Sistema')
    disp('     [2] Construir la Matriz Admitancia de Barra, Ybus')
    disp('     [3] Mostrar la Matriz Admitancia de Barra, Ybus')
    disp('     [7] Salir')
    disp('----------------------------------------------------')
    op= input(' Indique la operacion :');
    if op==1
        %   INTRODODUCIR DATOS DEL DIAGRAMA DE ADMITANCIAS
        disp('-----------------------------------------------------')
        disp(' Datos de Impedancias Reales de Lineas del Sistema')
        disp('-----------------------------------------------------')
        n_branch=input('     Numero de ramas :');
        n_bus=input('     Numero de Barras :');
        A=input('     Datos en [1] Impedancia [2] Admitancias');
        for s = 1:n_branch
            disp('-----------------------------------------------------')
            mm=[ 'Elemento ' int2str(s)];
            disp(mm)
            disp('-----------------------------------------------------') 
            bus_new=input(' Barra de Inicio :');  
            bus_old =input(' Barra Final :');   
            if A==1
                mm=[' Z[' int2str(bus_new) ',' int2str(bus_old) ']=' ];
                Zpq=input(mm); % Se introduce la Impedancia
                Ypq=1/Zpq;     % Se calcula la admitancia de cada elemento
            end
            if A==2
                mm=[' Y[' int2str(bus_new) ',' int2str(bus_old) ']=' ];
                Ypq=input(mm); % Se introduce la Impedancia
            end
            y(bus_new,bus_old)=Ypq; %asigna el valor a y
            y(bus_old,bus_new)=Ypq; %hace la matriz simetrica
        end % END FOR s
        sss=input(' Presione <ENTER>');
    end % END del IF OP ==1
    if op==2
        % Construir la Matriz Admitancia
        Y=zeros(n_bus);
        for k =1:n_branch
            for j = 1: n_branch
                suma=0; 
                if (j==k) %Calculo de la Admitancia propia Yjj
                    for m=1:n_branch % Efectua la busqueda de elementos conectados a la barra
                        if (m~=j)
                            suma=suma+y(j,m);
                        end
                    end % END FOR m
                    Y(j,j)= suma;
                end % Final del calculo de la admitancia propia
                if (j~=k) % Calculo de la admitancia de transferencia Yij
                    Y(k,j)=-y(k,j);
                end
            end % Final for j
        end % Final for k
    end % Final del IF OP = 2 
    if op ==3
        disp('-----------------------------------------------------')
        disp(' Matriz Admitancia de Barra ')
        disp('-----------------------------------------------------') 
        Ybus=Y
        sss=input(' Presione <ENTER>');
    end %Final del IF OP = 3
end % END MAIN WHILE
disp('   ')
disp('     Gracias por usar YbusMaker V 2.00')
disp('     Prof. Francisco M. Gonzalez-Longatt')
disp('     fglongatt@cantv.net, fglongatt@ieee.org')
clear sss A i j Y y  op  Ybus         bus_new  bus_old  k        m        mm       s        suma