Глава 4. Устойчивость, единственность, неравенство Гронуолла

Где стоит глава

От существования к чувствительности

Глава 8.3 дала гарантию: при липшицевости правой части оператор Пикара сжимает, а значит, процесс приближений сходится и решение единственно. Но у задачи Коши есть ещё один вопрос, не менее важный для приложений: насколько решение чувствительно к начальному условию? Если слегка сдвинуть $y_0$, как далеко разойдутся траектории? Ответ — устойчивость, или непрерывная зависимость от начальных данных, и инструмент для него — неравенство Гронуолла.

Эти три темы — единственность, устойчивость, Гронуолл — суть одна геометрия, прочитанная под тремя углами. Единственность — крайний случай: нулевое возмущение, траектории совпадают. Устойчивость — общий случай: возмущение ${\Delta }_0=y_0-{\tilde{y}}_0$ порождает расхождение, и это расхождение контролируемо — оно не превосходит геометрического множителя, умноженного на начальное возмущение. Гронуолл — общий принцип, стоящий за этим множителем: если величина растёт со скоростью, ограниченной кратным себя, то она ограничена геометрической (в пределе — экспоненциальной) оценкой.

Что эта глава утверждает — и что нет

Доказанное ядро — конкретные процессные оценки над ℚ. Расхождение двух Эйлеровых траекторий с разными начальными условиями на каждом шаге растёт не более чем в $(1+Lh)$ раз, а за $n$ шагов — не более чем в $(1+Lh{)}^n$; однородная итерация Гронуолла даёт ту же геометрическую форму $u_0(1+\alpha h{)}^n$, а при множителе не больше единицы — затухание.¹ Эти оценки — точные рациональные неравенства над ℚ, и их 0-аксиомный статус проверен Print Assumptions; переход к $\epsilon$-сходимости множителя к завершённому $e^{Lt}$ относится к L3-слою (как в Главе 8.3).

Чего глава не делает: она не предъявляет классический множитель $e^{Lt}$ как завершённую экспоненту, не утверждает безусловную оценку на бесконечном интервале и не доказывает глобальное продолжение решения. Множитель $(1+Lh{)}^n$ — рациональный процесс; его предел $e^{Lt}$ — роль-предел, а $\epsilon$-сходимость к нему наследует L3 (как и в Главе 8.3).

Три уровня. Первый — рациональные оценки расхождения и роста: доказаны над ℚ без аксиом. Второй — процессные чтения: расхождение как чувствительность, множитель как роль-роста, устойчивость как контролируемость. Третий — $e^{Lt}$ как завершённая экспонента и оценка на бесконечном интервале — роль-пределы.

E/R/R-каркас главы

Каркас в порождающем порядке; полный разбор — в § «Разбор в порождающем порядке». Опираемся на E/R/R-шапку ProcessODEStability.v. Глава вводит систему контроля расхождения решений.

Rules (Правила, $\leftrightarrow$ L5). Липшицевость даёт пошаговую оценку $|{\Delta }_{n+1}|\le (1+Lh)|{\Delta }_n|$, откуда геометрическая $|{\Delta }_n|\le (1+Lh{)}^n|{\Delta }_0|$; неравенство Гронуолла обобщает: рост, ограниченный кратным себя, влечёт геометрическую оценку $u_0(1+\alpha h{)}^n$; множитель $(1+Lh{)}^n$ есть стадия процесса для $e^{Lt}$.

Roles (Роли, $\leftrightarrow$ L4). Расхождение $|{\Delta }_n|$ — роль-чувствительность; множитель $(1+Lh{)}^n$ — роль-роста; константа Липшица $L$ — роль-скорость роста ошибки; возмущение ${\Delta }_0$ — роль-причина; устойчивость — роль-свойство (контролируемость расхождения); $e^{Lt}$ — роль-предел.

Elements (Элементы, $\leftrightarrow$ L1$+$P4). Рациональные расхождения $|{\Delta }_n|$, множители $(1+Lh{)}^n$, итерации Гронуолла $u_n$, конечный номер $n$. Под P4 завершённый $e^{Lt}$ и оценка на бесконечном интервале — вне элементов.

Единственность как нулевое возмущение

