Wert- und Zeitdiskrete Systeme

Wert- und Zeitdiskrete Systeme

Vorbemerkungen

Signale in kontinuierlicher und diskreter Zeit

kontinuierliche (konti.) Zeit

  • Zeit ist kontinuierliche Variable
  • Signal s(t)s(t) nimmt bestimmten Wert s(t)s^*(t^*) für beliebig kurze Zeitspanne an
  • Zwischen zwei beliebigen Zeitpunkte t1t_1 und t2t_2 liegen unendlich viele Zeitpunkt t1tt2t_1 \leq t \leq t_2
  • Werte könne kontinuierlich oder diskret sein
截屏2022-06-03 12.24.25 截屏2022-06-03 12.24.51
  • Kontinuierlich in Zeit und Wert \rightarrow analoges Signal

Diskrete Zeit

  • Diskrete Zeitpunkt tk,kZt_k, k \in \mathbb{Z}

    sk:=s(tk) s_k := s(t_k)

    wertdiskrete_systeme-deskrete_zeit

  • Zeitliche Anordnung der tkt_k ist beliebig, aber in viele Fällen äquidistant

    tk=kΔkZ t_k = k \cdot \Delta \quad k \in \mathbb{Z}
  • Wert können kontinuierlich oder diskret sein

    • Diskret in Zeit und Ort \rightarrow digitales Signal
    截屏2022-06-03 14.54.54 截屏2022-06-03 14.55.06
  • Signale können inhärent zeitdiskret sein, oder aus Abtastung kontinuierliche Signale entstehen.

Kategoriale und Kardinale Variablen

截屏2022-06-03 15.03.01

Kategoriale Variable

Nominal
  • The nominal scale is made up of pure labels.

    • The only meaningful question to ask is whether two variables have the same value: the nominal scale only allows to compare two values w.r.t. equivalence.
    • There is no meaningful transformation besides relabeling.
    • No empirical operation is permissible, i.e., there is no mathematical operation of nominal features that is also meaningful in the material world.
  • A typical example is the sex of a human.

    • The two possible values can be either written as “f” vs. “m,” “female” vs. “male”. The labels are different, but the meaning is the same.
  • Although nominal values are sometimes represented by digits, one must not interpret them as numbers.

    • For example, the postal codes used in Germany are digits, but there is no meaning in, e.g., adding two postal codes.
    • Similarly, nominal features do not have an ordering, i.e., the postal code 12345 is not “smaller” than the postal code 56789. Of course, most of the time there are options for how to introduce some kind of lexicographic sorting scheme, but this is purely artificial and has no meaning for the underlying objects. With respect to statistics, the permissible average is not the mean (since summa- tion is not allowed) or the median (since there is no ordering), but the mode, i.e., the most common value in the dataset.

Ordinal
  • The ordinal scale allows comparing values w.r.t. equivalence and rank.
    • Any transformation of the domain must preserve the order, which means that the transformation must be strictly increasing.
    • But there is still no way to add an offset to one value in order to obtain a new value or to take the difference between two values.
  • Example: school grades.
    • In the German grading system, the grade 1 (“excellent”) is better than 2 (“good”), which is better than 3 (“satisfactory”) and so on.
      • But quite surely the difference in a student’s skills is not the same between the grades 1 and 2 as between 2 and 3, although the “difference” in the grades is unity in both cases.
      • In addition, teachers often report the arithmetic mean of the grades in an exam, even though the arithmetic mean does not exist on the ordinal scale. In consequence, it is syntactically possible to compute the mean, even though the result, e.g., 2.47 has no place on the grading scale, other than it being “closer” to a 2 than a 3. The Anglo-Saxon grading system, which uses the letters “A” to “F”, is somewhat immune to this confusion.
  • The correct average involving an ordinal scale is obtained by the median.

Kardinale Variable

Interval
  • The interval scale allows adding an offset to one value to obtain a new one, or to calculate the difference between two values—hence the name.
  • However, the interval scale lacks a naturally defined zero. Values from the interval scale are typically represented using real numbers, which contains the symbol “0,” but this symbol has no special meaning and its position on the scale is arbitrary. For this reason, the scalar multiplication of two values from the interval scale is meaningless. Permissible transformations preserve the order, but may shift the position of the zero.
Verhältnis
  • The ratio scale has a well defined, non-arbitrary zero, and therefore allows calculating ratios of two values.

    • This implies that there is a scalar multiplication and that any transformation must preserve the zero.
  • Many features from the field of physics belong to this category and any transformation is merely a change of units.

Absolut
The absolute scale shares these properties, but is equipped with a natural unit and features of this scale can NOT be negative. In other words, features of the absolute scale represent counts of some quantities. Therefore, the only allowed transformation is the identity.

Wertdiskrete Systeme

Statische Systeme

