clear all;

TN1=10.92;
TN2=10;
KSA=1.6;
T1=8;
T2=6;
KSB=1.25;
T3=10;
Ta=1.4;
KP1=0.866;
% w_d2=0.0638;
% w=w_d2;
Tt=0;
syms w;

GSA=KSA/((1+T1*w*j)*(1+T2*w*j));
% GR1ohneKP1=((1+TN1*w*j)/(TN1*w*j));
GAH=(sin(Ta/2*w)/(Ta/2*w))*exp(-Ta/2*j*w);

GS1=GSA;

% G01_Stern=GR1ohneKP1*GS1*GTt;
% 
% G01_Stern_Betrag=abs(G01_Stern);
% KP1=1/G01_Stern_Betrag;



GR1=KP1*(1+TN1*w*j)/(TN1*w*j);

G01=GR1*GS1*GAH;

G01_Imag=imag(G01);

G01_Real=real(G01);

phi_0=atan(G01_Imag/G01_Real);

G01_Betrag=sqrt(G01_Real^2+G01_Imag^2);

GW1_Betrag=G01_Betrag/(1+2*G01_Betrag*cos(phi_0)+(G01_Betrag^2))^0.5;

phi_W1=atan(sin(phi_0)/(G01_Betrag+cos(phi_0)));
phi_W1=phi_W1+pi;

 

GW1=GW1_Betrag*exp(j*phi_W1);