Modelluntersuchungen betreffend die Stabilität - Institut für ...

Anhang
psi12=angle(z12)
psi21=angle(z21)
psi22=angle(z22)
% Berechnung der Komplementärwinkel
a11=(pi/2) - psi11
a12=(pi/2) - psi12
a21=(pi/2) - psi21
a22=(pi/2) - psi22
% Umrechnung rad-grad
a11g=a11 * 360 / (2*pi)
a12g=a12 * 360 / (2*pi)
a21g=a21 * 360 / (2*pi)
a22g=a22 * 360 / (2*pi)
% Berechnung des Betrages der komplexen Impedanzen
z11abs=abs(z11)
z22abs=abs(z22)
z12abs=abs(z12)
% Berechnung der Faktoren A1,A2,B zur Berechnung der Leistungen
A1=(u1*u1*sin(a11)) / z11abs
A2=(u2*u2*sin(a22)) / z22abs
B=(u1*u2) / z12abs
% Laufvariable des Polradwinkels
thetag= 0 : 5 : 180
thetar=thetag*pi/180;
% Berechnung der Leistungen
P1ges=3* ( A1 + (B*sin(thetar-a12)))
P2ges=3* (-A2 + (B*sin(thetar+a12)))
% Achsenbeschriftung und Plotkonfiguration
plot(thetag,P1ges,'-',thetag,P2ges,'--')
grid on
legend('In die Leitung hineingeschickte Leistung P1','Aus der Leitung herauskommende Leistung P2')
xlabel('Theta [°] (Winkel zwischen Polradspannung Ep und Netzspannung U2)')
ylabel('Wirkleistung [W]')
title('Übertragene Leistung als Fkt. des Polradwinkels')
% Maximalwert der übertagbaren Leistung
% p2max=max(P2ges);
p2max=3*(-A2+B);
P2MAX=p2max;
% Anzeigen der berechneten Werte am Bildschirm (Ep)
Teilung=p2max/10;
[a,erromsg]=sprintf('%5.2f',u1)
U1='Ep=';
V=' V';
k=strcat(U1,a,V);
text(60, Teilung*7,k)
% Anzeigen der berechneten Werte am Bildschirm (U2)
[u,erromsg]=sprintf('%5.2f',u2)
U2='U2=';
V=' V';
k=strcat(U2,u,V);
text(60, Teilung*6,k)
% Anzeigen der berechneten Werte am Bildschirm (Länge)
Teilung=p2max/10;
[s,erromsg]=sprintf('%5.2f',X)
L='Länge=';
KM=' km';
k=strcat(L,s,KM);
text(60, Teilung*5,k)
% Anzeigen der berechneten Werte am Bildschirm (Cbh)
[d,erromsg]=sprintf('%5.9f',cbh)
CBH='Cbh=';
F=' F';
k=strcat(CBH,d,F);
text(60, Teilung*4,k)
% Anzeigen der berechneten Werte am Bildschirm (RL)
[f,erromsg]=sprintf('%5.2f',RL)
RL='RL=';
O=' Ohm';
k=strcat(RL,f,O);
Diplomarbeit
90