Dyskretyzacja (matematyka)

W matematyce dyskretyzacja dotyczy procesu transformowania modeli i równań funkcji ciągłych na ich dyskretne odpowiedniki. Jest to zwykle pierwszy krok w procesie przygotowywania tych modeli (i równań) do ewaluacji numerycznej i implementacji na komputerach cyfrowych. Do przetwarzania na komputerze cyfrowym ponadto potrzebne jest wykonanie kwantyzacji.

Szczególnie istotne są tu dwie metody dyskretyzacji:

• dyskretyzacja Eulera (zob. metoda Eulera)
• dyskretyzacja ekstrapolatorem rzędu zerowego (ang. ZOH, Zero-order hold).

Rozwiązanie zdyskretyzowanego cząstkowego równania różniczkowego, uzyskane za pomocą metody elementów skończonych.

Dyskretyzacja związana jest także z matematyką dyskretną i jest ważną częścią (komputerowych) obliczeń ziarnistych (ang. granular computing) stosowanych w mechanice komputerowej. W tym kontekście dyskretyzacja odnosi się także do modyfikacji zmiennej w kategorii ziarnistości gdy agreguje się wiele zmiennych dyskretnych albo dokonuje się fuzji wielu kategorii dyskretnych.

Dyskretne równoważniki transmitancji operatorowej

W teorii sterowania, metodę projektowania układów dyskretnych polegająca na zaprojektowaniu kompensatora czasu ciągłego, a następnie zastąpieniu go równoważnikiem dyskretnym tak by można go zaimplementować w urządzeniu cyfowym nazywa się emulacją. Metoda ta jest bardzo szeroko używana przez inżynierów praktyków. Dyskretne równoważniki transmitancji operatorowej to transmitancje dyskretne, które aproksymują te same charakterystyki (w pewnym zakresie częstotliwości) jak dana transmitancja czasu ciągłego $G(s)\,$. Można w tym celu zastosować poniższe metody realizujące to zadanie:

• całkowanie numeryczne - w metodzie tej przeprowadza się całkowanie numeryczne równań różniczkowych opisujących wykonany projekt czasu ciągłego. Istnieje wiele technik pozwalających na całkowanie numeryczne w tym Metoda Eulera i techniki oparte na regułach prostokąta i trapezu.
• dyskretyzacja odpowiedzi impulsowej - w metodzie tej wyznacza się dla transmitancji ciągłej $G(s)\,$ ciągłej odpowiedz impulsową, którą następnie dyskretyzuje się. Ostatecznie dla dyskretnej odpowiedzi impulsowej wyznacza się transmitancję dyskretną $G(z)=Z[G(s)]\,$.
• przekształcenie zerowo-biegunowe - w metodzie tej porównuje się dziedzinę "s" oraz dziedzinę "z". Odpowiedź układu ciągłego z biegunem w pewnym punkcie $s=s_{0}\,$ w układzie spróbkowanym z okresem próbkowania $T\,$ reprezentowana jest przez odpowiedź układu dyskretnego z biegunem w punkcie $z=e^{s_{0}T}\,$. Ta własność może być wykorzystana do przekształcenia zer i biegunów, które aproksymują układ dyskretny.
• równoważność ekstrapolacji - metoda ta polega na pobieraniu próbek sygnału wejściowego, następnie ekstrapolacji pomiędzy próbkami do postaci aproksymacji sygnału i przesyłaniu tych aproksymacji przez transmitancję układu.

Całkowanie numeryczne

Całkowanie numeryczne jest zadaniem dość złożonym. Najbardziej elementarne techniki z tego zakresu to reguły o małej złożoności i ustalonym rozmiarze kroku. W metodzie tej daną transmitancję układu ciągłego G(s) zastępuje się przez równanie różniczkowe a następnie wyprowadza się równania różnicowe będące aproksymacją równań różniczkowych.

Niech dana będzie transmitancja integratora analogowego:

$G(s)={\frac {U(s)}{E(s)}}={\frac {1}{s}}\,$

gdzie $E(s)$ oraz $U(s)$ są odpowiednio transformatami wejścia i wyjścia integratora. Dla integratora tego można określić równoważne równanie różniczkowe

${\frac {du(t)dt}{e(t)}}=e(t)\,$

które można zapisać w postaci całkowej:

$u(t)=\int \limits _{0}^{t}e(\tau )d\tau \,$

