Відео

Muy buen video. Solo quiero hacer una correcion a la linea que actualiza delta V, la cual deberia ser asi: Vmo((find(3==Nodos(:,2))))=Vmo((find(3==Nodos(:,2))))+Delta_TV(pv+pq+1:length(Delta_TV));

buen video gracias

aqui esta el codigo, lo copie tal cual: %% clc % Limpia pantlla clearvars % Borra variables close all % Cierra ventanas %% CARGAR DATOS % Nombre del archivo [libro, dir] = uigetfile('.xlsread','archivo excel'); % Letrero espera wai waitbar(0,'Porfavor espere'); % Tipo -> 1>>Slack; 2>>PV; 3>>PQ % Datos odos % Nodos = [Bus_1 tipo pg Og Pd Ssh Vm Va baseKr Vmax Vmin] Nodos = readtable(LIBRO,'Nodos'); waitbar(0.5,wai); % Datos lienas % Lineas = [Ni Nf R X Bsh/2] Lineas = readtable(LIBRO, 'Lineas'); waitbar(1,wai); close(wai); NR=size(Lineas,1); % Numero de lineas NB=size(Nodos,1); % Numero de Nodos % Serie de lineas ZLienas=Lineas (:,3) + (sqrt(-1)*Lineas (:,4)); % Valores inicales Vm_o=Nodos(:,8); % Magnitud tensiones Va_o=Nodos(:,9); % Angulo de tensiones % Toleracia Tol=0.001; %% Construcción de la matriz nodal for k = 1:NR Ybij = -1/Zlinea(k); i = Lineas(k,1); j = Lineas(k,2); Ybus(i,j)=Ybij; Ybus(j,i)=Ybij; end % Elemntos diagonal for i = 1:NB for j = 1:NR if i ==Lineas(j,1) Ybus(1,i) = Ybus(i,i) + (1/Zlinea(j))+(sqrt(-1)*ineas(j,5)); end if i==Lineas(j,2) Ybus(i,i) = Ybus(i,i)+(1/Zlinea(j))+(sqrt(-1)*Lineas(j,5)); end end Ybus(i,i) = Ybus(i,i) + Nodos(1,7); end figure spy(Ybus) title('Estructura de la Y_(bus)') %% Construye el sistema de ecuaciónes de balance nodal % Contruye variables for k = 1:NB syms (['V'num2str(k)]); Vm(k,1) = eval(['V',num2str(k)]); syms (['Tets'num2str(k)]); Va(k,1) = eval(['Teta',num2str(k)]); end i=sqrt(-1); % Contruye ecuaciones de blace nodal for k1=1:NB for k2=1:NB T(k1,k2)=Va(k2)-Va(k1); syms(['V' num2str(k2)]); V(k1,k2) = eval(['V', num2str(k2)]); end end G_real(Ybus); B_imag(Ybus); Dp=(Nodos(:,3)-Nodos(:,5))-((diag((G_.*V_)*cos(T)).*Vm)+(diag((B_.*V_)*sin(T)).*Vm)); Dq=(Nodos(:,4)-Nodos(:,6))-((diag((G_.*V_)*cos(T)).*Vm)+(diag((B_.*V_)*sin(T)).*Vm)); %% Selecciona para los Dp (Pv,Pq) y para los Dq for h=1:NB if (Nodos (h,2)==1) % Slack Dp(h)=(inf); Dq(h)=(inf); Vm(h)=inf; Va(h)=inf; end if (Nodos(h,2)==2) %PV Dq(h)=(inf); Vm(h)=inf; end end Dp(find(inf==Dp))=[]; Dq(find(inf==Dq))=[]; Vm(find(inf==Vm))=[]; Va(find(inf==Va))=[]; X=[Va;Vm]; % Variables de estado desconocidas Y=[Dp;Dq]; % Sistema de ecuaciones no lineales %% Calculo de Jacobiano [jacobian(X,Y)] for k1=1:length(Y) for k2=1:length(X) J(k1,k2)=diff(Y(k1),X(k2)); end end %% INICIA METODO ITERATIVO for h = 1:NB evalc(['V' num2str(h) '= Vm_o(h)']); evalc(['Teta' num2str(h) '= Vm_o(h)']); end %Calcula delta P y delta O Delta=eval(Y); ITE=1; while (ITE<=100)&(sum((abs(Delta)>=Tol))>0) J_o=eval(-1*J); % Evalua Jacobiano Delt_T_V=inv(J_o)*Delta; % Calcula delta_Teta y Delta_V pv=length(find(2==Nodos(:,2))); % Numero de nodos pv pq=length(find(3==Nodos(:,2))); % Numero de nodos pq % Actualiza Delta_Teta Va_o((find(2<=Nodos(:,2))))=Va_o((find(2<=Nodos(:,2))))+Delt_T_V(1:pv+pq); % Actualiza delta_v Vm_o((find(3==Nodos(:,2))))=Vm_o((find(3==Nodos(:,2))))+Delt_T_V(pv+pq+1:length(Delt_T_V)); for h = 1:NB evalc(['V' num2str(h) ' = Vm_o(h)']); evalc(['Teta' num2str(h) ' = Va_o(h)']); end Delta=eval(Y); UTE=ITE+1; end %% Muestra reaultados disp('________________________________________') fprintf('Número de iteraciones: %d ',ITE) disp('________________________________________') disp(' ') disp('________________________________________') for k=1:length(Vm_o) fprintf(' V_%d = %.3f', k,Vm_o(k)) if mod(k,2)==0 %si k es par fprintf(' ') else fprintf(' :') end end disp('________________________________________') disp(' ') disp('________________________________________') disp(' Ángulo de las tensiones nodales(°)') for k=1:length(Va_o) if Va_o(k)<0 % Angulo negativo fpritnf(' V_%d = %.2f',k,Va_o(k)*180/pi) else fprintf(' V_%d = %.3f',k,Va_o(k)*180/pi) end if mod(k,2)==0 %si k es par fprintf(' ') else fprintf(' ;') end end %% Graficas % Magnitud tensiones subplot(1,2,1) bar(Vm_o) grid on title('Tension nodal') xlabel('nodo') ylabel('pu') % Potencia neta inyectada subplot(1,2,2) bar([pink,Qi]) grid on legend('P_i','Q_i') title(' Potencia neta inyectada') xlabel('nodo') ylabel('pu')

