примеры решения диффуров (ДУ) в математике
Часто одно лишь упоминание дифференциальных уравнений вызывает у студентов неприятное чувство. Почему так происходит? Чаще всего потому, что при изучении основ материала возникает пробел в знаниях, из-за которого дальнейшее изучение диффуров становиться просто пыткой. Ничего не понятно, что делать, как решать, с чего начать?
Однако мы постараемся вам показать, что диффуры – это не так сложно, как кажется.
Основные понятия теории дифференциальных уравнений
Со школы нам известны простейшие уравнения, в которых нужно найти неизвестную x. По сути дифференциальные уравнения лишь чуточку отличаются от них – вместо переменной х в них нужно найти функцию y(х), которая обратит уравнение в тождество.
Дифференциальные уравнения имеют огромное прикладное значение. Это не абстрактная математика, которая не имеет отношения к окружающему нас миру. С помощью дифференциальных уравнений описываются многие реальные природные процессы. Например, колебания струны, движение гармонического осциллятора, посредством дифференциальных уравнений в задачах механики находят скорость и ускорение тела. Также
Дифференциальное уравнение (ДУ) – это уравнение, содержащее производные функции y(х), саму функцию, независимые переменные и иные параметры в различных комбинациях.
Существует множество видов дифференциальных уравнений: обыкновенные дифференциальные уравнения, линейные и нелинейные, однородные и неоднородные, дифференциальные уравнения первого и высших порядков, дифуры в частных производных и так далее.
Решением дифференциального уравнения является функция, которая обращает его в тождество.
Существуют общие и частные решения ДУ.Общим решением ДУ является общее множество решений, обращающих уравнение в тождество. Частным решением дифференциального уравнения называется решение, удовлетворяющее дополнительным условиям, заданным изначально.
Порядок дифференциального уравнения определяется наивысшим порядком производных, входящих в него.
Решение уравнений
Обыкновенные дифференциальные уравнения
Обыкновенные дифференциальные уравнения – это уравнения, содержащие одну независимую переменную.
Рассмотрим простейшее обыкновенное дифференциальное уравнение первого порядка. Оно имеет вид:
Решить такое уравнение можно, просто проинтегрировав его правую часть.
Примеры таких уравнений:
Уравнения с разделяющимися переменными
В общем виде этот тип уравнений выглядит так:
Приведем пример:
Решая такое уравнение, нужно разделить переменные, приведя его к виду:
После этого останется проинтегрировать обе части и получить решение.
Математика
Линейные дифференциальные уравнения первого порядка
Такие уравнения имеют вид:
Здесь p(x) и q(x) – некоторые функции независимой переменной, а y=y(x) – искомая функция. Приведем пример такого уравнения:
Решая такое уравнение, чаще всего используют метод вариации произвольной постоянной либо представляют искомую функцию в виде произведения двух других функций y(x)=u(x)v(x).
Для решения таких уравнений необходима определенная подготовка и взять их “с наскока” будет довольно сложно.
Пример решения ДУ с разделяющимися переменными
Вот мы и рассмотрели простейшие типы ДУ. Теперь разберем решение одного из них. Пусть это будет уравнение с разделяющимися переменными.
Сначала перепишем производную в более привычном виде:
Затем разделим переменные, то есть в одной части уравнения соберем все «игреки», а в другой – «иксы»:
Теперь осталось проинтегрировать обе части:
Интегрируем и получаем общее решение данного уравнения:
Конечно, решение дифференциальных уравнений – своего рода искусство. Нужно уметь понимать, к какому типу относится уравнение, а также научиться видеть, какие преобразования нужно с ним совершить, чтобы привести к тому или иному виду, не говоря уже просто об умении дифференцировать и интегрировать. И чтобы преуспеть в решении ДУ, нужна практика (как и во всем). А если у Вас в данный момент нет времени разбираться с тем, как решаются дифференциальные уравнения или задача Коши встала как кость в горле или вы не знаете, как правильно оформить презентацию, обратитесь к нашим авторам. В сжатые сроки мы предоставим Вам готовое и подробное решение, разобраться в подробностях которого Вы сможете в любое удобное для Вас время. А пока предлагаем посмотреть видео на тему «Как решать дифференциальные уравнения»:
Однородные дифференциальные уравнения 1 порядка
Готовые ответы к примерам на однородные дифференциальные уравнения первого порядка ищут многие студенты (ДУ 1 порядка самые распространенные в обучении), далее Вы их сможете подробно разобрать. Но прежде чем перейти к рассмотрению примеров рекомендуем внимательно прочитать краткий теоретический материал.
Уравнения вида P(x,y)dx+Q(x,y)dy=0, где функции P(x,y) і Q(x,y) являются однородными функциями одного порядка называют однородным дифференциальным уравнением (ОДР).
1. Сначала нужно применить подстановку y=z*x, где z=z(x) – новая неизвестная функция (таким образом исходное уравнение сводится к дифференциальному уравнению с разделяющимися переменными.
2. Производная произведения равна y’=(z*x)’=z’*x+z*x’=z’*x+z или в дифференциалах dy=d(zx)=z*dx+x*dz.
3. Далее подставляем новую функцию у и ее производную y’ (или dy) в ДУ с разделяющимися переменными относительно x та z.
4. Решив дифференциальное уравнение с разделяющимися переменными, сделаем обратную замену y=z*x, поэтому z= y/х, и получим общее решение (общий интеграл) дифференциального уравнения.
5. Если задано начальное условие y(x0)=y0, то находим частное решение задачи Коши. В теории все звучит легко, однако на практике не у всех так весело получается решать дифференциальные уравнения. Поэтому для углубления знаний рассмотрим распространенные примеры. На легких задачах нет особо Вас научить, поэтому сразу перейдем к более сложным.
Пример 1. Решить дифференциальное уравнение
Решение: Делим правую сторону уравнения на переменную, которая стоит множителем возле производной. В результате придем к однородного дифференциального уравнения 0 порядка
И здесь многим пожалуй стало интересно, как определить порядок функции однородного уравнения?
Вопрос достаточно уместен, а ответ на него следующий:
в правую сторону подставляем вместо функции и аргумента значение t*x, t*y. При упрощении получают параметр «t» в определенном степени k, его и называют порядком уравнения. В нашем случае «t» сократится, что равносильно 0-м степени или нулевом порядке однородного уравнения.
Далее в правой стороне можем перейти к новой переменной y=zx; z=y/x .
При этом не забываем выразить производную «y» через производную новой переменной. По правилу части находим
Уравнения в дифференциалах примет вид
Совместные слагаемые в правой и левой части сокращаем и переходим к дифференциальному уравнению с разделенными переменными.
Проинтегрируем обе части ДУ
Для удобства дальнейших преобразований постоянную сразу вносим под логарифм
По свойствам логарифмов полученное логарифмическое уравнение эквивалентно следующему
Эта запись еще не решение (ответ), необходимо вернуться к выполненной замене переменных
Таким образом находят общее решение дифференциальных уравнений. Если Вы внимательно читали предыдущие уроки, то мы говорили, что схему вычисления уравнений с разделенными переменными Вы должны уметь применять свободно и такого рода уравнения придется вычислять для более сложных типов ДУ.
Пример 2. Найти интеграл дифференциального уравнения
Решение:Схема вычислений однородных и сводных к ним ДУ Вам тепер знакома. Переносим переменную в правую сторону уравнения, а также в числителе и знаменателе выносим x2, как общий множитель
Таким образом получим однородное ДУ нулевого порядка.
Следующим шагом вводим замену переменных z=y/x, y=z*x, о которой постоянно будем напоминать, чтобы Вы ее заучили
После этого ДУ записываем в дифференциалах
Далее преобразуем зависимость к дифференциальному уравнению с отделенными переменными
и интегрированием решаем его.
Интегралы несложные, остальные преобразования выполнены на основе свойств логарифма. Последнее действие включает экспонирования логарифма. Наконец возвращаемся к исходной замене и записываем решение дифференциального уравнения в форме
Константа «C» принимает любое значение. Все кто учится заочно имеют проблемы на экзаменах с данным типом уравнений, поэтому просьба внимательно посмотреть и запомнить схему вычислений.
Пример 3. Решить дифференциальное уравнение
Решение:Как следует из приведенной выше методики, дифференциальные уравнения такого типа решают методом введения новой переменной. Перепишем зависимость так, чтобы производная была без переменной
Далее по анализу правой части видим, что везде присутствует частка -ее и обозначаем за новую неизвестную
z=y/x, y=z*x.
Находим производную от y
С учетом замены первоначальное ДУ перепишем в виде
Одинаковые слагаемые упрощаем, а все получившие сводим к ДУ с отделенными переменными
Интегрированием обеих частей равенства
приходим к решению в виде логарифмов
Экспонируя зависимости находим общее решение дифференциального уравнения
которое после подстановки в него начальной замены переменных примет вид
Здесь С — постоянная, которую можно доопределить из условия Коши. Если не задана задача Коши то стала принимает произвольное действительное значение.
Вот и вся мудрость в исчислении однородных дифференциальных уравнений.
Автономные дифференциальные уравнения на прямой
Итак, давайте научимся решать какие-нибудь дифференциальные уравнения. Для начала — очень простые.В этой главе мы будем рассматривать дифференциальные уравнения вида
˙x=f(t,x),(2.1)
где x:D→R — неизвестная функция, D — связное подмножество прямой (вся прямая, луч, отрезок, полуинтервал, интервал), f:D×R→R — некоторая по меньшей мере непрерывная (а лучше бы гладкая, как мы увидим чуть позже) функция от двух переменных.Напомним, что решением уравнения (2.1) называется дифференцируемая функция φ, такая, что выполнено тождество
˙φ(t)=f(t,φ(t))∀t∈D
Обсудим для начала, как можно было бы находить значение функции φ, не пытаясь выписать решение в виде явной формулы. Оказывается, с помощью компьютера это довольно несложно сделать — правда, решение будет не точным, а приближённым. Обсуждение этого метода окажется полезным и для наших дальнейших аналитических построений.2.1Численное решение дифференциальных уравнений. Метод Эйлера
Пусть поставлена задача Коши: ˙x=f(t,x),x(t0)=x0.(2.2)Мы можем приблизительно решать её таким образом. Возьмём произвольную точку (t0,x0) расширенного фазового пространства. Интегральная кривая, проходящая через эту точку, имеет в ней касательную с угловым коэффициентом f(t0,x0). Касательная — это прямая, которая хорошо приближает график функции. Давайте на секундочку представим, что интегральная кривая в точности совпадает с касательной на некотором небольшом промежутке времени — начиная с момента t0 и заканчивая t0+Δt, где Δt — некоторое маленькое число. Иными словами, мы считаем, что на этом промежутке двигаемся с постоянной скоростью — той, которая была в момент времени t0, то есть f(t0,x0). В этом случае к моменту времени t0+Δt мы пройдём расстояние, равное f(t0,x0)⋅Δt и попадём в точку (t1,x1), задаваемую следующим образом: x1=x0+f(t0,x0)⋅Δtt1=t0+Δt Точка (t1,x1) лежит на касательной, проходящей через точку (t0,x0).Этот метод приближённого нахождения решений называется методом Эйлера. Он даёт представление о том, как можно использовать компьютер для исследования дифференциальных уравнений. На практике, впрочем, используются более сложные методы, хотя принцип их работы в целом очень схож.
На рис. 2.1 синим изображено истинное решение уравнения ˙x=t с начальным условием x(−3)=4, а красным, розовым, фиолетовым и зеленым изображены численные решения уравнения методом Эйлера 5, 10, 20 и 100 шагами соответственно. Заметим, что уже сто шагов дает достаточно хорошее приближение решения.
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py
ob.axes4x4()
ob.eulersplot(lambda t, x: t, -3, 4, 4, 5, color='red')
ob.eulersplot(lambda t, x: t, -3, 4, 4, 10, color='pink')
ob.eulersplot(lambda t, x: t, -3, 4, 4, 20, color='purple')
ob.eulersplot(lambda t, x: t, -3, 4, 4, 100, color='green')
ob.mplot(np.linspace(-4,4),lambda x: x**2 / 2 - 0.5,
color='steelblue')
Упражнение 1. Используя метод Эйлера (но не используя компьютер), найти решение дифференциального уравнения ˙x=x, удовлетворяюее начальному условию x(0)=1. Отсюда найти одну из известных формул для числа e.
2.2Аналитическое решение автономных дифференциальных уравнений на прямой
Вернёмся к аналитическому поиску решений.˙x=f(x)(2.3)
Рассмотрим задачу Коши для автономного дифференциального уравнения (2.3) с начальным условием x(t0)=x0. Пусть f(x0)≠0. В этом случае решение задаётся явной формулой (она называется формулой Барроу). Мы обсудим несколько способов её вывода и интерпретации.
2.2.1Геометрические соображения
В предыдущей главе мы обсуждали, что существует специальный класс дифференциальных уравнений, которые очень просто решаются: это уравнения вида ˙x=f(t), мгновенно сводящиеся к интегрированию (см. параграф 1.2.4). Мы будем называть такие уравнения простейшими, хотя это не общепринятый термин.Рассмотрим поля направлений двух уравнений: первое является простейшим, а второе автономным, см. рис. 2.2.
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py
plt.figure(figsize=(8, 4))
plt.subplot(121)
ob.axes4x4()
ob.normdirfield(np.arange(-4, 4, 0.5),
np.arange(-4, 4, 0.5),
lambda t, x: t,
color='red',
linewidth=1,
length=0.6)
plt.subplot(122)
ob.axes4x4()
ob.normdirfield(np.arange(-4, 4, 0.5),
np.arange(-4, 4, 0.5),
lambda t, x: x,
color='red',
linewidth=1,
length=0.6)
Рис. 2.2: Поля направлений для уравнения ˙x=t (слева) и ˙x=x (справа)
Это совсем разные уравнения, но их поля направлений обладают неким сходством:
они не меняются при сдвигах. Разница в том, что первое поле направлений не
меняется при Напомним, что задача отыскания решения дифференциального уравнения имеет простую геометрическую интерпретацию: нужно найти кривую, касающуюся в каждой своей точке соответствующего поля направлений. Вместе со сходством полей направлений это даёт надежду, что нам удастся придумать метод решения автономных уравнений, сводящий их к некоторым простейшим.
Оказывается, сделать это довольно просто: достаточно поменять роль осей и считать независимой переменной x, а неизвестной функцией — время. Ниже мы обсудим два способа реализации этого замысла.
2.2.2Механический подход
Решить дифференциальное уравнение — это значит научиться отвечать на вопрос о том, где окажется решение в произвольный момент времени t, если в момент времени t0 оно находилось в точке x0. В соответствии с выводами предыдущего пункта, поменяем роли переменных и зададим другой вопрос: сколько времени потребуется, чтобы добраться из точки x0 до некоторой другой точки x?Попробуем ответить на этот вопрос (хотя бы приближённо) с помощью аналога метода Эйлера (см. раздел 2.1). Предположим для определённости, что x>x0 и f(x0)>0 (вблизи x0 движение происходит вправо; обратный случай рассматривается полностью аналогично). Предположим также, что на всём отрезке [x0,x] функция f принимает только положительные значения (чуть ниже мы обсудим, что это вполне разумное предположение). Разобьем отрезок [x0,x] на n равных маленьких отрезочков длины Δx. Пусть xk=x0+k⋅Δx — концы наших отрезочков. Сколько времени нужно, чтобы попасть из точки xk в точку xk+1? Для этого нам придётся пройти расстояние, равное Δx. Мгновенная скорость движения в точке xk равна f(xk). Поскольку Δx мало, а функция f непрерывна, можно ожидать, что её значение не слишком сильно изменится, по крайней мере, пока мы находимся на том же отрезочке. Значит, можно считать (совершая некоторую ошибку, малую при малых Δx), что движение на всём отрезочке происходит с постоянной скоростью, равной f(xk). Тогда время движения вычисляется по школьной формуле: нужно расстояние Δx поделить на скорость f(xk). Обозначим вычисленное таким образом время прохождения k-го отрезочка через Δtk. Имеем:
Δtk=Δxf(xk)(2.4)
Пусть мы оказались в точке x в момент времени t. Тогда время прохождения всего отрезка от x0 до x равна t−t0 и получается сложением всех Δtk для k=0,…,n−1:t−t0≈n−1∑k=0Δtk=n−1∑k=01f(xk)Δx(2.5)
Ну-ка, что у нас тут в правой части? Это же интегральные суммы для функции 1f(x)! Равенство (2.5) является приближённым, но когда мы перейдём к пределу при Δx→0 (или, что то же самое, при n→∞), оно превратится в точное:t−t0=∫xx0dyf(y).(2.6)
Это соотношение и называется формулой Барроу. Его можно понимать как неявное выражение x через t. В некоторых ситуациях из него можно выразить функцию x(t) явно.2.2.3Аналитический подход
Приведём более формальный вывод формулы Барроу, опирающийся на математический анализ. Пусть функция x=φ(t) является решением уравнения (2.3) и удовлетворяет начальному условию φ(t0)=x0. Рассмотрим функцию t=ψ(x), обратную к функции x=φ(t). Рассмотрим произвольную точку (t1,x1), лежащую на графике решения: для неё выполняются соотношения x1=φ(t1), t1=ψ(x1) и ˙φ(t1)=f(x1) (поскольку φ является решением уравнения). Тогда по теореме о производной сложной функцииψ′(x1)=1˙φ(t1)=1f(x1)
Это равенство выполняется в любой точке x1. Значит, функция ψ является решением дифференциального уравненияψ′(x)=1f(x),
где x выступает в роли независимой переменной. Правая часть теперь не зависит от неизвестной функции и такое уравнение мы умеем решать:ψ(x)=∫xx0dyf(y)+t0
Вспоминая, что t=ψ(x) — обратная функция к решению x=φ(t), имеем:t−t0=∫φ(t)x0dyf(y)
Мы снова получили формулу Барроу.2.2.4Магия
Часто для вывода формулы Барроу используют такую символическую запись: ˙x=f(x)dxdt=f(x)dt=dxf(x)∫tt0dt=∫xx0dxf(x)(2.7)(2.8)(2.9)(2.10) Это может показаться некоторой магией — особенно загадочно выглядит уравнение (2. 9). Чуть позже мы дадим определение дифференциальной 1-формы, с помощью которого можно придать этим формулам аккуратный смысл, а пока обратим внимание, что уравнение (2.9) очень похоже на уравнение (2.4). В целом, эта формальная запись фактически повторяет наш вывод в параграфе 2.2.2. Пример 1. Решим уравнение ˙x=x с начальным условием x(t0)=x0. Пусть x=φ(t) — решение и t=φ−1(x) — обратная функция к решению. Имеем: (φ−1(x))′=1xφ−1(x)=ψ(x)ψ′(x)=1xψ(x)=∫dxx=ln|x|+Ct=ln|x|+Cx=±e−Cet=C1et Заметим, что если бы мы забыли модуль под логарифмом при интегрировании, то константа C1=e−C принимала бы только положительные значения. Но из-за модуля она может принимать и отрицательные значения.Заметим также, что в ходе преобразований (деления на x) мы «потеряли» решение x=0. Если в ответ подставить значение C1=0, получим как раз его. Таким образом, формула x(t)=Cet, C∈R даёт все известные нам решения. Мы пока не доказали, что других нет — на следующей лекции мы обсудим этот вопрос.
← Предыдущая глава Следующая глава →
Дифференциальное уравнение — Differential equation
Математическое уравнение с производными неизвестной функции
В математике дифференциальное уравнение — это уравнение, которое связывает одну или несколько функций и их производные . В приложениях функции обычно представляют физические величины, производные представляют скорость их изменения, а дифференциальное уравнение определяет взаимосвязь между ними. Такие отношения обычны; поэтому дифференциальные уравнения играют важную роль во многих дисциплинах, включая инженерию , физику , экономику и биологию .
В основном изучение дифференциальных уравнений состоит из изучения их решений (набора функций, удовлетворяющих каждому уравнению) и свойств их решений. Только простейшие дифференциальные уравнения разрешимы по явным формулам; однако многие свойства решений данного дифференциального уравнения могут быть определены без их точного вычисления. {п} \,}
для которой в следующем году Лейбниц получил решения, упростив ее.
Исторически проблема вибрирующей струны, такой как струна музыкального инструмента, изучалась Жаном ле Рондом Даламбером , Леонардом Эйлером , Даниэлем Бернулли и Жозефом-Луи Лагранжем . В 1746 году Даламбер открыл одномерное волновое уравнение , а через десять лет Эйлер открыл трехмерное волновое уравнение.
Уравнение Эйлера – Лагранжа было разработано в 1750-х годах Эйлером и Лагранжем в связи с их исследованиями проблемы таутохрон . Это проблема определения кривой, на которой взвешенная частица упадет в фиксированную точку за фиксированный промежуток времени, независимо от начальной точки. Лагранж решил эту проблему в 1755 году и отправил решение Эйлеру. Оба далее развили метод Лагранжа и применили его к механике , что привело к формулировке лагранжевой механики .
В 1822 году Фурье опубликовал свою работу о тепловом потоке в Théorie analytique de la chaleur (Аналитическая теория тепла), в которой он основывал свои рассуждения на законе охлаждения Ньютона , а именно, что поток тепла между двумя соседними молекулами пропорционален крайне малая разница их температур. В этой книге содержится предложение Фурье о его уравнении теплопроводности для кондуктивной диффузии тепла. Это уравнение в частных производных теперь преподается каждому студенту математической физики.
Пример
В классической механике движение тела описывается его положением и скоростью при изменении значения времени. Законы Ньютона позволяют динамически выражать эти переменные (с учетом положения, скорости, ускорения и различных сил, действующих на тело) в виде дифференциального уравнения для неизвестного положения тела как функции времени.
В некоторых случаях это дифференциальное уравнение (называемое уравнением движения ) может быть решено явно.
Примером моделирования реальной проблемы с использованием дифференциальных уравнений является определение скорости шара, падающего в воздухе, с учетом только силы тяжести и сопротивления воздуха. Ускорение мяча по направлению к земле — это ускорение свободного падения за вычетом замедления из-за сопротивления воздуха. Сила тяжести считается постоянной, а сопротивление воздуха можно моделировать пропорционально скорости мяча. Это означает, что ускорение мяча, которое является производной от его скорости, зависит от скорости (а скорость зависит от времени). Для определения скорости как функции времени необходимо решить дифференциальное уравнение и проверить его справедливость.
Типы
Дифференциальные уравнения можно разделить на несколько типов. Помимо описания свойств самого уравнения, эти классы дифференциальных уравнений могут помочь в выборе подхода к решению. Обычно используемые различия включают в себя то, является ли уравнение обычным или частным, линейным или нелинейным, однородным или неоднородным. Этот список далеко не исчерпывающий; существует множество других свойств и подклассов дифференциальных уравнений, которые могут быть очень полезны в определенных контекстах.
Обыкновенные дифференциальные уравнения
Обыкновенное дифференциальное уравнение ( ОДУ ) представляет собой уравнение , содержащее неизвестную функцию одной действительной или комплексной переменной х , его производных, а также некоторые заданные функции х . Неизвестная функция обычно представлена переменной (часто обозначаемой y ), которая, следовательно, зависит от x . Таким образом, x часто называют независимой переменной уравнения. Термин « обыкновенный » используется в отличие от термина « уравнение в частных производных» , которое может относиться к более чем одной независимой переменной.
Линейные дифференциальные уравнения — это дифференциальные уравнения, линейные относительно неизвестной функции и ее производных. Их теория хорошо разработана, и во многих случаях их решения можно выразить в терминах интегралов .
Большинство ОДУ, встречающихся в физике, являются линейными. Следовательно, большинство специальных функций могут быть определены как решения линейных дифференциальных уравнений (см. Голономная функция ).
Поскольку, как правило, решения дифференциального уравнения не могут быть выражены выражением в замкнутой форме , численные методы обычно используются для решения дифференциальных уравнений на компьютере.
Уравнения с частными производными
Парциальное дифференциальное уравнение ( ФДЭ ) представляет собой дифференциальное уравнение , которое содержит неизвестные многомерные функции и их частные производные . (В этом отличие от обычных дифференциальных уравнений , которые имеют дело с функциями одной переменной и их производными.) УЧП используются для формулирования задач, включающих функции нескольких переменных, и либо решаются в замкнутой форме, либо используются для создания соответствующего компьютера. модель .
PDE могут использоваться для описания широкого спектра явлений в природе, таких как звук , тепло , электростатика , электродинамика , поток жидкости , упругость или квантовая механика . Эти, казалось бы, различные физические явления можно формализовать аналогичным образом в терминах PDE. Так же, как обыкновенные дифференциальные уравнения часто моделируют одномерные динамические системы , уравнения в частных производных часто моделируют многомерные системы . Стохастические дифференциальные уравнения в частных производных обобщают уравнения в частных производных для моделирования случайности .
Нелинейные дифференциальные уравнения
Нелинейный дифференциальное уравнение представляет собой дифференциальное уравнение , которое не является линейным уравнением относительно неизвестной функции и ее производных (линейность или нелинейность в аргументах функции не рассматривается здесь). Существует очень мало методов точного решения нелинейных дифференциальных уравнений; те, которые известны, обычно зависят от уравнения, имеющего определенные симметрии . Нелинейные дифференциальные уравнения могут демонстрировать очень сложное поведение на длительных временных интервалах, характерное для хаоса . Даже фундаментальные вопросы существования, единственности и расширяемости решений нелинейных дифференциальных уравнений, а также корректности начальных и краевых задач для нелинейных уравнений в частных производных являются сложными задачами, и их разрешение в частных случаях считается значительным достижением в математике. теория (ср. существование и гладкость Навье – Стокса ). Однако, если дифференциальное уравнение является правильно сформулированным представлением значимого физического процесса, то можно ожидать, что оно имеет решение.
Линейные дифференциальные уравнения часто появляются как приближения к нелинейным уравнениям. Эти приближения действительны только при ограниченных условиях. Например, уравнение гармонического осциллятора является приближением к уравнению нелинейного маятника, которое справедливо для колебаний малой амплитуды (см. Ниже).
Порядок уравнения
Дифференциальные уравнения описываются их порядком, определяемым слагаемым с старшими производными . Уравнение, содержащее только первые производные, является дифференциальным уравнением первого порядка, уравнение , содержащее вторую производную, является дифференциальным уравнением второго порядка и т. Д. Дифференциальные уравнения, описывающие природные явления, почти всегда содержат только производные первого и второго порядка, но есть некоторые исключения, такие как уравнение тонкой пленки , которое является уравнением в частных производных четвертого порядка.
Примеры
В первой группе примеров u — неизвестная функция от x , а c и ω — константы, которые предполагается известными. Две широкие классификации как обыкновенных дифференциальных уравнений, так и уравнений в частных производных состоят из различения линейных и нелинейных дифференциальных уравнений, а также однородных дифференциальных уравнений и неоднородных . {n}}} + \ cdots + f_ {1} (x) {\ frac {dy} {dx}} + f_ {0} (x) y = g (x)}
такой, что
- y ( Икс 0 ) знак равно y 0 , y ′ ( Икс 0 ) знак равно y 0 ′ , y ″ ( Икс 0 ) знак равно y 0 ″ , … {\ displaystyle y (x_ {0}) = y_ {0}, y ‘(x_ {0}) = y’ _ {0}, y » (x_ {0}) = y » _ {0}, \ ldots}
Для любого ненулевого , если и непрерывны на некотором интервале, содержащем , единственно и существует. ж п ( Икс ) {\ displaystyle f_ {n} (x)} { ж 0 , ж 1 , … } {\ Displaystyle \ {е_ {0}, е_ {1}, \ ldots \}} грамм {\ displaystyle g} Икс 0 {\ displaystyle x_ {0}} y {\ displaystyle y}
Связанные понятия
- Дифференциальное уравнение задержки (DDE) представляет собой уравнение для функции одной переменной, обычно называют время , в котором производная функции в определенный момент времени дается в терминах значений функции в прежние времена.
- Интегро-дифференциальное уравнение (IDE) представляет собой уравнение , которое сочетает в себе аспекты дифференциального уравнения и к интегральному уравнению .
- Стохастическое дифференциальное уравнение (СД) является уравнением , в котором искомая величина является случайным процессом и уравнение включает некоторые известные случайные процессы, например, процесс Винер в случае уравнений диффузии.
- Стохастическое дифференциальное уравнение (SPDE) представляет собой уравнение , которое обобщает СДУ , чтобы включать в себя процессы шума пространства-времени, с приложениями в квантовой теории поля и статистической механики .
- Дифференциальное алгебраическое уравнение (ДАЭ) представляет собой дифференциальное уравнение , содержащее дифференциальные и алгебраические термины, приведенные в неявном виде.
Связь с разностными уравнениями
Теория дифференциальных уравнений тесно связана с теорией разностных уравнений , в которой координаты принимают только дискретные значения, а взаимосвязь включает значения неизвестной функции или функций и значения в близких координатах. Многие методы вычисления численных решений дифференциальных уравнений или изучения свойств дифференциальных уравнений включают аппроксимацию решения дифференциального уравнения решением соответствующего разностного уравнения.
Приложения
Изучение дифференциальных уравнений — обширная область чистой и прикладной математики , физики и техники . Все эти дисциплины связаны со свойствами дифференциальных уравнений различных типов. Чистая математика фокусируется на существовании и единственности решений, в то время как прикладная математика подчеркивает строгое обоснование методов приближения решений. Дифференциальные уравнения играют важную роль в моделировании практически всех физических, технических или биологических процессов, от движения небесных тел до конструкции мостов и взаимодействий между нейронами. Дифференциальные уравнения, такие как те, которые используются для решения реальных задач, не обязательно могут быть решаемы напрямую, т. Е. Не имеют решений в замкнутой форме . Вместо этого решения можно аппроксимировать численными методами .
Многие фундаментальные законы физики и химии можно сформулировать в виде дифференциальных уравнений. В биологии и экономике дифференциальные уравнения используются для моделирования поведения сложных систем. Математическая теория дифференциальных уравнений впервые развивалась вместе с науками, в которых эти уравнения возникли и где результаты нашли применение. Однако различные проблемы, иногда возникающие в совершенно разных областях науки, могут привести к идентичным дифференциальным уравнениям. Когда бы это ни происходило, математическую теорию, лежащую в основе уравнений, можно рассматривать как объединяющий принцип, лежащий в основе различных явлений. В качестве примера рассмотрим распространение света и звука в атмосфере и волн на поверхности пруда. Все они могут быть описаны одним и тем же уравнением в частных производных второго порядка , волновым уравнением , которое позволяет нам думать о свете и звуке как о формах волн, очень похожих на знакомые волны в воде. Теплопроводность, теория которой была разработана Джозефом Фурье , регулируется другим уравнением в частных производных второго порядка — уравнением теплопроводности . Оказывается, многие диффузионные процессы, хотя и кажутся разными, описываются одним и тем же уравнением; Блэка-Шоулза уравнение финансов, например, связанные с уравнением теплопроводности.
Количество дифференциальных уравнений, получивших название, в различных областях науки свидетельствует о важности данной темы. См. Список названных дифференциальных уравнений .
Программное обеспечение
Некоторые программы CAS могут решать дифференциальные уравнения. Стоит упомянуть эти программы CAS и их команды:
Смотрите также
Рекомендации
дальнейшее чтение
- Abbott, P .; Нил, Х. (2003). Научитесь исчислению . С. 266–277.
- Blanchard, P .; Девани, Р.Л . ; Холл, GR (2006). Дифференциальные уравнения . Томпсон.
- Boyce, W .; DiPrima, R .; Мид, Д. (2017). Элементарные дифференциальные уравнения и краевые задачи . Вайли.
- Коддингтон, EA; Левинсон, Н. (1955). Теория обыкновенных дифференциальных уравнений . Макгроу-Хилл.
- Инс, Э.Л. (1956). Обыкновенные дифференциальные уравнения . Дувр.
- Джонсон, В. (1913). Трактат об обыкновенных и дифференциальных уравнениях с частными производными . Джон Вили и сыновья. В исторической математической коллекции Мичиганского университета
- Полянин А.Д .; Зайцев, В. Ф. (2003). Справочник точных решений обыкновенных дифференциальных уравнений (2-е изд.). Бока-Ратон: Chapman & Hall / CRC Press. ISBN 1-58488-297-2 .
- Портер, Р.И. (1978). «XIX Дифференциальные уравнения». Дальнейший элементарный анализ .
- Тешл, Джеральд (2012). Обыкновенные дифференциальные уравнения и динамические системы . Провиденс : Американское математическое общество . ISBN 978-0-8218-8328-0 .
- Даниэль Цвиллинджер (12 мая 2014 г.). Справочник по дифференциальным уравнениям . Elsevier Science. ISBN 978-1-4832-6396-0 .
Внешние ссылки
Откуда берутся дифференциальные уравнения? // Владимир Побережный ≪ ∀ x, y, z
Математик Владимир Побережный об экспонентах, источниках дифференциальных уравнений и векторном пространстве функций.
Что такое дифференциальные уравнения? Это уравнения на какую-то неизвестную функцию или соотношения, которым должна удовлетворять эта функция и какие-то ее производные (если функция одной переменной, то просто производные, если функция многих переменных, то частные производные). Это обобщение наших обычных уравнений, например алгебраических. Мы сначала учим в школе линейные уравнения, их графики дают прямые на плоскости — бывают квадратичные, кубические и так далее. Это все алгебраические уравнения. Можно брать более сложные функции и более сложные уравнения, они дают какие-то более сложные графики. Объекты, которые они описывают, становятся более сложными, то есть линейные уравнения рисуют прямые, квадратичные — параболы, это все какие-то графики на плоскости или в более общем случае в большой размерности, какие-то поверхности в пространстве той или другой размерности. Поверхности или более сложные объекты, сделанные из поверхностей, — так называемые многообразия и так далее.
Дифференциальные уравнения — это следующий шаг. Уравнения, которые мы сейчас перечислили, задают в пространстве какие-то точки, подмножества точек. Уравнение задает множество точек на плоскости, и мы знаем, что эти точки выглядят как прямая. Это и есть график. Дифференциальные уравнения тоже задают какие-то подмножества, но они заданы уже в пространстве функций, то есть это соотношения, которым удовлетворяют функции. Решение дифференциального уравнения — это какой-то набор подмножества точек в пространстве функций. Пространство функций является бесконечномерным.
Возникает нужда в анализе: как это все устроено и почему мы вообще на это так смотрим? Такой взгляд действительно имеет вполне разумное содержание и смысл. Если мы рассматриваем линейные дифференциальные уравнения, то у нас возникает аналогия с обычными линейными уравнениями. Например, мы знаем, что линейные уравнения на плоскости — это прямая, в пространстве — какая-то гиперплоскость. То есть это какой-то плоский объект. Оказывается, что множество функций, удовлетворяющих линейному дифференциальному уравнению, устроено примерно так же, это в каком-то смысле плоскость, или прямая, или плоскость какой-то размерности, но уже в бесконечномерном пространстве функций (официально это называется векторным пространством). Множество решений линейного дифференциального уравнения образует векторное пространство во множестве всех функций.
Откуда берутся дифференциальные уравнения? Конечно, основной поставщик дифференциальных уравнений (это мы тоже со школы знаем) — это физика и механика. Законы Ньютона, например, ускорение материальной точки (силе, которая на нее действует). Но ускорение — это вторая производная. Вот у вас получилось дифференциальное уравнение (вторая производная координаты) равна какой-то силе . Свойство классической механики состоит в том, что, как правило, уравнения там второго порядка. Видимо, оттуда это возникло, причем, как принято у физиков (это не редкость), дифференциальные уравнения возникли чуть ли не раньше дифференциального исчисления, и решать их тоже (конечно, без построения общей теории) люди начали раньше, чем все эти понятия вообще были определены, и добивались каких-то успехов. Мы знаем, что введение основ дифференциального исчисления произошло как раз во времена Ньютона и Лейбница, то есть практически одновременно с законом Ньютона, в котором уже есть дифференцирование.
Физика не единственный источник этих уравнений. Практически любая околоестественная наука является таким источником. Например, в химии происходят какие-то реакции, скорость реакций зависит от количества и пропорций компонентов. Два вещества смешиваются и как-то превращаются в третье с какой-то скоростью, пропорциональной чему-то. Это дифференциальные уравнения. В биологии тоже есть дифференциальные уравнения.
Конечно, это не биология, а какой-то детский пример. Есть стандартная задача о размножении кроликов. У вас есть парочка кроликов, они с какой-то периодичностью рожают еще пару. У вас была пара кроликов, она родила — стало две пары. Каждая пара еще родила — стало четыре и так далее. Как устроен закон? Видно, что число растет очень быстро, это экспоненциальный рост. Здесь возникает очень интересный, но уже не совсем математический вопрос моделестроительства или адекватного построения модели. Вот мы хотим описать размножение кроликов. Если мы его описываем таким образом, то легко подсчитать, что если уравнение устроено так, что (это из физики идет такое стандартное обозначение; вообще производные функций обычно обозначаются , но если производная по времени, то ее удобно обозначать ) равняется , то есть скорость роста равна числу уже имеющихся пар. Такие уравнения мы умеем решать, это экспонента.
Эта модель, очевидно, не дает нам правильного приближения к жизни, на маленьких порядках немножко дает. С другой стороны, если бы все было в жизни устроено так, то кролики очень быстро бы захватили всю землю во много слоев, некуда было бы между ними наступить. Значит, надо как-то менять наше уравнение, подстраивать свойства модели под картинку, которую мы наблюдаем в жизни, и то, чему хотим быть адекватными. Например, чем больше кроликов, чем чаще они встречаются, тем больше вероятность, что у них возникнет какая-нибудь болезнь, которая будет заразной и будет передаваться от одного к другому, то есть надо вычесть какое-то слагаемое, пропорциональное частоте встреч. А как устроена частота встреч? Если кролики живут в каком-то лесу, каждый кролик занимает какое-то место, надо поделить площадь леса на площадь кроликов и так далее.
Стандартное, вполне обозримое и разумное приближение. Например, добавление в модель волков. У нас есть волки, есть кролики. Кролики как-то размножаются, и волки как-то размножаются. Кроликам для размножения нужен только лес и другие кролики, а волкам нужно что-то есть, им нужны, собственно, кролики. Поэтому скорость роста кроликов (), с одной стороны, равна числу пар (какому-то слагаемому ). С другой стороны, вычитается какое-то неудобство из-за перенаселенности, из-за ограниченности площади. С третьей стороны, вычитается какая-то пропорциональность числу волков, каждый волк кого-то съедает. А волки, в свою очередь, размножаются пропорционально своему имеющемуся числу (не как кролики, но все-таки), к тому же им надо что-то кушать, к тому же они тоже болеют. У нас получается набор, система уравнений. — это наши кролики, а , допустим, волки. Эти два уравнения должны выполняться одновременно, так модель усложняется и усложняется.
Даже в классической механике мы знаем, что если бросаем камень, то вблизи Земли у него ускорение постоянно . Но мы можем, например, добавлять сопротивление воздуха, оно уже зависит от скорости камня, то есть вторая производная будет не , а минус еще какое-то слагаемое, пропорциональное скорости . Например, падает дождевая капля. Во-первых, она падает из-за силы тяжести, во-вторых, тормозится воздухом, в-третьих, если воздух влажный, то она еще и конденсируется, растет, вбирает влажность из окружающего воздуха, то есть у нее меняется масса.
Можно строить разные модели, как-то их усложнять, исследовать те интересные вопросы, которые возникают почти в любом приложении, где как-то используется математика. Но математика ради математики здесь тоже имеется: дифференциальные уравнения — это очень большой отдельный разнообразный раздел со множеством вариаций. Он настолько большой, что даже практически не бывает конференций по дифференциальным уравнениям, потому что нужно более тонкое деление: качественная теория, асимптотические методы, интегрируемые системы, уравнения в частных производных и так далее. Это вполне большая развитая наука, продолжающая развиваться.
Какие основные свойства и характеристики есть у дифференциальных уравнений? Что можно о них сказать? Во-первых, краеугольный камень для обыкновенных дифференциальных уравнений для одной переменной (неважно, вещественной или комплексной, комплексной даже лучше, как всегда это устроено в анализе) — это теорема существования и единственности. Если у вас есть дифференциальное уравнение с достаточно разумными коэффициентами (эти слова формализуются разными способами, например гладкие) и есть начальные данные, то всегда есть локальное решение. Например, вы знаете, что ваш камень как-то падает, знаете, где он был в начальный момент времени и какая у него была в начальный момент времени скорость. После этого у него траектория считается по крайней мере локально, в окрестности этого положения.
Это очень сильный результат, опять-таки похожий на то, что у нас было с обычными уравнениями: мы знаем, что алгебраическое уравнение -того порядка имеет корней. В школе, конечно, учат, что бывает меньше, а потом если кто доучивается дальше, то учит, что нет, на самом деле столько же. Здесь есть аналогия: если уравнение -того порядка, то у него не решений, конечно, их бесконечно много, но множество решений параметризуется параметрами . Если есть уравнение второго порядка (наш камень), надо задать начальное положение и начальную скорость. И вообще, для уравнения -того порядка надо задать начальных данных, и тогда будет всегда существовать решение. Если уравнение линейное, то эти начальных данных — это просто его координаты в -мерном конечномерном векторном пространстве решений.
Это специфика обыкновенных уравнений от одной переменной, но при этом все-таки уравнение локально решается, то есть мы знаем, что решение существует, а вот найти его мы в явном виде можем не всегда. Мы можем использовать какие-то приближенные методы, как-то бороться, но гарантий, что мы напишем какое-то конечное выражение и оно будет решать наше уравнение, нет.
Это была деятельность XIX века, когда люди активно занимались этой областью и изучали уравнения математической физики, из этого возникла целая наука про классические многочленные специальные функции Лежандра, Лагерра, Чебышева. Это была попытка как-то решать уравнения, которые возникали при тогдашнем развитии науки. В явном и конечном виде решения не выписывались, но это совершенно не мешало заниматься их анализом: исследовать свойства, связи, асимптотики. Современная наука занимается более сложными уравнениями. Сейчас, например, вполне популярная деятельность — исследование уравнений Пенлеве. Это такие новые специальные функции — решения уравнений Пенлеве, сейчас занимаются их исследованиями, асимптотикой, связями, геометрическим смыслом, содержанием и так далее по аналогии с физикой XIX века.
Владимир Побережный, кандидат физико-математических наук, научный сотрудник международной лаборатории теории представлений и математической физики, доцент факультета математики НИУ ВШЭПостНаука
Однородные дифференциальные уравнения второго порядка с постоянными коэффициентами
- Главная
- Справочник
- Алгебра
- Однородные дифференциальные уравнения второго порядка с постоянными коэффициентами
Однородные линейные дифференциальные уравнения второго порядка с постоянными коэффициентами имеют вид
\[y» + py’ + q = 0,\]
где p и q — действительные числа. Рассмотрим на примерах, как решаются однородные дифференциальные уравнения второго порядка с постоянными коэффициентами.
Решение линейного однородного однородного дифференциального уравнения второго порядка зависит от корней характеристического уравнения. Характеристическое уравнение — это уравнение k²+pk+q=0.
1) Если корни характеристического уравнения — различные действительные числа:
\[{k_1} \in R,{k_2} \in R,{k_1} \ne {k_2},\]
то общее решение линейного однородного дифференциального уравнения второго порядка с постоянными коэффициентами имеет вид
\[y = {C_1}{e^{{k_1}x}} + {C_2}{e^{{k_2}x}}\]
2) Если корни характеристического уравнения — равные действительные числа
\[{k_1} \in R,{k_2} \in R,{k_1} = {k_2}\]
(например, при дискриминанте, равном нулю), то общее решение однородного дифференциального уравнения второго порядка есть
\[y = {e^{{k_1}x}}({C_1} + {C_2}x)\]
3) Если корни характеристического уравнения — комплексные числа
\[{k_{1,2}} = \alpha \pm \beta i\]
(например, при дискриминанте, равном отрицательному числу), то общее решение однородного дифференциального уравнения второго порядка записывается в виде
\[y = {e^{\alpha x}}({C_1}\cos \beta x + {C_2}\sin \beta x). \]
В вашем браузере отключен Javascript.Чтобы произвести расчеты, необходимо разрешить элементы ActiveX!
Источник
Больше интересного в телеграм @calcsboxДифференциальные уравнения численное решение — Справочник химика 21
Все проблемы, рассмотренные в этой главе, сводятся к решению обыкновенных дифференциальных уравнений. Мы уже замечали, что в некоторых случаях аналитическое решение невозможно, н решать задачу приходится численными методами. Существуют стандартные программы решения уравнений такого типа на вычислительных машинах. Тем не менее, знакомство с численными методами интегрирования уравнений полезно химику-технологу по двум важным причинам. Во-первых, вопреки распространенному мнению, вычислительная машина не умеет думать , и потому небезопасно давать ей задание, не имея понятия о том, как она его будет выполнять. Во-вторых, иногда возможно и даже желательно проводить вычисления вручную. Метод, который мы сейчас рассмотрим, применим к решению любой системы обыкновенных дифференциальных уравнений, включая уравнения, описывающие неизотермические процессы. Проиллюстрируем этот метод на примере одного уравнения и системы двух уравнений. [c.114]Общее решение этого дифференциального уравнения численно возможно двумя путями. Во-первых, рассмотрение баланса может проводиться с маленькими, но конечными Аи-значениями, причем получается [c.123]
Устойчивость решения. При решении дифференциальных уравнений численными методами помимо вопросов точности важную роль приобретают вопросы устойчивости решения. Под устойчивостью метода решения дифференциального уравнения понимается способность накопления и скорость роста ошибки интегрирования. Как уже отмечалось выше, при использовании формулы (12—17) ошибка вычислений накапливается в процессе интегрирования. [c.356]
Раздел V. 9. Относительно вопросов численного решения дифференциальных уравнений см. [c.118]
Численное интегрирование уравнения Эйлера. Как уже отмечалось выше, уравнение Эйлера (V, 133) обычно представляет собой нелинейное дифференциальное уравнение, аналитическое решение которого чаще всего найти нельзя. Кроме того, весьма существенным является то, что решение уравнения (V, 133) должно удовлетворять граничным условиям в двух точках экстремали, которые в простейшем случае имеют вид соотношений (V, 135). [c.227]
Рассмотренная неустойчивость решения является серьезным препятствием при решении дифференциальных уравнений численными методами, когда невольно приходится ограничиваться конечной точностью представления чисел, в результате чего погрешность решения может достигать значительной величины. Следует отметить, что если при решении одного дифференциального уравнения первого порядка еще можно предусмотреть, некоторые методы. устранения неустойчивости, то при интегрировании систем дифференциальных уравнений задача обеспечения устойчивости решения становится весьма серьезной и иногда даже непреодолимой на пути получения решения оптимальной задачи. [c.231]
Многие процессы с распределенными параметрами, которые на первый взгляд нельзя представить как многостадийные из-за непрерывности изменения величин, определяющих их состояние и управление (например, реактор вытеснения), могут быть описаны как предельный случай многостадийного процесса, если в качестве отдельной стадии принять достаточно малый элемент, аналогично тому, как при решении дифференциальных уравнений численными методами используется их конечно-разностная форма. [c.258]
В частных производных не сводятся к обыкновенным дифференциальным уравнениям. Поэтому решения следует получать другими методами, например численными или методом локальной [c.184]
Задание 114. Решите в численном виде приведенные ниже дифференциальные уравнения найдите решения при нескольких значениях независимой переменной и нанесите найденные точки на диаграмму. Сравните результаты численного и аналитического метода (в тех случаях, когда уравнение можно решить аналитически). [c.223]
После подстановки соотношений (2.54), (2.55) в уравнения (2.56), (2.57) последние преобразуются аналогично тому, как это сделано в примере 2.1, к системе обыкновенных нелинейных дифференциальных уравнений. Их решения можно получить одним из численных методов. После чего производится расчет по соотношениям (2.58), (2.59). [c.95]
Постоянные интегрирования Су (г = 1, 2,. .8 / — О, 1, 2) из (9.4), (9.5) необходимо определять из указанных краевых условий и условий непрерывности (9.8) и (9.10) при 1, которые дадут систему 24 линейных алгебраических уравнений. Численное решение указанной системы уравнении затруднено, так как матрица данной системы плохо обусловлена. Это связано с тем, что коэффициенты системы дифференциальных уравне- [c.57]
B. Э. M илн, Численное решение дифференциальных уравнений, ИЛ, 1955. [c.118]
Исследовать внутреннюю диффузию нри конечной скорости адсорбции гораздо труднее, поскольку мы сразу же сталкиваемся с нелинейными дифференциальными уравнениями. Общий метод, описанный в конце предыдущего раздела, можно применить к решению уравнений с кинетическими зависимостями типа (VI.20). Получить какие-либо общие результаты здесь, однако, трудно, вследствие большого числа параметров, входящих в кинетическую зависимость, и необходимости численного интегрирования. [c.141]
Следовательно, расчет реактора сводится к решению системы двух уравнений с двумя неизвестными а и Т. Для проточного реактора полного перемешивания это будет система алгебраических уравнений. В остальных случаях получается система дифференциальных уравнений. Как правило, для решения необходимо использовать численные или графические методы. Ниже будет рассмотрено несколько примеров расчета неизотермических реакторов. [c.332]
Дифференциальные уравнения в частных производных получаются в тех случаях, когда рассматривается одновременное изменение более чем двух переменных. Для этих уравнений справедливы те же соображения, какие были высказаны по поводу обыкновенных дифференциальных уравнений. Для полного математического описания физической проблемы, помимо самого дифференциального уравнения, необходимы еще дополнительные указания начальные условия, из которых определяются константы, возникающие при интегрировании обыкновенных дифференциальных уравнений, а также начальные и граничные условия, из которых находятся параметры, полученные при точном решении дифференциальных уравнений в частных производных. (Разумеется, начальные и граничные условия в равной мере необходимы и при численных методах.—Прим. ред.) [c.385]
Чтобы представить дифференциальные уравнения в форме, пригодной для решения на цифровых вычислительных машинах, следует их аппроксимировать и заменить конечно-разностными алгебраическими уравнениями. Численная модель состоит из полученной системы уравнений и построения численного алгоритма их решения. Численные модели основных процессов фильтрации пластовых флюидов обсуждаются в следующих параграфах. [c.381]
При моделировании на ЦВМ получается совокупность чисел, отражающих конечный результат протекания процесса. Картину же изменения внутренних связей между физико-химическими величинами в ходе решения получить нельзя. Причиной этого является сам принцип дискретности работы цифровой машины и вытекающая отсюда при решении необходимость предварительного преобразования дифференциального уравнения методами численного анализа. Естественно, что это в некоторой степени обесценивает результаты моделирования на ЦВМ. Однако возможность получения значительного объема числового материала при моделировании различных вариантов частично компенсирует [c.11]
Аналитические решения рассмотренных выше дифференциальных уравнений первого и второго порядка известны лишь для частных случаев с единичными простыми реакциями в изотермических условиях. Поэтому для интегрирования их в настоящее время используются в основном численные методы и решение производится на ЭЦВМ [6, 7, 68]. [c.43]
Методика отыскания численных значений вероятностных характеристик по экспериментально найденным распределениям общеизвестна и детально описана во многих руководствах по математической статистике, например в работах [74, 80]. Поэтому, опуская непосредственно вычисление указанных характеристик, установим лишь связь между ними и числами Пекле. Эта связь определяет т из решения дифференциального уравнения диффузионной модели, составленного применительно к изменению концентрации [c.49]
При численном решении уравнения (4.127) при условиях (4.128) вначале интегрировалось обыкновенное дифференциальное уравнение [c.202]
В университете штата Канзас (где преподает автор—доп. ред.) в начале семестра одна неделя отводится ознакомлению студентов с математическими методами, примерно в объеме, соответствующем объему главы XII этой книги. Сюда относится знакомство с типами дифференциальных уравнений, часто встречающимися в учении о химической кинетике, и методами численного интегрирования. Приближенные методы расчета находят широкое применение, так как экономят время и труд, а точность получаемых решений обычно вполне соответствует точности исходных экспериментальных данных. Применение указанных методов в тексте сохраняет элементарный характер изложения, принятый нами для настоящей книги. Точные решения, как правило, настолько сложны, что их использование могло бы оттолкнуть начинающего и затруднило бы понимание основных идей. [c.10]
В некоторых случаях решение дифференциального уравнения в частных производных может быть сведено к решению обыкновенных дифференциальных уравнений. Уравнение в декартовых координатах приводит к обыкновенным линейным дифференциальным уравнениям второго порядка с постоянными коэффициентами, решение которых выражается в виде показательных или тригонометрических функций. Цилиндрические координаты ведут к обыкновенным дифференциальным уравнениям, решение которых имеет вид бесконечных рядов, называемых функциями Бесселя. Метод решения дифференциального уравнения в частных производных может быть пояснен примером в декартовых координатах, поскольку свойства тригонометрических функций, возможно, лучше известны, чем свойства функций Бесселя. Ниже будут показаны как аналитическое, так и численное решения. [c.247]
Пример VIU-2. Применим к предыдущему примеру численный метод решения дифференциальных уравнений в частных производных, описанный в главе XII (стр. 399). [c.250]
Матрица решения системы обыкновенных дифференциальных уравнений численным методом Рунге—Ку гга на интервале от х1 до х2 с переменным шагом, при минимальном числе шагов п, причем правые части уравнений в символьной форме задаются в векторе D, а начальные условия — в векторе V (только для Math ad Professional) Матрица решения системы обыкновенных дифференциальных уравнений методом Рунге—Кутта на интервале от х1 до х2 при фиксированном числе шагов п, причем правые части уравнений записаны в символьном векторе D, а начальные условия — в векторе v [c.453]
Весьма затруднительно дать краткое изложение методов решений дифференциальных уравнений в частных производных . Более подробные сведения можно найти в литературе . Пример точного решения дифференциального уравнения в частных производных приведен выше (см. стр. 246). Нахождение точных решений таких уравнений часто довольно трудно. В таких случаях необходимо прибегать к численным методам. Существует большое число методов для решения дифференциальных [c.385]
Численное решение дифференциальных уравнений [c.395]
Если величины Ь, с н а не зависят от , то уравнения (137а и (1376) сводятся к системе обыкновенных дифференциальных уравнений. Их решение называется автомодельным. В таком случае профили скорости в различных точках X отличаются только масштабными множителями. Таким же образом ведут себя и профили температуры. Большое число автомодельных решений было рассчитано численно [72—75]. [c.112]
Последняя строка в этой таблице представляет табличную форму написи функции/I ( ), удовлетворяющей уравнению (3). Изложенный метод расчета принципиально может быть применен для всех типов дифференциальных уравнений, но точность его невысока. В настоящее время проводятся работы 1 , цель которых повысить точность II скорость решения дифференциальных уравнений численными методами. [c.37]
Хотя многократное рассеяние происходит по тем же законам, что и однократное, расчет интенсивности света, прошедшего через плотное облако, представляет значительные математические трудности. Эта проблема исследовалась в разных направлениях и известна как проблема переноса излучения (Чандрасекар. Точное решение получено только для весьма идеализированных условий, в основном для изотропного рассеяния и для случаев точечного источника света и сферических рассеивающих систем, а также плоского источника и плоскопараллельных систем. Практическая важность этой оптической проблемы и аналогичной проблемы рассеяния нейтронов плотными материалами способствовала разработке нескольких приближенных теорий. Можно получить решение некоторых задач, используя упрощенные методы расчета индикатрисы однократного рассеяния. Чу и Черчилль предложили шестипоточную модель, в которой индикатриса рассеяния одной сферой представлена в виде суммы шести компонент направленной вперед, направленной назад и четырех равных боковых. Интегродиф-ференвд1альное уравнение, описывающее интенсивность излучения в рассеивающей среде, сводится при этом к системе обычных дифференциальных уравнений, и решения, пригодные для численных расчетов, могут быть получены для различных геометрических конфигураций источника света и рассеивающей системы. В некоторых случаях можно использовать двухпоточную модель, в которой боковые компоненты приравниваются нулю. Опубликованы такие расчеты для многократного рассеяния плотной суспензии, имеющей частично отражающие границы. Экспериментально исследовано прохождение света через многократно рассеивающие суспензии частиц латекса и изучено влияние расстояния между частицами на многократное рассеяние 2. Согласно выводам авторов, слой плотного гидрозоля толщиной в несколько миллиметров может применяться для моделирования рассеяния в грубодисперсных атмосферных облаках с размерами порядка нескольких километров. [c.128]
Уравнение (8.7) необходимо использовать дважды один раз — для молекул А и другой раз — для молекул В. Аналитическое решение указанного уравнения пока не найдено, исключая особые случаи, которые, по счастливой случайности, представляют большой практический интерес. В общем с пользой могут быть применены цифровые вычислительные машины, дающие численные решения [12, 72]. Поскольку дифференциальные уравнения нелинейны, решение для модели Данквертса непосредственно не получено и по существу к нему обычно и не стремятся. Приближенные аналитические решения были найдены Ван-Кревеленом и Хофтицером [98] при использовании пленочной модели. [c.355]
Именно это обстоятельство, т. е. необходимость выполиения гранпч1п11х условий, заданных в различных точках экстремали, зачастую и осложняет получение численного решения. Для того чтобы попять, какие при этом возникают трудности, рассмотрим простейший метод численного интегрирования дифференциальных уравнений, используемый для выполнения расчетов на вычислительных мап]пнах. [c.215]
Важный вопрос о соответствии значений констант скоростп реакций эксперпментальным данным вынесен в этой главе в упражнения. Сделано так потому, что, с одной стороны, этот вопрос относится скорее к области чистой, чем прикладной кинетики, и, с другой стороны, его решаюш,ее значение для всей проблемы расчета химических реакторов не вызывает сомнений. Если кинетические зависимости изображаются прямыми линиями, как на логарифмическом графике для реакции первого порядка в упражнении У.2, то оценка точности найденных значений констант скорости реакций может быть получена из отклонения экспериментальных данных от прямой линии, наилучшим образом оиисываюш ей ход процесса. Если дифференциальные уравнения, описывающие систему реакций, должны с самого начала интегрироваться численно, то провести оценку значений констант скорости и их точности значительно труднее. В простейших случаях уравнения можно решать с помощью аналоговой вычислительной машины, где константы скорости представляются переменными сопротивлениями. Эти сопротивления можно изменять вручную, пока не будет достигнуто наилучшее возможное соответствие между расчетными и экспериментальными данными. Если решение проводится на цифровой вычислительной машине, следует использовать метод проб и ошибок. Предположим, [c.116]
В работе Амундсона, Коста и Рудда (см. библиографию на стр. 305) показано, что модель ячеек идеального смешения с N = PJ2 дает хорошее приближение к решению не только простого дифференциального уравнения, но и системы нелинейных уравнений для степени полноты реакции и температуры при Р = Р а. Это позволяет искать решение с помош ью алгебраических, а не дифференциальных уравнений. Полученные значения переменных у выхода реактора Г (1) и (1) можно затем использовать в качестве начальных условий при интегрировании дифференциальных уравнений в обратном направлении (от выхода к входу). Так как в этом направлении интегрирование численно устойчиво, можно найти путем итераций точное решение дифференциальных уравнений. [c.297]
Решение дифференциальных уравнений для двухмерного зернистого слоя представляет значительные трудности. В работе [128] получено численное решение с учетом экзотермической реакции в слое с сильным тепловьш эффектом, однако расчетная разница температур фаз не превышает 2°С при максимальной разности температур слоя и стенки трубы 52 °С.. Определение коэффициентов теплопроводности в зернистом слое на основе двухфазной модели [44] дало результаты на 4% выше, чем для квазигомогенной модели, в интервале Re, = 40 — 500. [c.170]
Математическая модель объекта, характеризуемого не очень сложными дифференциальными уравнениями, часто может быть реализована на аналоговой вычислительной маншне. 2 d \ theta = sin (t + 0.3} / 3 = -cos (t + 0,2) + K`
Мы проинтегрировали по θ слева и по t справа.
Вот график нашего решения, взяв K = 2:
Типичный график решения для примера 2 DE: `theta (t) = root (3) (- 3cos (t + 0.2) +6)`.
Решение дифференциального уравнения
Из приведенных выше примеров мы видим, что решение DE означает нахождение уравнение без производных, удовлетворяющее заданной DE.Решение дифференциального уравнения всегда включает одно или несколько интеграции шага.
Важно уметь идентифицировать тип DE , с которым мы имеем дело, прежде чем пытаться Найди решение.
Определения
DE первого порядка: Содержит только первые производные
DE второго порядка: Содержит вторые производные (и возможно также первые производные)
Степень: наивысшая мощность из наивысшая производная , встречающаяся в DE.7-5лет = 3`
Это DE имеет порядок 2 (самая высокая производная вторая производная ) и степень 4 (степень старшей производной 4.)
Общие и частные решения
Когда мы впервые выполнили интеграцию, мы получили общий раствор (с постоянной K ).
Мы получили частное решение заменой известных значения для x и y .Эти известные условия называется граничными условиями (или начальных условия ).
Это та же концепция, что и при решении дифференциальных уравнений — сначала найдите общее решение, а затем замените заданные числа, чтобы найти частные решения.
Рассмотрим несколько примеров ДУ первого порядка и первой степени.
Пример 4
а. Найдите общее решение для дифференциала уравнение
`dy + 7x dx = 0`
г. 2) / 2 + Bx + C` .
Сейчас
y (0) = 3 дает C = 3.
и
y ‘ (2) = 6 дает A = 6
(на самом деле y » = 6 для любого значения x в этой задаче, поскольку нет члена x )
Наконец,
y ‘ (1) = 4 дает B = -2.
Итак, конкретное решение этого вопроса:
y = 3 x 2 — 2 x + 3
Проверка решения путем дифференцирования и подстановки начальных условий:
y ‘= 6 x — 2
y ‘ (1) = 6 (1) — 2 = 4
y ‘= 6
y » = 0
Наше решение правильное.
Пример 7
После решения дифференциала уравнение,
`(dy) / (dx) ln x-y / x = 0`
(мы увидим, как решить эту DE в следующих раздел Разделение переменных), получаем результат
`y = c ln x`
Приняли ли мы правильное общее решение?
Ответ
Теперь, если `y = c ln x`, то` (dy) / (dx) = c / x`
[См. Производную логарифмической функции, если вы не знаете этого.)
Так
`» LHS «= (dy) / (dx) ln x-y / x`
`= (c / x) ln x — ((c ln x)) / x`
`= 0`
`=» RHS «`
Делаем вывод, что у нас есть правильное решение.
DE второго порядка
Мы включили сюда еще два примера, чтобы дать вам представление о DE второго порядка. Позже в этой главе мы увидим, как решать такие линейные DE второго порядка.
Пример 8
Общее решение второго порядка DE
y ‘+ a 2 y = 0
это
`y = A cos ax + B sin ax`
Пример 9
Общее решение второго порядка DE
y ‘- 3 y ‘ + 2 y = 0
это
y = Ae 2 x + Be x
Если у нас есть следующие граничные условия:
y (0) = 4, y ‘ (0) = 5
, то конкретное решение дает:
y = e 2 x + 3 e x
Теперь мы рассмотрим несколько примеров с использованием DE второго порядка, где нам дается окончательный ответ, и нам нужно проверить, является ли это правильным решением. (m_2x)`
Решение линейных дифференциальных уравнений первого порядка
Возможно, вам сначала захочется прочитать о дифференциальных уравнениях и разделении переменных!
Дифференциальное уравнение — это уравнение с функцией и одной или несколькими производными:
Пример: уравнение с функцией y и ее производная dy dx
Здесь мы рассмотрим решение специального класса дифференциальных уравнений, называемых линейных дифференциальных уравнений первого порядка
Первый орден
Они «Первого Ордена», когда их всего dy dx , а не d 2 y dx 2 или d 3 y dx 3 и т. Д.
Линейный
Дифференциальное уравнение первого порядка является линейным , если его можно сделать так:
dy dx + Р (х) у = Q (х)
Где P (x) и Q (x) — функции от x.
Для ее решения есть специальный метод:
- Мы изобретаем две новые функции от x, называем их u и v и говорим, что y = uv .
- Затем мы решаем, чтобы найти и , а затем найти v , убрать, и все готово!
И мы также используем производную от y = uv (см. Производные правила (правило продукта)):
dy dx = u дв dx + v du dx
ступеньки
Вот пошаговый метод их решения:
Давайте посмотрим на примере, чтобы увидеть:
Пример 1: Решите это:
dy dx — л х = 1
Во-первых, это линейно? Да, так как в форме
dy dx + P (x) y = Q (x)
, где P (x) = — 1 х и Q (x) = 1
Итак, давайте выполним шаги:
Шаг 1: Подставляем y = uv и . dy dx = u дв dx + v du dx
Так это: dy dx — л х = 1
Становится этим: u дв dx + v du dx — УФ х = 1
Шаг 2: Разложите на множители детали v
Фактор v : u дв dx + v ( du dx — u х ) = 1
Шаг 3. Положите член v равным нулю
v член равен нулю: du dx — u х = 0
Итак: du dx знак равно u х
Шаг 4: Решите, используя разделение переменных, чтобы найти u
Отдельные переменные: du u знак равно dx х
Поставьте знак интеграла: ∫ du u = ∫ dx х
Интегрировать: ln (u) = ln (x) + C
Сделайте C = ln (k): ln (u) = ln (x) + ln (k)
Итак: u = kx
Шаг 5: подставьте u обратно в уравнение на шаге 2
(помните, что термин v равен 0, поэтому его можно игнорировать): kx дв dx = 1
Шаг 6: Решите это, чтобы найти v
Отдельные переменные: k dv = dx х
Поставить знак интеграла: ∫ k дв. = ∫ dx х
Интегрировать: kv = ln (x) + C
Сделайте C = ln (c): kv = ln (x) + ln (c)
И так: kv = ln (cx)
И так: v = 1 к ln (сх)
Шаг 7: Подставьте в y = uv , чтобы найти решение исходного уравнения.
y = uv: y = kx 1 к ln (сх)
Упростить: y = x ln (cx)
И он производит это прекрасное семейство кривых:
y = x ln (cx) для различных значений c
Что означают эти кривые? Они являются решением уравнения dy dx — л х = 1
Другими словами:
В любом месте на любой из этих кривых
наклон минус л х равно 1
Давайте проверим несколько точек на c = 0. 6 кривая:
Расчет по графику (до 1 знака после запятой):
Точка | х | y | Уклон ( dy dx ) | dy dx — л х |
---|---|---|---|---|
А | 0,6 | -0,6 | 0 | 0 — −0.6 0,6 = 0 + 1 = 1 |
B | 1,6 | 0 | 1 | 1 — 0 1,6 = 1 — 0 = 1 |
С | 2,5 | 1 | 1,4 | 1,4 — 1 2,5 = 1,4 — 0,4 = 1 |
Почему бы не проверить несколько пунктов самостоятельно? Здесь вы можете построить кривую.
Может, вам поможет еще один пример? Может, посложнее?
Пример 2: Решите это:
dy dx — 3 года х = х
Во-первых, это линейно? Да, так как в форме
dy dx + P (x) y = Q (x)
, где P (x) = — 3 х и Q (x) = x
Итак, давайте выполним шаги:
Шаг 1: Подставляем y = uv и . dy dx = u дв dx + v du dx
Так это: dy dx — 3 года х = х
Становится этим: u дв dx + v du dx — 3uv х = х
Шаг 2: Разложите на множители детали v
Фактор v : u дв dx + v ( du dx — 3u х ) = х
Шаг 3. Положите член v равным нулю
v член = ноль: du dx — 3u х = 0
Итак: du dx знак равно 3u х
Шаг 4: Решите, используя разделение переменных, чтобы найти u
Отдельные переменные: du u = 3 dx х
Поставьте знак интеграла: ∫ du u = 3 ∫ dx х
Интегрировать: ln (u) = 3 ln (x) + C
Сделайте C = −ln (k): ln (u) + ln (k) = 3ln (x)
Тогда: uk = x 3
Итак: u = x 3 к
Шаг 5: подставьте u обратно в уравнение на шаге 2
(помните, что термин v равен 0, поэтому его можно игнорировать) 🙁 x 3 к ) дв dx = х
Шаг 6: Решите это, чтобы найти v
Отдельные переменные: dv = k x −2 dx
Поставьте знак интеграла: ∫ dv = ∫ k x −2 dx
Интегрировать: v = −k x −1 + D
Шаг 7: Подставьте в y = uv , чтобы найти решение исходного уравнения.
у = УФ: у = x 3 к (−k x −1 + D)
Упростить: y = −x 2 + D к х 3
Заменить D / k одной константой c : y = c х 3 — х 2
И он производит это прекрасное семейство кривых:
у = с
x 3 — x 2 для различных значений c
И еще один пример, на этот раз еще на сложнее :
Пример 3: Решите это:
dy dx + 2xy = −2x 3
Во-первых, это линейно? Да, так как в форме
dy dx + P (x) y = Q (x)
, где P (x) = 2x и Q (x) = −2x 3
Итак, давайте выполним шаги:
Шаг 1: Подставляем y = uv и . dy dx = u дв dx + v du dx
Так это: dy dx + 2xy = −2x 3
Становится этим: u дв dx + v du dx + 2xuv = −2x 3
Шаг 2: Разложите на множители детали v
Фактор v : u дв dx + v ( du dx + 2xu ) = −2x 3
Шаг 3. Положите член v равным нулю
v член = ноль: du dx + 2xu = 0
Шаг 4: Решите, используя разделение переменных, чтобы найти u
Отдельные переменные: du u = −2x dx
Поставить знак интеграла: ∫ du u = −2 ∫ x dx
Интегрировать: ln (u) = −x 2 + C
Сделайте C = −ln (k): ln (u) + ln (k) = −x 2
Тогда: uk = e −x 2
Итак: u = e −x 2 к
Шаг 5: подставьте u обратно в уравнение на шаге 2
(помните, что термин v равен 0, поэтому его можно игнорировать) 🙁 e −x 2 к ) дв dx = −2x 3
Шаг 6: Решите это, чтобы найти v
Отдельные переменные: dv = −2k x 3 e x 2 dx
Поставьте знак интеграла: ∫ дв. = ∫ −2k x 3 e x 2 dx
Интегрировать: v = о нет! это трудно!
Посмотрим… мы можем интегрировать по частям … где написано:
∫ RS dx = R ∫ S dx — ∫ R ‘(∫ S dx) dx
(Боковое примечание: здесь мы используем R и S, использование u и v может сбивать с толку, поскольку они уже означают что-то другое.)
Выбор R и S очень важен, это лучший выбор, который мы нашли:
Итак, вперед:
Первое вытягивание k: v = k ∫ −2x 3 e x 2 dx
R = −x 2 и S = 2x e x 2 : v = k ∫ (−x 2 ) (2xe x 2 ) dx
Теперь интегрировать по частям: v = kR ∫ S dx — k ∫ R ‘(∫ S dx) dx
Положим R = −x 2 и S = 2x e x 2
А также R ‘= −2x и ∫ S dx = e x 2
Таким образом, получается: v = −kx 2 ∫ 2x e x 2 dx — k ∫ −2x (e x 2 ) dx
Теперь интегрируйте: v = −kx 2 e x 2 + k e x 2 + D
Упростить: v = ke x 2 (1 − x 2 ) + D
Шаг 7: Подставьте в y = uv , чтобы найти решение исходного уравнения.
у = УФ: у = e −x 2 к (ke x 2 (1 − x 2 ) + D)
Упростить: y = 1 — x 2 + ( D к ) д — x 2
Заменить D / k на одну константу c : y = 1 — x 2 + c е — x 2
И мы получаем красивое семейство кривых:
y = 1 — x 2 +
c
e — x 2 для различных значений c
Дифференциальные уравнения — Линейные уравнения
Показать мобильное уведомление Показать все заметки Скрыть все заметкиПохоже, вы используете устройство с «узкой» шириной экрана ( i.е. вы, вероятно, пользуетесь мобильным телефоном). Из-за особенностей математики на этом сайте лучше всего просматривать в ландшафтном режиме. Если ваше устройство не находится в альбомном режиме, многие уравнения будут отображаться сбоку от вашего устройства (должна быть возможность прокручивать, чтобы увидеть их), а некоторые элементы меню будут обрезаны из-за узкой ширины экрана.
Раздел 2-1: Линейные дифференциальные уравнения
Первым частным случаем дифференциальных уравнений первого порядка, который мы рассмотрим, является линейное дифференциальное уравнение первого порядка.В этом случае, в отличие от большинства случаев первого порядка, которые мы рассмотрим, мы действительно можем вывести формулу для общего решения. Общее решение выводится ниже. Однако мы бы посоветовали вам не запоминать саму формулу. Вместо того, чтобы запоминать формулу, вы должны запомнить и понять процесс, который я собираюсь использовать для получения формулы. На самом деле с большинством проблем легче справиться, используя процесс, а не формулу.
Итак, давайте посмотрим, как решить линейное дифференциальное уравнение первого порядка.Помните, когда мы проходим этот процесс, цель состоит в том, чтобы прийти к решению в форме \ (y = y \ left (t \ right) \). Иногда легко упустить из виду цель, когда мы проходим через этот процесс впервые.
Чтобы решить линейное дифференциальное уравнение первого порядка, мы ДОЛЖНЫ начать с дифференциального уравнения в форме, показанной ниже. Если дифференциальное уравнение не в такой форме, то процесс, который мы собираемся использовать, не будет работать.
\ [\ begin {уравнение} \ frac {{dy}} {{dt}} + p \ left (t \ right) y = g \ left (t \ right) \ label {eq: eq1} \ end {уравнение} \]Где и \ (p (t) \), и \ (g (t) \) — непрерывные функции.Напомним, что быстрое и грязное определение непрерывной функции состоит в том, что функция будет непрерывной при условии, что вы можете рисовать график слева направо, даже не беря в руки карандаш / ручку. Другими словами, функция является непрерывной, если в ней нет дыр и разрывов.
Теперь мы собираемся предположить, что где-то в мире существует некоторая магическая функция, \ (\ mu \ left (t \ right) \), называемая интегрирующим множителем . На этом этапе не беспокойтесь о том, что это за функция и откуда она взялась.Мы выясним, что такое \ (\ mu \ left (t \ right) \), когда у нас будет формула для общего решения.
Итак, теперь, когда мы предположили существование \ (\ mu \ left (t \ right) \), умножаем все в \ (\ eqref {eq: eq1} \) на \ (\ mu \ left (t \ right) \). Это даст.
\ [\ begin {уравнение} \ mu \ left (t \ right) \ frac {{dy}} {{dt}} + \ mu \ left (t \ right) p \ left (t \ right) y = \ mu \ left (t \ right) g \ left (t \ right) \ label {eq: eq2} \ end {уравнение} \]Вот здесь и вступает в игру магия \ (\ mu \ left (t \ right) \).Мы собираемся предположить, что что бы ни было \ (\ mu \ left (t \ right) \), оно будет удовлетворять следующему.
\ [\ begin {уравнение} \ mu \ left (t \ right) p \ left (t \ right) = \ mu ‘\ left (t \ right) \ label {eq: eq3} \ end {уравнение} \]Снова не беспокойтесь о том, как мы можем найти \ (\ mu \ left (t \ right) \), который будет удовлетворять \ (\ eqref {eq: eq3} \). Как мы увидим, при условии непрерывности \ (p (t) \) мы можем его найти. Итак, подставив \ (\ eqref {eq: eq3} \), мы приходим к.
\ [\ begin {уравнение} \ mu \ left (t \ right) \ frac {{dy}} {{dt}} + \ mu ‘\ left (t \ right) y = \ mu \ left (t \ right) g \ left (t \ right) \ label {eq: eq4} \ end {Equation} \]На этом этапе мы должны признать, что левая часть \ (\ eqref {eq: eq4} \) — не что иное, как следующее правило продукта.\ prime} \, dt}} = \ int {{\ mu \ left (t \ right) g \ left (t \ right) \, dt}} \] \ [\ begin {уравнение} \ mu \ left (t \ right) y \ left (t \ right) + c = \ int {{\ mu \ left (t \ right) g \ left (t \ right) \, dt}} \ label {eq: eq6} \ end {уравнение} \]
Обратите внимание, что сюда включена константа интегрирования \ (c \) из левой части интегрирования. Это жизненно важно. Если его не указывать, вы каждый раз будете получать неправильный ответ.
Последний шаг — это некоторая алгебра, которую нужно решить для решения, \ (y (t) \).
\ [\ begin {align *} \ mu \ left (t \ right) y \ left (t \ right) & = \ int {{\ mu \ left (t \ right) g \ left (t \ right) \, dt}} — c \\ y \ left (t \ right) & = \ frac {{\ int {{\ mu \ left (t \ right) g \ left (t \ right) \, dt}} — c} } {{\ mu \ left (t \ right)}} \ end {align *} \]Теперь, с точки зрения обозначений, мы знаем, что постоянная интегрирования \ (c \) — это неизвестная константа, и поэтому, чтобы облегчить нашу жизнь, мы включим знак минус перед ней в константу и вместо этого будем использовать плюс. .Это НЕ повлияет на окончательный ответ решения. Итак, с этим изменением у нас есть.
\ [\ begin {уравнение} y \ left (t \ right) = \ frac {{\ int {{\ mu \ left (t \ right) g \ left (t \ right) \, dt}} + c}} {{\ mu \ left (t \ right)}} \ label {eq: eq7} \ end {уравнение} \]Опять же, изменение знака константы не повлияет на наш ответ. Если вы решите сохранить знак минус, вы получите то же значение \ (c \), что и мы, за исключением того, что у него будет противоположный знак. Подключив \ (c \), мы получим точно такой же ответ.
С константами интеграции в этом разделе очень много шуток, так что вам нужно к этому привыкнуть. Когда мы делаем это, мы всегда будем стараться очень ясно дать понять, что происходит, и попытаться оправдать, почему мы сделали то, что мы сделали.
Итак, теперь, когда у нас есть общее решение для \ (\ eqref {eq: eq1} \), нам нужно вернуться и определить, что это за магическая функция \ (\ mu \ left (t \ right) \). . На самом деле это более простой процесс, чем вы думаете.\ prime} = p \ left (t \ right) \]
Как и в случае с процессом, прежде всего, нам нужно объединить обе стороны, чтобы получить.
\ [\ begin {align *} \ ln \ mu \ left (t \ right) + k & = \ int {{p \ left (t \ right) \, dt}} \\ \ ln \ mu \ left (t \ right) & = \ int {{p \ left (t \ right) \, dt}} + k \ end {align *} \]Вы заметите, что константа интегрирования с левой стороны, \ (k \), была перемещена в правую часть и снова поглощена знаком минус, как мы это делали ранее. {\ int {{p \ left (t \ right) \, dt}}}} \ label {eq: eq10} \ end {формула} \]
Теперь реальность такова, что \ (\ eqref {eq: eq9} \) не так полезен, как может показаться.Часто проще просто выполнить процесс, который привел нас к \ (\ eqref {eq: eq9} \), чем использовать формулу. Мы не будем использовать эту формулу ни в одном из наших примеров. Нам нужно будет регулярно использовать \ (\ eqref {eq: eq10} \), так как эту формулу проще использовать, чем процесс ее получения.
Процесс решения
Процесс решения линейного дифференциального уравнения первого порядка выглядит следующим образом.
- Приведите дифференциальное уравнение в правильную начальную форму, \ (\ eqref {eq: eq1} \).
- Найдите интегрирующий коэффициент \ (\ mu \ left (t \ right) \), используя \ (\ eqref {eq: eq10} \).
- Умножьте все в дифференциальном уравнении на \ (\ mu \ left (t \ right) \) и убедитесь, что левая часть становится правилом произведения \ (\ left ({\ mu \ left (t \ right) y \ left ( t \ right)} \ right) ‘\) и напишите это как таковое.
- Интегрируйте обе стороны, убедитесь, что вы правильно применили постоянную интеграции.
- Найдите решение \ (y (t) \).
Давайте поработаем пару примеров. Начнем с решения дифференциального уравнения, полученного нами в разделе «Поле направления».
Пример 1 Найдите решение следующего дифференциального уравнения. \ [\ frac {{dv}} {{dt}} = 9,8–0,196 об. \] Показать решениеВо-первых, нам нужно получить дифференциальное уравнение в правильной форме.
\ [\ frac {{dv}} {{dt}} + 0,196v = 9.{- 0,196т}} \]Из решения этого примера мы теперь можем понять, почему постоянная интеграции так важна в этом процессе. Без него в этом случае мы получили бы одно постоянное решение \ (v (t) = 50 \). Используя постоянную интегрирования, мы получаем бесконечно много решений, по одному для каждого значения \ (c \).
Вернувшись в раздел поля направления, где мы впервые вывели дифференциальное уравнение, используемое в последнем примере, мы использовали поле направления, чтобы помочь нам наметить некоторые решения.Посмотрим, правильно ли мы их поняли. Чтобы набросать некоторые решения, все, что нам нужно сделать, это выбрать разные значения \ (c \), чтобы получить решение. Некоторые из них показаны на графике ниже.
Итак, похоже, мы неплохо сделали набросок графиков в секции поля направлений.
Теперь вспомните из раздела «Определения», что начальные условия позволят нам сосредоточиться на конкретном решении. В решениях дифференциальных уравнений первого порядка (не только линейных, как мы увидим) будет одна неизвестная константа, поэтому нам понадобится ровно одно начальное условие, чтобы найти значение этой константы и, следовательно, найти решение, к которому мы пришли.Начальное условие для дифференциальных уравнений первого порядка будет иметь вид
\ [y \ left ({{t_0}} \ right) = {y_0} \]Напомним также, что дифференциальное уравнение с достаточным количеством начальных условий называется задачей начального значения (IVP).
Пример 2 Решите следующую IVP. \ [\ frac {{dv}} {{dt}} = 9,8–0,196v \ hspace {0,25 дюйма} v \ left (0 \ right) = 48 \] Показать решениеЧтобы найти решение IVP, мы должны сначала найти общее решение дифференциального уравнения, а затем использовать начальное условие для определения точного решения, которое мы ищем.{- 0,196т}} \]
Теперь, чтобы найти решение, которое мы ищем, нам нужно определить значение \ (c \), которое даст нам решение, которое мы ищем. Для этого мы просто подставляем начальное условие, которое даст нам уравнение, которое мы можем решить для \ (c \). {- 0.{\ ln \, \, \ sec \ left (x \ right)}} = \ sec \ left (x \ right) \]
А интеграл сделать можно? Если нет, перепишите касательную обратно в синусы и косинусы, а затем используйте простую замену. Обратите внимание, что мы могли опустить столбцы абсолютного значения на секансе из-за ограничений на \ (x \). Фактически, это причина ограничений на \ (x \). Отметим также, что есть две формы ответа на этот интеграл. Они эквивалентны, как показано ниже. Что вы используете — это действительно вопрос предпочтений.{\ ln f \ left (x \ right)}} = f \ left (x \ right) \ label {eq: eq11} \ end {уравнение} \]
Это важный факт, о котором вы всегда должны помнить при возникновении подобных проблем. Мы захотим максимально упростить интегрирующий коэффициент во всех случаях, и этот факт поможет в этом упрощении.
Вернемся к примеру. Умножьте интегрирующий коэффициент на дифференциальное уравнение и убедитесь, что левая часть соответствует правилу произведения. Также обратите внимание, что мы умножаем интегрирующий коэффициент на переписанное дифференциальное уравнение, а НЕ на исходное дифференциальное уравнение.2} \ left (x \ right) \, dx}} \\ \ sec \ left (x \ right) y \ left (x \ right) & = — \ frac {1} {2} \ cos \ left ({ 2x} \ right) — \ tan \ left (x \ right) + c \ end {align *} \]
Обратите внимание на использование тригонометрической формулы \ (\ sin \ left ({2 \ theta} \ right) = 2 \ sin \ theta \ cos \ theta \), которая упростила интеграл. Затем найдите решение.
\ [\ begin {align *} y \ left (x \ right) & = — \ frac {1} {2} \ cos \ left (x \ right) \ cos \ left ({2x} \ right) — \ cos \ left (x \ right) \ tan \ left (x \ right) + c \ cos \ left (x \ right) \\ & = — \ frac {1} {2} \ cos \ left (x \ right) \ cos \ left ({2x} \ right) — \ sin \ left (x \ right) + c \ cos \ left (x \ right) \ end {align *} \]Наконец, примените начальное условие, чтобы найти значение \ (c \). {\ frac {t } {2}}} \]
Теперь, когда у нас есть решение, давайте посмотрим на долгосрочное поведение ( i.е. \ (t \ to \ infty \)) решения. Первые два члена решения останутся конечными при всех значениях \ (t \). Это последний член, который будет определять поведение решения. Экспонента всегда будет стремиться к бесконечности как \ (t \ to \ infty \), однако в зависимости от знака коэффициента \ (c \) (да, мы уже нашли его, но для простоты этого обсуждения мы продолжим называть это \ (c \)). В следующей таблице показано долгосрочное поведение решения для всех значений \ (c \).
Диапазон \ (c \) | Поведение решения как \ (t \ to \ infty \) |
---|---|
\ (с \) <0 | \ (y \ left (t \ right) \ to — \ infty \) |
\ (с \) = 0 | \ (y \ left (t \ right) \) остается конечным |
\ (c \)> 0 | \ (y \ left (t \ right) \ to \ infty \) |
Это поведение также можно увидеть на следующем графике некоторых решений.
Теперь, поскольку мы знаем, как \ (c \) относится к \ (y_ {0} \), мы можем связать поведение решения с \ (y_ {0} \). В следующей таблице показано поведение решения в терминах \ (y_ {0} \) вместо \ (c \).
Диапазон \ (y_ {0} \) | Поведение решения как \ (t \ to \ infty \) |
---|---|
\ ({y_0} <- \ frac {{24}} {{37}} \) | \ (y \ left (t \ right) \ to — \ infty \) |
\ ({y_0} = — \ frac {{24}} {{37}} \) | \ (y \ left (t \ right) \) остается конечным |
\ ({y_0}> — \ frac {{24}} {{37}} \) | \ (y \ left (t \ right) \ to \ infty \) |
Обратите внимание, что для \ ({y_0} = — \ frac {{24}} {{37}} \) решение останется конечным.Так бывает не всегда.
Исследование долгосрочного поведения решений иногда бывает важнее, чем само решение. Предположим, что указанный выше раствор дал температуру в металлическом бруске. В этом случае нам нужно решение (я), которое останется конечным в долгосрочной перспективе. Благодаря этому исследованию у нас теперь будет значение начального условия, которое даст нам это решение, и, что более важно, значения начального условия, которых нам нужно было бы избежать, чтобы мы не расплавили стержень.
4.1 Основы дифференциальных уравнений — Исчисление Том 2
Цели обучения
- 4.1.1 Определите порядок дифференциального уравнения.
- 4.1.2 Объясните, что подразумевается под решением дифференциального уравнения.
- 4.1.3. Различают общее решение и частное решение дифференциального уравнения.
- 4.1.4 Определите проблему начального значения.
- 4.1.5 Определите, является ли данная функция решением дифференциального уравнения или задачей с начальным значением.
Исчисление — это математика изменений, а скорость изменения выражается производными. Таким образом, один из наиболее распространенных способов использования исчисления — это составить уравнение, содержащее неизвестную функцию y = f (x) y = f (x) и ее производную, известное как дифференциальное уравнение . Решение таких уравнений часто дает информацию о том, как меняются количества, и часто дает представление о том, как и почему происходят изменения.
Методы решения дифференциальных уравнений могут принимать различные формы, включая прямое решение, использование графиков или компьютерные вычисления.Мы представляем основные идеи в этой главе и более подробно опишем их позже в курсе. В этом разделе мы изучаем, что такое дифференциальные уравнения, как проверять их решения, некоторые методы, которые используются для их решения, а также некоторые примеры общих и полезных уравнений.
Общие дифференциальные уравнения
Рассмотрим уравнение y ′ = 3×2, y ′ = 3×2, которое является примером дифференциального уравнения, поскольку оно включает производную. Между переменными xx и y существует связь: yy: y — неизвестная функция от x.Икс. Кроме того, левая часть уравнения — это производная от y.y. Следовательно, мы можем интерпретировать это уравнение следующим образом: начните с некоторой функции y = f (x) y = f (x) и возьмите ее производную. Ответ должен быть равен 3×2,3×2. У какой функции есть производная, равная 3×2? 3×2? Одна из таких функций — y = x3, y = x3, поэтому эта функция считается решением дифференциального уравнения.
Определение
Дифференциальное уравнение — это уравнение, включающее неизвестную функцию y = f (x) y = f (x) и одну или несколько ее производных.Решение дифференциального уравнения — это функция y = f (x) y = f (x), которая удовлетворяет дифференциальному уравнению, когда ff и его производные подставляются в уравнение.
Некоторые примеры дифференциальных уравнений и их решений приведены в таблице 4.1.
Уравнение | Решение |
---|---|
y ′ = 2xy ′ = 2x | у = х2у = х2 |
y ′ + 3y = 6x + 11y ′ + 3y = 6x + 11 | y = e − 3x + 2x + 3y = e − 3x + 2x + 3 |
y ′ ′ — 3y ′ + 2y = 24e − 2xy ′ ′ — 3y ′ + 2y = 24e − 2x | y = 3ex − 4e2x + 2e − 2xy = 3ex − 4e2x + 2e − 2x |
Таблица 4.1 Примеры дифференциальных уравнений и их решений
Обратите внимание, что решение дифференциального уравнения не обязательно является уникальным, прежде всего потому, что производная константы равна нулю. Например, y = x2 + 4y = x2 + 4 также является решением первого дифференциального уравнения в таблице 4.1. Мы вернемся к этой идее чуть позже в этом разделе. А пока давайте сосредоточимся на том, что означает, что функция является решением дифференциального уравнения.
Пример 4.1
Проверка решений дифференциальных уравнений
Убедитесь, что функция y = e − 3x + 2x + 3y = e − 3x + 2x + 3 является решением дифференциального уравнения y ′ + 3y = 6x + 11.у ‘+ 3у = 6х + 11.
Решение
Чтобы проверить решение, мы сначала вычисляем y′y ′, используя цепное правило для производных. Это дает y ′ = — 3e − 3x + 2.y ′ = — 3e − 3x + 2. Затем подставляем yy и y′y ′ в левую часть дифференциального уравнения:
(−3e − 2x + 2) +3 (e − 2x + 2x + 3). (- 3e − 2x + 2) +3 (e − 2x + 2x + 3).Результирующее выражение можно упростить, предварительно распространив его, исключив круглые скобки, что даст
−3e − 2x + 2 + 3e − 2x + 6x + 9. − 3e − 2x + 2 + 3e − 2x + 6x + 9.Объединение одинаковых членов приводит к выражению 6x + 11,6x + 11, которое равно правой части дифференциального уравнения.Этот результат подтверждает, что y = e − 3x + 2x + 3y = e − 3x + 2x + 3 является решением дифференциального уравнения.
КПП 4.1
Убедитесь, что y = 2e3x − 2x − 2y = 2e3x − 2x − 2 является решением дифференциального уравнения y′ − 3y = 6x + 4.y′ − 3y = 6x + 4.
Удобно определять характеристики дифференциальных уравнений, которые упрощают их обсуждение и категоризацию. Самая основная характеристика дифференциального уравнения — это его порядок.
Определение
Порядок дифференциального уравнения — это наивысший порядок любой производной неизвестной функции, которая появляется в уравнении.
Пример 4.2
Определение порядка дифференциального уравнения
Каков порядок каждого из следующих дифференциальных уравнений?
- y′ − 4y = x2−3x + 4y′ − 4y = x2−3x + 4
- x2y ‴ −3xy ″ + xy′ − 3y = sinxx2y ‴ −3xy ″ + xy′ − 3y = sinx
- 4xy (4) −6x2y ″ + 12x4y = x3−3×2 + 4x − 124xy (4) −6x2y ″ + 12x4y = x3−3×2 + 4x − 12
Решение
- Старшая производная в уравнении — y ′, y ′, поэтому порядок равен 1.1.
- Старшая производная в уравнении — y ‴, y ‴, поэтому порядок равен 3.3.
- Старшая производная в уравнении — y (4), y (4), поэтому порядок равен 4,4.
Контрольно-пропускной пункт 4.2
Каков порядок следующего дифференциального уравнения?
(x4−3x) y (5) — (3×2 + 1) y ′ + 3y = sinxcosx (x4−3x) y (5) — (3×2 + 1) y ′ + 3y = sinxcosxОбщие и частные решения
Мы уже отмечали, что дифференциальное уравнение y ′ = 2xy ′ = 2x имеет по крайней мере два решения: y = x2y = x2 и y = x2 + 4.y = x2 + 4. Единственная разница между этими двумя решениями — это последний член, который является константой.Что, если последний член — другая константа? Будет ли это выражение по-прежнему решением дифференциального уравнения? Фактически, любая функция вида y = x2 + C, y = x2 + C, где CC представляет любую константу, также является решением. Причина в том, что производная x2 + Cx2 + C равна 2x, 2x, независимо от значения C.C. Можно показать, что любое решение этого дифференциального уравнения должно иметь вид y = x2 + C.y = x2 + C. Это пример общего решения дифференциального уравнения. График некоторых из этих решений представлен на рисунке 4.2. ( Примечание : на этом графике мы использовали четные целые значения для CC в диапазоне от -4-4 до 4,4. Фактически, нет никаких ограничений на значение C; C; оно может быть целым или нет.)
Рис. 4.2 Семейство решений дифференциального уравнения y ′ = 2x.y ′ = 2x.В этом примере мы можем выбрать любое решение, какое захотим; например, y = x2−3y = x2−3 является членом семейства решений этого дифференциального уравнения. Это называется частным решением дифференциального уравнения. Конкретное решение часто можно однозначно идентифицировать, если нам предоставить дополнительную информацию о проблеме.
Пример 4.3
Поиск конкретного решения
Найти частное решение дифференциального уравнения y ′ = 2xy ′ = 2x, проходящее через точку (2,7). (2,7).
Решение
Любая функция вида y = x2 + Cy = x2 + C является решением этого дифференциального уравнения. Чтобы определить значение C, C, мы подставляем значения x = 2x = 2 и y = 7y = 7 в это уравнение и решаем относительно C: C:
y = x2 + C7 = 22 + C = 4 + CC = 3. y = x2 + C7 = 22 + C = 4 + CC = 3.Следовательно, частное решение, проходящее через точку (2,7) (2,7), имеет вид y = x2 + 3.у = х2 + 3.
КПП 4.3
Найдите частное решение дифференциального уравнения
, проходящее через точку (1,7), (1,7), при условии, что y = 2×2 + 3x + Cy = 2×2 + 3x + C, является общим решением дифференциального уравнения.
Задачи начального значения
Обычно данное дифференциальное уравнение имеет бесконечное число решений, поэтому естественно спросить, какое из них мы хотим использовать. Чтобы выбрать одно решение, необходима дополнительная информация. Некоторая конкретная информация, которая может быть полезна, — это начальное значение, которое представляет собой упорядоченную пару, которая используется для поиска конкретного решения.
Дифференциальное уравнение вместе с одним или несколькими начальными значениями называется задачей начальных значений. Общее правило состоит в том, что количество начальных значений, необходимых для задачи с начальными значениями, равно порядку дифференциального уравнения. Например, если у нас есть дифференциальное уравнение y ′ = 2x, y ′ = 2x, тогда y (3) = 7y (3) = 7 является начальным значением, и вместе эти уравнения образуют задачу с начальным значением. Дифференциальное уравнение y ″ −3y ′ + 2y = 4exy ″ −3y ′ + 2y = 4ex второго порядка, поэтому нам нужны два начальных значения.В задачах с начальным значением порядка больше единицы одно и то же значение следует использовать для независимой переменной. Примером начальных значений для этого уравнения второго порядка может быть y (0) = 2y (0) = 2 и y ′ (0) = — 1.y ′ (0) = — 1. Эти два начальных значения вместе с дифференциальным уравнением образуют начальную задачу. Эти проблемы названы так потому, что часто независимой переменной в неизвестной функции является t, t, которая представляет время. Таким образом, значение t = 0t = 0 представляет начало проблемы.
Пример 4.4
Проверка решения проблемы с начальным значением
Проверить, что функция y = 2e − 2t + ety = 2e − 2t + et является решением задачи начального значения
y ′ + 2y = 3et, y (0) = 3. y ′ + 2y = 3et, y (0) = 3.Решение
Чтобы функция удовлетворяла задаче с начальным значением, она должна удовлетворять как дифференциальному уравнению, так и начальному условию. Чтобы показать, что yy удовлетворяет дифференциальному уравнению, мы начнем с вычисления y′.y ′. Это дает y ′ = — 4e − 2t + et.y ′ = — 4e − 2t + et. Затем подставляем yy и y′y ′ в левую часть дифференциального уравнения и упрощаем:
y ′ + 2y = (- 4e − 2t + et) +2 (2e − 2t + et) = — 4e − 2t + et + 4e − 2t + 2et = 3et.y ′ + 2y = (- 4e − 2t + et ) +2 (2e − 2t + et) = — 4e − 2t + et + 4e − 2t + 2et = 3et.Это равно правой части дифференциального уравнения, поэтому y = 2e − 2t + ety = 2e − 2t + et решает дифференциальное уравнение. Затем мы вычисляем y (0): y (0):
y (0) = 2e − 2 (0) + e0 = 2 + 1 = 3. y (0) = 2e − 2 (0) + e0 = 2 + 1 = 3.Этот результат подтверждает исходное значение. Следовательно, данная функция удовлетворяет начальной задаче.
КПП 4.4
Убедитесь, что y = 3e2t + 4sinty = 3e2t + 4sint является решением задачи начального значения
y′ − 2y = 4cost − 8sint, y (0) = 3. y′ − 2y = 4cost − 8sint, y (0) = 3.В примере 4.4 задача начального значения состояла из двух частей. Первой частью было дифференциальное уравнение y ′ + 2y = 3ex, y ′ + 2y = 3ex, а второй частью было начальное значение y (0) = 3.y (0) = 3. Эти два уравнения вместе образуют начальную задачу.
То же и в целом. Задача начальной стоимости будет состоять из двух частей: дифференциального уравнения и начального условия.Дифференциальное уравнение имеет семейство решений, а начальное условие определяет значение C.C. Семейство решений дифференциального уравнения из примера 4.4 имеет вид y = 2e − 2t + Cet.y = 2e − 2t + Cet. Это семейство решений показано на рис. 4.3, с помеченным частным решением y = 2e − 2t + ety = 2e − 2t + et.
Рис. 4.3. Семейство решений дифференциального уравнения y ′ + 2y = 3et.y ′ + 2y = 3et. Отмечено частное решение y = 2e − 2t + ety = 2e − 2t + et.Пример 4.5
Решение задачи с начальным значением
Решите следующую задачу с начальным значением:
y ′ = 3ex + x2−4, y (0) = 5.y ′ = 3ex + x2−4, y (0) = 5.Решение
Первым шагом в решении этой начальной задачи является поиск общего семейства решений. Для этого находим первообразную обеих частей дифференциального уравнения
∫y′dx = ∫ (3ex + x2−4) dx, ∫y′dx = ∫ (3ex + x2−4) dx,а именно
y + C1 = 3ex + 13×3−4x + C2.y + C1 = 3ex + 13×3−4x + C2.4,1
Мы можем интегрировать обе стороны, потому что термин y появляется сам по себе. Обратите внимание на две константы интегрирования: C1C1 и C2.C2. Решение уравнения 4.1 для yy дает
y = 3ex + 13×3−4x + C2 − C1.y = 3ex + 13×3−4x + C2 − C1.Поскольку C1C1 и C2C2 являются константами, C2-C1C2-C1 также является константой. Следовательно, мы можем определить C = C2 − C1, C = C2 − C1, что приводит к уравнению
y = 3ex + 13×3−4x + C.y = 3ex + 13×3−4x + C.Далее мы определяем значение C.C. Для этого мы подставляем x = 0x = 0 и y = 5y = 5 в уравнение 4.1 и решаем относительно C: C:
5 = 3e0 + 1303−4 (0) + C5 = 3 + CC = 2,5 = 3e0 + 1303−4 (0) + C5 = 3 + CC = 2.Теперь мы подставляем значение C = 2C = 2 в уравнение 4.1. Решение начальной задачи: y = 3ex + 13×3−4x + 2.y = 3ex + 13×3−4x + 2.
Анализ
Разница между общим решением и частным решением состоит в том, что общее решение включает в себя семейство функций, определенных явно или неявно, от независимой переменной. Начальное значение или значения определяют, какое конкретное решение в семействе решений удовлетворяет желаемым условиям.
КПП 4.5
Решите задачу начального значения
y ′ = x2−4x + 3−6ex, y (0) = 8.y ′ = x2−4x + 3−6ex, y (0) = 8.В физических и инженерных приложениях мы часто рассматриваем силы, действующие на объект, и используем эту информацию, чтобы понять результирующее движение, которое может произойти. Например, если мы начнем с объекта на поверхности Земли, первичной силой, действующей на этот объект, будет гравитация. Физики и инженеры могут использовать эту информацию вместе со вторым законом движения Ньютона (в форме уравнения F = ma, F = ma, где FF представляет силу, mm представляет массу, а aa представляет ускорение), чтобы вывести уравнение, которое может быть решено. .
Рис. 4.4. При падении бейсбольного мяча в воздухе единственной силой, действующей на него, является сила тяжести (без учета сопротивления воздуха).
На рисунке 4.4 мы предполагаем, что единственная сила, действующая на бейсбольный мяч, — это сила тяжести. Это предположение игнорирует сопротивление воздуха. (Сила, создаваемая сопротивлением воздуха, рассматривается в более позднем обсуждении.) Ускорение свободного падения на поверхности Земли, g, g, составляет примерно 9,8 м / с и 2,9,8 м / с2. Мы вводим систему отсчета, в которой поверхность Земли находится на высоте 0 метров.Пусть v (t) v (t) представляет скорость объекта в метрах в секунду. Если v (t)> 0, v (t)> 0, мяч поднимается, а если v (t) <0, v (t) <0, мяч падает (рисунок 4.5).
Рис. 4.5. Возможные скорости восходящего / падающего бейсбольного мяча.
Наша цель — найти скорость v (t) v (t) в любой момент времени t.t. Для этого мы создаем задачу с начальным значением. Предположим, масса мяча равна m, m, где мм измеряется в килограммах. Мы используем второй закон Ньютона, который гласит, что сила, действующая на объект, равна его массе, умноженной на его ускорение (F = ma).(F = ma). Ускорение — это производная скорости, поэтому a (t) = v ′ (t) .a (t) = v ′ (t). Следовательно, сила, действующая на бейсбольный мяч, определяется выражением F = mv ′ (t) .F = mv ′ (t). Однако эта сила должна быть равна силе тяжести, действующей на объект, которая (опять же с использованием второго закона Ньютона) задается формулами Fg = −mg, Fg = −mg, поскольку эта сила действует в направлении вниз. Следовательно, мы получаем уравнение F = Fg, F = Fg, которое превращается в mv ′ (t) = — mg.mv ′ (t) = — mg. Разделив обе части уравнения на миллиметры, получим уравнение
v ′ (t) = — g.v ′ (t) = — g.Обратите внимание, что это дифференциальное уравнение остается неизменным независимо от массы объекта.
Теперь нам нужно начальное значение. Поскольку мы решаем для скорости, в контексте задачи имеет смысл предположить, что мы знаем начальную скорость или скорость в момент времени t = 0. t = 0. Это обозначается v (0) = v0.v (0) = v0.
Пример 4.6
Скорость движущегося бейсбольного мяча
Бейсбольный мяч подбрасывается вверх с высоты 33 метра над поверхностью Земли с начальной скоростью 10 м / с, 10 м / с, и единственная сила, действующая на него, — это сила тяжести.Мяч имеет массу 0,15 кг 0,15 кг на поверхности Земли.
- Найдите скорость v (t) v (t) бейсбольного мяча в момент времени t.t.
- Какова его скорость через 22 секунды?
Решение
- Из предыдущего обсуждения дифференциальное уравнение, которое применяется в этой ситуации, равно
v ′ (t) = — g, v ′ (t) = — g,
где g = 9,8 м / с2. g = 9,8 м / с2. Начальное условие: v (0) = v0, v (0) = v0, где v0 = 10 м / с. V0 = 10 м / с. Следовательно, начальная задача v ′ (t) = — 9.8 м / с2, v (0) = 10 м / с. V ′ (t) = — 9,8 м / с2, v (0) = 10 м / с.
Первым шагом в решении этой начальной задачи является взятие первообразной обеих частей дифференциального уравнения. Это дает
V ′ (t) dt = ∫ − 9.8dtv (t) = — 9.8t + C.∫v ′ (t) dt = ∫ − 9.8dtv (t) = — 9.8t + C.
Следующий шаг — найти C.C. Для этого подставляем t = 0t = 0 и v (0) = 10: v (0) = 10:
v (t) = — 9,8t + Cv (0) = — 9,8 (0) + C10 = C. v (t) = — 9,8t + Cv (0) = — 9,8 (0) + C10 = C.
Следовательно, C = 10C = 10 и функция скорости определяется как v (t) = — 9,8t + 10.v (t) = — 9.8т + 10. - Чтобы найти скорость через 22 секунды, подставьте t = 2t = 2 в v (t) .v (t).
v (t) = — 9,8t + 10v (2) = — 9,8 (2) + 10v (2) = — 9,6.v (t) = — 9,8t + 10v (2) = — 9,8 (2) + 10v ( 2) = — 9,6.
Единицы измерения скорости — метры в секунду. Поскольку ответ отрицательный, объект падает со скоростью 9,6 м / с. 9,6 м / с.
КПП 4.6
Предположим, что скала падает с высоты 100-100 метров, и единственная сила, действующая на нее, — это сила тяжести. Найдите уравнение для скорости v (t) v (t) как функции времени, измеряемой в метрах в секунду.
Естественный вопрос, который следует задать после решения проблемы такого типа, — насколько высоко объект будет находиться над поверхностью Земли в данный момент времени. Пусть s (t) s (t) обозначает высоту объекта над поверхностью Земли, измеряемую в метрах. Поскольку скорость является производной от положения (в данном случае от высоты), это предположение дает уравнение s ′ (t) = v (t) .s ′ (t) = v (t). Необходимо начальное значение; в этом случае подходит начальная высота объекта. Пусть начальная высота задается уравнением s (0) = s0.s (0) = s0. Вместе эти предположения дают задачу начального значения
. s ′ (t) = v (t), s (0) = s0.s ′ (t) = v (t), s (0) = s0.Если функция скорости известна, то можно также найти функцию положения.
Пример 4.7
Высота движущегося бейсбольного мяча
Бейсбольный мяч подбрасывается вверх с высоты 33 метра над поверхностью Земли с начальной скоростью 10 м / с, 10 м / с, и единственная сила, действующая на него, — это сила тяжести. У поверхности Земли шар имеет массу 0,150,15 кг.
- Найдите положение s (t) s (t) бейсбольного мяча в момент времени t.t.
- Какова его высота через 22 секунды?
Решение
- Мы уже знаем, что функция скорости для этой задачи равна v (t) = — 9,8t + 10.v (t) = — 9,8t + 10. Начальная высота бейсбольного мяча составляет 33 метра, поэтому s0 = 3.s0 = 3. Следовательно, задача с начальным значением для этого примера —
Чтобы решить задачу с начальным значением, мы сначала находим первообразные:
∫s ′ (t) dt = ∫ − 9.8t + 10dts (t) = — 4.9t2 + 10t + C.∫s ′ (t) dt = ∫ − 9.8t + 10dts (t) = — 4.9t2 + 10t + C.
Затем мы подставляем t = 0t = 0 и решаем относительно C: C:
s (t) = — 4.9t2 + 10t + Cs (0) = — 4.9 (0) 2 + 10 (0) + C3 = Cs (t) = — 4.9t2 + 10t + Cs (0) = — 4.9 (0 ) 2 + 10 (0) + С3 = С.
Следовательно, функция положения равна s (t) = — 4.9t2 + 10t + 3.s (t) = — 4.9t2 + 10t + 3. - Высота бейсбольного мяча после 2s2s определяется выражением s (2): s (2):
s (2) = — 4,9 (2) 2 + 10 (2) + 3 = −4,9 (4) + 23 = 3,4. s (2) = — 4,9 (2) 2 + 10 (2) + 3 = −4,9 (4) + 23 = 3,4.
Следовательно, через 22 секунды бейсбольный мяч оказывается на высоте 3,43,4 метра над поверхностью Земли.Стоит отметить, что масса мяча полностью уменьшилась в процессе решения задачи.
Раздел 4.1 Упражнения
Определите порядок следующих дифференциальных уравнений.
2.(y ′) 2 = y ′ + 2y (y ′) 2 = y ′ + 2y
3.y ‴ + y ″ y ′ = 3x2y ‴ + y ″ y ′ = 3×2
4.y ′ = y ″ + 3t2y ′ = y ″ + 3t2
6.dydx + d2ydx2 = 3x4dydx + d2ydx2 = 3×4
7.(dydt) 2 + 8dydt + 3y = 4t (dydt) 2 + 8dydt + 3y = 4t
Убедитесь, что следующие функции являются решениями данного дифференциального уравнения.
8.y = x33y = x33 решает y ′ = x2y ′ = x2
9.y = 2e − x + x − 1y = 2e − x + x − 1 решает y ′ = x − yy ′ = x − y
10.y = e3x − ex2y = e3x − ex2 решает y ′ = 3y + exy ′ = 3y + ex
11.y = 11 − xy = 11 − x решает y ′ = y2y ′ = y2
12.y = ex2 / 2y = ex2 / 2 решает y ′ = xyy ′ = xy
13.y = 4 + lnxy = 4 + lnx решает xy ′ = 1xy ′ = 1
14.y = 3 − x + xlnxy = 3 − x + xlnx решает y ′ = lnxy ′ = lnx
15.y = 2ex − x − 1y = 2ex − x − 1 решает y ′ = y + xy ′ = y + x
16.y = ex + sinx2 − cosx2y = ex + sinx2 − cosx2 решает y ′ = cosx + yy ′ = cosx + y
17.y = πe − cosxy = πe − cosx решает y ′ = ysinxy ′ = ysinx
Проверьте следующие общие решения и найдите частное решение.
18.Найдите частное решение дифференциального уравнения y ′ = 4x2y ′ = 4×2, которое проходит через (−3, −30), (- 3, −30), учитывая, что y = C + 4x33y = C + 4×33 является общим решением .
19.Найдите частное решение дифференциального уравнения y ′ = 3x3y ′ = 3×3, которое проходит через (1,4.75), (1,4.75), учитывая, что y = C + 3x44y = C + 3×44 является общим решением.
20.Найдите частное решение дифференциального уравнения y ′ = 3x2yy ′ = 3x2y, которое проходит через (0,12), (0,12), учитывая, что y = Cex3y = Cex3 является общим решением.
21.Найдите частное решение дифференциального уравнения y ′ = 2xyy ′ = 2xy, которое проходит через (0,12), (0,12), учитывая, что y = Cex2y = Cex2 является общим решением.
22.Найдите частное решение дифференциального уравнения y ′ = (2xy) 2y ′ = (2xy) 2, которое проходит через (1, −12), (1, −12), учитывая, что y = −3C + 4x3y = −3C + 4×3 — общее решение.
23.Найдите частное решение дифференциального уравнения y′x2 = yy′x2 = y, которое проходит через (1,2e), (1,2e), учитывая, что y = Ce − 1 / xy = Ce − 1 / x является общее решение.
24.Найдите частное решение дифференциального уравнения 8dxdt = −2cos (2t) −cos (4t) 8dxdt = −2cos (2t) −cos (4t), которое проходит через (π, π), (π, π), при условии, что x = C − 18sin (2t) −132sin (4t) x = C − 18sin (2t) −132sin (4t) — общее решение.
25.Найдите частное решение дифференциального уравнения dudt = tanududt = tanu, которое проходит через (1, π2), (1, π2), учитывая, что u = sin − 1 (eC + t) u = sin − 1 (eC + t ) — общее решение.
26.Найдите частное решение дифференциального уравнения dydt = e (t + y) dydt = e (t + y), которое проходит через (1,0), (1,0), учитывая, что y = −ln (C − et ) y = −ln (C − et) — общее решение.
27.Найдите частное решение дифференциального уравнения y ′ (1 − x2) = 1 + yy ′ (1 − x2) = 1 + y, которое проходит через (0, −2), (0, −2), учитывая, что y = Cx + 11 − x − 1y = Cx + 11 − x − 1 — общее решение.
Найдите общее решение дифференциального уравнения для следующих задач.
29.y ′ = lnx + tanxy ′ = lnx + tanx
30.y ′ = sinxecosxy ′ = sinxecosx
32.y ′ = sin − 1 (2x) y ′ = sin − 1 (2x)
34.x ′ = cotht + lnt + 3t2x ′ = cotht + lnt + 3t2
Решите следующие задачи с начальным значением, начиная с y (0) = 1y (0) = 1 и y (0) = — 1.y (0) = — 1. Нарисуйте оба решения на одном графике.
Решите следующие задачи с начальным значением, начиная с y0 = 10.y0 = 10. В какое время yy увеличится до 100100 или упадет до 1? 1?
Напомним, что семейство решений включает решения дифференциального уравнения, которые отличаются на константу.Для следующих задач используйте свой калькулятор, чтобы построить график семейства решений данного дифференциального уравнения. Используйте начальные условия от y (t = 0) = — 10y (t = 0) = — 10 до y (t = 0) = 10y (t = 0) = 10 с увеличением на 2,2. Есть ли критическая точка, в которой поведение решения начинает меняться?
51.[T] y ′ = x + yy ′ = x + y ( Подсказка: y = Cex − x − 1y = Cex − x − 1 — общее решение)
52.[T] y ′ = xlnx + sinxy ′ = xlnx + sinx
53.Найдите общее решение для описания скорости шара массой 1 фунт1 фунт, который подбрасывается вверх со скоростью aa ft / sec.
54.В предыдущей задаче, если начальная скорость мяча, брошенного в воздух, равна a = 25a = 25 фут / с, запишите частное решение для скорости мяча. Решите, чтобы найти время, когда мяч упадет на землю.
55.Вы подбрасываете два объекта разной массы m1m1 и m2m2 вверх в воздух с одинаковой начальной скоростью aa ft / s. Какая разница в их скорости через 11 секунд?
56.[T] Вы бросаете шар массой 11 кг вверх со скоростью a = 25a = 25 м / с на Марс, где сила тяжести g = −3.711g = −3,711 м / с 2 . Воспользуйтесь калькулятором, чтобы приблизительно определить, насколько дольше мяч находится в воздухе на Марсе, чем на Земле, где g = -9,8 м / с2g = -9,8 м / с2.
57.[T] Для предыдущей задачи воспользуйтесь калькулятором, чтобы приблизительно определить, насколько выше поднялся шар на Марсе, где g = -9,8 м / с2g = -9,8 м / с2.
58.[T] Автомобиль на автостраде ускоряется в соответствии с a = 15cos (πt), a = 15cos (πt), где tt измеряется в часах. Задайте и решите дифференциальное уравнение, чтобы определить скорость автомобиля, если его начальная скорость составляет 5151 миль в час.Какова скорость водителя после 4040 минут езды?
59.[T] Для автомобиля в предыдущей задаче найдите выражение для расстояния, которое автомобиль проехал за время t, t, предполагая, что начальное расстояние равно 0,0. Сколько времени нужно автомобилю, чтобы проехать 100100 миль? Округлите ответ до часов и минут.
60.[T] Для предыдущей задачи найдите общее расстояние, пройденное за первый час.
61.Подставьте y = Be3ty = Be3t в y′ − y = 8e3ty′ − y = 8e3t, чтобы найти конкретное решение.
62.Подставьте y = acos (2t) + bsin (2t) y = acos (2t) + bsin (2t) в y ′ + y = 4sin (2t) y ′ + y = 4sin (2t), чтобы найти конкретное решение.
63.Замените y = a + bt + ct2y = a + bt + ct2 на y ′ + y = 1 + t2y ′ + y = 1 + t2, чтобы найти конкретное решение.
64.Замените y = aetcost + betsinty = aetcost + betsint на y ′ = 2etcosty ′ = 2etcost, чтобы найти конкретное решение.
65.Решите y ′ = ekty ′ = ekt с начальным условием y (0) = 0y (0) = 0 и решите y ′ = 1y ′ = 1 с тем же начальным условием. Что вы замечаете, когда kk приближается к 0,0?
Методы решения дифференциальных уравненийСуществует несколько различных способов решения дифференциальных уравнений, которые я перечислю в приблизительном порядке их популярности. Я также классифицирую их иначе, чем в учебниках. Узнай или найди .Конечно! Уже решено очень много дифференциальных уравнений. Некоторые из них вы узнаете, а другие можете найти. Это , безусловно, наиболее распространенный способ, которым ученые или математики «решают» дифференциальные уравнения. Таким же образом некоторые (нечисловые) компьютерные программы решают дифференциальные уравнения. Замена . Часто дифференциальное уравнение можно упростить, заменив ту или иную переменную. Это может превратить ее в решение, которое уже решено (см. Выше) или которое может быть решено одним из других методов.(Пакеты программного обеспечения тоже это делают.) Эта категория решений включает в себя ряд методов, которым вы научитесь на втором курсе курса математики. Угадай и попробуй . Другой очень распространенный метод решения дифференциальных уравнений: угадайте, каким может быть решение, подставьте его и, если это не решение или не полное решение, измените предположение, пока не получите полное решение. Это используется часто — чаще, чем можно предположить, читая книги и статьи, где процесс обычно кажется довольно элегантным.Во многих случаях вы знаете что-то об изучаемой системе, что дает вам ключ к разгадке. Опыт, конечно, тоже помогает. Однако ниже мы увидим, что угадать иногда легко. Измените более простое решение . Если вы знаете решение уравнения, которое является упрощенной версией того, с которым вы столкнулись, попробуйте изменить решение более простого уравнения, чтобы преобразовать его в решение более сложного. Преобразование .Некоторые дифференциальные уравнения легче решать при математическом преобразовании. Это основное использование преобразований Лапласа. Численное решение. Если все вышеперечисленное терпит неудачу, то алгоритм, обычно реализуемый на компьютере, может решить эту проблему явно, вычисляя производные как отношения. Обычно это крайний метод по двум причинам. Во-первых, он дает вам решение только для одного конкретного набора граничных условий и параметров, тогда как все вышеперечисленное дает вам общие решения.Во-вторых, он имеет ограниченную точность: числовые дервиативы по своей сути зашумлены. Интеграция . Этот метод элегантен, но часто бывает трудным (или невозможным). Иногда можно умножить уравнение на коэффициент интегрирования, чтобы сделать интегрирование возможным. Особые типы . Это расплывчатое название включает специальные методы, которые работают с определенными типами уравнений. Это тоже касается обучения на высших курсах математики. Аналоговое решение. Некоторые дифференциальные уравнения легко решаются аналоговыми компьютерами. Они очень быстрые и поэтому подходят для задач управления «в реальном времени». Их недостатки — ограниченная точность и то, что аналоговые компьютеры сейчас редкость. Ниже мы показываем два примера решения общих уравнений. Они просты, потому что имеют только постоянные коэффициенты, но именно с ними вы столкнетесь на первом курсе физики. Эти уравнения можно решить несколькими способами, описанными выше, но мы проиллюстрируем только два метода. Пример 1: Экспоненциальный рост и спад
Кстати, на этом стоит остановиться и отметить, что дифференциальные уравнения почти всегда являются только приближениями. Невозможно иметь систему, описываемую этим уравнением.Например, популяция любого вида не может расти экспоненциально. Процитируем только одно ограничение: как только организмы занимают твердую сферу, радиус которой увеличивается со скоростью света, любой дальнейший рост не может быть экспоненциальным. (Об этом стоит помнить, когда политики становятся одержимыми достижением роста чего угодно, но особенно населения.)
Пример 2: Простое гармоническое движение
В направлении x второй закон Ньютона говорит нам, что F = ma = m.d 2 x / dt 2 , а здесь сила — kx. Это дает нам дифференциальное уравнение: где x — это смещение массы m от равновесия в момент времени t, а k — жесткость пружины, к которой прикреплена масса.
Это хорошее время, чтобы использовать метод Guess и попробовать метод решения . Нам нужно решение, которое постоянно колеблется и обладает тем свойством, что его вторая производная пропорциональна самой себе, но отрицательна. Все это делает функция синуса. Теперь мы не можем записать x = sin t по причинам размерности: аргумент синусоидальной функции не может иметь размерностей: он задается в радианах (что является соотношением или числом).Можем написать
Однако sin (ωt) — это число, и нам нужно, чтобы длина имела те же размеры, что и x, поэтому возможное решение: Когда мы описывали простое гармоническое движение, мы назвали A амплитудой : функция синуса изменяется от -1 до +1, поэтому движение изменяется от -A до + A.Однако есть проблема с этим предлагаемым решением: у него x = 0 при t = 0. Что было бы нормально, если бы я дал ему толчок, чтобы запустить его из состояния покоя, но что, если я отпущу массу из состояния покоя в какой-то точке из равновесия? В последнем случае мне понадобится x = A cos (ωt). Общее решение должно учитывать эти и любые другие начальные условия. Вместо этого мы пишем:
Итак, вернемся к рассмотрению φ. Если мы начнем движение (t = 0) с v = 0 при x = A, тогда φ должен быть 90 °: у нас есть функция cos вместо синуса. В качестве альтернативы, если мы начнем с максимальной (положительной) скорости при x = 0, тогда нам понадобится φ = 0. Мы приводим примеры этих случаев на странице фона для колебаний.Однако мы могли бы начать с любой комбинации начального смещения x = x 0 и v = v 0 . Итак, для общего случая (x 0 ≠ 0, v 0 ≠ 0) мы можем заменить, чтобы получить Мы можем решить их в терминах A и φ, сначала разделив два уравнения, затем возведя их в квадрат и сложив. Итак, для этих данных начальных условий мы можем найти комбинацию констант A и φ, так что это общее решение. Сколько граничных условий? В нашем первом примере нам нужно было найти только одну константу интегрирования, поэтому нужно было найти только одно начальное условие (или другое граничное условие).Второй пример — уравнение второго порядка, требующее двух интегрирований или двух граничных условий. Здесь мы могли бы указать два из начального смещения, скорость и ускорение, или некоторые другие два параметра. Затухающие и вынужденные колебания
Физически этот термин соответствует силе, пропорциональной скорости. Что мы можем догадаться о решении и как нам изменить полученное выше решение, чтобы оно удовлетворяло нашему новому дифференциальному уравнению? Опять же, мы можем использовать наши знания о физической системе: когда мы прикладываем силу, направление которой противоположно направлению скорости, мы замедляем ее.Таким образом, мы ожидаем одного из двух возможных ответов: либо он должен колебаться, при этом величина колебаний постепенно уменьшается с течением времени, либо (если демпфирование достаточно велико), он может замедлиться до остановки, даже не колеблясь. Это подсказывает нам возможность решения вида x = A e −βt sin (ωt + φ). Мы уже можем это попробовать. Но это не совсем решение. А что, если демпфирующая сила замедляет вибрацию? Почему бы не попробовать (ω + δω) вместо ω = k / m и посмотреть, дает ли это решение для подходящего значения δω? Давайте добавим еще одно усложнение: давайте начнем встряхивать частицу с дополнительной осциллирующей силой, скажем F = F 0 sin Ωt.Это дает нам новое дифференциальное уравнение:
Уравнения в частных производных: волновое уравнение
∂y / ∂x. Думайте об этом как dy / dx в заданное постоянное время t. Представьте, что вы делаете снимок (время постоянно): на изображении в момент времени t, это наклон формы y (x) в момент фотографии. ∂y / ∂t. Думайте об этом как dy / dt в данной позиции x. Это просто скорость в направлении y в определенной точке x на струне. (Кстати, не скорость волны). y = A sin (kx — ωt), поэтому Два нижних графика — вторые производные по тем же переменным:
∂y 2 / ∂t 2 = — ω 2 A sin (kx — ωt), что является ускорением точка на веревочке. Они имеют важное физическое значение: первый определяет кривизну струны. Если ∂y 2 / ∂x 2 = 0, то наклон постоянный, значит, прямой. Это означает, что натяжение T действует в противоположных направлениях на противоположных концах, не создавая чистой силы. Однако, если сегмент изогнут (∂y 2 / ∂x 2 ≠ 0), на него действует сила. Для постоянной кривизны на небольшой длине L чистая сила пропорциональна L. Нам известно ускорение, поэтому мы можем применить второй закон Ньютона. Масса сегмента равна мкл, где μ — масса на единицу длины μ. Запись закона Ньютона в виде a = F / m дает:
Оглядываясь назад на наши выражения для двух вторых производных, мы видим, что они являются нашими постоянными, умноженными на нашу исходную функцию y = A sin (kx — ωt). Это означает, что y = A sin (kx — ωt) является решением волнового уравнения при условии, что T / μ = ω 2 / к 2 .Мы также видели в Волны я, что ω / k — скорость волны, v. Который, наконец, связывает скорость волны с физическими свойствами T и μ струны: Скорость волны больше, если струна натянута сильнее, и меньше, если струна имеет большую массу на единицу длины. |
исчисление — Есть ли причина, по которой мы так редко можем решать дифференциальные уравнения?
Формулируя математические задачи, мы обычно формулируем их в терминах элементарных функций, но, безусловно, вычислимых функций, потому что это единственные, которые мы умеем записывать в конечном пространстве!
Поскольку наш мозг может только явным образом концептуализировать вычислимые функции, у нас есть врожденное предубеждение думать об этих функциях и придавать им центральное место в мире чисел.Когда вы читаете о диагонализации, возникает соблазн подумать: «Невычислимые числа — такая сложная задача для построения! Конечно, они должны быть редкостью!» Но реальность такова, что вычислимые функции — это виды, находящиеся под бесконечной угрозой исчезновения! Таких функций всего $ \ aleph_0 $, но не менее $ \ mathfrak c $ невычислимых функций.
Есть много способов перейти от вычислимого числа / функции к невычислимым (хорошо известный пример — диагонализация, а другим — проблема остановки), но я был бы удивлен, если бы кто-нибудь мог назвать «естественную» проблему который начинается с невычислимой функции / числа и чье решение вычислимо (и под «естественным» я подразумеваю такое, которое специально не для этого придумано).
Уравнение определяет пересечение двух функций. Если взять две произвольные функции, каковы шансы, что эти функции пересекутся по одной из бесконечно маловероятных счетных функций? Вот почему математики удивляются, когда результат имеет красивую замкнутую форму. Обычно этим свойством обладают только тривиальные задачи.
Конечно, есть задача решить, что такое «элементарная функция» или «аналитическое решение». Ответ: «Это не имеет значения.«Эти вопросы совершенно неуместны. Выберите любой конечный набор задач, который вам нравится. Давайте дадим имена решениям этих проблем, независимо от того, вычислимы они или нет. Теперь мы значительно расширили сферу« элементарных функций » Мы даже сделали что-то удивительное … мы добавили несколько невычислимых функций, которые действительно должны увеличить нашу способность решать проблемы, не так ли? Что ж, если вам не очень повезет, я бы поставил против.
Произвольная невычислимая функция — мусор.Это менее чем бесполезно. Хотя это решение бесконечно многих проблем и, таким образом, расширяет ваши возможности по написанию «решений закрытой формы» по крайней мере в $ \ aleph_0 $, я готов поспорить, что это также , а не решение (или удаленно актуальное ) до большей бесконечности проблем. Эти проблемы требуют различных невычислимых функций, отличных от тех, которые вы назвали.
Хорошо, хорошо … Я позволю тебе обмануть. Я позволю вам открыть панель инструментов и добавить еще несколько функций.Я не сказал, сколько вы можете добавить раньше, только то, что их должно было быть конечных . Вы могли бы добавить функции googol, мне все равно. На этот раз я буду очень щедрым. Я позволю вам добавить бесконечных невычислимых функций, до $ \ aleph_0 $ из них!
Конечно, теперь мы можем записать хорошие «алгебраические» решения для большинства проблем, учитывая, что мы увеличили наш набор инструментов в бесконечность! Но, к сожалению, нет. Наша бесконечность недостаточно велика.Независимо от того, насколько вы умны в выборе невычислимых функций, все равно будет бесконечно много проблем, для решения которых потребуется невычислимая функция, которую вы, , не выбрали .
Знаете что? Я чувствую себя щедрым. Мне плохо, потому что вы хотите, чтобы математика была красивой и красивой, а пока она выглядит как гигантский беспорядок. Мы пытались навести порядок, называя множество решений, у которых раньше не было названий. И пока мы остаемся ниже порога $ \ aleph_0 $, мы можем назначать конечные имена нашему «Набору инструментов расширенных алгебраических функций».Я сделаю тебе последнее одолжение. Я позволю вам добавить столько невычислимых функций, сколько вы хотите ! Это должно решить эту проблему раз и навсегда, верно?
Ну нет.