Das Entfernen struktureller Singularitäten mittels Pantelides Algorithmus

Das Entfernen struktureller Singularitäten mittels Pantelides Algorithmus • Diese Vorlesung stellt ein Verfahren vor, welches dazu verwendet werden kann, strukturelle Singularitäten in systematischer und algorithmischer Weise aus einem Modell zu entfernen. Das Verfahren wird Pantelides Algorithm genannt. • Beim Algorithmus von Pantelides handelt es sich um ein symbolisches Indexreduktionsverfahren. 3. November, 2004 Anfang Präsentation

Übersicht • Structurelle Singularitäten und der Strukturdigraph • Pantelides Algorithmus 3. November, 2004 Anfang Präsentation

Strukturelle Singularitäten: Ein Beispiel I Wir stellen ein Modell unter Verwendung der Ströme, Spannungen und Potentiale auf. Die Maschengleichungen werden daher ignoriert. v 0 i. R v 3 i L 2 i. L 1 i. C v 2 I 1 v 1 I 3 Wir haben 7 Netzwerkkomponenten plus die Erde, somit 2 7 + 1 = 15 Gleichungen. Dazu kommen vier Knoten, die zu 3 zusätzlichen Gleichungen führen. Somit erwarten wir 18 Gleichungen in 18 Unbekannten. Die Spannungen werden bei passiven Komponenten in die gleiche Richtung positiv normiert wie die Ströme. Bei aktiven Komponenten (Quellen) ist es umgekehrt. 3. November, 2004 Anfang Präsentation

Strukturelle Singularitäten: Ein Beispiel II I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 /dt i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 3. November, 2004 02 I 2 03 I 3 04 u. R 05 i. R 06 u. L 1 07 di. L 1 /dt 08 u. L 2 09 I 1 di. L 2 /dt 10 i. C 11 du. C /dt 12 v 0 13 v 1 14 v 2 15 v 3 16 u 1 17 u 2 18 1: 2: 3: 4: 5: 6: 7: 8: 01 u 3 Anfang Präsentation

Strukturelle Singularitäten: Ein Beispiel III 01 02 03 I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 /dt i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 3. November, 2004 I 2 03 I 3 u. R 05 i. R 06 u. L 1 07 di. L 1 /dt 08 u. L 2 09 02 04 04 I 1 di. L 2 /dt 10 i. C 11 du. C /dt 12 v 0 13 v 1 14 v 2 15 v 3 16 u 1 17 u 2 18 1: 2: 3: 4: 5: 6: 7: 8: 01 u 3 Anfang Präsentation 13 14 15 16 17 18

Strukturelle Singularitäten: Ein Beispiel IV 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 /dt i. C = C · du. C /dt v 0 = 0 3. November, 2004 03 I 3 u. R i. R u. L 1 07 di. L 1 /dt 08 u. L 2 09 di. L 2 /dt 10 i. C 11 du. C /dt 12 v 0 v 1 14 v 2 15 v 3 16 u 1 17 Beschränkungsgleichung I 2 06 05 02 05 I 1 04 04 01 13 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 01 02 03 u 2 18 u 3 Anfang Präsentation Alle Verbindungen sind blau 13 14 15 16 17 18

Das Färben des Strukturdigraphen • Der Algorithmus zum Färben des Strukturdigraphen ist völlig äquivalent zum bisher angewandten Verfahren der Kausalisierung von Gleichungen. • Eine Implementierung des Verfahrens unter Verwendung eines Computerprogramms wird vermutlich den Digraphen vorziehen, da dieser direkt auf Datenstrukturen gängiger Programmiersprachen abgebildet werden kann. • Für das menschliche Auge ist das Färben der Gleichungen vermutlich lesbarer. Darum wird in der Vorlesung in Zukunft dem Färben der Gleichungen wieder ein Vorzug eingeräumt. • Das vertikale Sortieren kann gleichzeitig durch Umnummerieren der Gleichungen erfolgen. 3. November, 2004 Anfang Präsentation

Der Algorithmus von Pantelides I • Wenn eine Beschränkungsgleichung gefunden wurde, muss diese abgeleitet werden. • Beim Algorithmus von Pantelides wird die abgeleitete Beschränkungsgleichung dem Gleichungssystem zugefügt. • Somit hat das Gleichungssystem nun eine überzählige Gleichung. • Um die Anzahl von Gleichungen und Unbekannten wieder auszugleichen, wird ein mit der Beschränkungsgleichung verbundener Integrator eliminiert. 3. November, 2004 Anfang Präsentation

