Instr Cw6 2016 odblokowany

WYDZIAA ETI PG
Katedra Systemów Elektróniki Mórskiej
Technika obliczeniowa i symulacyjna
- laboratorium
Automatyczne konstruowanie równań liniowych
opisujących obwód elektryczny
oraz ich numeryczne rozwiązywanie.
MATLAB i SPICE jako narzędzia do obliczania prądów, napięć
i mocy oraz elementów zastępczych w obwodzie elektrycznym
(Materiały dó ćwiczenia nr 6)
Opracówał: Czesław Stefanski
Gdansk 2015
C. Stefański, Materiały do ćw 6 Lab. TOiS 2/12 Ostatnia modyfikacja: maj 2016
1. Wstęp
Niniejsze ćwiczenie jest kolejnym, które jest poświęcone wdrożeniu umiejętności obliczania prądów i napięć w ob-
wodzie przy wykorzystaniu narzędzia MATLAB z jego algorytmami numerycznego rozwiązywania układów równań linio-
wych oraz nabyciu obycia w korzystaniu z symulacyjnego narzędzia Micro5CAP (jedna z odmian SPICE a)1. Jednak tym
razem będziemy starali się tak zorganizować proces badania obwodu, by jeszcze bardziej odciążyć użytkownika od ręcz-
nego udziału w tym badaniu w tym celu spróbujemy zmusić komputer, by nie tylko rozwiązywał układy równań
opisujące obwód, ale by je również sam tworzył! W tej sytuacji warto przypomnieć, że przy tworzeniu algorytmów
analizy automatycznej (komputerowej) wygodnie jest trzymać się zasady, że każdy dwójnik ma swój prąd i swoje
napięcie, przeciwnie względem siebie zastrzałkowane (czyli stosujemy tzw. strzałkowanie skojarzone). Przez to algo-
rytmy stają się bardziej jednorodne i dzięki temu zwykle są szybsze mimo większej niż przy strzałkowaniu chaotycznym
liczby niewiadomych.
Przypominamy, że na rysunkach 1 mamy do czynienia z tzw. skojarzonym strzałkowaniem prądów i napięć polega
ono na tym, że zwrot spadku napięcia na każdym z elementów obwodu jest przeciwny do kierunku przepływającego
przez ten element prądu.
W prostym obwodzie elektrycznym, na przykład składającym się z podobnej liczby elementów jak w obwodzie
przedstawionym na rys. 1a, poszukiwanie takich podstawowych wielkości elektrycznych jak natężenie prądu (prąd) 5؊�,
czy też różnica potencjałów (napięcie) 5�
�, może być względnie łatwym zadaniem, również dlatego, że są to liniowe ob-
wody prądu stałego. Jednak ze wzrostem liczby i różnorodności elementów tworzących strukturę obwodu stopniowo
zmniejsza się możliwość czy też opłacalność ręcznego (analitycznego) poszukiwania prądów i napięć gałęziowych.
u u
1 2
v�
u�
i i
1 2
i i J
R
R i 3 i 4 5
2 3
1 2
i i
0 1
u u u u u u u
0 2 3 1 3 4 5
E E
0 1
R R
3 4
i
5
R R
2 3
t
Rys. 1a. Prosty obwód prądu stałego Rys. 1b. Przykład strzałkowania skojarzonego
2. Poszukiwanie prądów, napięć i mocy w rozbudowanym obwodzie elektrycznym
Dla ustalenia uwagi przyjmiemy, że naszym zadaniem będzie dokonanie w miarę automatycznej i wszechstronnej
analizy obwodu z rysunku 1b, w tym wyznaczenie prądów i napięć, a także mocy na każdym z jego elementów. Ponadto
porównamy wyniki naszej analizy z wynikami symulacji przeprowadzonej w Micro5Cap-ie, a także sprawdzimy bilanso-
wanie się mocy.
Jednak tym razem wykorzystamy opis obwodu elementarnego poziomu, tzw. opis tableau, którego tworzenie jest
w miarę łatwe do oprogramowania.
Pozostałe pojęcia, algorytmy i narzędzia, wykorzystywane w tym ćwiczeniu, zostały już opisane w materiale do
wcześniejszych ćwiczeń. W szczególności wykorzystamy macierz incydencji 5�4� grafu obwodu oraz tzw. transformację
węzłową 5�|� = 5�h�5�{� �" 5�}�.
6
m�
3
.
4
9
1 -1 1 0 0 0 0 0 0
2
k
.
.
j� .
l
"A = [0 1 0 -1 1 0 -1 0 0]
7 0 0 0 0 0 -1 1 1 -1
5
0 0 -1 1 0 1 0 0 1
1
8
i
Rys. 2. Przykładowy graf prądowy i jego zredukowana macierz incydencji
1
Narzędzie to, po narysowaniu obwodu w jego oknie graficznym przez użytkownika, potrafi automatycznie, acz bez ujawniania ich użytkownikowi,
utworzyć odpowiednie równania opisujące obwód, rozwiązać je (podobnymi jak MATLAB algorytmami) i udostępnić rozwiązania między innymi we
wspomnianym oknie z obwodem.
C. Stefański, Materiały do ćw 6 Lab. TOiS 3/12 Ostatnia modyfikacja: maj 2016
Graf z rysunku 2, przy założeniu, że węzeł i obwodu jest węzłem odniesienia, stanowił podstawę do zapisania jego
zredukowanej macierzy incydencji 5�h� = [5؂�5�Ś�5�ć�] (dalej nazywanej krótko macierzą incydencji) według następującej umowy:
5�N�5�X�5�Y� = 1 w przypadku, gdy prąd w gałęzi o numerze 5�Y� wypływa z węzła 5�X� ,
5�N�5�X�5�Y� = -1 w przypadku, gdy prąd w gałęzi o numerze 5�Y� wpływa do węzła 5�X� ,
5�N�5�X�5�Y� = 0 w przypadku, gdy gałąz o numerze 5�Y� nie łączy się z węzłem 5�X� .
Metoda tableau opisu obwodu
Dla obwodu o 5�O� gałęziach oraz (5�[� + 1) węzłach [dla obwodu z rysunku 1b mamy 5�O� = 5 i 5�[� = 2] węzły numerujemy
tak, by węzeł odniesienia (masy) miał numer zero, a pozostałe węzły numery od 1 do 5�[�. Wtedy napięcia i prądy w takich
gałęziach można przedstawić w formie następujących wektorów kolumnowych:
( )
5�|� = col 5�H�1, 5�H�2, & , 5�H�5�O� ,
( )
5�p� = col 5�<�1, 5�<�2, & , 5�<�5�O� , (1)
( )
5�}� = col 5�I�1, 5�I�2, & , 5�I�5�[� ,
gdzie 5�}�, to wektor potencjałów węzłowych czyli potencjałów 5�[� węzłów w odniesieniu do potencjału węzła odniesienia
(czyli węzła o numerze zero).
Prądy i napięcia we wszystkich gałęziach obwodu złożonego z 5�O� gałęzi (na przykład o postaci jak na rysunku 1b) oraz
5�[� węzłów niezależnych i węzła odniesienia, podlegają prawom Kirchhoffa, które można przedstawić za pomocą nastę-
pujących zapisów macierzowych
5�h� �" 5�p� = 5��
(2)
5�|� = 5�h�5�{� �" 5�}�
gdzie 5�h�5�{� to transponowana macierz incydencji.
Do pełnego opisu obwodu brakuje równań opisujących zachowanie się elementów; dla oporników są to równania prawa
Ohma, dla niezależnych idealnych zródeł napięciowych równania postaci 5�H�5�8� = 5�8�, itd. Takich równań powinno być łącz-
nie 5�O�. W zapisie macierzowym można im nadać formę
5�t� �" 5�p� + 5�u� �" 5�|� = 5�l�5ؙ�5�"�, (3)
( )
gdzie 5�t� = [5�Z�5�V�5�W�]5�O��5�O�, 5�u� = [5�[�5�V�5�W�]5�O��5�O�, 5�l�5ؙ�5�"� = col 5�R�5�e�1, 5�R�5�e�2, & , 5�R�5�e�5�O� .
Na przykład, gdy element o numerze 5�X� jest opornikiem o oporze 5�E�5�X�, to 5�Z�5�X�5�X� = 5�E�5�X�, 5�[�5�X�5�X� = -1, 5�R�5�e�5�X� = 0, co odpo-
wiada (pod warunkiem, że 5�Z�5�X�5�W� = 5�[�5�X�5�W� = 0 dla 5�W� `" 5�X�) zapisowi prawa Ohma dla tego elementu w postaci
5�E�5�X� �" 5�<�5�X� - 5�H�5�X� = 0.
Podobnie, gdy element o numerze 5�Y� jest zródłem prądowym o SPM 5�=�5�Y�, to 5�Z�5�Y�5�Y� = 1, 5�[�5�Y�5�Y� = 0, 5�R�5�e�5�X� = 5�=�5�Y�, co odpowiada
(pod warunkiem, że 5�Z�5�Y�5�W� = 5�[�5�Y�5�W� = 0 dla 5�W� `" 5�Y�) zapisowi definicyjnego równania dla tego elementu w postaci
1 �" 5�<�5�Y� - 0 �" 5�H�5�Y� = 5�=�5�Y�
(założono, że strzałki 5�<�5�Y� oraz 5�=�5�Y� są zgodne).
Powyższe równania macierzowe dają się zapisać w postaci
5�h� 5�� 5�� 5�p� 5��
[ ] �" [ ] = [ ] . (4)
5�� 5�h�5�G� -5�� 5�}� 5��
5�t� 5�� 5�u� 5�|� 5�l�5ؙ�5�"�
Macierz kwadratowa z ostatniego równania nosi nazwę macierzy tableau.
Powyższe równanie może być łatwo przekształcone do postaci
[5�h� 5�� ] �" [5�p�] = [5�� ] . (5)
5�t� 5�u��"5�h�5�G� 5�}� 5�l�5ؙ�5�"�
Równanie (5) jest macierzowym zapisem układu (5�O�+5�[�) równań liniowych z (5�O�+5�[�) niewiadomymi (na marginesie
zauważmy, że, przy naszych oznaczeniach, 5�O�+5�[� byłoby liczbą niezależnych cykli w napięciowym grafie obwodu, co ozna-
cza, że ten układ równań jest równoliczny z układem równań, jaki uzyskalibyśmy stosując metodę prądów obwodowych
(w szczególności oczkowych)).
Rozwiązaniem układu równań (5) są potencjały wszystkich (niezależnych) 5�[� węzłów obwodu oraz prądy dwójników
obwodu. Przy oznaczeniach 5�z� = [5�� ], 5�� = [5�p�], 5�G� = [5�h� 5�� ] rozwiązanie 5�� możemy uzyskać jedną instrukcją MA-
5�l�5ؙ�5�"� 5�}� 5�t� 5�u��"5�h�5�G�
TLAB-a:
X=S\T
(6)
Bilans mocy w obwodzie
W poprzednich ćwiczeniach stwierdziliśmy, że bilans mocy w obwodzie układamy przyjmując, że iloczyn prądu i napięcia przy przeciwnej orien-
tacji napięcia na dwójniku i płynącego przez ten dwójnik prądu oznacza moc traconą na tym dwójniku; zawsze na rezystorach moc tak liczona jest
dodatnia, zaś bardzo często liczona w ten sam sposób moc na zródłach jest ujemna; ujemną moc traconą należy interpretować (po zmianie znaku)
jako dodatnią moc dostarczoną.
Spójrzmy na zagadnienie bilansu najpierw ogólnie, pamiętając, że dla prądów i napięć gałęzi mieliśmy zależności
(patrz wzory (2)):
5�h� �" 5�p� = 5��, 5�|� = 5�h�5�{� �" 5�}�.
C. Stefański, Materiały do ćw 6 Lab. TOiS 4/12 Ostatnia modyfikacja: maj 2016
Wykorzystamy je w następującym rachunku:
"5�C� = 5�|�5�G� �" 5�p� = 5�h�5�{� �" 5�}� �" 5�p� = 5�}�5�{� �" 5�h� �" 5�p� = 5�}�5�{� �" 5�h� �" 5�p� = 5�}�5�{� �" 5�� = 0 . (7)
( )5�G� ( ) ( )
Oznacza to, że suma mocy traconych na wszystkich elementach w obwodzie o 5�O� dwójnikach wynosi zero (jak zresztą
wynika z twierdzenia Tellegena).
Bilans mocy winien być oczywiście zerowy (suma po wszystkich elementach mocy traconych na elementach = zero):
tracona
""5�C�dwójników = 0. (8)
Koncepcja zródeł zastępczych Thevenina i Nortona
Jak pamiętamy z wcześniejszego opracowania, jeżeli nasze potrzeby ograniczają się do znalezienia prądu (napięcia, mocy) tylko w kilku zgrupo-
wanych elementach (w skrajnym przypadku w jednym elemencie) obwodu, to całą pozostałą część struktury możemy zastąpić tzw. zródłem Thevenina
lub zródłem Nortona. Przyjmując, tym razem, że wybranym ważnym elementem jest idealne zródło prądowe 5�=�5 możemy,
zgodnie z powyższą koncepcją, wieloelementowy obwód zredukować do jednej z przedstawionych na rysunku 3a,b po-
staci.
) c)
a b
)
(x)
k k
q
Obwód liniowy
J I
5
RN
J
5
zastępowany
5
RT
5
zródłem
U U
q
U
Thevenina
JN
lub Nortona
ET
i i
(y)
Rys. 3. Obwód zastępczy Thevenina a), obwód zastępczy Nortona b), schemat ogólny do wyznaczania parame-
trów obwodu zastępczego c).
Poza parametrami idealnego zródła prądowego 5�=�5, parametry pozostałych elementów z rysunków 3a,b muszą zo-
stać dopiero obliczone na tej podstawie, że na zródłe 5�=�5 winien wystąpić taki sam spadek napięcia 5�H�, jaki występował
na tym elemencie w ramach obwodu z rysunku 1b.
Wprowadzone przez Thevenina oraz Nortona definicje dają nam następujące zapisy, które odnoszą się do ry-
sunku 1b i rysunków 3a,b:
5�8�5�G�
5�8�5�G� = 5�H�5| 5�=�5=0 , 5�=�5�A� = 5�<�5| 5�=�5 zastąpione zwarciem , 5�E�5�G� = 5�E�5�A� = . (9a)
5�=�5�A�
Inaczej sprawę ujmując, napięcie Thevenina otrzymujemy wtedy, kiedy w obwodzie z rysunku 1b wyzerujemy ele-
ment 5�=�5, co odpowiada rozwarciu zacisków obwodu zastępowanego zródłem Thevenina. W sytuacji analizy metodą ta-
bleau wymaga to tylko zmiany wartości elementu 5�=�5 bez zmiany kształtu równania (5).
Prąd Nortona otrzymamy na podstawie rozwiązania równań (5) po zastąpieniu prądowego zródła 5�=�5 zródłem napię-
ciowym 5�8�5 = 0. Co prawda nie wymaga to tym razem od nas przekonstruowania równań typu (5) opisujących obwód
(bowiem algorytm uczyni to automatycznie), ale dane wejściowe trzeba będzie skorygować.
Na szczęście, jak już wiemy, do wyznaczenia parametrów zródeł Thevenina i Nortona można podejść bardziej jed-
nolicie, mianowicie parametry te można obliczyć z ogólniejszych zależności, przy oznaczeniach z rysunku 3c:
2 2 1
5�<� �"15�H� -15�<� �"25�H� 5�H� �"15�<� -15�H� �"25�<� 5�H� -25�H�
5�8�5�G� = , 5�=�5�A� = , 5�E�5�G� = 5�E�5�A� = (9b)
2 2 2
5�<� -15�<� 5�H� -15�H� 5�<� -15�<�
Wynikiem analizy tableau naszego obwodu są prądy i napięcia elementowe (a także potencjały węzłowe), zatem
1
wzór (9b) może być wykorzystany bezpośrednio. Na przykład można przyjąć, że (15�H�, 5�<� ) są odpowiednio napięciem
2
i prądem zródła 5�=�5 w obwodzie bez modyfikacji, zaś (25�H� , 5�<� ) to odpowiednio napięcie i prąd zródła 5�=�5, gdy wartość 5�=�5 =
0 :
2 2
5�H� �"15�<� 5�H� -15�H�
2
5�8�5�G� = 5�H�, 5�=�5�A� = , 5�E�5�G� = 5�E�5�A� = . (9c)
2 1
5�H� -15�H� 5�<�
Znowu do wyznaczenia zródeł Thevenina lub Nortona, gdy korzystamy z wzorów (9b, c) wystarczy odpowiednio
manipulować obciążeniem , a nie strukturą całego obwodu.
Policzmy jeszcze moce w obwodach 3a-b, starając się trzymać w konwencji, że korzystamy z wyliczonych prądów
i napięć elementowych na dopełniającym zródło zastępcze dwójniku 5�=�5.
Moce (liczone jako tracone (!!!)) w obwodzie z rysunku 3a
obwodu, q=1 lub 2
q-ty wariant obciążenia
I
I
C. Stefański, Materiały do ćw 6 Lab. TOiS 5/12 Ostatnia modyfikacja: maj 2016
tracona 2
5�C�rezystorów = 5�=�5 �" 5�E�5�G�,
(13)
tracona
5�C�zródeł = -5�=�5 �" 5�8�5�G� + 5�=�5 �" 5�H�5.
Moce (liczone jako tracone (!!!)) w obwodzie z rysunku 3b
2
5�H�5
tracona
5�C�rezystorów = ,
5�E�5�A�
(14)
tracona
5�C�zródeł = -5�=�5�A� �" 5�H�5 + 5�=�5 �" 5�H�5.
Warto zastanowić się, czy (i dlaczego) suma mocy traconych na wszystkich, oprócz idealnego zródła prądowego 5�=�5,
elementach obwodu 1b jest równa sumie mocy traconych na wszystkich, oprócz idealnego zródła prądowego 5�=�5, ele-
mentach obwodu z rysunku 3a, jak i 3b.
3. Wspomagane programem MATLAB rozwiązywanie układu równań liniowych
Przypomnimy teraz krótko to, co zostało bardziej szczegółowo omówione w instrukcji do poprzedniego ćwiczenia, a co dotyczyło rozwiązywania
układu równań liniowych wspomaganego narzędziami typu MATLAB (lub mu podobnymi, jak OCTAVE, czy SciLab). Układ równań można zapisywać w
formie klamrowej (15):
5�N�115�e�1 + 5�N�125�e�2 + �" 5�N�15�[�5�e�5�[� = 5�O�1
5�N�215�e�1 + 5�N�225�e�2 + �" 5�N�25�[�5�e�5�[� = 5�O�2
{ (15)
& & & & & & & & & & .
5�N�5�[�15�e�1 + 5�N�5�[�25�e�2 + �" 5�N�5�[�5�[�5�e�5�[� = 5�O�5�[�
lub w postaci macierzowej
5�N�11 5�N�12 5�N�13 & . 5�N�15�A� 5�e�1 5�O�1
5�N�21 5�N�22 5�N�23 & . 5�N�25�A� 5�e�2 5�O�2
5�N�31 5�N�32 5�N�33 & 5�N�35�A� �" 5�e�3 = 5�O�3 . (16)
& & & & & &
&
[5�N�5�A�1 5�N�5�A�2 5�N�5�A�3 & 5�N�5�A�5�A�] [5�e�5�A�]
[5�O�5�A�]
Skrócony zapis (przy notacji bez nawiasów kwadratowych) równania (16) wygląda następująco:
5�h� " 5ؙ� = 5؃� (17)
W tym przypadku rozwiązanie jest w naszym zasięgu nawet bez wspomagania komputerowego, przy czym często korzystamy z przejrzystych
zapisów metody Cramera
5��5�� 5��5�� 5��5�Ń� 5��5��
5ؙ�5�� = , 5ؙ�5�� = , 5ؙ�5�Ń� = , 5ؙ�5�� = , 5�� `" 5�� . (18)
5�� 5�� 5�� 5��
gdzie " jest wyznacznikiem macierzy 5�h�, zaś każdy z "5�[� to wyznacznik macierzy powstałej z 5�h� przez zastąpienie w niej wektorem 5؃� kolumny o numerze
5�[�-tym .
Współcześnie, spośród wielu opracowanych metod numerycznego rozwiązywania układów równań liniowych, duże znaczenie mają różne wa-
rianty metody eliminacji Gaussa. Algorytm tej metody dla równania (16) jest następujący (eliminacja niewiadomych w wyrażeniu (16) przebiega
następująco):
�� W pierwszym kroku od wierszy 2, 3, & . , 5�A� odejmujemy wiersz pierwszy pomnożony przez taki współczynnik, który zapewnia jako
wynik odejmowania nowy wiersz nie zawierający już niewiadomej 5�e�1 Po pierwszym kroku (krok 5�`� = 1; ang.: step) otrzymujemy
.
5�N�11 5�N�12 5�N�13 & 5�N�15�A� 5�O�1
5�e�1
(1) (1) (1)
0 5�N�22 5�N�23 & 5�N�25�A� 5�O�2
5�e�2 (1)
(1) (1) (1)
0 5�N�32 5�N�33 & 5�N�35�A� 5�O�3
5�e�3 (1)
�" = (19a)
& & & & & & &
(1) (1) (1)
(1)
[ 0 5�N�5�A�2 5�N�5�A�3 & 5�N�5�A�5�A� [5�e�5�A�]
] 5�O�5�A�
[ ]
gdzie
5�N�5�V�1
(1) (1)
5�N�5�V�5�X� = 5�N�5�V�5�X� - 5�Y�5�V�1 �" 5�N�15�X� , 5�O�5�V� = 5�O�5�V� - 5�Y�5�V�1 �" 5�O�1, przy 5�Y�5�V�1 = oraz 5�V�, 5�X� = 2, 3, & , 5�A�
5�N�11
�� W drugim kroku od wierszy 3, 4, & . , 5�A� odejmujemy wiersz drugi pomnożony przez taki współczynnik, który zapewnia jako wynik
odejmowania nowy wiersz nie zawierający już niewiadomej 5�e�2 Po drugim kroku (krok 5�`� = 2) otrzymujemy
.
5�N�11 5�N�12 5�N�13 & 5�N�15�A� 5�O�1
5�e�1
(1) (1) (1)
0 5�N�22 5�N�23 & 5�N�25�A� 5�O�2
5�e�2 (1)
(2) (2)
0 0 5�N�33 & 5�N�35�A� 5�O�3
5�e�3 (2)
�" = (19b)
& & & & & & &
(2) (2)
(2)
[ 0 0 5�N�5�A�3 & 5�N�5�A�5�A� [5�e�5�A�]
] 5�O�5�A�
[ ]
gdzie
(1
(2) (1) (1 (2) (1) ( 5�N�5�V�2)
5�N�5�V�5�X� = 5�N�5�V�5�X� - 5�Y�5�V�2 �" 5�N�25�X�) , 5�O�5�V� = 5�O�5�V� - 5�Y�5�V�2 �" 5�O�21), przy 5�Y�5�V�2 = oraz 5�V�, 5�X� = 3, 4, & , 5�A�
(1
5�N�22)
C. Stefański, Materiały do ćw 6 Lab. TOiS 6/12 Ostatnia modyfikacja: maj 2016
�� Powyższe postępowanie stosujemy do otrzymanego w drugim kroku równania (20b) i kolejno w następnych krokach do następnych
równań ze zwiększającą się ilością zer pod przekątną główną przetwarzanej macierzy; po (N-1) krokach otrzymujemy
5�N�11 5�N�12 5�N�13 & 5�N�15�A� 5�e�1 5�O�1
(1) (1) (1)
0 5�N�22 5�N�23 & 5�N�25�A� 5�O�2
5�e�2 (1)
(2) (2)
5�e�3 (2)
0 0 5�N�33 & 5�N�35�A� 5�O�3
�" = (20)
& & & & & & &
(5�A�-1)
(5�A�-1)
0 0 0 & 5�N�5�A�5�A� [5�e�5�A�]
[ ]
[5�O�5�[� ]
gdzie
(5�`� (5�`�-1) (5�`�-1) (5�`�) (5�`�-1) (
5�N�5�V�5�X�) = 5�N�5�V�5�X� - 5�Y�5�V�5�`� �" 5�N�5�`�5�X� , 5�O�5�V� = 5�O�5�V� - 5�Y�5�V�5�`� �" 5�O�5�`�5�`�-1) ,
(5�`�-1)
5�N�5�V�5�`�
(0
przy 5�Y�5�V�5�`� = , 5�`� = 1, 2, & , 5�A�-1, 5�V� = 5�`�+1, 5�`�+2, & , 5�A�, 5�X� = 5�V�-1, 5�V�, & 5�A� (przy czym 5�N�5�]�5�^�) = 5�N�5�]�5�^�).
(5�`�-1)
5�N�5�`�5�`�
Przykład obliczeń według tego algorytmu był podany w instrukcji do poprzedniego ćwiczenia warto z nim ponownie zapoznać się.
Poniżej przypominamy też pewien wariant metody Gaussa noszący nazwę metody LU. Polega ona na przyjęciu lewostronnej macierzy ze wzoru
(20) jako macierzy 5�|� (trójkątnej górnej) i utworzeniu dodatkowej macierzy 5�s� (trójkątnej dolnej) zbudowanej ze współczynników 5�Y�5�V�5�`�. Przykładowo dla
przypadku macierzy o wymiarach 4x4 (5�A� = 4) otrzymujemy
5�N�11 5�N�12 5�N�13 5�N�14
1 0 0 0
(1) (1) (1)
0 5�N�22 5�N�23 5�N�24 5�Y�21 1 0 0
5�|� = , 5�s� = (21)
(2) (2)
0 0 5�N�33 5�N�34 5�Y�31 5�Y�32 1 0
(3)
0 0 0 5�N�44 [5�Y�41 5�Y�42 5�Y�43 1]
[ ]
Macierze te, dla 5�A� e" 2, mają następujące pożyteczne właściwości
| | | | | |
5�s� " 5�|� = 5�h�, 5�s� = 1, 5�h� = 5�|� (22)
Wzory (23) i (24) wskazują na korzystniejszy od tradycyjnego sposób obliczania wyznacznika macierzy 5�h�; dodajmy, sposób wykorzystywany w
operacji det(A) MATLAB-a.
Poszukiwanie wektora rozwiązań 5ؙ� możemy teraz przedstawić w postaci
5�s� " 5�|� " 5ؙ� = 5؃� ,
(23)
5ؚ� = 5�s�-5�� " 5؃� ,
5ؙ� = 5�|�-5�� " 5ؚ� ,
gdzie 5ؚ� to pomocniczy wektor kolumnowy w procedurze (23).
Eliminację Gaussa prowadzącą do otrzymania macierzy LU można także zastosować nie do oryginalnego układu równań (15), tylko do układu
równań z przestawionymi wierszami. Najbardziej oczywista potrzeba przestawienia wierszy może w ramach algorytmu eliminacji wynikać z potrzeby
ominięcia sytuacji dzielenia przez zero, ale możemy także w ten sposób zwiększać dokładność obliczeń.
Wspomnianą optymalizację dokładności obliczeń w odniesieniu do układu równań (17) otrzymujemy w MATLAB-ie po zastosowaniu do dekom-
pozycji funkcji lu MATLAB-a z argumentami wyjściowymi 5�s�, 5�|�, 5�w�:
[L,U,P] = lu(A). (24)
Funkcja ta w miejsce macierzy 5�h� zwraca nam trzy macierze, które możemy podstawić do równania (17). Po podstawieniu, zapis tego równania przyj-
muje postać
5�w� " 5�s� " 5�|� " 5ؙ� = 5؃� . (25)
Zawierająca informację o przestawieniach wierszy macierz 5�w� ma następujące właściwości
5�w� = 5�w�-5��, | |
5�w� = 1 . (26)
Stąd rozwiązanie równania (17) możemy w tym przypadku uzyskać na podstawie wyrażenia
5�s� " 5�|� " 5ؙ� = 5؃�5�ę� (27)
gdzie
5؃�5�ę� = 5�w� " 5؃�.
Sposób zamiany kolejności wierszy daje się odczytać z macierzy 5�w�. Ściśle rzecz biorąc, niezerowy wyraz 5�]�5�V�5�W� = 1 tej macierzy informuje nas, że
nowym 5�V�-tym równaniem w zapisie (15) ma być stare 5�W�-te równanie tego zapisu.
Realizując poleceniem lu(A)opisaną powyżej procedurę stosujemy optymalizowaną dekompozycję [L, U, P] w sposób jawny, ale oprócz
tego dekompozycja taka jest automatycznie wykonywana niejawnie w ramach następujących poleceń Matlaba:
�� w przypadku równania postaci (16), gdy je rozwiązujemy poleceniem
x=A\b (28)
�� w przypadku obliczania wyznacznika macierzy A
DA = det(A) (29)
5. Pliki MATLABA do wspomagania analizy obwodów
Przyjęto podobny do stosowanego w programach typu SPICE sposób zadawania struktury i wartości elementów
obwodu czyli strukturę danych o obwodzie. Sposób ten w literaturze bywa nazywany tablicą połączeń.
C. Stefański, Materiały do ćw 6 Lab. TOiS 7/12 Ostatnia modyfikacja: maj 2016
Zakładamy, że wiersz danych ma postać tekstową o następującej strukturze:
ReNe nv1 nv2 val
gdzie: Re to rodzaj elementu (na razie jedno z {R, G, E, J}), Ne to nr elementu (wymaga się ciągłej numeracji elementów
poczynając od 1); nv1 i nv2 to pierwszy i drugi węzeł, między które załączono element (wymaga się ciągłej numeracji
węzłów poczynając od zera (węzeł odniesienia)), wreszcie val oznacza wartość liczbową elementu; dopuszcza się krot-
ności {p,n,u,m,k,M,G,T}.
Na jeden dwójnik przypada jeden wiersz tablicy połączeń. Przy pisaniu algorytmów założono, że nv1 oznacza węzeł
o umownie wyższym potencjale (umowny kierunek prądu od nv1 ku nv2).
Aby nie komplikować oprogramowania przyjęto, że liczba węzłów i gałęzi jest podawana na początku danych (pro-
fesjonalne programy takiego wymagania nie mają).
% cramer.m
function [x]=cramer(A,b)
%Metoda Cramera rozwiązywania układów równań liniowych
%Zapis układu równań A*x=y
N=rank(A);
for k=1:N
K=A(:,k); % Kolumna o numerze k w macierzy A na starcie
A(:,k)=b; % Macierze A , A ,& ,A
1 2 m
DK(k)=det(A); % Wyznaczniki D , D ,& ,D
1 2 m
A(:,k)=K; %Powrót do startowej postaci macierzy A
end
D=det(A); % Wyznacznik główny
x=DK./D; % Tablicowe dzielenie wyznaczników
%Koniec cramer.m
% dekompozycja.m
function [LA,UA]=dekompozycja(A)
% Dekompozycja (rozkład) macierzy A na dwie macierze
% trójkątną-dolną LA i trójkątną- górną UA
N=rank(A); % Określa liczbę niezależnych równań
el=eye(N); % Generuje macierz z "jedynkami" na głównej przekątnej
a=zeros(N); % Generuje macierz wypełnioną zerami
a=A; % Przepisuje elementy macierzy
for s=1:(N-1) % Litera s oznacza numer kroku w metodzie Gaussa
% Litery i oraz k oznaczają numery wiersza i kolumny
n=s+1;
for i=n:N
el(i,s)=a(i,s)/a(s,s)
for k=s:N
a(i,k)=a(i,k)-el(i,s)*a(s,k);
end
end
end
LA=el; % Macierz trójkątna dolna
UA=a; % Macierz trójkątna górna
% Koniec dekompozycja.m
% DaneOobwodzie.m
%Do Cw4
%***Program wspomagający wpisywanie tabeli połączeń elementów***
%
% Nową tabelę połączeń realizujemy poprzez edycję niniejszej
% tabeli połączeń
%
% Dane o obwodzie
TablicaPolaczen=[
"* to jest proba-test";
"E1 1 0 1";
"R2 2 1 1";
"* to jest proba-test";
"R3 2 0 1";
"*koniec";
"R4 3 2 1";
"R5 3 0 1";
"*koniec";
C. Stefański, Materiały do ćw 6 Lab. TOiS 8/12 Ostatnia modyfikacja: maj 2016
"J6 3 4 1";
"R7 4 0 1";
"R8 4 5 1";
"R9 5 0 1";
"R10 5 6 1";
"E11 6 0 1";
"koniec"];
%
% UzupelnijMacierzIncydencji.m
function UzupelnijMacierzIncydencji(A,d)
global A;
if (d(2)>0)
A(d(2),d(1))=1;
endif
if (d(3)>0)
A(d(3),d(1))=-1;
endif
endfunction
% UzupelnijWpisyEWmacierzachNiExc.m
function UzupelnijWpisyEWmacierzachNiExc(N,Exc,d)
global N;
global Exc;
if (d(1)>0)
N(d(1),d(1))=1;
Exc(d(1),1)=d(4);
else
print("4E. Blad danych\n");
close;
endif
endfunction
% UzupelnijWpisyGWmacierzachMiN.m
function UzupelnijWpisyGWmacierzachMiN(M,N,d)
global M;
global N;
if (d(1)>0)
M(d(1),d(1))=-1;
N(d(1),d(1))=d(4);
else
print("4G. Blad danych\n");
close;
endif
endfunction
% UzupelnijWpisyJWmacierzachMiExc.m
function UzupelnijWpisyJWmacierzachMiExc(M,Exc,d)
global M;
global Exc;
if (d(1)>0)
M(d(1),d(1))=1;
Exc(d(1),1)=d(4);
else
print("4J. Blad danych\n");
close;
endif
endfunction
% UzupelnijWpisyRWmacierzachMiN.m
function UzupelnijWpisyRWmacierzachMiN(M,N,d)
global M;
C. Stefański, Materiały do ćw 6 Lab. TOiS 9/12 Ostatnia modyfikacja: maj 2016
global N;
if (d(1)>0)
M(d(1),d(1))=d(4);
N(d(1),d(1))=-1;
else
print("4R. Blad danych\n");
close;
endif
endfunction
% EJRG.m
% Pojedynczy wiersz z danymi oznaczamy WDanych
%
% Zakładamy, że wiersz danych ma postać tekstową o następującej strukturze:
%
% ReNe nv1 nv2 val
% gdzie
% Re to rodzaj elementu (na razie jedno z {R, G, E, J})
% Ne to nr elementu; wymaga się ciągłej numeracji elementów poczynając od 1;
% nv1 i nv2 to pierwszy i drugi węzeł, między które załączono element;
% wymaga się ciągłej numeracji węzłów poczynając od zera (węzeł odniesienia)
% val oznacza wartość liczbową elementu; dopuszcza się krotności {p,n,u,m,k,M,G,T}
global A;
global M;
global N;
global Exc;
LE=input("Liczba elementow wynosi (numeracja od 1 do LE): LE=");
NN=input("Wezel o najwyzszym numerze to (numeracja od 0 do NN): NN=");
A=zeros(NN,LE);
M=zeros(LE,LE);
N=M;
Exc=zeros(LE,1);
T=zeros(LE+NN,LE+NN);
S=zeros(LE+NN,1);
RElem=0;
i=1;
k=0;
DaneOobwodzie;
while (abs(RElem-5)>0.1)
if (i>LE)
break;
else
k=k+1;
endif;
% printf("%d %d \n",i,k);
WDanych=TablicaPolaczen(k,:);
zp=WDanych;
X = strsplit(zp);
if (zp(1,1)=="R")
RElem=1;
elseif (zp(1,1)=="G")
RElem=2;
elseif (zp(1,1)=="E")
RElem=3;
elseif (zp(1,1)=="J")
RElem=4;
elseif (zp(1,1:6)=="koniec")
RElem=5;
elseif (zp(1,1)=="*")
RElem=6;
else
printf ("1. Nierozpoznany element\n");
close;
endif;
C. Stefański, Materiały do ćw 6 Lab. TOiS 10/12 Ostatnia modyfikacja: maj 2016
if (and(RElem==5,i printf("2a. Blad danych; nacisnij Ctrl+C\n");
keyboard;
endif;
if(and(size(X)<4, RElem<5))
printf ("2b. Blad danych\n");
close;
endif;
if (RElem<6)
b=strjoin(X(1,1));
a=b;
a(1)=[];
d(1)=str2num(a);
d(2)=str2num(strjoin(X(1,2)));
d(3)=str2num(strjoin(X(1,3)));
c = strrep(strrep (X(1,4), "m", "e-3"),"k","e3");
c = strrep(strrep (c, "u", "e-6"),"M","e6");
c = strrep(strrep (c, "n", "e-9"),"G","e9");
c = strrep(strrep (c, "p", "e-12"),"G","e12");
d(4)=str2num(strjoin(c));
i=i+1;
if(or(eq(d(2),d(3)), gt(d(2),NN), gt(d(3),NN), gt(d(1),LE)))
printf ("3. Blad danych\n");
close;
endif;
if (RElem==1)
UzupelnijWpisyRWmacierzachMiN(M,N,d);
endif;
if (RElem==2)
UzupelnijWpisyGWmacierzachMiN(M,N,d);
endif;
if (RElem==3)
UzupelnijWpisyEWmacierzachNiExc(N,Exc,d);
endif;
if (RElem==4)
UzupelnijWpisyJWmacierzachMiExc(M,Exc,d);
endif;
if (RElem <= 4)
UzupelnijMacierzIncydencji(A,d(1:3));
endif;
endif;
endwhile;
T(1:NN,1:LE)=A;
T(NN+1:NN+LE,1:LE)=M;
T(NN+1:NN+LE,LE+1:LE+NN)=N*A';
S(NN+1:NN+LE,1)=Exc;
x=T\S;
I=x(1:LE,1);
Vn=x(LE+1:LE+NN,1);
U=A'*Vn;
I',U',Vn'
%
6. Symulacja obwodu elektrycznego za pomocą programu SPICE (Micro5CAP)
Z wielu wariantów SPICE a w laboratorium jest do dyspozycji program Micro5CAP w wersji 10 lub 11 (dalej będziemy
na oznaczenie tego narzędzia używać skrótu MC).
C. Stefański, Materiały do ćw 6 Lab. TOiS 11/12 Ostatnia modyfikacja: maj 2016
W edytorze programu MC tworzymy wirtualny obwód złożony: ze zródeł niezależnych, zródeł sterowanych, rezy-
storów, cewek indukcyjnych i kondensatorów i w miarę potrzeby z innych elementów pochodzących z banku elementów
SPICE a. SPICE w zewnętrznym oglądzie raczej symuluje właściwości obwodu niż wspomaga poszukiwanie rozwiązań
równań opisujących taki obwód, choć w ukryciu przed użytkownikiem naprawdę te równania automatycznie tworzy
i jedno- lub wielokrotnie te równania rozwią-
.define N 30
zuje.
Obwód prądu stałego z rys. 1b, wprowa-
1
2
dzony do MC, wygląda jak na rys. 4.
R2
Edytor MC nie pozwala w wygodny sposób
2.22*N
określać kierunki odniesienia dla przepływu prą-
J5
E1
R3 R4
dów w gałęziach oraz nadawać numery porząd- 1.11*N
3.33/N 4.44/N
DC 5.55/N
kowe węzłów. We wprowadzonym do edytora
Rys. 4. Schemat obwodu z rys. 1b generowany przez program
programu obwodzie można poprzez kliknięcie
myszką w obrębie wybranego elementu zmienić Micro5CAP
w otwierającym się wtedy okienku wartość licz-
bową parametru (np. E1, J5, R2, R3). Do równoczesnego i sprawnego nadawania wartości wielu elementom zdefinio-
wano stałą symboliczną N (linia .define N 30). Od wartości tej stałej uzależnione są tu SEM zródła E1, SPM zródła
J5, oporności oporników R2, R3 i R4.
Jeżeli z kilku wariantów badania obwodu DC (prądu stałego) wybierzemy analizę Dynamic DC (Analysisą�Dy-
namic DC& ) otrzymamy napięcia węzłowe (V1, V2), czyli potencjały każdego z ponumerowanych węzłów względem
napięcia odniesienia (masa, ziemia) oraz prądy gałęziowe (I1, I2, & , I5). Na ekran wartości tych napięć i prądów można
wywołać wybierając odpowiednie ikonki z paska narzędziowego programu1.
Oprócz dotychczas pokazanych sposobów wyznaczania parametrów zródeł zastępczych Thevenina lub Nortona
można też posłużyć się następującymi podejściami, sugerowanymi rysunkiem 5, dość wygodnymi, gdy korzystamy z pro-
gramu MC. Mianowicie analiza Dynamic DC przy 5�=�5�a�5�R�5�`�5�a�=0 pozwala wyznaczyć 5�8�5�G� bezpośrednio, zaś przy 5�=�5�a�5�R�5�`�5�a�=-1 A licz-
bowa wartość napięcia 5�H�5�=�5�a�5�R�5�`�5�a� różni się liczbowo o 5�E�5�G� od liczbowej wartości 5�8�5�G�; podobnie analiza Dynamic DC przy
5�8�5�a�5�R�5�`�5�a� = 0 pozwala wyznaczyć 5�=�5�A� bezpośrednio, zaś przy 5�8�5�a�5�R�5�`�5�a� = -1 V liczbowa wartość prądu 5�<�5�8�5�a�5�R�5�`�5�a� różni się liczbowo
o 5�:�5�A� od liczbowej wartości 5�=�5�A�.
(a)
I =J �E
Etest N-G
N test
R �J E =0 �� I =J
T test test Etest N
G �E
N test
(a)
U =E �J D�I =-G �D�E
Jtest T-R Etest N test
T test
R J
T N
J =0 �� U =E D�E =-1V �� D�I =G �1V
test Jtest T test Etest N
G E
N test
J D�U =-R �D�J
test Jtest T test
E
T
D�J =-1A �� D�U =R �1A
test Jtest T
(b)
(b)
Rys. 5. Alternatywne sposoby wyznaczania parametrów zródeł zastępczych
Kontynuując myśl sugerowaną rysunkiem 5 możemy pobudzić badany obwód (ważne, by był to układ zbudowany
tylko z elementów E, J, R, ZS, WO) ze żródła 5�=�5�a�5�R�5�`�5�a� o przebiegu okresowym w postaci fali prostokątnej unipolarnej, o mi-
nimalnej wartości -1 A (minus jeden amper) i wartości maksymalnej 0 A. W przedziałach czasu, w których przebieg 5�=�5�a�5�R�5�`�5�a�
przyjmuje wartość zero, napięcie na zródle 5�=�5�a�5�R�5�`�5�a� wyniesie 5�8�5�G� (będzie szukanym napięciem rozwarcia), zaś zmiana war-
tości napięcia na zródle 5�=�5�a�5�R�5�`�5�a� przy zmianie prądu na wartość minimalną (-1 A) liczbowo będzie równa rezystancji Theve-
nina. Badanie to w Micro"Cap-ie przeprowadzić można stosując analizę TRANSIENT. Odpowiednie ilustracje rysunku 6
przybliżają ten pomysł dla przypadku, gdy szukamy zródła Thevenina dla układu bez R4 z punktu widzenia zacisków (2)
i (0).
.define N 30
2
R2
1
Jtest
2.22*N
J5
E1
R3
1.11*N
3.33/N
DC 5.55/N
1
W MC11 ikonki te wyglądają, jak pokazano tutaj .
C. Stefański, Materiały do ćw 6 Lab. TOiS 12/12 Ostatnia modyfikacja: maj 2016
Micro-Cap 10 Evaluation Version
ZnalezienieThevenin.cir
0.25
0.00
-0.25
-0.50
-0.75
-1.00
0.00u 5.00u 10.00u 15.00u 20.00u 25.00u
I(Jtest) (A)
T (Secs)
180.00m
320n,34.906822m
150.00m
.define N 30
120.00m
34.906822m
90.00m
110.815308m
R2
33.3
Jtest
2.22*N
60.00m
J5
E1
R3
30.00m
0.00u 5.00u 10.00u 15.00u 20.00u 25.00u
1.11*N
v(2) (V)
3.33/N
DC 5.55/N
T (Secs)
Rys. 6. Obwód, okienko nastaw analizy TRANSIENT i jej rezultaty dla wyznaczenia 5�y�5�{� i 5�l�5�{�
Na kolejnym rysunku porównano uzyskane w programie Micro"Cap napięcie na R4=4.4/N omów w obwodzie
oryginalnym (po lewej) i w obwodzie z zastępczym zródłem Thevenina (po prawej).
34.906822m
RT
33.3 19.96099m 19.96099m
2.22*N 110.815308m
ET
R4
34.906822m
1.11*N
3.33/N 4.44/N 4.44/N
DC 5.55/N
Rys. 7. Napięcie na R4 w badanym obwodzie i w obwodzie ze zródłem Thevenina
Sugeruje się by Czytelnik zaproponował analogiczną procedurę dla znalezienia 5�:�5�A� i 5�=�5�A� .
Pytania kontrolne
1. Na czym polega tzw. metoda LU rozwiązywania układu równań liniowych?
2. Co to jest macierz incydencji obwodu? Narysuj graf o zredukowanej macierzy incydencji
1 0 -1 1 0
]
5�4� = [-1 1 0 0 0
0 -1 1 0 -1
3. Podać szczegóły tworzenia tzw. zródła Thevenina dla wieloelementowego obwodu liniowego.
4. Podać szczegóły tworzenia tzw. zródła Nortona dla wieloelementowego obwodu liniowego.
5. Podać szczegóły metody rozwiązywania układu równań liniowych związaną z nazwiskiem Cramera.
6. W jaki sposób wprowadzamy dane badanego obwodu do programu Micro5CAP?
7. Na czym polega tzw. dekompozycja LU macierzy współczynników obwodu?
8. Opisz, jaką ma składnię, jak jest interpretowane i do czego można wykorzystać polecenie .define
9. Objaśnij poszczególne równania macierzowe: 5�h� �" 5�p� = 5��, 5�|� = 5�h�5�{� �" 5�}�, 5�t� �" 5�p� + 5�u� �" 5�|� = 5�l�5ؙ�5�"�. W szczególności
podaj również rozmiary macierzy i wektorów, gdy wiesz, że dotyczą one obwodu złożonego tylko z dwójników
(w liczbie 5�O�) połączonych w węzłach (węzłów jest 5�[�).
10. Co to jest tablica połączeń? Objaśnij szczegóły wiersza takiej tablicy na przykładzie elementu wybranego spo-
śród {E, J, R, G}.

Wyszukiwarka

Podobne podstrony:
InstrCw7 16 odblokowany
InstrukcjaDoCw5zLTOiS 16 odblokowany
Instr Cw2 16CS odblokowany
cw6 arkusz obliczeniowy przyklad
Scenariusz 16 Rowerem do szkoły
r 1 nr 16 1386694464
16 narrator
16 MISJA
Fakty nieznane , bo niebyłe Nasz Dziennik, 2011 03 16
990904 16

więcej podobnych podstron