Ein-/Ausgang: Zufallsvariable uku_k (Eingang) und yky_k (Ausgang), kN0k \in \mathbb{N}_0

wertdiskrete_systeme-statistische_systeme.drawio

uku_k und yky_k sind wertdiskret, wobei o.B.d.A

uk{1,2,,p}yk{1,2,,M} \begin{array}{l} u_{k} \in\{1,2, \cdots, p\} \\ y_{k} \in\{1,2, \ldots, M\} \end{array}

Stochastische Abhängigkeit yky_k von uku_k:

P(yk=iuk=j)j{1,,p},i{1,,m} P\left(y_{k}=i \mid u_{k}=j\right) \qquad j \in\{1, \cdots, p\}, i \in\{1, \ldots, m\}

Anordnung der Wahrscheinlichkeit in Matrix AkA_k:

Ak=(P(yk=1uk=1)P(yk=Muk=1)P(yk=1uk=P)P(yk=Muk=P)) \mathbf{A}_{k}=\left(\begin{array}{ccc} P\left(y_{k}=1 \mid u_{k}=1\right) & \cdots & P\left(y_{k}=M \mid u_{k}=1\right) \\ \vdots & & \vdots \\ P\left(y_{k}=1 \mid u_{k}=P\right) & \cdots & P\left(y_{k}=M \mid u_{k}=P\right) \end{array}\right)
  • Elemente 0\geq 0

  • Zeilensumme =1= 1

  • Auftrittswahrscheinlichkeit als Vektoren:

    ηku=(P(uk=1)P(uk=2)P(uk=P))ηky=(P(yk=1)P(yk=2)P(yk=M)) \eta_{k}^{u}=\left(\begin{array}{c} P\left(u_{k}=1\right) \\ P\left(u_{k}=2\right) \\ \vdots \\ P\left(u_{k}=P\right) \end{array}\right) \qquad \eta_{k}^{y}=\left(\begin{array}{c} P\left(y_{k}=1\right) \\ P\left(y_{k}=2\right) \\ \vdots \\ P\left(y_{k}=M\right) \end{array}\right)

Berechnung von ηky\eta_k^y aus ηku\eta_k^u (in Vektor-Matrix-Form):

ηky=Akηku \eta_{k}^{y}=\mathbf{A}_{k}^{\top} \eta_{k}^{u}
Details
P(yk=i)=j=1PP(yk=i,uk=j)=j=1pP(yk=iuk=j)P(uk=j) \begin{aligned} P\left(y_{k}=i\right) &=\sum_{j=1}^{P} P\left(y_{k}=i, u_{k}=j\right) \\\\ &=\sum_{j=1}^{p} P\left(y_{k}=i \mid u_{k}=j\right) \cdot P\left(u_{k}=j\right) \end{aligned}

Spezialfall: uk=ju_k = j^* ist bekannt, also

P(uk=j)=1P(uk=j)=0j=1,M,jj \begin{array}{l} P\left(u_{k}=j^{*}\right)=1 \\ P\left(u_{k}=j\right)=0 \quad j=1, \cdots M, j \neq j^{*} \end{array} P(yk=i)=j=1pp(yk=iuk=j)P(uk=j)=P(yk=iuk=j) \begin{aligned} \Rightarrow \quad P\left(y_{k}=i\right) &=\sum_{j=1}^{p} p\left(y_{k}=i \mid u_{k}=j\right) P\left(u_{k}=j\right) \\ &=P\left(y_{k}=i \mid u_{k}=j^{*}\right) \end{aligned}

In Vektor-Matrix-Form:

ηky=Ak(j,:)die jte Zeile von Ak=(P(yk=1uk=j)P(yk=Muk=j)) \eta_{k}^{y}={\underbrace{\mathbf{A}_{k}\left(j^{*}, :\right)}_{\text{die } j^*-\text{te Zeile von } A_k}}^\top=\left(P\left(y_{k}=1 \mid u_{k}=j^{*}\right) \cdots P\left(y_{k}=M \mid u_{k}=j^{*}\right)\right)^{\top}

Dynamische Systeme

  • Der aktuellen Ausgang yky_k ist abhängig von
    • dem aktuellen Eingang uku_k
    • dem aktuellen Zustand xkx_k
  • Aufteilung des dynamischen Systems in zwei Teile
    • Systemabbildung (dynamischer Teil): beschreibt zeitliche Entwicklung des Zustands xkx_k
    • Messabbildung (statischer Teil): beschreibt die Abbildung des Ausgang yky_k von Zustand xkx_k (und evtl. von aktuellem Eingang uku_k)

Systemabbildung

  • Zufallsvariable xk,kN0x_k, k \in \mathbb{N}_0 mit xk{1,2,,N}x_k \in \{1, 2, \dots, N\}

  • Entwicklung des Zustands xkx_k bescrhieben ducrch

    P(xk+1=ixk,,x1,x0,uk) P(x_{k+1}=i | x_k, \dots, x_1, x_0, u_k)

    (uku_k oft explizit forgelassen)

