Глава 5. Экспонента, линейные уравнения, применения

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

Замыкание Части VIII

Часть VIII прошла путь от онтологии к гарантиям: Глава 8.1 переопределила решение задачи Коши как процесс приближений, Глава 8.2 предъявила метод Эйлера на сетке, Глава 8.3 доказала, что оператор Пикара есть сжатие на пространстве грид-функций (сжимающий механизм и геометрическое сближение итератов), а Глава 8.4 разобрала устойчивость и Гронуолл. Настоящая, заключительная глава делает три вещи: показывает процесс на канонических уравнениях (экспонента, линейные, системы), доводит дело до применений (динамика цены) и синтезирует всю Часть, обозначая переход к уравнениям в частных производных (Часть IX).

Все примеры подчиняются единой онтологии, выстроенной в предыдущих главах: решение — процесс рациональных приближений; правая часть — закон скорости; завершённые функции $e^t$, $\sin$, $\cos$ — роль-пределы, а не предзаданные трансцендентности. Глава не вводит нового аппарата; она показывает, что построенный язык работает на конкретных уравнениях и применениях.

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

Доказанное ядро — конкретные процессы над ℚ: экспоненциальная траектория $y^{\prime }=y$ как лестница $(1+h{)}^n{y}_0$, затухание $y^{\prime }=-\alpha y$, общее линейное $y^{\prime }=Ly$, векторный шаг для систем и гармонический осциллятор $y»=-y$, дискретная динамика цены.¹ Все стадии — конечные рациональные числа, проверяемые точной арифметикой.

Чего глава не делает: она не предъявляет $e^t$, $\sin$, $\cos$ как завершённые трансцендентности, не доказывает точное сохранение энергии осциллятора и не строит глобальных во времени решений. Эти объекты — роль-пределы, и $\epsilon$-сходимость к ним наследует L3 (как в Главах 8.3–8.4).

Экспонента как процесс

Канонический пример — уравнение $y^{\prime }=y$, $y(0)=1$, чьё решение $y(t)=e^t$. Метод Эйлера для него особенно прозрачен: поскольку наклон равен значению, каждый шаг умножает на $(1+h)$, и за $n$ шагов получается $(1+h{)}^n$.² Замечательно, что этот процесс совпадает с однородной итерацией Гронуолла из Главы 8.4: экспоненциальный рост и Гронуолл-оценка — буквально одна и та же рекурсия.³

Достигая момента $t=1$ за $N$ шагов размера $1/N$, получаем лестницу $(1+1/N{)}^N$ — $e$-процесс, разобранный в Главе 8.2: рациональные стадии $2,9/4,64/27,\dots$, ограниченные снизу двойкой для всех $N$ и начинающие классический процесс приближения к $e$ (завершённое $e$ — роль-предел).⁴ Для произвольного времени аналогично: $e^t$ — роль-предел процесса $(1+t/N{)}^N$, чьи стадии суть рациональные числа.⁵

$e$ и $e^t$ — не завершённые трансцендентности, а роль-пределы рациональных процессов: лестниц $(1+1/N{)}^N$ и $(1+t/N{)}^N$, порождённых уравнением $y^{\prime }=y$. Спросить <<чему равно $e^t$>> значит спросить <<какой процесс его задаёт>>; ответ — Эйлерова траектория уравнения роста. Актуальна стадия, $e^t$ — направление.

Линейные уравнения: рост и затухание

Уравнение $y^{\prime }=Ly$ — общее линейное однородное — ведёт себя как масштабированная экспонента: шаг Эйлера умножает на $(1+Lh)$, и знак $L$ решает судьбу процесса. При $L>0$ — рост (множитель больше единицы), при $L<0$ — затухание.

Затухание — $y^{\prime }=-\alpha y$ при $\alpha >0$ — особенно важно для приложений (остывание, распад, релаксация).⁶ Здесь множитель шага есть $1-\alpha h$, и при $0\le 1-\alpha h\le 1$ процесс не растёт, а убывает — ровно случай затухающего Гронуолла из Главы 8.4: $u_n\le u_0$. Решение $y(t)=y_0{e}^{-\alpha t}$ снова есть роль-предел рационального процесса $y_0(1-\alpha h{)}^n$, и завершённая убывающая экспонента не строится как объект.

Линейное уравнение задаёт роль-темп: константа $L$ (или $-\alpha$) — скорость, с которой процесс умножает себя на шаге. Рост и затухание — два режима одного правила; в обоих решение есть процесс множителей, а завершённая экспонента — роль-предел.

Системы и гармонический осциллятор

До сих пор речь шла об одном уравнении. Но многие задачи — это системы: пара (или $n$-ка) величин, эволюционирующих совместно. Векторный шаг Эйлера переносит всё сказанное на пары: состояние $(y,v)$ обновляется покомпонентно, $(y,v)\mapsto (y+h{f}_1,v+h{f}_2)$.⁷

