Quelltext der Seite Systemidentifikation

'''Systemidentifikation''' (auch '''Systemidentifizierung''') ist die theoretische und/oder experimentelle Ermittlung der quantitativen Abhängigkeit der Ausgangs- von den Eingangsgrößen eines Systems. Dazu wird das System mit definierten Testsignalen (Sprung, Impuls, Rampe o.&nbsp;Ä.) angeregt und der Ausgang aufgezeichnet. Die zur mathematischen Auswertung angewandten Verfahren können deterministisch oder stochastisch sein.

== Theoretische Systemidentifikation ==
Bei der theoretischen Systemidentifikation erfolgt die Modellbildung auf der Grundlage von Bilanzgleichungen unter Berücksichtigung von [[Erhaltungssatz|Erhaltungssätzen]]. Das Ergebnis ist ein die Beziehung zwischen Ein- und Ausgangsgrößen beschreibendes [[Differentialgleichung]]ssystem. Im Fall eines [[Lineares zeitinvariantes System|linearen, zeitinvarianten Systems]] gilt:

:<math> y^{(n)} + \ldots + a_{1}y^{(1)} + a_{0}y = b_{m}u^{(m)} + \ldots + b_{1}u^{(1)} + b_{0}u </math>

und, da in diesem Fall die [[Laplace-Transformation]] durchführbar ist, gilt für die [[Übertragungsfunktion]]

:<math> G(s) = \frac{b_{m}s^{m} + \ldots + b_{1}s + b_{0}}{s^{n} + a_{n-1}s^{n-1} + \ldots + a_{1}s + a_{0}}</math>.

Sind alle Koeffizienten <math>a_n</math> und <math>b_m</math> bekannt, ist die Identifizierungsaufgabe gelöst. Ansonsten müssen die unbekannten Koeffizienten durch die experimentelle Systemidentifikation bestimmt werden.

