clc;
clear;

Z_ref = 1.20822;        
U_DC = 11760.6;
R = 0.00556;
L = 0.1207;
T_vsi = 330e-9;              %%T_vsi ... konst

wc_hz = 200;                 %%Grenzfrequenz in Hz
wc_rad = wc_hz*2*pi;         %%Grenzfrequenz in rad/s ... Überschwingen 3%
                              
G = -(Z_ref*U_DC)/(R*sqrt(3))*tf([0 1],[(L/R) 1]);  %%ÜF - Trafo
G_vsi = tf([0 1],[T_vsi 1]);                        %%ÜF - Stromrichter
G_s = G*G_vsi;                                      %%ÜF - Trafo+Stromrichter
L1_s = G_s*tf([0 1],[1 0]);                         %%offender Kreis mit Integrator

G_jwc = -(Z_ref*U_DC)/(R*sqrt(3))*(1/(1i*wc_rad*(L/R) + 1)); 
G_vsi_jwc = 1/(1 + 1i*wc_rad*T_vsi);
G_s_jwc = G_jwc * G_vsi_jwc;          %%Strecke bei Grenzfrequenz
L1_jwc = 1/(1i*wc_rad)*G_s_jwc;       %%offene Strecke bei Grenzfrequenz mit Integrator

T_Ri = tan((-180+70)*pi/180-(angle(L1_jwc)))/wc_rad;   %%Zeitkonstante PI-Regler

L_jwc = L1_jwc * (1+1i*wc_rad*T_Ri);  %%Offene Strecke mit bekannten Parametern zur Berechnung der Verstärkung

V_Ri = 1/abs(L_jwc);    %%Verstärkung PI-Regler

Ri = V_Ri * tf([T_Ri 1],[1 0]); %%Ausgelegter PI-Regler

sisotool(G_s,Ri)