clear all;
f = 60;
w = 2*pi*f;
%
% FONTE
%
Vl    = 230E3;
Vf    = Vl/sqrt(3);
Scc1f = 1E9*exp(1i*80*pi/180);
Scc3f = 1.5E9*exp(1i*80*pi/180);
Z1    = Vl^2/conj(Scc3f);
Z0    = 3*Vl^2/conj(Scc1f) - 2*Z1;
R0TH  = real(Z0);
L0TH  = imag(Z0)/w;
R1TH  = real(Z1);
L1TH  = imag(Z1)/w;
%
% CARGA
%
P = 48E6;
Q = 18E6;
%
% TPC
%
VPTP = 30E3;
VSTP = 115;
C1   = 9660E-12;
C2   = 64400E-12;
L    = 1/(w^2*(C1+C2));
%
%
%
% satura = [1.41421356E-02 3.37618619E-02;5.36733089E-02 3.37618619E-01;1.31694552E-01 1.60556410E+00;1.75046597E-01 1.87565899E+00;1.89134128E-01 2.25079079E+00;3.41310866E-01 2.62592259E+00;5.61072569E-01 2.92602803E+00;9.75998771E-01 3.00105439E+00;9.43968011E+01 3.47747177E+00]
sim('TC_TP.slx')