Conf MATLAB clear all clc disp(['OBTENCIN DE Ybus POR TRANSFORMACIONES SINGULARES']) f=input('Dame el numero de impedancias= '); disp(['f=' num2str(f)]); c=input('dame el numero de nodos= '); disp(['c=' num2str(c)]); disp(['FORMACIN DE LA MATRIZ DE INDICENCIA']) disp([' Si el elemento sale coloca -1']) disp([' Si el elemento entra coloca 1']) disp([' Si no hay elemento coloca 0']) for m= 1:f for n= 1:c sprintf('El elemento A[%d %d.]es', m,n) valor1=input(' ') disp(['valor2=' num2str(valor1)]); A(m,n)=valor1 n=n+1; end m= m+1; end disp(['La matriz de indicencia es']) A B=A'; disp(['La matriz transpuesta es']) B disp(['FORMACIN DE LA MATRIZ DE ADMITANCIAS PRIMITIVA']) P=eye(f); for k= 1:f j=k; sprintf('Dame el valor de la impedancia Z%d.', k) valor3=input(' ') valor4=1/valor3; disp(['Matriz Admitancias Primitiva=' num2str(valor4)]); P(k,j)=valor4 k= k+1; end disp(['La matriz de admitancias primitiva es']) P Y=B*P*A; disp(['La matriz Ybus es']) Y disp(['CREACIN DE LA MATRIZ Zbus']) I=eye(c) for k=1:c Simulacin de disturbios en Sistemas Elctricos de Potencia. for j=1:k if(k==j) diag=Y(k,j); for t=1:c Y(k,t)=Y(k,t)/diag; I(k,t)=I(k,t)/diag; end else if(Y(k,j)~=0) diag=Y(k,j); for t=1:c Y(k,t)=(-1*diag*Y(j,t))+Y(k,t); I(k,t)=(-1*diag*I(j,t))+I(k,t); end end end end end for k=1:c for j=(k+1):c if(Y(k,j)~=0) diag=Y(k,j); for t=1:c Y(k,t)=(-1*diag*Y(j,t))+Y(k,t); I(k,t)=(-1*diag*I(j,t))+I(k,t); end end end end disp(['LA MATRIZ ZBUS ES:']) Zbus=I; Zbus disp(['FUENTES DE CORRIENTE EN LA RED:']) disp(['Dame los elementos de la matriz de corrientes en coordenadasrectangulares']) disp(['Indica si la fuente entra (+) o si sale (-) al nodocorrespondiente']) for l=1:c sprintf('El elemento de corriente del nodo %d. es',l) valor6= input(' ') D(l,1) = valor6 end Vbus=Zbus*D; disp(['LA MATRIZ Vbus ES:']) Vbus