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

Anhang
text(60, Teilung*3,k)
% Anzeigen der berechneten Werte am Bildschirm (XL)
[g,erromsg]=sprintf('%5.2f',XL)
XL='XL=';
O=' Ohm';
k=strcat(XL,g,O);
text(60, Teilung*2,k)
% Anzeigen der berechneten Werte am Bildschirm (P2max)
[h,erromsg]=sprintf('%5.2f',P2MAX)
P2MAX='P2max=';
W=' Watt';
k=strcat(P2MAX,h,W);
text(60, Teilung,k)
hold on;
8.1.4 Energieübertragung mit Spannungsregler und Kompensation
(Blindleistungsbetrachtung nur für die Leitung und Transformator)
function EUMRCBL (X);
%EUMRCBL ... Energie Übertragung mit Regler und Blindleistungsbetrachtung nur für die Leitung mit Kompensation der
Längsreaktenz
format long e ;
% Berechnung von Omega
w=2*pi*50;
% Widerstandsbelag der Freileitung
Rstrich=3e-002;
% Induktivitätsbelag der Freileitung
XLstrich=3e-001;
% Berechnung der Betriebsreaktanz XL und der Induktivität L
XL=(XLstrich*X);
L=(XL/w);
% Halber wert der Betriebskapazität, laut Ifea Modell
cbh=L/(2*274^2);
% Berechnung des Betriebswiderstandes RL
RL=Rstrich*X;
% Berechnung der Suszeptanzen der Freileitung ,wobei der Index p für Pi-Ersatzschaltung steht
y1p=w*cbh*i;
y2p=y1p;
% Berechnung der Betriebsimpedanz der Freileitung
z3p=RL+XL*i;
z1p=inv(y1p);
z2p=z1p;
% Korrektur der Impedanzen für Feileitungslängen größer 300km
z3pk=((sinh(sqrt(2*z3p/z2p)))/(sqrt(2*z3p/z2p)))*z3p;
z2pk=((sqrt(z3p/(2*z2p)))/(tanh(sqrt(z3p/(2*z2p)))))*z2p;
z1pk=z2pk;
% Umrechnung Impedanz-Admittanz
y1pk=inv(z1pk);
y2pk=inv(z2pk);
y3pk=inv(z3pk);
% Vereinfachung-Faktorisierung (Beginn der Umrechnung Pi-Ersatzschaltung in die T-Ersatzschaltung)
k1=y1pk*y2pk;
k2=y2pk*y3pk;
k3=y1pk*y3pk;
k=k1+k2+k3;
% Berechnung der Ersatzimpedanzen der T-Ersatzschaltung
z1t=y1pk/k;
z2t=y2pk/k;
z3t=y3pk/k,
% Parameter der Trafo und Synchronreaktanz
zd=117.9*i;
zt=5.74*i;
%Parameter des Serienkondensators
zc=-XL*i
% Zusammenfassung der Impedanzen der T-Ersatzschaltung
z1g=z1t+zd+zt;
Diplomarbeit
91