A alguien le paso el codigo y la hoja excel ?

Podrías compartir el código y el archivo excel, por favor

bonita voz Luisa fernanda. lo as heredado de esa madre v maravillosa y talentosa gloria un saludo

Seria bueno q comparta el codigo :/

Buen video me podrías pasar el código por favor mi correo es fuentesjonathan66@gmail.com

enokdj62@gmail.com me puede pasar el código por favor

me gustaria tener el codigo me lo puedes pasar al correo danieldelgado20.drdt@gmail.com

Buen aporte .Podrías compartir el código edwin.lenovo95@gmail.com gracias

Hola, si es posible que me puedas enviar el código por fa: facundo.bravo1@gmail.com

esta genial, me gustaria tener el codigo de matlab, muchas gracias de antemano technicaloption_@hotmail.com

muy bueno el video y muy bien explicado será que me podrías pasar el flujo mas las hojas excel

Hola. Me parece que hay un error al calcular el Ybus. No se debería considerar las impedancias de los bancos de condensadores o reactores conectados a las barras. He visto programas en donde no se le considera. Si alguien podría aclarar el tema se lo agradecería mucho. Gracias!

Muchas gracias por el apoyo ingeniero

Ingeniero me puede enviar el código por favor kevinyed@hotmail.com

muy bueno me podria pasar el codigo porfavor mi correo dayanajimena333@gmail.com, gracias

Buen aporte ingeniero, me podría pasar el código por favor, es un muy buen aporte: ubeflores@gmail.com

