%% Faktoren

%% Xi Berechung
tau_P=r_i*pi/p;
aN=2*pi*p/N;
q=N1/(3*2*p);

Xi_s=sin(1*pi/2*Ws/tau_P);
Xi_z=sin(1*q*aN/2)/(q*sin(aN/2));
Xi_1=Xi_s*Xi_z;
Xi_Dd=Xi_1;
Xi_Dq=Xi_1;


%% Widerstände umrechnen auf Ständerseite bezogen
rf_Strich=3/2*c_d1/C1*4/pi*w1^2*Xi_1^2/wF^2*rf;
rDd_Strich=3/2*(w1^2*Xi_1^2)/(wDd^2*Xi_Dd^2)*rDd;
rDq_Strich=3/2*(w1^2*Xi_1^2)/(wDq^2*Xi_Dq^2)*rDq;


%% Streuinduktivitäten umrechnen auf Ständerseite bezogen
lc_Strich=3/2*c_d1/C1*4/pi*w1^2*Xi_1^2/wF^2*l_c;
lf_Strich=3/2*c_d1/C1*4/pi*w1^2*Xi_1^2/wF^2*l_f;
lDd_Strich=3/2*(w1^2*Xi_1^2)/(wDd^2*Xi_Dd^2)*l_Dd;
lDq_Strich=3/2*(w1^2*Xi_1^2)/(wDq^2*Xi_Dq^2)*l_Dq;

lfs_Strich=lf_Strich+lc_Strich;
lDds_Strich=lDd_Strich+lc_Strich;
lDqs_Strich=lDq_Strich+lc_Strich;

ld=l_1s+l_hd;
lq=l_1s+l_hq;
lf=lfs_Strich+l_hd;
lDd=lDds_Strich+l_hd;
lDq=lDqs_Strich+l_hq;



%% Umrechnungsfaktoren zurück
Fak_zu_If=Xi_1*w1/wF*c_d1/C1*6/pi;
Fak_zu_IDd=3/2*(w1*Xi_1)/(wDd*Xi_Dd);
Fak_zu_IDq=3/2*(w1*Xi_1)/(wDq*Xi_Dq);

Fak_zu_Uf_Strich=Xi_1*w1/wF;




%% Matrix-Elemente (Flussverkettung-Gleichungssystem => I="a-Matrix" *phi)
var_1=ld*(lDd*lf-l_hd^2)-(lDd+lf-2*l_hd)*l_hd^2;
var_2=lDq*lq-l_hq^2;

a11=(lDd*lf-l_hd^2)/var_1;
a12=0;
a13=(-(lDd-l_hd)*l_hd)/var_1;
a14=(-(lf-l_hd)*l_hd)/var_1;
a15=0;

a21=0;
a22=l_hq^2/(var_2*lq)+1/lq;
a23=0;
a24=0;
a25=-l_hq/var_2;

a31=(-(lDd-l_hd)*l_hd)/var_1;
a32=0;
a33=((lDd^2-2*lDd*l_hd+l_hd^2)*l_hd^2)/(var_1*(lDd*lf-l_hd^2))+l_hd^2/((lDd*lf-l_hd^2)*lf)+1/lf;
a34=((lDd-l_hd)*(lf-l_hd)*l_hd^2)/(var_1*(lDd*lf-l_hd^2))-l_hd/(lDd*lf-l_hd^2);
a35=0;

a41=(-(lf-l_hd)*l_hd)/var_1;
a42=0;
a43=(-(ld-l_hd)*l_hd)/var_1;
a44=(ld*lf-l_hd^2)/var_1;
a45=0;

a51=0;
a52=-l_hq/var_2;
a53=0;
a54=0;
a55=lq/var_2;