Majorizace ODR: porovnání řešení

Majorizace systémů ODR: Porovnání řešení podle pravých stran

Při analýze dynamických systémů je přesné řešení obyčejných diferenciálních rovnic často nedosažitelné. Věta o majorizaci umožňuje odhadnout chování "složitého" systému pomocí řešení "jednoduchého" se známými vlastnostmi. Pokud pravá strana a počáteční podmínky jednoho systému jsou po složkách větší než druhého, pak i trajektorie řešení zachovávají tuto nerovnost.

Pro vektory $x=(x^1,\ldots,x^m)$, $y=(y^1,\ldots,y^m) \in \mathbb{R}^m$ označme $x \le y$, pokud $x^i \le y^i$ pro všechna $i=1,\dots,m$. Kvádr $[x,y]:=\{z\in\mathbb{R}^m\mid x\le z\le y\}$.

Uvažujme oblast $D \subset \mathbb{R}^m$ takovou, že pro $x,y\in D$, $x\le y$ platí $[x,y]\subset D$. Funkce $f_1,f_2\in C(I\times D,\mathbb{R}^m)$ se spojitými jakobiány $\frac{\partial f_i}{\partial x}\in C(I\times D,\mathbb{R}^{m\times m})$, kde $I=[t_0,t_1]$.

Dva Cauchyho úlohy:

$\dot x = f_i(t,x)$, $x(t_0)=\hat x_i$, $i=1,2$.

Formulace hlavní věty

Nechť $f_1(t,x)\le f_2(t,x)$ pro všechna $(t,x)\in I\times D$ a $\hat x_1\le\hat x_2$. Existují diagonální matice $U_i(t,x)$ z nezáporných prvků a matice $B_i(t,x)$ s nezápornými prvky:

$\frac{\partial f_i}{\partial x}=U_i + B_i$, $i=1,2$.

Tvrzení: Pokud $x_i(t)$ jsou řešení Cauchyho úloh na $I$, pak $x_1(t)\le x_2(t)$ pro všechna $t\in I$.

Tato vlastnost umožňuje majorizovat neřešitelný systém jednoduchým, jehož vlastnosti jsou známé.

Struktura důkazu

Zavedeme homotopii $y=y(t,s)$ pro $s\in[0,1]$:

$y_t=(1-s)f_1(t,y)+s f_2(t,y)$, $y(t_0,s)=(1-s)\hat x_1 + s\hat x_2$.

Pro $s=0$ dostaneme $x_1(t)$, pro $s=1$ — $x_2(t)$. Stačí ukázat, že $s\mapsto y(t,s)$ je neklesající pro pevné $t$, tj. $y_s\ge 0$.

Derivování podle parametru

Derivujeme podle $s$:

$$\begin{aligned} y_{st} &= f_2 - f_1 + (u + b)y_s, \\ u(t,s) &= sU_2(t, y(t,s)) + (1-s)U_1(t, y(t,s)), \\ b(t,s) &= sB_2(t, y(t,s)) + (1-s)B_1(t, y(t,s)), \\ y_s|_{t=t_0} &= \hat x_2 - \hat x_1 \ge 0.\end{aligned}$$

Oblast monotonie

$M=A=\{t\in I\mid y_s(\xi,s)\ge0,\ \xi\in[t_0,t],\ s\in[0,1]\}$. Množina je neprázdná ($t_0\in A$), uzavřená. Nechť $\tau=\sup A$.

Předpokládejme $\tau<t_1$. V integrálním tvaru:

$$y_s(t,s)=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)+\int_\tau^t e^{\int_\xi^t u\,d\lambda}(f_2-f_1+by_s)\,d\xi=:\Phi(y_s).$$

Metoda komprese

Podle principu kompresních zobrazení se rovnice $\Phi(y_s)=y_s$ řeší pro malé $t-\tau>0$. Nulová aproximace $Y_0=e^{\int_\tau^t u\,d\lambda}y_s(\tau,s)\ge0$. Indukčně $Y_{k+1}=\Phi(Y_k)\ge0$, protože operátor $\Phi$ zachovává nezápornost. Limitní $y_s\ge0$ na $[\tau,\tau+\varepsilon]\times[0,1]$, spor s definicí $\tau$.

Aplikace: Gronwallova-Bellmanova nerovnost

Pro nezápornou $w(t)\le C+\int_{t_0}^t v(\xi)w(\xi)\,d\xi$, $v\ge0$.

Pomocná $\phi(t)=C+\int_{t_0}^t v w\,d\xi$, $w\le\phi$, $\dot\phi=vw\le v\phi$.

Systém 1: $\dot x=f_1(t,x):=vx-(v\phi-\dot\phi)$, $x(t_0)=C$. Závorky jsou nezáporné.

Majorizující: $\dot x=f_2(t,x):=vx$, $x(t_0)=C$, řešení $x(t)=Ce^{\int_{t_0}^tv\,d\xi}$.

Podle věty $w\le\phi\le Ce^{\int_{t_0}^tv\,d\xi}$.

Co je důležité:

• Majorizace funguje při $f_1\le f_2$ a $\hat x_1\le\hat x_2$.
• Jakobiány se rozkládají na diagonální $U$ (nezápornou) a $B$ (nezápornou).
• Důkaz pomocí homotopie a metody komprese.
• Přímá souvislost s Gronwallovou nerovností pro integrální odhady.
• Aplikace: odhad řešení nelineárních ODR jednoduchými exponenciálními funkcemi.

Závěr

Věta poskytuje nástroj pro kvalitativní analýzu dynamiky bez explicitního integrování. Podmínky na jakobiány zajišťují monotonii homotopie, což je klíčové pro zachování nerovnosti.

— Editorial Team