Método Newton-Rapshon no Matlab para um sistema de Quatro barras

%-------Universidade Federal de Sergipe ------------------%-------Disciplina: Sistemas Elétricos de Potência--------- Elabore uma rotina em Matlab para resolução do fluxo de carga pelo método de Newton-Raphson para o sistema abaixo. Aluno: %-------Dados de entrada do Sistema elétrico com Quatro barras----------------------YBarra = [3-12i, (-2+8i), (-1+4i), 0;(-2+8i), (3.666-14.64i),-0.666+2.64i,-1+4i;-1+4i,- 0.666+2.64i,3.66614.64i,-2+8i;0 ,-1+4i , -2+8i,3-12i]; G=[3,-2,-1,0;-2,3.666,-0.666,-1;-1,-0.666,3.666,-2;0,-1,-2,3]; %Matriz condutância %Matriz Suspectância B=[-12,8,4,0;8,-14.64,2.64,4;4,2.64,-14.64,8;0,4,8,-12]; nBarras=4; % número de barras V=[1.06;1;1;1]; % Tensão nas barras W=[0;0;0;0]; % fase das Tensões nas barras % A ordem é 2, 3, 4 em relação as barras Pprog=[-0.5;-0.4;-0.3]; Qprog=[-0.2;-0.3;-0.1]; % A ordem é 2, 3, 4 em relação as barras MagYbar=[12.37, 8.25, 4.123, 0; 8.25, 15.09, 2.72, 4.123; 4.123, 2.72, 15.09, 8.25; 0, 4.123, 8.25, 12.37]; FasYbar=(pi/180)*[-75.9, 104, 104, 0; 104, -75.9, 104, 104;104, 104, -75.9, 104; 0, 104, 104, -75.9]; % Teste de Convergência for L= 2:4 Acumulator=0; for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L,C)*cos(FasYbar(L,C)+W(C)-W(L))+ Acumulator; end end Pcal(L) = (V(L)^2)*G(L,L) + Acumulator; end DPQ(1)= Pprog(1)-Pcal(2); % valores de potência que serão adicionados DPQ(2)= Pprog(2)-Pcal(3); % valores de potência que serão adicionados DPQ(3)= Pprog(3)-Pcal(4); % valores de potência que serão adicionados for L= 2:4 Acumulator=0; for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L,C)*sin(FasYbar(L,C)+W(C)-W(L))+ Acumulator; end end Qcal(L) = (-1)*(V(L)^2)*B(L,L) - Acumulator; end DPQ(4)= Qprog(1)-Qcal(2); % valores de potência que serão adicionados DPQ(5)= Qprog(2)-Qcal(3); % valores de potência que serão adicionados DPQ(6)= Qprog(3)-Qcal(4); % valores de potência que serão adicionados fprintf('\n\niteração V1 th1 V2 th2 V3 th3 V4 th4 \n'); fprintf('%6i %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f\n',... 0,V(1),W(1),V(2),W(2),V(3),W(3),V(4),W(4)); i)-2*(V(i)^2)*B(i.a2)= N(a2+1.l)+N(k.j)+W(j)-W(i)). for k=2:4 for l=1:4 if l~=k N(k. end end end %Calculando LH matriz que contém as derivadas parciais da Potência ativa em relação as Tensões for i= 2:4 for j=2:4 if i~=j LH(i.2)=0.2)=0.i) = N(i.k)= -N(k.j)= -(abs(V(j)*V(i))*MagYbar(i.4)=0.i)+2*(V(i)^2)*G(i. N(2. end end end %Calculando UH matriz que contém as derivadas parciais da Potência reativa em relação as Tensões for i= 2:4 for j=2:4 if i~=j UH(i.M(2.4)=0. for k=2:4 for l=1:4 if l~=k M(k.j).k)= -M(k.d+1).j))*sin(FasYbar(i. end end .1)=0. N(4. N(3. M(3.d2+1).k). end end end M(1.i). end end end N(1.i).%Primeira Iteração %|-----------| %| M | LH | %|-----------| %| N | UH | %|-----------| %Calculando M a matriz que contém as derivadas parciais da Potência ativa em relação as fases for i=2:4 for j=1:4 if i~=j M(i. else LH(i. end end end for a2=1:3 for d2=1:3 Jacob(d2+3. else UH(i.l)+M(k. end end end %Calculando N.a)= M(a+1.j)=-N(i. M(4. a matriz que contém as derivadas parciais da Potência reativa em relação as Fases for i=2:4 for j=1:4 if i~=j N(i.k).j))*cos(FasYbar(i.1)=0. end end end %Obtenção da Matriz do Jacobiano for a=1:6 for d=1:6 if (d<=3 && a<=3) Jacob(d.j)= -abs((V(i)*V(j))*MagYbar(i.3)=0.j).3)=0.j)= M(i.j)+W(j)-W(i)).i)= -M(i. C)*sin(FasYbar(L.4)..a4+3)= UH(a4+1. % valores de potência para serem adicionados for L= 2:4 Acumulator=0. % Matriz inversa do Jacobiano TV=Jacob2*DPQ'. LH(2.V(3).3).Acumulator.L) + Acumulator. %Segunda Iteração %|-----------| %| M | LH | %|-----------| %| N | UH | %|-----------| %Calculando M a matriz que contém as derivadas parciais em relação as fases for i=2:4 for j=1:4 if i~=j M(i.j)+W(j)-W(i)). end end Pcal(L) = (V(L)^2)*G(L.6f %10.M(2.2). LH(4.L) . end end end M(1.6f %10.2)=0. end end Qcal(L) = (-1)*(V(L)^2)*B(L.5)= Jacob(3. end DPQ(4)= Qprog(1)-Qcal(2).5)= Jacob(1. 3 e 4 W(2)=W(2)+TV(1).4)= Jacob(2.C)+W(C)-W(L))+ Acumulator. 3 e 4 e tensões 2.d4+1).6)= LH(2. %atualização dos valores das fases W(3)=W(3)+TV(2). V(4)=V(4)*( 1 + TV(6)). % Vetor de 6 posições 3 fases e 3 tensões fase 2. % Teste de Convergência para Primeira Iteração for L= 2:4 Acumulator=0.W(1). %atualização dos valores das fases W(4)=W(4)+TV(3).6f %10.6f\n'. for a4=1:3 for d4=1:3 Jacob(d4+3.V(2).1)=0.j)= -(abs(V(j)*V(i))*MagYbar(i. LH(4.4).5)= Jacob(2.V(1). M(3. % valores de potência para serem adicionados DPQ(6)= Qprog(3)-Qcal(4).4). LH(2.6)= Jacob(2..2).3). end end Jacob2=inv(Jacob).W(2). % valores de potência para serem adicionados fprintf('%6i %10.2).4)=0. LH(4.W(3). for k=2:4 for l=1:4 if l~=k .Jacob(1.C)*cos(FasYbar(L.W(4)). DPQ(2)= Pprog(2)-Pcal(3).j))*sin(FasYbar(i. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L. % valores de potência para serem adicionados DPQ(5)= Qprog(2)-Qcal(3).3). for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L. M(4. end % valores de potência para serem adicionados DPQ(1)= Pprog(1)-Pcal(2).. LH(3.6f %10.6f %10. % valores de potência para serem adicionados DPQ(3)= Pprog(3)-Pcal(4). LH(3.4)= Jacob(1. %atualização dos valores das fases V(2)=V(2)*( 1 + TV(4)).6f %10.3)=0.4)= Jacob(3. V(3)=V(4)*( 1 + TV(5)).6f %10.V(4). LH(3.C)+W(C)-W(L))+ Acumulator.6)= Jacob(3. 1. k)= -N(k.5)= Jacob(1. LH(2.2).d+1).d4+1). end end end %Calculando LH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j LH(i. N(2. LH(3. end end end for a2=1:3 for d2=1:3 Jacob(d2+3.j)=-N(i. end end end %Obtenção da Matriz do Jacobiano for a=1:6 for d=1:6 if (d<=3 && a<=3) Jacob(d. end end Jacob(1. LH(4.i).k).M(k.i)-2*(V(i)^2)*B(i.d2+1).l)+N(k.4)= Jacob(1.i). end end end %Calculando N matriz que contém as derivadas parciais em relação as Tensões for i=2:4 for j=1:4 if i~=j N(i.6)= Jacob(2.2).4)= Jacob(2.k).3).i) = N(i.2)=0.5)= Jacob(2. LH(2. end end end N(1. end end . LH(3.2). else UH(i.4). LH(4.l)+M(k. end end end %Calculando UH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j UH(i.j)= -abs((V(i)*V(j))*MagYbar(i.3).a2)= N(a2+1. for a4=1:3 for d4=1:3 Jacob(d4+3.5)= Jacob(3.a)= M(a+1.j)+W(j)-W(i)). LH(4.i)= -M(i. else LH(i.1)=0.6)= LH(2.k)= -M(k. N(4.j))*cos(FasYbar(i.4).3)=0. LH(3.6)= Jacob(3.i)+2*(V(i)^2)*G(i.4)= Jacob(3. for k=2:4 for l=1:4 if l~=k N(k.4)=0.j). N(3.a4+3)= UH(a4+1.3).j)= M(i.j).4). k). W(2)=W(2)+TV(1).j)+W(j)-W(i)). % valores de potência para serem adicionados for L= 2:4 Acumulator=0.L) + Acumulator.V(3).W(1). end DPQ(4)= Qprog(1)-Qcal(2). W(4)=W(4)+TV(3). %Terceira Iteração %|-----------| %| M | LH | %|-----------| %| N | UH | %|-----------| %Calculando M a matriz que contém as derivadas parciais em relação as fases for i=2:4 for j=1:4 if i~=j M(i.. M(3.W(2). % valores de potência para serem adicionados %10. % valores de potência para serem adicionados DPQ(6)= Qprog(3)-Qcal(4).6f %10.M(2.6f %10. M(4. end end Qcal(L) = (-1)*(V(L)^2)*B(L.j))*sin(FasYbar(i. end DPQ(1)= Pprog(1)-Pcal(2). N(2.4)=0.4)=0. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L.6f %10.2)=0. W(3)=W(3)+TV(2).V(4).1)=0. % valores de potência para serem adicionados DPQ(3)= Pprog(3)-Pcal(4).W(3). N(3.j))*cos(FasYbar(i.L) .Jacob.l)+M(k. end end end M(1.6f\n'.k)= -M(k. TV=Jacob2*DPQ'.1)=0.V(1).3 e 4 e tensões 2. % Até aqui está OK Vetor de 6 posições 3 fases e 3 tensões fase2. for k=2:4 for l=1:4 if l~=k .3)=0.6f %10.C)*cos(FasYbar(L.6f %10. % TV. fprintf('%6i 2.V(2).j)= -(abs(V(j)*V(i))*MagYbar(i.3 e 4 %atualização dos valores das fases %atualização dos valores das fases %atualização dos valores das fases V(2)=V(2)*( 1 + TV(4)). N(4. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L.C)*sin(FasYbar(L.C)+W(C)-W(L))+ Acumulator. % Teste de Convergência p/ Segunda Iteração for L= 2:4 Acumulator=0. V(3)=V(4)*( 1 + TV(5)).Acumulator. Jacob2=inv(Jacob). V(4)=V(4)*( 1 + TV(6)).3)=0. end end Pcal(L) = (V(L)^2)*G(L.6f %10.2)=0. % valores de potência para serem adicionados DPQ(5)= Qprog(2)-Qcal(3).j)= -abs((V(i)*V(j))*MagYbar(i. % valores de potência para serem adicionados DPQ(2)= Pprog(2)-Pcal(3).C)+W(C)-W(L))+ Acumulator. end end end N(1. for k=2:4 for l=1:4 if l~=k M(k..6f %10.W(4)).j)+W(j)-W(i)).. end end end %Calculando N matriz que contém as derivadas parciais em relação as Tensões for i=2:4 for j=1:4 if i~=j N(i. a2)= N(a2+1. %atualização dos valores das fases V(2)=V(2)*( 1 + TV(4)).4).5)= Jacob(1. LH(3. end end Jacob(1.C)+W(C)-W(L))+ Acumulator.4)= Jacob(3. 3 e 4 e tensões 2.d2+1).k).6)= Jacob(3. %atualização dos valores das fases W(3)=W(3)+TV(2).a)= M(a+1.2). %atualização dos valores das fases W(4)=W(4)+TV(3).3).i)= -M(i. else UH(i. 3 e 4 TV. LH(4.5)= Jacob(3. % Vetor de 6 posições 3 fases e 3 tensões fase 2. end end end for a2=1:3 for d2=1:3 Jacob(d2+3. else LH(i.4).3).4)= Jacob(2. end end Jacob2=inv(Jacob).4)= Jacob(1.N(k.j). V(3)=V(4)*( 1 + TV(5)). % Teste de Convergência para Terceira Iteração for L= 2:4 Acumulator=0. LH(3. W(2)=W(2)+TV(1).i)+2*(V(i)^2)*G(i.i)-2*(V(i)^2)*B(i.5)= Jacob(2.2). LH(4.6)= Jacob(2. end end end %Calculando UH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j UH(i.i) = N(i. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L.j)=-N(i. LH(2. end end end %Obtenção da Matriz do Jacobiano for a=1:6 for d=1:6 if (d<=3 && a<=3) Jacob(d. LH(3.l)+N(k. end end end %Calculando LH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j LH(i. for a4=1:3 for d4=1:3 Jacob(d4+3.d+1).a4+3)= UH(a4+1.k)= -N(k. LH(2.i). LH(4. % Inversa do Jacobiano TV=Jacob2*DPQ'.C)*cos(FasYbar(L.i).3). V(4)=V(4)*( 1 + TV(6)).2).j).d4+1).j)= M(i. end .4).6)= LH(2. 4)=0.C)*sin(FasYbar(L.4)=0..3)=0.L) .6f %10. end end end %Calculando UH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 . N(4. Pprog(1)-Pcal(2).W(2).k). %Quarta Iteração %|-----------| %| M | LH | %|-----------| %| N | UH | %|-----------| %Calculando M a matriz que contém as derivadas parciais em relação as fases for i=2:4 for j=1:4 if i~=j M(i. M(3.2)=0..6f %10. end end Qcal(L) = (-1)*(V(L)^2)*B(L.end Pcal(L) end DPQ(1)= DPQ(2)= DPQ(3)= = (V(L)^2)*G(L. end end end %Calculando N matriz que contém as derivadas parciais em relação as Tensões for i=2:4 for j=1:4 if i~=j N(i. end end end N(1. end end end M(1.l)+N(k.j)= -abs((V(i)*V(j))*MagYbar(i.1)=0. N(2.3)=0. % valores de potência para serem adicionados fprintf('%6i %10.i).6f\n'.L) + Acumulator.W(3).6f %10. % valores de potência para serem adicionados DPQ(5)= Qprog(2)-Qcal(3).Acumulator.j)= -(abs(V(j)*V(i))*MagYbar(i. end end end %Calculando LH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j LH(i.2)=0.j)=-N(i.M(2.6f %10.j).V(3). for k=2:4 for l=1:4 if l~=k M(k. N(3. Pprog(2)-Pcal(3).k)= -N(k.i) = N(i.i)+2*(V(i)^2)*G(i. % valores de potência para serem adicionados DPQ(6)= Qprog(3)-Qcal(4). 3.l)+M(k.V(2).6f %10. % valores de potência para serem adicionados % valores de potência para serem adicionados % valores de potência para serem adicionados for L= 2:4 Acumulator=0.j))*sin(FasYbar(i.j))*cos(FasYbar(i.6f %10.k). else LH(i.W(4))..C)+W(C)-W(L))+ Acumulator. end DPQ(4)= Qprog(1)-Qcal(2). M(4. for k=2:4 for l=1:4 if l~=k N(k.W(1).k)= -M(k.1)=0.V(1).j)+W(j)-W(i)). Pprog(3)-Pcal(4).V(4).j)+W(j)-W(i)).6f %10. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L. L) + Acumulator. end DPQ(1)= Pprog(1)-Pcal(2).C)+W(C)-W(L))+ Acumulator.i). % valores de potência para serem adicionados DPQ(3)= Pprog(3)-Pcal(4).5)= Jacob(1. end end Jacob(1.3). TV=Jacob2*DPQ'. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L. % TV.2). W(4)=W(4)+TV(3). end DPQ(4)= Qprog(1)-Qcal(2).5)= Jacob(3. end end end %Obtenção da Matriz do Jacobiano for a=1:6 for d=1:6 if (d<=3 && a<=3) Jacob(d. LH(4.4)= Jacob(3.C)+W(C)-W(L))+ Acumulator.i)= -M(i.j)= M(i. LH(4.d+1).5)= Jacob(2.d2+1). Jacob2=inv(Jacob).2).a4+3)= UH(a4+1.a)= M(a+1. % Teste de Convergência p/ Quarta Iteração for L= 2:4 Acumulator=0. V(4)=V(4)*( 1 + TV(6)). end end Qcal(L) = (-1)*(V(L)^2)*B(L. % valores de potência para serem adicionados DPQ(2)= Pprog(2)-Pcal(3).Acumulator. LH(3.4)= Jacob(2.if i~=j UH(i. % valores de potência para serem adicionados for L= 2:4 Acumulator=0.C)*sin(FasYbar(L.3).4).4). LH(2. else UH(i. LH(4. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L.3 e 4 e tensões 2.j).3). end end Pcal(L) = (V(L)^2)*G(L.L) . V(3)=V(4)*( 1 + TV(5)). % Até aqui está OK Vetor de 6 posições 3 fases e 3 tensões fase2.a2)= N(a2+1. LH(2.i)-2*(V(i)^2)*B(i.3 e 4 %atualização dos valores das fases %atualização dos valores das fases %atualização dos valores das fases V(2)=V(2)*( 1 + TV(4)). W(2)=W(2)+TV(1). LH(3.6)= LH(2. LH(3.6)= Jacob(3. % valores de potência para serem adicionados .4)= Jacob(1. end end Jacob.4).6)= Jacob(2. end end end for a2=1:3 for d2=1:3 Jacob(d2+3.d4+1). W(3)=W(3)+TV(2). for a4=1:3 for d4=1:3 Jacob(d4+3.C)*cos(FasYbar(L.2). 4)=0.i)-2*(V(i)^2)*B(i.i)+2*(V(i)^2)*G(i. end end end %Calculando UH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j UH(i.i) = N(i.i).M(2. for k=2:4 for l=1:4 if l~=k M(k. end end end %10. M(3.j)= M(i.6f 4.i).k)= -N(k.W(2).3)=0.k).V(2).6f %10. M(4.W(4)). %10. end end end %Calculando LH matriz que contém as derivadas parciais em relação as Tensões for i= 2:4 for j=2:4 if i~=j LH(i.V(1).6f %10. end end end N(1. for k=2:4 for l=1:4 if l~=k N(k.6f %Quinta Iteração %|-----------| %| M | LH | %|-----------| %| N | UH | %|-----------| %Calculando M a matriz que contém as derivadas parciais em relação as fases for i=2:4 for j=1:4 if i~=j M(i.V(4).a)= M(a+1.j))*cos(FasYbar(i.6f %10.i)= -M(i.j))*sin(FasYbar(i. % valores de potência para serem adicionados % valores de potência para serem adicionados DPQ(6)= Qprog(3)-Qcal(4).j).W(1).6f\n'. end end end %Calculando N matriz que contém as derivadas parciais em relação as Tensões for i=2:4 for j=1:4 if i~=j N(i.k). fprintf('%6i %10.l)+N(k..6f %10.DPQ(5)= Qprog(2)-Qcal(3). else UH(i.j). N(4.j)= -abs((V(i)*V(j))*MagYbar(i. .W(3).l)+M(k.j)+W(j)-W(i)).2)=0.j)+W(j)-W(i)). end end end %Obtenção da Matriz do Jacobiano for a=1:6 for d=1:6 if (d<=3 && a<=3) Jacob(d.j)=-N(i.j)= -(abs(V(j)*V(i))*MagYbar(i.1)=0.6f %10. end end end M(1. N(2. else LH(i..3)=0..d+1).V(3). N(3.4)=0.k)= -M(k.1)=0.2)=0. C)+W(C)-W(L))+ Acumulator. end DPQ(1)= Pprog(1)-Pcal(2).1)*cos(W(1))) + V(2)*( G(1.4). for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L. end end Qcal(L) = (-1)*(V(L)^2)*B(L.3).C)*cos(FasYbar(L. LH(4.6f %10. % valores de potência para serem adicionados % valores de potência para serem adicionados DPQ(2)= Pprog(2)-Pcal(3). 5. TV=Jacob2*DPQ'.5)= Jacob(3.V(1).W(4)).6)= Jacob(3. % Este Vetor Apresenta as tensões nas barras corrigidas W.d2+1). LH(3. for C = 1:4 if C~=L Acumulator = abs(V(L)*V(C))* MagYbar(L. .2). % valores de potência para serem adicionados DPQ(5)= Qprog(2)-Qcal(3).2)*cos(W(1)-W(2)) )+ V(3)*(G(1.C)*sin(FasYbar(L..4). Jacob2=inv(Jacob).V(2).6f\n'.2)*sin(W(1)-W(2))-B(1.L) . LH(2.4).W(3).for a2=1:3 for d2=1:3 Jacob(d2+3.V(4). Q1=V(1)*( V(1)*(-B(1. % valores de potência para serem adicionados %Finalmente os resultados das iterações V.6f %10.4)= Jacob(1. %Potência Ativa e Reativa na Barra Oscilante P1=V(1)*( V(1)*(G(1. % Teste de Convergência p/ Quinta Iteração for L= 2:4 Acumulator=0.6f %10. % Até aqui está OK Vetor de 6 posições 3 fases e 3 tensões fase2.3)*sin(W(1)-W(3)) ) ).3 e 4 %atualização dos valores das fases %atualização dos valores das fases %atualização dos valores das fases V(2)=V(2)*( 1 + TV(4)).W(1). LH(4.5)= Jacob(2. end end Jacob. LH(3.L) + Acumulator..6f %10. W(3)=W(3)+TV(2).V(3).3).a2)= N(a2+1.2)*sin(W(1)-W(2)) )+ V(3)*(G(1.2).3). for a4=1:3 for d4=1:3 Jacob(d4+3.d4+1). V(4)=V(4)*( 1 + TV(6)). end end Jacob(1. V(3)=V(4)*( 1 + TV(5)).6f %10.2).3)*sin(W(1)-W(3))-B(1. W(4)=W(4)+TV(3).W(2). % valores de potência para serem adicionados DPQ(6)= Qprog(3)-Qcal(4).3)*cos(W(1)-W(3))+B(1.Acumulator. % Este Vetor Apresenta as fases das tensões nas barras corrigidas em graus fprintf('%6i %10. LH(4. % TV. end DPQ(4)= Qprog(1)-Qcal(2).1)*cos(W(1))) + V(2)*( G(1. LH(3. DPQ(3)= Pprog(3)-Pcal(4). W(2)=W(2)+TV(1).6)= LH(2.6f %10. end end Pcal(L) = (V(L)^2)*G(L. % valores de potência para serem adicionados for L= 2:4 Acumulator=0.3)*cos(W(1)-W(3)) ) ).C)+W(C)-W(L))+ Acumulator.4)= Jacob(2.4)= Jacob(3. LH(2.6)= Jacob(2.2)*cos(W(1)-W(2))+B(1.5)= Jacob(1.3 e 4 e tensões 2..6f %10.a4+3)= UH(a4+1. 7586 .2776 Q1 = 0.P1 = 1.