Единственность решения — уже доказанный факт (Глава 8.3): у сжимающего оператора Пикара не более одной неподвижной точки на сетке, и эта единственность установлена конструктивно, через конечный максимум, без закона исключённого третьего.² В контексте устойчивости у неё появляется наглядное прочтение: единственность есть устойчивость при нулевом возмущении. Если начальные условия совпадают, ${\Delta }_0=0$, то оценка расхождения $|{\Delta }_n|\le (1+Lh{)}^n\cdot 0=0$ заставляет траектории совпадать на каждом шаге. Одно и то же начало по одному и тому же закону порождает одну и ту же траекторию.

Единственность — роль-жёсткость, прочитанная как предельный случай устойчивости: нулевое возмущение даёт нулевое расхождение. Решение определено своим началом однозначно — не потому, что <<так постулировано>>, а потому, что процесс сжатия не оставляет места второй траектории.

Расхождение траекторий: непрерывная зависимость

Возьмём две Эйлеровы траектории $y$ и $\tilde{y}$ одного уравнения, стартующие из разных точек, и обозначим их расхождение ${\Delta }_n=y_n-{\tilde{y}}_n$. Один шаг Эйлера меняет расхождение предсказуемо: новые значения отличаются на старое расхождение плюс разность наклонов, а липшицевость ограничивает эту разность,

$$|{\Delta }_{n+1}|\le |{\Delta }_n|+hL|{\Delta }_n|=(1+Lh)|{\Delta }_n|.$$

Каждый шаг растягивает расхождение не более чем в $(1+Lh)$ раз.³ Накапливая за $n$ шагов, получаем непрерывную зависимость от начальных данных:

$$|{\Delta }_n|\le (1+Lh{)}^n|{\Delta }_0|.$$

Расхождение в любой момент ограничено начальным возмущением, умноженным на геометрический множитель.⁴ На конечном интервале множитель конечен, и потому малое возмущение даёт малое расхождение — это и есть содержательная устойчивость над ℚ — дискретная процессная форма непрерывной зависимости (малое начальное расхождение даёт малое расхождение на любом конечном числе шагов).

Расхождение $|{\Delta }_n|$ — роль-чувствительность: мера того, как решение реагирует на сдвиг начала. Множитель $(1+Lh{)}^n$ — роль-роста этой чувствительности. Устойчивость — не отсутствие расхождения, а его контролируемость: конечный множитель на конечном интервале.

Неравенство Гронуолла

Геометрический множитель $(1+Lh{)}^n$ — частный случай общего принципа, который в теории ОДУ называют неравенством Гронуолла. Принцип таков: если величина растёт со скоростью, ограниченной кратным самой себя (плюс, возможно, внешний приток), то она ограничена геометрической оценкой. В дискретной процессной форме рекуррентный шаг Гронуолла — это

$$u_{n+1}=u_n+h(\alpha u_n+\beta )=(1+\alpha h)u_n+\beta h,$$

где $\alpha$ — скорость самовоздействия, $\beta$ — внешний приток.⁵ В однородном случае ($\beta =0$) итерация решается явно и даёт в точности геометрический рост:

$$u_n=u_0(1+\alpha h{)}^n.$$

Это и есть процессная форма Гронуолла, и при $\alpha =L$, $u_0=|{\Delta }_0|$ она совпадает с оценкой расхождения предыдущего раздела — устойчивость есть Гронуолл для расхождения.⁶ А если множитель не превосходит единицы ($r=1+\alpha h$ в диапазоне $0\le r\le 1$ — затухающий режим), итерация не растёт, а убывает: $u_n\le u_0$.⁷

Гронуолл — роль-правило контроля: оно превращает локальную информацию о скорости роста (<<не быстрее, чем кратное себя>>) в глобальную оценку всей траектории. Геометрический множитель $(1+\alpha h{)}^n$ — роль-роста; устойчивость и затухание — два его режима ($\alpha$ положительно или отрицательно.)

Множитель роста и $e^{Lt}${e^(Lt)} как роль-предел

Классическая теория устойчивости пишет оценку непрерывной зависимости в экспоненциальной форме:

$$|y(t)-\tilde{y}(t)|\le |y_0-{\tilde{y}}_0|e^{Lt}.$$

