Нелинейные дифференциальные уравнения, дискретность пространства-времени и эпсилон произведение
Всё содержание этой статьи является следствием решения задачи уровня первого года институтского курса математического анализа.
Здесь мы вводим новое произведение для решетчатых функций (мы будем называть его эпсилон-умножением), что дает нам возможность увидеть интересную связь между интегрально-дифференциальными уравнениями (в том числе и нелинейными) и рекуррентными соотношениями. Это позволяет с необычного ракурса взглянуть на некоторые методы их решения.
Но, что мне показалось особенно интересным, это то, что новое произведение также позволяет нам «порассуждать» о таких фундаментальных вещах, как непрерывность пространства-времени. Это рассуждение, конечно, нужно рассматривать, скорее, как интеллектуальное упражнение.
Мне показалось, что этот математический ребус или, если хотите, небольшое математическое путешествие на уровне знаний первого-второго курса технического ВУЗа может заинтересовать читателей Хабра, интересующихся математикой.
Замечание
- вы не увидите никаких ссылок и, вообще, всё это результат размышлений отдельно взятого человека
- возможно (и даже скорее всего), всё это и не ново, и, возможно, некоторые термины, которые я использую не соответствуют математической традиции, но так как это не научная публикация, здесь нет претензии на новизну и всё это довольно просто, то это терминологическое отличие, на мой взгляд, не является чем-то действительно важным
Краткое изложение
В этой статье речь будет идти о двух типах функций (определенных в области действительных чисел): аналитических и решетчатых.
Замечание 1Точнее речь будет идти об аналитических справа от нуля функциях, которые мы будем обозначать заглавными буквами и о решетчатых функциях, определенных на множестве целых неотрицательных чисел (), для обозначения которых мы будем использовать строчные буквы . Множество таких решетчатых функций мы будем обозначать как . Точные определения этих функций будут приведены ниже в данной статье
Замечание 2Иногда, в примерах мы будем расширять наш подход и на область комплексных чисел. Строго говоря, это должно быть обосновано. Но я опускаю это обоснование (как и доказательства некоторых других утверждений), чтобы не перегружать и так длинную статью
Мы установим взаимно однозначное соответствие между функциями из этих двух классов. Мы также введем новые операции для решетчатых функций. Эти операции будут аналогичны «обычным» операциям над «обычными» функциями. Мы будем отличать эти новые операции добавлением приставки «эпсилон». Так «обычной» производной будет соответствовать эпсилон-производная, обычному произведению — эпсилон-произведение, обычному интегрированию — эпсилон-интегрирование и т.д. При этом, если операция остается неизменной (как в случае суммирования и вычитание) название меняться не будет.
Одна из основных целей введения этих операций проста — сохранить свойства производной произведения функций:
Что в случае решетчатых функций () будет выглядеть как
где — левый оператор эпсилон-производной
— эпсилон-умножение
Это позволит нам (по аналогии, например, с преобразованием Лапласа) ввести новое преобразование, которое мы будем называть эпсилон-преобразованием. Это даст нам возможность преобразовывать в общем случае нелинейные интегрально-дифференциальные уравнения в рекуррентные выражения и решать их (численно или аналитически), соответственно, рекуррентными методами.
ЗамечаниеТакже справедливо и обратное: мы получаем возможность решать рекуррентные соотношения методами интегрально-дифференциальных уравнений
Важным фактом при этом является то, что при стремлении шага решетчатой функции к нулю (), все наши эпсилон-функции, эпсилон-операции и рекуррентные уравнения стремятся к «обычным» функциям (аналитическим справа от нуля), обычным операциям между ними и обычным интегрально-дифференциальным уравнениям.
Но тогда, это наводит на фундаментальный вопрос:, а не является ли эпсилон-умножение «истинным» умножением в нашем физическом мире (по крайней мере там, где речь идет о пространстве или о времени)?
Так, например, говоря о времени, при достаточно малом шаге (например, порядка Планковского времени сек) эпсилон-уравнение Шредингера c для временных процессов с характерным временем (или в области частот с частотами ) будет давать тот же результат, что и обычное уравнение Шредингера, но при этом, естественным образом вводится квантование времени.
В принципе, аналогичное рассуждение возможно и для пространства.
— множество всех функций определенных выше
— иногда, для обозначения функции для простоты будем писать просто
— заглавными буквами мы будем обозначать «обычные» функции. Все подобные функции в данной статье будут аналитическими справа от нуля функциями (определение будет дано ниже в статье)
— шаг решетчатой функции: . В зависимости от контекста может быть или обозначением самой решетчатой функции или значением этой функции в точке
— левый оператор, обозначающий производную дифференцируемой функции:
— левый оператор, обозначающий эпсилон-производную решетчатой функции:
— левый оператор, введенный для удобства записи некоторых выражений, и обозначающий производную дифференцируемой функции, применяемый только к функции (или к производной функции любого порядка) в произведении функций, например,
— единичный левый оператор для решетчатой функции:
— эпсилон-производная решетчатой функции:
— левый оператор, введенный для удобства записи некоторых выражений, и обозначающий эпсилон-производную решетчатой функции, применяемый только к решетчатой функции (или эпсилон-производной функции любого порядка) в произведении (обычном) решетчатых функций, например,
— эпсилон произведение
— -ая эпсилон-степень функции . Так, например,
— -ая эпсилон-производная функции . Так, например,
— эпсилон-преобразование (прямое)
Если , то и являются эпсилон-сопряженными, при этом является эпсилон-образом и — эпсилон-прообразом
— обратное эпсилон-преобразование
- Мы будем рассматривать решетчатые функции , …, определенные на множестве , где , может быть любым положительным действительным числом. Поэтому иногда в примерах, в тех случаях, когда мы не используем стремление , мы принимаем для простоты
- Назовем эпсилон дифференцированием следующий левый оператор:
, где — шаг решетчатой функции - Введем новое произведение, которое будем называть эпсилон-произведением и обозначать как .
Здесь три эквивалентные формулы для эпсилон-произведения:
где — биномиальный коэффициент Ньютона: , ()
Его основные свойства: - Вводим эпсилон-ряд Тейлора. Он будет аналогичен обычному ряду Тейлора, но произведение, возведение в степень и производные должны быть заменены на эпсилон-произведение, эпсилон-степень и эпсилон-производные:
- Введение эпсилон-ряда Тейлора позволяет нам ввести эпсилон-функции, аналогичные аналитическим функциям (например, эпсилон-тригонометрические функции). Для этого потребуем, чтобы эпсилон-производные в нуле (всех порядков) для этих эпсилон-функций равнялись производным в нуле (справа) для обычных функций. Для обозначения эпсилон-функции мы будем добавлять приставку . Так, например, для эпсилон-экспоненты получим:
- Введем понятие эпсилон преобразования. Суть этого преобразования заключается в замене всех операций (производная, умножение, возведение в степень) на эпсилон-операции, а всех функций на эпсилон-функции. При данном преобразовании тождества будут сохраняться. Так, например, для эпсилон-тригонометрических функций все тождества будут аналогичны обычным тригонометрическим тождествам. Например,
Если эпсилон преобразование обозначить через , то этот факт можно записать, как
Эпсилон-функцию мы также будем называть эпсилон-образом. Так эпсилон-косинус является эпсилон-образом косинуса. При этом сам косинус является эпсилон-прообразом для эпсилон-косинуса. При этом эти две функции (например, эпсилон-косинус и косинус) являются эпсилон-сопряженными.
- Соответственно, обратное эпсилон-преобразование — это преобразование обратное (прямому) эпсилон-преобразованию. Будем обозначать его, как .
- Таким образом, через прямые и обратные эпсилон-преобразования мы устанавливаем соответствие между дифференциальными уравнениями (ДУ) (в том числе и нелинейными) и рекуррентными соотношениями (РС).
Теперь для решения нелинейных дифференциальных уравнений мы можем применить следующую схему:
ДУ РС (решение РС) (решение ДУ)
То есть, мы делаем эпсилон-преобразование исходного ДУ, получаем РC, решаем его и делаем обратное эпсилон-преобразование, получая таким образом решение дифференциального уравнения. Важно понимать, что таким образом будут найдены лишь решения, являющиеся функциями, аналитическими справа от нуля. Это тот же подход, что применяется при решении, например, линейных дифференциальных уравнений с помощью преобразований Лапласа.
- Аналогично, рекуррентные уравнения можно привести к дифференциальным и решать их соответствующими методами. В этом случае мы имеем следующую схему:
РC ДУ (решение ДУ) (решение РC)
- Также легко можно ввести и дополнительные операции, и преобразования, например, эпсилон-интегрирование и эпсилон-преобразование Лапласа
- В статье приводятся три примера
- решение линейного дифференциального уравнения гармонических колебаний
- решение нелинейного дифференциального уравнения (Бесселя)
- решение линейного рекуррентного уравнения (формула для чисел Фибоначчи)
- При стремлении (шаг решетчатой функции) эпсилон-функции и эпсилон-операторы стремятся к своим прообразам. Так, например, при стремлении эпсилон произведение стремится к обычному произведению, эпсилон-производная — к обычной производной, … эпсилон-косинус — к обычному косинусу, эпсилон-синус к обычному синусу… Это значит, что обратное эпсилон-преобразование может быть выполнено обычным стремлением . Этот факт позволяет нам рассматривать дифференциально-интегральные уравнения (в том числе и нелинейные) с решениями аналитическими справа от нуля, как предельный случай рекуррентных уравнений с шагом .
- Это наводит на мысль о том, а не является ли эпсилон-умножение «истинным» умножением (во всяком случае в некоторых случаях) в нашем физическом мире? Так, например, рассмотрим одномерное уравнение Шредингера:
Это уравнение предполагает непрерывность пространства времени и возможность бесконечно малой величины и (что заложено в том факте, что мы используем дифференцирование), что в случае квантовой механики подразумевает бесконечную энергию, что невозможно.
Но давайте рассмотрим эпсилон-преобразование этого уравнения (заменив производные на эпсилон-производные, функции на решетчатые эпсилон-функции, произведение — на эпсилон-произведение, …). При достаточно малом мы получим решения очень близкие к тем, которые мы бы получили при решении изначального уравнения Шредингера. При этом естественным образом вводится дискретность пространства-времени, что устраняет проблему бесконечной энергии.
- Это в свою очередь ставит вопрос о том, как можно проверить, а какое же умножение является «истинным» в нашем физическом мире, обычное или эпсилон-умножение? Возможно ли поставить мыслительный или физический эксперимент для ответа на этот вопрос?
Далее каждый из этих пунктов рассматривается более детально
Решетчатые функции и эпсилон-производная
Решетчатая функция и множество
Решетчатые функции — это функции, определенные на дискретном множестве действительных чисел с постоянным шагом. Мы будем рассматривать только функции, определенные на множестве , где , , где — целое неотрицательное число. При этом . Будем обозначать множество таких функций через . Все решетчатые функции, которые мы будем рассматривать ниже будут принадлежать этому множеству .
Пример решетчатой функции
Функция , где :
Эпсилон-производная
Будем называть эпсилон-производной функции в следующую величину:
Соответственно, функцию, значение которой в для любого равно эпсилон-производной функции в , мы будем называть эпсилон-производной функции и обозначать как (или просто ).
Очевидно, что функция и поэтому по отношению к ней можно еще раз применить эпсилон дифференцирование, и т.д. Мы будем обозначать символом -ю -производную функции. По индукции легко доказывается, что
Эпсилон-интегрирование
Процедуру, обратную эпсилон-дифференцированию, назовем соответственно эпсилон-интегрированием. Функцию , значение которой в для любого
будем называть эпсилон-первообразной функции .
Эпсилон производная произведения
К сожалению, эпсилон-производная не обладает свойствами обычной производной. Так, например, очевидно, что в общем случае
Часть отличает эту «производную» от обычной производной произведения дифференцируемых функций. В случае решения численными методами нелинейных дифференциальных уравнений это дает дополнительную погрешность в вычислениях.
ЗамечаниеЭто дополнительное слагаемое приводит также, в частности, к тому, что мы не можем применить обычную формулу Тейлора для разложения решетчатой функции в ряд
Теперь давайте представим, что у нас не было бы этого дополнительного члена. Чтобы нам это могло дать?
По крайней мере, мы могли бы применять некоторые методы решения дифференциальных уравнений для решения рекуррентных уравнений (уравнений с функциями из множества ). Например, мы имели бы разложение в ряд Тейлора, мы смогли бы пользоваться интегрированием (эпсилон-интегрированием) по частям, применять преобразования (эпсилон-преобразование) Лапласа. Но в действительности все даже интереснее, и этому и посвящена данная статья.
Но как избавиться от этого дополнительного элемента? Мы можем создать, например, другое умножение (эпсилон-умножение).
Базовая задача
Это именно та задача, следствием решения которой и является вся эта статья.
Мы хотим найти такую формулу для эпсилон-произведения , чтобы удовлетворялись следующие требования:
- выполнялись бы свойства обычного умножения, такие, как коммутативность, ассоциативность, дистрибутивность, существование единицы (единственной) и нулевого элемента (единственного)
- для эпсилон-производной эпсилон-произведения решетчатых функций мы должны иметь ту же формулу, что и для производной произведения (обычных) дифференцируемых функций, то есть
- при стремлении (шаг решетчатой функции) эпсилон-произведение должно стремиться к обычному произведение , где и соответственно функции, к которым стремятся решетчатые функции и при , то есть:
Эпсилон-умножение
Всем этим свойствам удовлетворяет функция , определяемая как (все три формулы эквивалентны):
,
где — биномиальный коэффициент Ньютона: , ()
ЗамечаниеЧтобы не перегружать статью мы не будем приводить доказательства, но первые два свойства доказываются легко, последнее — не так очевидно
Чтобы немного освоиться с новым произведением, рассмотрим несколько примеров.
Таким образом получаем, что эпсилон-аналог (далее будет введено понятие эпсилон-образа) функции является функция
Так, например,
Очевидно, что равна нулю для
ЗамечаниеВ последней формуле пропущено , что не является проблемой, т.к. может быть любым (), в том числе и 1. Конечно, мы не должны забывать про , если исследуется зависимость от него, например, при
Эпсилон-ряд Тейлора
Теперь мы имеем полный аналог разложения Тейлора:
Будем называть этот ряд эпсилон-рядом Tейлора, а само разложение — эпсилон-разложением в ряд Tэйлора.
Эпсилон преобразование
Определение. (Эпсилон преобразование)
Пусть ряд Tэйлора функции в правой окрестности нуля сходится к самой функции. Тогда функцию мы будем называть аналитической функцией в нуле справа. Тогда функции и будем называть эпсилон-сопряженными, если для любого целого или нулевого -ая производная справа функции в нуле равняется -ой эпсилон-производной функции в нуле.
Функцию тогда будем называть эпсилон-образ или в случае, если эта функция имеет название (например, cos), то будем называть просто эпсилон-'название функции' (например, эпсилон-косинус) и обозначать как ''обозначение функции'' (cos). Функцию будем называть эпсилон-прообраз функции .