Der Algorithmus von Pantelides II dx dt unbekannt x bekannt, da Zustandsvariable dx dt x unbekannt Eine zusätzliche Unbekannte wurde durch die Elimination des Integrators geschaffen. x und dx sind nun algebraische Variablen, für die Gleichungen gefunden werden müssen. 3. November, 2004 dx unbekannt Anfang Präsentation

Der Algorithmus von Pantelides III • Beim Ableiten der Beschränkungsgleichung kann es geschehen, dass zusätzliche neue Variablen erzeugt werden, z. B. v dv, wobei v eine algebraische Variable ist. • Nachdem v bereits blau war (sonst wäre es ja keine Beschränkungsgleichung), existiert eine andere Gleichung, die v ermittelt. • Diese Gleichung muss nun ebenfalls abgeleitet werden. • Das Ableiten zusätzlicher Gleichungen hört erst dann auf, wenn keine neuen Variablen mehr erzeugt werden. 3. November, 2004 Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel I eliminierter Integrator 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 /dt i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 d. I 1 + di. C + di. L 2 + d. I 3 = 0 neu eingeführte Variabeln 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 3. November, 2004 Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel II 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 /dt i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel III 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 neu eingeführte Variable 23: d. I 2 = df 2(t)/dt Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel IV 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 23: d. I 2 = df 2(t)/dt Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel V 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 23: d. I 2 = df 2(t)/dt Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel VI 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 23: d. I 2 = df 2(t)/dt Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel VII 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 23: d. I 2 = df 2(t)/dt Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel VIII 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 Es findet sich ein algebraisch gekoppeltes System mit 7 Gleichungen in 7 Unbekannten. 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 23: d. I 2 = df 2(t)/dt Wahl 3. November, 2004 Anfang Präsentation

Der Pantelides Algorithmus: Ein Beispiel IX 1: 2: 3: 4: 5: 6: 7: 8: I 1 = f 1(t) I 2 = f 2(t) I 3 = f 3(t) u. R = R · i R u. L 1 = L 1 · di. L 1 /dt u. L 2 = L 2 · di. L 2 i. C = C · du. C /dt v 0 = 0 9: u 1 = v 0 – v 1 10: u 2 = v 3 – v 2 11: u 3 = v 0 – v 1 12: u. R = v 3 – v 0 13: u. L 1 = v 2 – v 0 14: u. L 2 = v 1 – v 3 15: u. C = v 1 – v 2 16: i. C = i. L 1 + I 2 17: i. R = i. L 2 + I 2 18: I 1 + i. C + i. L 2 + I 3 = 0 19: d. I 1 + di. C + di. L 2 + d. I 3 = 0 3. November, 2004 20: d. I 1 = df 1(t)/dt 21: d. I 3 = df 3(t)/dt 22: di. C = di. L 1 /dt + d. I 2 23: d. I 2 = df 2(t)/dt Anfang Präsentation

Zusammenfassung I • Zunächst findet man einen vollständigen Satz a-kausaler Algebrodifferentialgleichungen. • Auf diesen Satz wendet man den Färbealgorithmus von Tarjan an. • Falls sich eine Gleichung findet, die völlig blau gefärbt ist, ist das System strukturell singulär. • Das strukturell singuläre System wird mittels Anwendung des Algorithmus von Pantelides regulär gemacht. • Es mag nötig sein, den Pantelides Algorithmus mehrfach anzuwenden. 3. November, 2004 Anfang Präsentation

Zusammenfassung II • Auf das nunmehr reguläre Algebrodifferentialgleichungssystem wendet man wiederum den Färbealgorithmus von Tarjan an. • Falls der Algorithmus ins Stocken kommt, hat man es mit einem algebraisch gekoppelten System zu tun. Nach der Anwendung des Pantelides Algorithmus zur Indexreduktion eines strukturell singulären Systems treten algebraische Schleifen häufig auf. • Dieses System muss nun zunächst weiterverarbeitet werden. Das Aufschneideverfahren, welches bereits vorgestellt wurde, ist ein mögliches Verfahren, um mit solchen algebraisch gekoppelten Systemen umzugehen. 3. November, 2004 Anfang Präsentation

Referenzen • Cellier, F. E. and H. Elmqvist (1993), “Automated formula manipulation supports object-oriented continuous-system modeling, ” IEEE Control Systems, 13(2), pp. 28 -38. • Pantelides, C. C. (1988), “The consistent initialization of differential-algebraic systems, ” SIAM Journal Scientific Statistical Computation, 9(2), pp. 213 -231. 3. November, 2004 Anfang Präsentation