Векторная форма позволяет свести уравнение высшего порядка к системе первого: $y»=g(t,y,y^{\prime })$ становится парой $y^{\prime }=v$, $v^{\prime }=g(t,y,v)$ — введение скорости как самостоятельной роли расширяет пространство состояний и понижает порядок.⁸ Классический пример — гармонический осциллятор $y»=-y$, то есть пара $y^{\prime }=v$, $v^{\prime }=-y$. Его точное решение — $(\cos t,\sin t)$, и эта траектория — роль-предел векторного Эйлерова процесса; $\sin$ и $\cos$ как завершённые трансцендентности — не Элементы, а роль-пределы, ровно как $e$.

Осциллятор обнажает и честную сторону численного метода: явный Эйлер не сохраняет энергию. Величина $E=y^2+v^2$, классически постоянная, на сетке дрейфует вверх: за два шага с $y_0=1,v_0=0,h=1/10$ она вырастает с $1$ до $1,0201$.⁹ Это не дефект изложения, а свойство явного метода, честно вычисленное: точное сохранение энергии принадлежит непрерывному потоку осциллятора (или специально сохраняющим структуру, симплектическим, схемам), тогда как явный Эйлер на конечной сетке даёт систематический рациональный дрейф.

Система — совместный векторный процесс; скорость в редукции порядка — самостоятельная роль, а не вспомогательный символ. Энергия осциллятора — роль-инвариант непрерывного потока (и сохраняющих структуру схем); её дрейф на явной сетке — честный рациональный Элемент. $\cos$ и $\sin$ — роль-пределы векторного процесса, как $e$ — роль-предел скалярного.

Применения: динамика цены

Язык ОДУ-процессов выходит за пределы чистой математики. В прикладных областях уравнение изменения часто задаёт дискретную динамику с самого начала — и тогда процессная форма не приближение, а собственно модель. Простой пример — динамика цены: цена обновляется на шаге внешним воздействием (приток спроса, новость, поток заявок),

$$p_{n+1}=p_n+f_n,$$

а вся траектория — процесс по списку воздействий.¹⁰ Конкретно: цена $100$ под воздействиями $1,-\frac{1}{2},\frac{3}{4},-\frac{1}{4},\frac{1}{2}$ проходит стадии $100,101,\frac{201}{2},\frac{405}{4},101$ и приходит в $\frac{203}{2}=101,5$ — точная рациональная траектория.¹¹ Файл вводит и проверку устойчивости через липшицев коэффициент — аналогичное липшицево условие (коэффициент сжатия меньше единицы), родственное условию $LT<1$ сжатия Пикара (Глава 8.3), прочитанное как признак устойчивой динамики.

В прикладной динамике процесс — не приближение к завершённой кривой, а сама модель: цена есть роль-процесс, обновляемый воздействиями. Устойчивость управляется липшицевым условием того же типа, что и сходимость Пикара — единый E/R/R-механизм работает и в чистой, и в прикладной задаче.

Синтез Части VIII и переход к уравнениям в частных производных

Единая E/R/R-структура Части

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

Rules (Правила, $\leftrightarrow$ L5). В основании — само уравнение $y^{\prime }=f(t,y)$ (закон скорости) и правило порождения приближений (шаг Эйлера или оператор Пикара). Над ними — гарантия: при липшицевости итерация сжимает (Глава 8.3), откуда Коши-сходимость процесса и контроль расхождения по Гронуоллу (Глава 8.4). Универсальный слой — L1–L5, причём L3 входит лишь там, где роль-предел читается как завершённый объект.

Roles (Роли, $\leftrightarrow$ L4). Правая часть $f$ — роль-закон-скорости; приближение $y_n$ — роль-стадия; неподвижная точка оператора — роль-решение; константа Липшица — роль-темп; а завершённые $e^t$, $\sin$, $\cos$ и решение на бесконечном интервале — роль-пределы.

Elements (Элементы, $\leftrightarrow$ L1$+$P4). На каждой стадии — конечные рациональные данные: значения $y_n$, узлы $t_n$, множители $(1+Lh{)}^n$, номер $n$. Под P4 завершённые функции-решения и трансцендентности не суть элементы: актуальны конечные стадии процессов.

Компонент	Что фиксирует	E/R/R
уравнение $y^{\prime }=f$; шаг Эйлера / оператор Пикара; сжатие $\Rightarrow$ Коши, Гронуолл	КАК порождается и контролируется процесс решения	Rules ($L5$)
закон скорости $f$; стадия $y_n$; решение; темп $L$; $e^t,\sin ,\cos$ как предел	ЗАЧЕМ значимы носители: роли в процессе	Roles ($L4$)
$y_n,t_n,(1+Lh{)}^n\in Q$; номер $n\in N$	ЧТО есть на каждой стадии (конечно, P4)	Elements ($L1$, P4)

Что доказано и где граница