Wiele reguł opiera się na właściwej sobie metodzie aproksymacji składnika powiększania pola (pod krzywą funkcją, która w powyższym wzorze podlega całkowaniu). Należą do nich:

• reguła prostokątna wprzód
• reguła prostokatna wstecz
• reguła trapezu.

W regule prostokatnej wprzód obszar aproksymuje się przez prostokąt wyznaczany wprzód od chwili kT do chwili kT+T i bierze jako amplitudę prostokąta wartość napotkaną w kT. Szerokość takiego prostokąta wynosi T. Można więc zapisać równanie w pierwszej aproksymacji:

$u_{1}(kT+T)=u_{1}(kT)+Te(kT)\,$

gdzie wyrażenie $u_{1}(kT)\,$ reprezentuje obszar pod całkowaną krzywą e(t) w przedziale od t = 0 do t = kT. Po zastosowaniu transformaty Z do powyższej zależności otrzymuje się:

$G_{F}(z)={\frac {U_{1}(z)}{E(z)}}={\frac {T}{z-1}}={\frac {1}{\frac {z-1}{T}}}\,$.

W regule prostokatnej wstecz obszar aproksymuje się przez prostokąt wyznaczany wstecz od chwili kT do kT-T i bierze jako amplitudę prostokąta wartość napotkaną w kT. Szerokość takiego prostokąta wynosi T. Można więc zapisać równanie w pierwszej aproksymacji:

$u_{2}(kT)=u_{2}(kT-T)+Te(kT)\,$

Po zastosowaniu transformaty Z do powyższej zależności otrzymuje się:

$H_{B}(z)={\frac {U_{2}(z)}{E(z)}}={\frac {zT}{z-1}}={\frac {1}{{\frac {1}{T}}{\frac {z-1}{z}}}}\,$

W regule trapezu obszar aproksymuje się przez pole trapezu umieszczonego pod całkowaną krzywą. Równanie aproksymacji ma wówcas postać :

$u_{3}(kT+T)=u_{3}(kT)+{\frac {T}{2}}[e(kT)+e(kT+1)]\,$

Po zastosowaniu transformaty Z do powyższej zależności otrzymuje się:

$H_{T}(z)={\frac {U_{3}(z)}{E(z)}}={\frac {T}{2}}{\frac {z+1}{z-1}}={\frac {1}{{\frac {2}{T}}{\frac {z-1}{z+1}}}}\,$

Metoda reguły trapezu jest również znana jako metoda Tustina lub pod nazwą transformacji biliniowej (bilinearnej). Metoda projektowania wykorzystująca tę regułę polega na tym, że daną transmitancję ciągłą, $G(s)\,$, równoważną transmitancja dyskretnej wyznacza się przez podstawienie:

$G_{T}(z)=G(s)|_{s={\frac {2}{T}}{\frac {z-1}{z+1}}}\,$

Porównując transmitancje operatorowe z trzema aproksymacjami dyskretnymi można zauważyć, że transmitancję dyskretną można uzyskać bezpośrednio z transformaty operatorowej podstawiając za zmienną zespoloną "s" jej aproksymatę.

Każda z powyższych aproksymacji może być potraktowana jako przekształcenie płaszczyzny s na płaszczyznę z.

Dyskretyzacja modelu układu liniowego w przestrzeni stanów

Dyskretyzacja stosowana jest też przy transformacji ciągłych równań różniczkowych do dyskretnych równań różnicowych, odpowiednich dla analizy numerycznej.

Następujący model zmiennych stanu czasu ciągłego

${\dot {\mathbf {x} }}(t)=\mathbf {A} \mathbf {x} (t)+\mathbf {B} \mathbf {u} (t)+\mathbf {w} (t)$

$\mathbf {y} (t)=\mathbf {C} \mathbf {x} (t)+\mathbf {D} \mathbf {u} (t)+\mathbf {v} (t)$

gdzie $v\,$ i $w\,$ to żródła ciągłego szumu białego o zerowej średniej z kowariancjami

$\mathbf {w} (t)\sim N(0,\mathbf {Q} )$

$\mathbf {v} (t)\sim N(0,\mathbf {R} )$

można zdyskretyzować, przyjmując ekstrapolator rzędu zerowego dla wejścia $u\,$ i ciągłe całkowanie dla szumu $v\,$, do postaci:

$\mathbf {x} [k+1]=\mathbf {A} _{d}\mathbf {x} [k]+\mathbf {B} _{d}\mathbf {u} [k]+\mathbf {w} [k]$

