O Sistema Elétrico de Potência no Matlab

clear all; close all;  clc; 

  

%%dados globais 

R = 25; X = 125; bsh = 4e-6;  Vbase = 500e3; 

Sbase = 100e6; 

  

Tensao_E1= 500e3; 

Angulo_E1=0; 

E1=Tensao_E1*exp(deg2rad(Angulo_E1)*j); 

Tensao_E2=501.6e3; 

Angulo_E2=-4.5739; 

E2=Tensao_E2*exp(deg2rad(Angulo_E2)*j); 

  

E1_pu=E1/Vbase; 

E2_pu=E2/Vbase; 

  

Zbase = (Vbase^2)/Sbase; 

Z = sqrt((R^2)+(X^2)); %impedância 

  

R_pu = R/Zbase; 

X_pu = X/Zbase; 

Z_pu = Z/Zbase; %impedância 

  

Ykm = (R_pu-j*X_pu)/(Z_pu^2); 

Izkm=(E1_pu-E2_pu)*Ykm; 

  

Zsh=1/bsh; Zsh_pu=Zsh/Zbase; bsh_pu=1/Zsh_pu; 

  

Ikm_sh=E1_pu*j*bsh_pu; 

Imk_sh=E2_pu*j*bsh_pu; 

  

Ikm=Izkm+Ikm_sh; 

Imk=(-1*Izkm)+Imk_sh; 

  

Pkm=(abs(E1_pu^2)*real(Ykm))-

(abs(E1_pu)*abs(E2_pu)*real(Ykm)*cos(imag(E1_pu)-imag(E2_pu)))(abs(E1_pu)*abs(E2_pu)*imag(Ykm)*sin(imag(E1_pu)-imag(E2_pu))); 

Qkm=-(abs(E1_pu^2)*(bsh_pu+imag(Ykm)))(abs(E1_pu)*abs(E2_pu)*real(Ykm)*sin(imag(E1_pu)imag(E2_pu)))+(abs(E1_pu)*abs(E2_pu)*imag(Ykm)*cos(imag(E1_pu)imag(E2_pu))); 

    vm=0.95:0.1:1.05; for i=1:length(vm) 

    Pkm(i)=(abs(E1_pu^2)*real(Ykm))-

(abs(E1_pu)*vm(i)*real(Ykm)*cos(imag(E1_pu)-imag(E2_pu)))(abs(E1_pu)*vm(i)*imag(Ykm)*sin(imag(E1_pu)-imag(E2_pu)));     Qkm(i)=-(abs(E1_pu^2)*(bsh_pu+imag(Ykm)))(abs(E1_pu)*vm(i)*real(Ykm)*sin(imag(E1_pu)-

imag(E2_pu)))+(abs(E1_pu)*vm(i)*imag(Ykm)*cos(imag(E1_pu)-imag(E2_pu))); end figure(2) hold on 

plot(vm,Pkm); %plot de potência ativa com ângulo variável plot(vm,Qkm);%plot de potência ativa com amplitude variável 

  th=-5:0.1:5; th=th*pi/180; for i=1:length(th) 

    Pkm(i)=(abs(E1_pu^2)*real(Ykm))-

(abs(E1_pu)*abs(E2_pu)*real(Ykm)*cos(imag(E1_pu)-th(i)))(abs(E1_pu)*abs(E2_pu)*imag(Ykm)*sin(imag(E1_pu)-th(i))); 

    Qkm(i)=-(abs(E1_pu^2)*(bsh_pu+imag(Ykm)))(abs(E1_pu)*abs(E2_pu)*real(Ykm)*sin(imag(E1_pu)-

th(i)))+(abs(E1_pu)*abs(E2_pu)*imag(Ykm)*cos(imag(E1_pu)-th(i))); end figure(1) hold on 

plot(th,Pkm); %plot de potência ativa com ângulo variável plot(th,Qkm);%plot de potência ativa com amplitude variável 

[pic 1]

[pic 2]

 

%%limpar clear all; close all; clc;   

  

  

%% parametros globais 

Sbase = 100e6; %pot base 

Vbase = 500e3; %tensao base 

R = 25; 

X = 125; 

E1 = 1; %e gerador 1 pu Sbase_pu = 1; cont = 1; 

  

%% encontrando valores puimpedancias 

Zbase = (Vbase^2)/Sbase; 

Rfase_pu = R/Zbase; 

XLfase_pu = X/Zbase;