Definition

Bei xkx_k handelt es sich um eine Markov-Ketter (erster Ordnung), falls gilt

P(xk+1=ixk,,x1,x0,uk)=P(xk+1=ixk,uk) P\left(x_{k+1}=i \mid x_{k}, \ldots, x_{1}, x_{0}, u_{k}\right)=P\left(x_{k+1}=i \mid x_{k}, u_{k}\right)
  • Die zukünftige Entwicklung xk+1x_{k+1} ist bedingt unabhängig von vergangen Zuständen xk1,,x1,x0x_{k-1}, \dots, x_1, x_0, falls aktueller Zustand xkx_k bekannt ist

  • Vereinfachte Übergangswahrscheinlichkeit

    P(xk+1=jxk=i) P(x_{k+1} = j| x_k = i)

Definition

Eine Markov-Kette wird als Zeithomogen oder allg. als zeitinvariant bezeichnet, falls die Übergangswahrscheinlichkeit nicht von Zeitindex abhängen, d.h. es gilt

P(xk+1=jxk=i)=A(i,j) P\left(x_{k+1}=j \mid x_{k}=i\right)=\mathbf{A}(i, j)

Übergangsmatrix (zeithomogen):

A=(A(1,1)A(1,2)A(1,N)A(2,1)A(2,2)A(2,N)A(N,1)A(N,2)A(N,N)) \mathbf{A}=\left(\begin{array}{cccc} A(1,1) & A(1,2) & \ldots & A(1, N) \\\\ A(2,1) & A(2,2) & \cdots & A(2, N) \\\\ \vdots & \vdots & & \vdots \\\\ A(N, 1) & A(N, 2) & \cdots & A(N, N) \end{array}\right)

Definition

Eine quadratische Matrix A\mathbf{A} heißt Markov-Matrix, falls

  • Alle Elemente nicht-negative sind

    A(i,j)0 fu¨i,j1,,N A(i, j) \geq 0 \quad \text{ für } i, j \in \\{1, \dots, N\\}
  • Die Zeilensumme gleich 1

    i=1NA(i,j)=1fu¨i1,,N \sum_{i=1}^{N} A(i, j)=1 \quad \text{für } i \in \\{1, \dots, N\\}

Graphische Darstellung einer Markov-Kette:

z.B. N=2,xk1,2N=2, x_k \in \\{1, 2\\}

wertdiskrete_systeme-Markov_Kette.drawio

Messabbildung

  • Zustand typischerweise NICHT direkt verfügbar (latente Variable)

  • Messabbildung vom Zustand xkx_k und dem aktuelle Eingang uku_k auf aktuelle Ausgang yky_k

    P(yk=jxk=i,uk=m) P\left(y_{k}=j \mid x_{k}=i, u_{k}=m\right)
    • uku_k oft explizit forgelassen
  • Zeithomogen (allg. zeitinvariant)

    P(yk=jxk=i)=B(i,j) P\left(y_{k}=j|x_{k}=i\right)=B(i, j)
  • Messe-/Beobachtungsmatrix

    B=[B(1,1)B(1,M)B(N,1)B(N,M)] \mathbf{B}=\left[\begin{array}{ccc} B(1,1) & \cdots & B(1, M) \\ \vdots & & \vdots \\ B(N, 1) & \cdots & B(N, M) \end{array}\right]

Gesamtes Dynamisches System

Hidden Markov Model

  • Zustand

    • Wert xk,k=1,2,x_k, k=1,2,\dots
    • Verteilung ηkx,k=1,2,\eta_k^x, k=1,2,\dots
  • Initialer Zustand

    • Wert x0x_0
    • Verteilung η0x\eta_0^x
  • Eingänge

    • Werte uk,k=0,1,u_k, k=0,1,\dots
    • Verteilung ηku,k=0,1,\eta_k^u,k=0,1,\dots
  • Ausgänge

    • Werte yk,k=0,1,y_k, k=0,1,\dots
    • Verteilung ηky,k=0,1,\eta_k^y,k=0,1,\dots
  • Systemabbildung Ak\mathbf{A}_k

  • Messabbildung Bk\mathbf{B}_k

Graphische Darstellung

  • Ausgerollte zeitliche Abhängigkeit der Zufallsvariablen

    Markot-Kette (ausgerollte Darstellung)

    Markot-Kette (ausgerollte Darstellung)

  • Rekursive Darstellung der zeitliche Abbildung der Zufallsvariablen

    Markot-Kette (rekursive Darstellung)

    Markot-Kette (rekursive Darstellung)

  • Betont Übergange und Wahrscheinlichkeit

    Markot-Kette (betont Übergange und Wahrscheinlichkeit)

    Markot-Kette (betont Übergange und Wahrscheinlichkeit)