Calculo de Ybus o Ybarra

by

 

Esta función permite el calculo de Ybus de un sistema eléctrico de potencia.

Calc_Ybarra.m
function Ybar=Calc_Ybarra()
%Calc_Ybarra  Calcula el Ybarra de un circuito de potencia
%             por el mtodo de simple inspeccin

%   Author(s): Kelman Belloso
%   Revision : 1.0 
%   Date     : 2013/02/08 20:41:30 
format compact

%*************************************************************************
%Solicitud de datos para el clculo del sistema
Num_nodos=input('Digite el nmero de nodos: ');
Num_ramas=factorial(Num_nodos)/(factorial(2)*factorial(Num_nodos-2));
Num_ramas=Num_nodos+Num_ramas;
Ybar=zeros(Num_nodos);

for m=1:1:Num_nodos
    fprintf('Digite la admitancia de los generadores y%d0',m);
    Yadm(1,m)=input(': ');
end
n2=Num_ramas-Num_nodos;
n=n2;
k=1;
while n~=0
    for m=1:1:Num_nodos
        if k>=m
        else
            fprintf('Digite la admitancia y%d%d',k,m);
            f=input(': ');
            Yadm2{n2-n+1}=['y',num2str(k),num2str(m),num2str(f)];
            n=n-1;
        end
    end
    k=k+1;
end
for m=1:1:Num_nodos
    for n=1:1:Num_nodos 
        if m==n
            Ybar(m,n)=Ybar(m,n)+Yadm(1,m);
            for k=1:1:length(Yadm2)
                if 1==strcmp(num2str(m),Yadm2{k}(2)) || 1==strcmp(num2str(m),Yadm2{k}(3))
                    Ybar(m,n)=Ybar(m,n)+eval(Yadm2{k}(4:end));
                end
            end
        else
            for k=1:1:length(Yadm2)
                if 1==strcmp(strcat(num2str(m),num2str(n)),Yadm2{k}(2:3)) || 1==strcmp(strcat(num2str(n),num2str(m)),Yadm2{k}(2:3))
                    Ybar(m,n)=Ybar(m,n)-eval(Yadm2{k}(4:end));
                end
            end
        end
    end
end
fprintf('\n');
Ybar;

Contact us