NEWTON RAPHSON ROTIINA GENERICA exemplo 09 barras 1

Rotina genérica pelo método de Newton-Raphson

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%APS2 de Analise de Sistemas de Potencia 1%%%%
%%%Realizar a Analise de Fluxo de Carga %%%%
%%% do Sistema proposto com 18 Barras %%%%
%%%Aluno: Sidnei Pereira 22/06/2010 %%%%
%%%Engenharia Eletrica 7° Periodo %%%%
%%% Método: Newton Raphson %%%%
%%%programa: fluxo9Barras.m versao: 1.0 %%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all clc %Dados dos ramos de conexão (linhas e trafos)[pu]
% | Da | Para | R | X | Bshunt |
% | Barra| Barra| | | | d_l = [ 1 4 0.0 0.0576 0.0; % Matriz dados_linhas
2 7 0.0 0.0625 0.0;
3 9 0.0 0.0586 0.0;
4 5 0.01 0.085 0.088;
4 6 0.017 0.0920 0.079;
5 7 0.032 0.161 0.153;
6 9 0.039 0.1700 0.179;
7 8 0.0085 0.072 0.0745;
8 9 0.0119 0.1008 0.1045;];
% Monta Matriz admitancia de barras Y fb = d_l(:,1); % Da barra numero tb = d_l(:,2); % Para a barra numero r = d_l(:,3); % Resistencia, R x = d_l(:,4); % Reatancia, X b = d_l(:,5); % Admitancia B Shunt z = r + i*x; % Matrix Z y = 1./z; % Inverte cada elemento de Z b = i*b; % Muda Bshunt para imaginario nbarra = max(max(fb),max(tb)); % identifica numero de barras nramos = length(fb); % numero de ramos
Ybarra = zeros(nbarra,nbarra); % Inicialisa YBarra
% Completa Elementos fora da diagonal for k=1:nramos
Ybarra(fb(k),tb(k)) = -y(k);
Ybarra(tb(k),fb(k)) = Ybarra(fb(k),tb(k)); end % Completa os elementos da Diagonal for m=1:nbarra for n=1:nramos if fb(n) == m | tb(n) == m
Ybarra(m,m) = Ybarra(m,m) + y(n) + b(n); end end end Ybarra; % Matrix Admitancia de Barra
%---------------------------------------------------------------------
%Dados fornecidos das barras
% Tipos: 1- Referencia; 2-PV(tensão controlada); 3-PQ(barra de carga)
% |Barra|Tipo| V | delta | PGi | QGi | PDi | QDi | Qmin | Qmax | dados_barra = [ 1 1 1.04 0 0.0 0 0 0 0 0;
2 2 1.025 0 1.63 0 0 0 0.05025 0.08375;
3 2 1.025 0 0.85 0 0 0 -0.0817 -0.1362;
4 3 1.0 0 0.0 0 0 0