Tworzenie skryptów i współpraca z plikami danych 33
Instrukcja pomocnicza return powoduje opuszczenie danego skryptu lub funkcji, a następnie powrót do jego (jej) wywołania (normalnie funkcja jest opuszczana po jej wykonaniu).
Dana jest zespolona macierz admitancji węzłowych
45 - j 45 -25 + j 25 0
-25 + j25 0
50 - j 50 -25 + y'25
- 25 + j25 25 - j25
(3.1)
oraz zespolony wektor prądów węzłowych
I = I.+jIb =
0
0
(3.2)
0.1-yl
Należy we własnym katalogu zapamiętać część rzeczywistą G macierzy admitancji w pliku G.dat oraz część urojoną B - w pliku B.dat. Podobnie należy postąpić z wektorem prądów węzłowych. Następnie utworzyć skrypt o nazwie uyi.m, realizujący następujące działania:
• czytanie z dysku macierzy G oraz B,
• tworzenie zespolonej macierzy admitancji węzłowych Y = G + jB,
• czytanie z dysku wektorów I„ oraz I/„
• tworzenie zespolonego wektora prądów węzłowych I = I„ + jI/„
• rozwiązanie układu równań liniowych Y U = /,
• wydzielenie części rzeczywistych i urojonych napięć węzłowych U = Uu + jUh
• zapisanie na dysku wartości obliczonych napięć węzłowych w plikach Ua.dat oraz Ub.dat
Zadanie 2
Napisać skrypt tworzący zespoloną macierz admitancji węzłowych na podstawie wczytanych z klawiatury numerów węzłów początkowych i końcowych oraz admitancji gałęzi liniowego obwodu elektrycznego pokazanego na rysunku 3.1. Utworzyć macierz wp - numerów węzłów początkowych gałęzi, wk - numerów węzłów końcowych gałęzi, g - konduktancji gałęzi, b - susceptancji gałęzi.