$\mathbf {y} [k]=\mathbf {C} _{d}\mathbf {x} [k]+\mathbf {D} _{d}\mathbf {u} [k]+\mathbf {v} [k]$

z kowariancjami

$\mathbf {w} [k]\sim N(0,\mathbf {Q} _{d})$

$\mathbf {v} [k]\sim N(0,\mathbf {R} _{d})$

gdzie:

$\mathbf {A} _{d}=e^{\mathbf {A} T}={\mathcal {L}}^{-1}\{(s\mathbf {I} -\mathbf {A} )^{-1}\}_{t=T}$

$\mathbf {B} _{d}=\left(\int _{\tau =0}^{T}e^{\mathbf {A} \tau }d\tau \right)\mathbf {B} =\mathbf {A} ^{-1}(\mathbf {A} _{d}-I)\mathbf {B}$, jeśli $\mathbf {A}$ jest nieosobliwa

$\mathbf {C} _{d}=\mathbf {C}$

$\mathbf {D} _{d}=\mathbf {D}$

$\mathbf {Q} _{d}=\int _{\tau =0}^{T}e^{\mathbf {A} \tau }\mathbf {Q} e^{\mathbf {A} ^{T}\tau }d\tau$

$\mathbf {R} _{d}=\mathbf {R}$

a $T\,$ jest czasem próbkowania.

Zręczne wyliczenie $Ad\,$ i $Bd\,$ w jednym kroku można wykonać korzystając z następującej własności:

$\mathbf {e} ^{\mathbf {\begin{bmatrix}\mathbf {A} &\mathbf {B} \\\mathbf {0} &\mathbf {0} \end{bmatrix}} T}={\begin{bmatrix}\mathbf {M_{11}} &\mathbf {M_{12}} \\\mathbf {0} &\mathbf {I} \end{bmatrix}}$

i wówczas mając:

$\mathbf {A} _{d}=M_{11}$

$\mathbf {B} _{d}=M_{12}$

Dyskretyzacja szumu procesu

Numeryczna ewaluacja $\mathbf {Q} _{d}$ jest nieco bardziej złożona z uwagi na całkę eksponenty macierzy. Można ją, jednakże, wyliczyć poprzez skonstruowanie najpierw macierzy a następnie wyliczenie na komputerze jej eksponenty:

$\mathbf {F} ={\begin{bmatrix}-\mathbf {A} &\mathbf {Q} \\\mathbf {0} &\mathbf {A} ^{T}\end{bmatrix}}T$

$\mathbf {G} =e^{\mathbf {F} }={\begin{bmatrix}\dots &\mathbf {A} _{d}^{-1}\mathbf {Q} _{d}\\\mathbf {0} &\mathbf {A} _{d}^{T}\end{bmatrix}}.$

Zdyskretyzowany szum procesu jest wówczas wyliczany poprzez przemnożenie transponowanej dolnej, prawej partycji macierzy G z górną, prawą partycją macierzy G:

$\mathbf {Q} _{d}=(\mathbf {A} _{d}^{T})^{T}(\mathbf {A} _{d}^{-1}\mathbf {Q} _{d}).$

Wyprowadzenie

Rozpoczynając z modelem ciągłym

$\mathbf {\dot {x}} (t)=\mathbf {A} \mathbf {x} (t)+\mathbf {B} \mathbf {u} (t)$

wiadomo, że eksponenta macierzy jest następująca:

${\frac {d}{dt}}e^{\mathbf {A} t}=\mathbf {A} e^{\mathbf {A} t}=e^{\mathbf {A} t}\mathbf {A}$

i przez wcześniejsze przemnożenie modelu uzyskuje się:

$e^{-\mathbf {A} t}\mathbf {\dot {x}} (t)=e^{-\mathbf {A} t}\mathbf {A} \mathbf {x} (t)+e^{-\mathbf {A} t}\mathbf {B} \mathbf {u} (t)$

co zapisać można jako

${\frac {d}{dt}}(e^{-\mathbf {A} t}\mathbf {x} (t))=e^{-\mathbf {A} t}\mathbf {B} \mathbf {u} (t)$

a następnie całkując:

$e^{-\mathbf {A} t}\mathbf {x} (t)-e^{0}\mathbf {x} (0)=\int _{0}^{t}e^{-\mathbf {A} \tau }\mathbf {B} \mathbf {u} (\tau )d\tau$