Откуда здесь экспонента? Она — роль-предел нашего геометрического множителя. Если взять шаг $h=t/n$ (то есть достичь момента $t$ за $n$ шагов), множитель станет $(1+Lt/n{)}^n$ — а это в точности стадия процесса для $e^{Lt}$, того самого $e$-процесса, что разбирался в Главах 8.2 и 8.5.⁸ Так классическая форма $e^{Lt}$ оказывается не предзаданной экспонентой, а роль-пределом рационального процесса $(1+Lt/n{)}^n$ — лестницы множителей, каждый из которых рационален и вычислим.

Доказано конкретное и процессное: пошаговая и геометрическая оценки расхождения, их совпадение с однородным Гронуоллом, тождество множителя со стадией $e^{Lt}$-процесса. $\epsilon$-сходимость множителя к завершённому $e^{Lt}$ и безусловная оценка на бесконечном интервале — роль-предел и граница, наследующие L3 (как в Главе 8.3).

Классический множитель $e^{Lt}$ — роль-предел процесса роста, а не предзаданная экспонента. Оценка устойчивости $|\Delta (t)|\le |{\Delta }_0|e^{Lt}$ читается как направление, к которому стремится рациональная оценка $|{\Delta }_n|\le (1+Lh{)}^n|{\Delta }_0|$; актуальна конечная оценка, $e^{Lt}$ — её роль-предел.

E/R/R-разбор: расхождение как Роль-чувствительность

Постановка

Разберём построенную систему — контроль расхождения решений — в терминах E/R/R, читая E/R/R-шапку ProcessODEStability.v.⁹ Это онтологическое осмысление и чтение шапки. Аббревиатура задаёт эпистемический порядок (Elements $\to$ Roles $\to$ Rules); разбор ведём в порождающем.

Разбор в порождающем порядке

Rules (Правила, $\leftrightarrow$ L5). В основании — пошаговая оценка расхождения $|{\Delta }_{n+1}|\le (1+Lh)|{\Delta }_n|$ из липшицевости. Над ней — геометрическая оценка $|{\Delta }_n|\le (1+Lh{)}^n|{\Delta }_0|$ и её обобщение — неравенство Гронуолла ($u_n=u_0(1+\alpha h{)}^n$ в однородном случае, затухание при множителе $\le 1$). Венчает — тождество множителя $(1+Lt/n{)}^n$ со стадией $e^{Lt}$-процесса. Универсальный слой включает L3 лишь там, где множитель читается как $\epsilon$-сходимость к $e^{Lt}$; сами рациональные оценки его не требуют.

Roles (Роли, $\leftrightarrow$ L4). Расхождение $|{\Delta }_n|$ — роль-чувствительность (как решение реагирует на сдвиг начала). Множитель $(1+Lh{)}^n$ — роль-роста. Константа Липшица $L$ — роль-скорость роста ошибки. Возмущение ${\Delta }_0$ — роль-причина. Устойчивость — роль-свойство (контролируемость расхождения конечным множителем). $e^{Lt}$ — роль-предел процесса роста.

Elements (Элементы, $\leftrightarrow$ L1$+$P4). Носители конечны: рациональные расхождения $|{\Delta }_n|\in Q$, множители $(1+Lh{)}^n\in Q$, итерации Гронуолла $u_n\in Q$, константы $L,h\in Q$, номер $n\in N$. Под P4 завершённый $e^{Lt}$ и оценка на бесконечном интервале не суть элементы: актуальны конечные оценки на конечных интервалах.

Компонент	Что фиксирует	E/R/R
$|{\Delta }_{n+1}|\le (1+Lh)|{\Delta }_n|$; Гронуолл $u_0(1+\alpha h{)}^n$; множитель$=$стадия $e^{Lt}$	КАК контролируется рост расхождения	Rules ($L5$)
расхождение $|{\Delta }_n|$; множитель $(1+Lh{)}^n$; $L$; возмущение ${\Delta }_0$; устойчивость; $e^{Lt}$	ЗАЧЕМ значимы носители: роли в контроле	Roles ($L4$)
$|{\Delta }_n|,(1+Lh{)}^n,u_n\in Q$; $L,h\in Q$; $n\in N$	ЧТО есть на каждой стадии (конечно, P4)	Elements ($L1$, P4)