{ Подведём баланс Части. Доказано над ℚ без аксиом: метод Эйлера и его конкретные траектории; оператор Пикара как сжатие на грид-функциях ($Q^{N+1}$, sup-метрика) с геометрическим зазором итератов; единственность решения на сетке через конечный максимум; пошаговый и геометрический контроль расхождения; однородный Гронуолл и затухание; экспонента, линейные уравнения, осциллятор и динамика цены как конкретные рациональные процессы. Граница (L3 / роль-предел): $\epsilon$-сходимость итераций и множителей к завершённому пределу; завершённые функции-решения, $e^t$, $\sin$, $\cos$ как объекты; глобальное во времени продолжение, blow-up, максимальный интервал. Эта граница — не пробел, а честная демаркация: то, что над ℚ есть процесс, не выдаётся за завершённый объект.}

Переход к уравнениям в частных производных

Часть VIII решала обыкновенные уравнения — эволюцию во времени. Следующий шаг — уравнения в частных производных (Часть IX): эволюция функции, зависящей и от времени, и от пространства. Процессная онтология переносится без изменения принципа: пространство дискретизуется сеткой узлов, время — шагами, а решение — процесс на пространственно-временной сетке. Уравнение теплопроводности станет процессом сглаживания, волновое — процессом распространения, а интеграл и производные — конечными разностями над ℚ. То, что в Части VIII было лестницей чисел $(1+1/N{)}^N$, в Части IX станет лестницей сеточных функций; но онтология — решение есть процесс, а не завершённый объект — остаётся той же.

Решение дифференциального уравнения — обыкновенного или в частных производных — есть процесс приближений, чьи стадии конечны и рациональны, а завершённая функция-решение — роль-предел. Это и есть итог Части VIII: уравнения изменения решаются не извлечением готового объекта из функционального пространства, а разворачиванием Коши-процесса, гарантированного сжатием и контролируемого Гронуоллом.

---

---

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

Навигация: ← Глава 4. Устойчивость, единственность, неравенство Гронуолла · Глава 1. Решение как процесс на пространственно-временной сетке — Часть IX →

Footnotes

1. Опорные файлы — ProcessODEExamples.v (экспонента exp_rhs, затухание decay_rhs, линейное), ProcessExpProcess.v ($e$-процесс, батч C), ProcessODESystem.v (системы и осциллятор, батч A), PriceDynamicsODE.v (динамика цены, каталог stdlib/trading/). Новые файлы — 0 аксиом (Print Assumptions). ↩

2. exp_rhs$=\lambda y.y$ и exp_euler в ProcessODEExamples.v: exp_euler_succ даёт $y_{n+1}=(1+h)y_n$. exp_rhs_lipschitz — константа Липшица $1$. ↩

3. exp_euler_is_gronwall в ProcessODEExamples.v: exp_euler y0 h n == gronwall_iterate 1 0 h y0 n. Экспонента есть Гронуолл с $\alpha =1$, $\beta =0$. ↩

4. exp1_process, exp1_is_euler (тождество <<$e$-процесс$=$траектория Эйлера для $y^{\prime }=y$ при $t=1$>>), exp1_lower ($\ge 2$ для всех $N$, Бернулли) в ProcessExpProcess.v, 0 аксиом. ↩

5. expt_process в ProcessExpProcess.v: стадия $N$ есть $(1+t/(N+1){)}^{N+1}$; конкретно $‘exp{t}_{p}rocess‘00=1$ (значение $e^0$), $‘exp{t}_{p}rocess‘10=2$ (первая стадия $e^1$). ↩

6. decay_rhs$=\lambda y.-\alpha y$ и decay_rhs_lipschitz (константа $\alpha$) в ProcessODEExamples.v. ↩

7. sys_euler в ProcessODESystem.v (батч A): векторный Эйлер для системы $y^{\prime }=f_1$, $v^{\prime }=f_2$; законы шага sys_euler_fst_S, sys_euler_snd_S. 0 аксиом. ↩

8. order2_euler в ProcessODESystem.v; order2_position_step: в редуцированной системе позиция растёт на $h\cdot v$, то есть уравнение $y^{\prime }=v$ выполнено. 0 аксиом. ↩

9. sho_euler (осциллятор $y^{\prime \prime }=-y$) в ProcessODESystem.v: sho_energy2 ($E_2=10201/10000$) и sho_energy_drifts ($1<E_2$). Точное сохранение энергии — роль-предел; дрейф на сетке — точный рациональный факт. ↩

10. price_step$=\lambda pf.p+f$ и price_trajectory в PriceDynamicsODE.v (stdlib/trading/). Это Эйлеров шаг с внешним притоком — частный случай правой части $f(t,p)$. ↩

11. traj_final в PriceDynamicsODE.v: price_trajectory p0_ex forces_ex == 203#2; промежуточные шаги — traj_step1–traj_step4. Точная арифметика над ℚ. ↩