esta genial, me gustaria tener el codigo de matlab, muchas gracias de antemano garymena777@gmail.com

Vaya, con todo lo explicado, ya me corrió el código, chevere capelaezhoy!

hola, soy estudiante de ing. Electrica, me ayudo mucho tu codigo y tu explicacion. sera posible que puedas compartir el codigo a mi correo jfmolinap@gmail.com... muchas gracias

GENIAL!! soy estudiante de electricidad me ayudaría mucho si me ayuda con el código de matlab, he tratado de seguir los pasos pero me he perdido... Muchas gracias de antemano ktsu98@gmail.com

Excelente aporte Ing. Carlos, favor si me puedes copiar el codigo al correo saul_fiee@hotmail.com

Excelente video, me puede pasar el código. daniel109035@gmail.com gracias

Excelente video. No le voy a pedir el código no se preocupe. Ando mas bien peleando para que por fsolve me logre andar. Saludos

genial me agrado. me podria mandar el codigo por fa. gracias.. mi correo: emersonchileno27@gmail.com

Muy bueno el tutorial, el mejor que he visto, me funciona muy bien el código. Gracias

Excelente video, Ing. Saldarriaga. Muchísimas gracias!!!

Que pasa con mi matriz jacobiana si sobre paso la potencia de mi carga en este caso con la reactiva?? Si es en un nodo PV??

Felicidades Ingeniero, es de muy buena ayuda... me puedes enviar el codigo por favor.. te lo agradesco anticipadamente: joserolando_9@hotmail.com

buen video amigo facilitanos el codigo si se puede

Felicidades Ingeniero, es de muy buena ayuda... me puedes enviar el codigo por favor.. te lo agradesco anticipadamente: misaeldanny10@gmail.com

Amigo podras pasarme el archivo xfavor

Donde esta la variable V_?

buen video ,,,

Saludos tu vídeo es muy interesante estoy iniciando el estudio de sistemas eléctricos de potencia seria posible que me envíes el código de matlab lo eh realizado siguiendo tus pasos pero al parecer estoy fallando en alguna parte gracias

Hola, saludos cordiales, tengo este problema en la ejecución, me podrían ayudar a solucionarlo. Warning: Matrix is singular to working precision. Delt_T_V=inv(J_o)*Delta;%Calcula delta_Teta y delta_V

no q q tal

sólo tiene que formular la hoja de cálculo en Excel, por favor enviarme la hoja de cálculo...

No tiene un ejemplo solo usando el Metodo de Gauss

Buenas me ayudariais a crecer ¿? estoy empezando en you tube y tengo varios tutoriales de programacion (C++, MatLab, Wireshark,... ) si quereis aprender a programar de forma facil y sencilla visita mi canal y suscribete saludos

ingeniero carlos muy buenas tardes un favor como se crea un proyecto en matlab??

Usted puede enviarme EL ARCHIVO MATLAB ?

Muy buen aporte.

buenisimo el video estoy viendo la materia de programacion y este video aclaro muchas dudas gracias por tomarse el tiempo de explicarnos

JAJAJAJAJAJA POBRES PENDEJOS QUIEREN EL CODIGO DE ESTE GRAN TRABAJO, CHINGUENLE PUTOS, NO ES NADA SENCILLO, QUE LES HACE PENSAR QUE SE LOS VAN A PASAR ASI DE FACIL, ES COMO REGALAR HORAS Y HORAS DE TRABAJO POR IMBECILES. ESO LO HACE UN VERDADERO INGENIERO Y NO SANGANOS COMO LO EXIGEN POR ESTE MEDIO.

Amigos existe una libreria muy poderosa para el calculo de flujos de potencia en MATLAB esta libreria se llama MATPOWER y se puede descargar gratuitamente de www.pserc.cornell.edu/matpower hay que descargar la libreria y hacerla correr desde matlab en la misma paguina existen manuales. Esta libreria es la que usa en el video para el calculo de flujos de potencia.

me fue de mucha ayuda, gracias

muy educativo , gracias por subir el video