Проверка сформированности и что даёт разбор

{ Система сформирована корректно: каждый компонент попадает ровно в одну E/R/R-категорию, и нет самоотнесения (P1). Существенно, что устойчивость взята как контролируемость расхождения конечным множителем, а не как утверждение о завершённом поведении на бесконечном интервале. Множитель $(1+Lh{)}^n$ — Элемент (рациональное число на каждой стадии), а его предел $e^{Lt}$ — роль-предел; смешать их значило бы поставить завершённую экспоненту на место рационального множителя.}

Что даёт разбор. Он переводит результаты главы на язык структуры: расхождение — Роль-чувствительность, множитель роста — Роль-роста, устойчивость — Роль-свойство контролируемости, а сами числа $|{\Delta }_n|$, $(1+Lh{)}^n$, $u_n$ — Элементы. Единственность вписывается как предельный случай (нулевое возмущение), а Гронуолл — как общее Правило, частным случаем которого служит оценка расхождения. Брать $e^{Lt}$ как завершённую экспоненту значит ставить роль-предел на место Элемента — то же смешение, что <<число как объект>> (Часть II) или <<решение как готовая функция>> (Глава 8.1).

Заключительная глава Части соберёт нити воедино: экспонента как процесс (решение $y^{\prime }=y$), линейные уравнения и затухание, системы и гармонический осциллятор, применения — и подведёт итог онтологии <<решение есть процесс>>, обозначив переход к уравнениям в частных производных.

---

---

Часть: Часть VIII. Обыкновенные дифференциальные уравнения · Том: «Математика»

Навигация: ← Глава 3. Пикар–Линделёф: оператор-сжатие на пространстве функций · Глава 5. Экспонента, линейные уравнения, применения →

Footnotes

1. Опорные файлы — ProcessODEStability.v (непрерывная зависимость: euler_step_diff, euler_diff_bound; 0 аксиом, Print Assumptions <>) и ProcessGronwall.v (gronwall_step, gronwall_iterate, gronwall_beta_zero: $u_0(1+\alpha h{)}^n$; gronwall_decay_factor: затухание при множителе $\le 1$). Единственность — picard_unique (ProcessPicardOperator.v, Глава 8.3). ↩

2. picard_unique в ProcessPicardOperator.v: если две грид-функции неподвижны под оператором Пикара, они совпадают во всех узлах. 0 аксиом (Глава 8.3). ↩

3. euler_step_diff в ProcessODEStability.v: $|{\Delta }_{k+1}|\le (1+Lh)|{\Delta }_k|$, через неравенство треугольника и Липшица. 0 аксиом. ↩

4. euler_diff_bound в ProcessODEStability.v: $|{\Delta }_n|\le (1+Lh{)}^n|{\Delta }_0|$, индукцией по $n$ через euler_step_diff. 0 аксиом. ↩

5. gronwall_step и gronwall_iterate в ProcessGronwall.v; gronwall_process — соответствующий RealProcess. ↩

6. gronwall_beta_zero в ProcessGronwall.v: gronwall_iterate alpha 0 h u0 n == u0 * Qpow (1+alpha*h) n. 0 аксиом. При $\alpha =L$ это правая часть euler_diff_bound. ↩

7. gronwall_decay_factor в ProcessGronwall.v: для $0\le r\le 1$ выполнено $u_0{r}^n\le u_0$. 0 аксиом. ↩

8. euler_growth_is_expt_process в ProcessODEStability.v: множитель $(1+s/(k+1){)}^{k+1}$ есть стадия $‘exp{t}_{p}rocess‘sk$ из ProcessExpProcess.v (батч C). 0 аксиом. При $s=Lt$ это $(1+Lt/n{)}^n$. ↩

9. Шапка размечает прямо: Rules — Липшиц $\Rightarrow$ $|{\Delta }_{n+1}|\le (1+Lh)|{\Delta }_n|$ $\Rightarrow$ $|{\Delta }_n|\le (1+Lh{)}^n|{\Delta }_0|$; Roles — $|{\Delta }_n|$ роль-расхождение, $(1+Lh{)}^n$ роль-роста, $e^{Lt}$ роль-предел; Elements — рациональные $|{\Delta }_n|$, $(1+Lh{)}^n$, конечный $n$. ↩