== Experimentelle Systemidentifikation ==
[[Bild:SFP_Systemidentifikation.png|mini|Signalflussplan der experimentellen Systemidentifikation]]
Das System wird mit geeigneten [[Testsignal]]en ([[Heaviside-Funktion|Sprung]], [[Impulsfunktion|Impuls]], [[Rampenfunktion|Rampe]] o.&nbsp;a.) <math>u(t)</math> angeregt. Diese Signale werden außerdem einem mathematischen Modell, welches freie Parameter hat, zugeführt. Das Modell ist aus einer vorhergehenden theoretischen Prozessidentifikation bekannt. Das Modell kann entweder im Zeitbereich oder im [[Fourier-Analysis#Anwendungen|Frequenzbereich]] vorliegen. Aus beiden Ausgangssignalen (System und Modell) wird die Abweichung (Differenz) berechnet und von einem Gütekriterium in Form eines [[Funktional]]s bewertet. Das Ergebnis der Bewertung wird von einem Algorithmus benutzt, um die Parameter des Modells anzupassen. Dieser Prozess wird solange wiederholt, bis die gewünschte Güte erreicht ist.

Die [[Iterationsverfahren|iterative Anpassung]] der Modellparameter kann durch Unterstützung entsprechender Softwarewerkzeuge verkürzt werden.

== Wendetangentenverfahren ==
[[Bild:Tutg.svg|mini|Sprungantwort einer Regelstrecke mit Wendetangente zur Bestimmung von <math>T_u</math> und <math>T_g</math>.]]
Ein System mit Ausgleich und ohne Überschwingen hat einen Wendepunkt in der [[Sprungantwort]]. Es tritt bei Systemen mit der Reihenschaltung mehrerer Verzögerungselemente ([[PT1-Glied]]er) auf. Durch Anlegen der Tangente an den Wendepunkt können die Verzugszeit <math>T_u</math> und die Ausgleichszeit <math>T_g</math> bestimmt werden. Ziel ist, aus diesen experimentell ermittelten Werten des Systems die Zeitkonstanten der Übertragungsfunktion des Modells zu ermitteln. Die Modellfunktion muss ebenfalls ein System mit Ausgleich beschreiben und die Sprungantwortfunktion des Modells in analytischer Form vorliegen. 

=== Prinzip ===
Aus der [[Übertragungsfunktion]] des Modells
:<math>G(s) = \frac{x_a(s)}{x_e(s)}</math>
und der Sprungfunktion im [[Bildbereich]]
:<math>x_e(s) = \frac{K}{s}</math>
erhält man durch inverse [[Laplace-Transformation]] die Sprungantwortfunktion des Modells
:<math>x_a(t) = L^{-1}\{G(s) \cdot x_e(s)\}</math>.
Der Anstieg <math>m(t)</math> der Sprungantwortfunktion des Modells zum Zeitpunkt <math>t</math> ist
:<math>m(t) = \frac{dx_a(t)}{dt}</math>.
Der Wendezeitpunkt <math>t_w</math> errechnet sich aus der Bedingung
:<math>\frac{d^2x_a(t)}{dt^2}=0</math>.
Mit diesen Daten ist
:<math>x_T(t) = m(t_w) \cdot (t-t_w)+x_a(t_w)</math>
die Gleichung der Wendetangente des Modells. Weiterhin gilt für den Durchtrittspunkt der Wendetangente durch die Zeitachse
:<math>x_T(t_0)=0 = m(t_w) \cdot (t_0-t_w)+x_a(t_w)</math>
und durch den stationären Wert
:<math>x_T(t_1) = K = m(t_w) \cdot (t_1-t_w)+x_a(t_w)</math>.
Aus der Skizze ist zu erkennen, dass <math>m(t_w) = \frac{K}{T_g}</math>, <math>t_0 = T_u</math> und <math>t_1 = T_u+T_g</math> gilt.
Damit lautet die Beziehung zwischen den Parametern des Systems <math>T_u</math>, <math>T_g</math> und den Eigenschaften des Modells <math>t_w</math> und <math>x_a(t_w)</math>
:<math>T_u = t_w-T_g \cdot \frac{x_a(t_w)}{K}</math>.
Diese Beziehungen sind als Tabellen und Nomogramme für konkrete Modelle verfügbar.

=== Modell-Übertragungsfunktionen ===
Mit dem beschriebenen Verfahren können nur Funktionen mit zwei Zeitkonstanten oder bei <math>n</math> gleichen Zeitkonstanten die Anzahl <math>n</math> und die [[Zeitkonstante]] bestimmt werden. Als Übertragungsfunktionen werden unter anderem eingesetzt:
*Zwei unterschiedliche Zeitkonstanten <math>T_1</math> und <math>T_2</math>
:<math>G(s) = \frac{K}{(1+T_1 \cdot s) \cdot (1+T_2\cdot s)}</math>
*<math>n</math> gleiche Zeitkonstanten <math>T_1</math>
:<math>G(s) = \frac{K}{(1+T_1 \cdot s)^n}</math>
*Zwei gleiche Zeitkonstanten <math>T_1</math> und eine Zeitkonstante <math>T_2</math> 
:<math>G(s) = \frac{K}{(1+T_1 \cdot s)^2\cdot (1+T_2\cdot s)}</math>
*Bei großen Totzeiten oder Zeitverzögerungen geht auch
:<math>G(s) = \frac{K \cdot e^{-T_t \cdot s}}{(1+T_1\cdot s)}</math>

=== Beispiel ===
[[Bild:Ident_pt3.png|mini|Vergleich gemessener und identifizierter Sprungantwort.]]

Für den Fall <math>n</math> gleiche Zeitkonstanten <math>T</math>
:<math>G(s) = \frac{K}{(1+T\cdot s)^n}</math>
lautet die Sprungantwortfunktion
:<math>x_a(t) = K \cdot \left\{1-\left[\sum_{m=0}^{n-1}\frac{(\frac{t}{T})^m}{m!} \right]\cdot e^{-t/T}\right\}</math>
Der Anstieg der Wendetangente ist
:<math>m(t)=\frac{dx_a(t)}{dt}=\left( \sum_{m=0}^{n-1}\frac{t^m}{T^{m+1}\cdot m!}-\sum_{m=1}^{n-1}\frac{m\cdot t^{m-1}}{T^m\cdot m!} \right)\cdot e^{-\frac{t}{T}} </math>

und aus
:<math>\frac{d^2x_a(t)}{dt^2} = \left( -\sum_{m=2}^{n-1} \frac{m\cdot (m-1)\cdot t^{m-2}}{T^m\cdot m!} +2\cdot \sum_{m=1}^{n-1} \frac{m\cdot t^{m-1}}{T^{m+1}\cdot m!} -\sum_{m=0}^{n-1}\frac{t^m}{T^{m+2}\cdot m!} \right) \dot e^{-\frac{t}{T}} </math>
folgt mit
:<math>\frac{d^2x_a(t)}{dt^2}=0</math>
nach umindizieren der Summen
:<math>{-\sum_{m=0}^{n-3}\frac{(\frac{t}{T})^m}{m!}+2\cdot \sum_{m=0}^{n-3}\frac{(\frac{t}{T})^m}{m!}+2\cdot \frac{(\frac{t}{T})^{n-2}}{(n-2)!}-\sum_{m=0}^{n-3}\frac{(\frac{t}{T})^{m-1}}{(m-1)!}-\frac{(\frac{t}{T})^{n-2}}{(n-2)!}-\frac{(\frac{t}{T})^{n-1}}{(n-1)!}=0}</math>
und
:<math>\frac{(\frac{t}{T})^{n-2}}{(n-2)!}=\frac{(\frac{t}{T})^{n-1}}{(n-1)!}</math>
die Wendezeit
:<math>t_w = (n-1)\cdot T</math>.
Der Anstieg im Wendepunkt ist
:<math>m(t_w) = K \cdot \frac{e^{-(n-1)}\cdot (n-1)^{n-2}}{T\cdot (n-2)!}</math>
und
:<math>x_a(t_w) = K \cdot \left\{1-\left[\sum_{m=0}^{n-1}\frac{(n-1)^m}{m!} \right]\cdot e^{-(n-1)}\right\}</math>

Die Wendetangentenkonstruktion liefert mit <math>T_g=\frac{K}{m(t_w)}</math> die numerisch auszuwertenden Beziehungen
:<math>\frac{T_g}{T} = \frac{(n-1)!}{(n-1)^{n-1}}\cdot e^{n-1}</math>
und mit <math>T_u = t_w-T_g\cdot \frac{x_a(t_w)}{K}</math>
:<math>\frac{T_u}{T} = n-1-\frac{(n-1)!}{(n-1)^{n-1}}\cdot\left[e^{n-1}-\sum_{m=0}^{n-1}\frac{(n-1)^m}{m!}\right]</math>.

Das [[Scilab]]-Script

 N=10;
 printf("\n");
 printf(" n     |Tg/T  |Tu/T  |Tg/Tu\n");
 printf("----------------------------\n");
 for n=2:N,
  su=0;
  for i=0:n-1,su=su+(n-1)^i/factorial(i);end;
  fa=factorial(n-1)*exp(n-1)/(n-1)^(n-1);
  fu=n-1-factorial(n-1)*(exp(n-1)-su)/(n-1)^(n-1);
  printf(" %-5d |%-5.3f |%-5.3f |%-5.3f\n",n,fa,fu,fa/fu);
 end;
 printf("----------------------------\n");

erzeugt eine Tabelle der Werte nach den oben angegebenen Beziehungen.

 n     |Tg/T  |Tu/T  |Tg/Tu
 ----------------------------
 2     |2.718 |0.282 |9.649
 3     |3.695 |0.805 |4.587
 4     |4.463 |1.425 |3.131
 5     |5.119 |2.100 |2.437
 6     |5.699 |2.811 |2.027
 7     |6.226 |3.549 |1.754
 8     |6.711 |4.307 |1.558
 9     |7.164 |5.081 |1.410
 10    |7.590 |5.869 |1.293
 ----------------------------

Aus der gemessenen Sprungantwort wurden die Werte <math>T_u=1{,}030</math>, <math>T_g=5{,}183</math> und <math>T_g/T_u=5{,}031</math> mit einem Scilab-Script numerisch ermittelt. Aus der Tabelle folgt <math>n=3</math>, <math>T_g/T=3{,}695</math> und <math>T_u/T=0{,}805</math>. Weil die aus <math>T_g</math> und <math>T_u</math> berechneten Zeitkonstanten unterschiedlich sind, wurde die mittlere Zeitkonstante zu <math>T=1{,}341</math> berechnet. Damit gilt für die Übertragungsfunktion des Modells
:<math>G(s)=\frac{1}{(1+1{,}341s)^3}</math>.

== Anwendungen ==
{{Überarbeiten}}
Ein wesentliches Hilfsmittel der Systemidentifikation ist die [[lineare Regression]]sanalyse. Man setzt hierbei als [[funktionale Abhängigkeit]] eine Linearkombination willkürlich gewählter Ansatzfunktionen an. Jede Ansatzfunktion ist arithmetischer Ausdruck der verursachenden physikalischen Größen. Die verursachte physikalische Größe wird errechnet, indem man jede Ansatzfunktion mit einem zunächst unbekannten Koeffizienten multipliziert und zum Ergebnis addiert. Die Koeffizienten werden dergestalt bestimmt, dass die mittlere quadratische Abweichung des gemessenen vom berechneten Ergebnis minimal wird. Dies bedeutet: die partielle Ableitung der mittleren quadratischen Abweichung nach jedem einzelnen Koeffizienten muss null sein. Daraus ergibt sich ein lineares inhomogenes Gleichungssystem zur Bestimmung der Koeffizienten.
Die Matrix des Gleichungssystems besteht aus Produkten aus je zwei Ansatzfunktionen, gemittelt über alle Messungen. Die rechte Seite des Gleichungssystems besteht aus Produkten der verursachten Größe mit je einer Ansatzfunktion, gemittelt über alle Messungen. Bei der schrittweise durchgeführten linearen Regressionsanalyse wird iterativ bestimmt, welche Reihenglieder den meisten und welche den geringsten Einfluss auf die Genauigkeit haben, und die Reihenglieder ohne nennenswerten Einfluss werden weggelassen.

Anstelle der schrittweise erfolgenden linearen Regressionsanalyse kann bei vielen Problemstellungen alternativ ein mehrlagiges [[Perzeptron]] (engl. [[multi-layer perceptron]], MLP) verwendet werden, was häufig mit dem Oberbegriff [[Künstliches neuronales Netz|Neuronale Netze]] bezeichnet wird. 

Die Systemidentifikation kommt beispielsweise in der [[Strömungsmechanik]] zum Einsatz, sei es um Widerstand und Auftrieb eines Profils zu berechnen oder sei es, um [[manövrieren]]de Schiffe numerisch zu simulieren. Ein anderes Anwendungsgebiet ist die [[Schwingung]]stechnik, wo mit Übertragungsfunktionen (engl.: RAO = response amplification operator) berechnet wird, mit welcher Vergrößerung und Phasenverschiebung ein schwingfähiges System auf die einzelnen Frequenzen der Schwingungsursache reagiert.

In der Luftfahrt wird die Systemidentifizierung beispielsweise eingesetzt, um die aerodynamischen Parameter zu ermitteln, die aus analytischen Verfahren und Windkanalversuchen häufig nur ungenau bekannt sind. Hier ist die „Quad-M“-Systematik gebräuchlich, welche als die wesentlichen Teilaufgaben der Systemidentifizierung die Manöver (d.&nbsp;h. die Systemanregung), die Messung der Systemreaktion, das mathematische Modell des Systems sowie die Methoden zur Parameterschätzung angibt.

== Literatur ==
* {{Literatur|Autor=Heinz Unbehauen|Titel=Regelungstechnik I|Verlag=Vieweg|Jahr=1997|ISBN=3-528-83332-7}}
* {{Literatur|Autor=Holger Lutz, [[Wolfgang Wendt]]|Titel=Taschenbuch der Regelungstechnik mit MATLAB und Simulink|Auflage=12.|Verlag=Europa-Lehrmittel|Jahr=2021|ISBN=978-3-8085-5870-6}}
* {{Literatur|Autor=Lennart Ljung|Titel=System Identification: theory for the User|Auflage=2.|Verlag=Prentice Hall|Ort=Upper Saddle River|Jahr=2006|ISBN=978-0-13-656695-3}}
* {{Literatur|Autor=Oliver Nelles|Titel=Nonlinear System Identification|Auflage=1.|Verlag=Springer|Jahr=2000|ISBN=978-3540673699}}
* {{Literatur|Autor=[[Rolf Isermann]], Marco Münchhof|Titel=Identification of Dynamic Systems – An Introduction with Applications|Auflage=1.|Verlag=Springer|Jahr=2010|ISBN=978-3540788782}}
* {{Literatur|Autor=Ravindra V. Jategaonkar |Titel=Flight Vehicle System Identification: A Time-Domain Methodology|Auflage=2.|Verlag=American Institute of Aeronautics and Astronautics, Inc.|Jahr=2015|ISBN=978-1624102783}}

[[Kategorie:Systemtheorie (Kybernetik)]]
[[Kategorie:Regelungstheorie]]