$\mathbf {x} (t)=e^{\mathbf {A} t}\mathbf {x} (0)+\int _{0}^{t}e^{\mathbf {A} (t-\tau )}\mathbf {B} \mathbf {u} (\tau )d\tau$

co jest rozwiązaniem analitycznym dla modelu ciągłego.

Teraz należy zdyskretyzować powyższe wyrażenie. Można przyjąć, że $u\,$ jest stała podczas każdego kroku czasowego.

$\mathbf {x} [k]\ {\stackrel {\mathrm {def} }{=}}\ \mathbf {x} (kT)$

$\mathbf {x} [k]=e^{\mathbf {A} kT}\mathbf {x} (0)+\int _{0}^{kT}e^{\mathbf {A} (kT-\tau )}\mathbf {B} \mathbf {u} (\tau )d\tau$

$\mathbf {x} [k+1]=e^{\mathbf {A} (k+1)T}\mathbf {x} (0)+\int _{0}^{(k+1)T}e^{\mathbf {A} ((k+1)T-\tau )}\mathbf {B} \mathbf {u} (\tau )d\tau$

$\mathbf {x} [k+1]=e^{\mathbf {A} T}\left[e^{\mathbf {A} kT}\mathbf {x} (0)+\int _{0}^{kT}e^{\mathbf {A} (kT-\tau )}\mathbf {B} \mathbf {u} (\tau )d\tau \right]+\int _{kT}^{(k+1)T}e^{\mathbf {A} (kT+T-\tau )}\mathbf {B} \mathbf {u} (\tau )d\tau$

Wyrażenie w nawiasie można zapisać jako $\mathbf {x} [k]$ a drugie wyrażenie można uprościć przez podstawienie $v=kT+T-\tau$. Ponadto można przyjąć, że $\mathbf {u}$ jest stałe podczas całkowania, co z koleii daje:

$\mathbf {x} [k+1]=e^{\mathbf {A} T}\mathbf {x} [k]+\left(\int _{0}^{T}e^{\mathbf {A} v}dv\right)\mathbf {B} \mathbf {u} [k]=e^{\mathbf {A} T}\mathbf {x} [k]+A^{-1}\left(e^{\mathbf {A} T}-I\right)\mathbf {B} \mathbf {u} [k]$

co stanowi dokładne rozwiązanie dyskretyzowanego problemu.

Aproksymacje

Dokładna dyskretyzacja czasami może być trudna z uwagi na dużą eksponentę macierzy i związane z tym operacje całkowania. Znacznie łatwiej wyliczyć, w oparciu o nią, przybliżony model dyskretny dla małych kroków czasowych $e^{\mathbf {A} T}\approx \mathbf {I} +\mathbf {A} T$. Przybliżone rozwiązanie przyjmuje wówczas postać:

$\mathbf {x} [k+1]\approx (\mathbf {I} +\mathbf {A} T)\mathbf {x} [k]+(\mathbf {I} T+{\frac {1}{2}}\mathbf {A} T^{2})\mathbf {B} \mathbf {u} [k]$

co można dalej aproksymować jeśli ${\frac {1}{2}}\mathbf {A} T^{2}$ jest małe; co daje:

$\mathbf {x} [k+1]\approx (\mathbf {I} +\mathbf {A} T)\mathbf {x} [k]+T\mathbf {B} \mathbf {u} [k]$

Inne możliwe aproksymacje to: $e^{\mathbf {A} T}\approx \left(\mathbf {I} -\mathbf {A} T\right)^{-1}$ i $e^{\mathbf {A} T}\approx \left(\mathbf {I} +{\frac {1}{2}}\mathbf {A} T\right)\left(\mathbf {I} -{\frac {1}{2}}\mathbf {A} T\right)^{-1}$. Każda z nich ma inne własności związane ze stabilnością. Ostatnia znana jest jako transformacja Tustina (transformacja bilinearna) i zachowuje stabilność lub odpowiednio niestabilność układu czasu ciągłego.

Dyskretyzacja własności ciągłych

W statystyce i w uczeniu maszynowym termin dyskretyzacja odnosi się do procesu konwersji ciągłych własności lub zmiennych na zdyskretyzowane lub nominalne własności. Może to być użyteczne przy tworzeniu masowych funkcji prawdopodobieństwa.

Zobacz też

• przestrzeń dyskretna
• przestrzeń czasowa