Метод характеристик при решение задачи коши для уравнений гиперболического типа
Стерлитамакский филиал Башкирский государственный университет
NovaInfo58, с. 11-15
Опубликовано 25 января 2017
Раздел: Физико-математические науки
Просмотров за месяц: 83
CC BY-NC
Аннотация
В статье рассматривается решение задачи Коши для уравнения гиперболического типа. Продемонстрировано решение данного уравнения методом характеристик.
Ключевые слова
Текст научной работы
Многие задачи физики, в частности механики, приводят к исследованию дифференциальных уравнений в частных производных второго порядка. Так, например, при изучении различных видов волн: звуковых, электромагнитных и других колебательных явлений приходят к волновому уравнению
где u=u(x,y,z,t), a — скорость распространения волны в данной среде. В одномерном случае это уравнение примет вид
которое является уравнением вынужденных колебаний однородной струны [1, 12].
В одномерном случае рассмотрим уравнение струны [2, 26]:
Задача Коши: Найти решение u(x,y) данного уравнения, удовлетворяющее начальным условиям:
Задача Коши для уравнения струны является математической моделью физической задачи о колебаниях настолько большой струны, что влияние ее концов уже не сказывается на колебаниях других точек струны. По этой причине в этой задаче отсутствуют граничные условия.
Приведем уравнение (1) к каноническому виду. Для этого составим уравнение характеристик
где A=0, 2B=e y , C=-1. Вычислим D=B^2-AC=\frac
. Следовательно, уравнение (1) является уравнением гиперболического типа.
Подставляя в уравнение характеристик наши значения, получим:
Основные типы уравнений математической физики
Основные типы уравнений
К основным уравнениям математической физики относятся следующие дифференциальные уравнения в частных производных второго порядка.
1. Волновое уравнение:
.
Это уравнение является простейшим уравнением гиперболического типа. К его исследованию приводит изучение процессов поперечных колебаний струны, продольных колебаний стержня, электрических колебаний в проводах и т. д.
2. Уравнение теплопроводности, или уравнение Фурье:
.
Это уравнение является простейшим уравнением параболического типа. К его исследованию приводит рассмотрение процессов распространения тепла, фильтрации жидкости и газа в пористой среде, изучение некоторых вопросов теории вероятностей и т. д.
3. Уравнение Лапласа:
.
Это уравнение относится к простейшим уравнениям эллиптического типа. К его исследованию приводит изучение задач об электрических и магнитных полях, о стационарном тепловом состоянии, задач гидродинамики и т. д.
В выписанных уравнениях искомая функция u зависит от двух переменных t, x или x, y. Рассматриваются также уравнения и для функций с большим числом переменных. Например, волновое уравнение с тремя независимыми переменными имеет вид
,
и уравнение Лапласа
.
Уравнение колебаний струны.
Формулировка краевой задачи
В математической физике струной называют гибкую упругую нить. Пусть струна в начальный момент времени расположена на отрезке 0≤x≤l оси Ox. Предположим, что ее концы закреплены в точках x=0 и x=l. Если струну отклонить от первоначального положения, а потом предоставить самой себе или придать ее точкам некоторую скорость, то точки струны будут совершать движение. Задача заключается в определении формы струны в любой момент времени и в определении закона движения каждой точки струны в зависимости от времени.
Если предположить, что движение точек струны происходит перпендикулярно оси Ox и в одной плоскости, то процесс колебания струны описывается одной функцией u(x,t), которая определяет величину перемещения точки струны с абсциссой x в момент t.
Доказано, что при отсутствии внешней силы функция u(x,t) должна удовлетворять дифференциальному уравнению в частных производных второго порядка
.
Для полного определения движения струны одного уравнения недостаточно. Искомая функция u(x,t) должна удовлетворять граничным условиям, указывающим, что делается на концах струны (при x=0 и x=l), и начальным условиям, описывающим состояние струны в начальный момент (t=0). Совокупность граничных и начальных условий называется краевыми условиями.
Пусть, например, концы струны при x=0 и x=l неподвижны. Тогда при любом t должны выполняться равенства
Это – граничные условия для рассматриваемой задачи. В начальный момент t=0 струна имеет определенную форму, которую мы ей придали. Пусть эта форма определяется функцией f(x), т. е.
Далее в начальный момент должна быть задана скорость в каждой точке струны, которая определяется функцией φ(x), т. е.
.
Эти два условия называются начальными условиями.
Колебания бесконечной струны.
Формула Даламбера решения задачи Коши
для волнового уравнения
Прежде чем решать задачу о колебаниях закрепленной струны, рассмотрим более простую задачу – о колебаниях бесконечной струны. Если представить очень длинную струну, то ясно, что на колебания, возникающие в ее средней части, концы струны не будут оказывать заметного влияния.
Рассматривая свободные колебания, мы должны решить однородное уравнение
при начальных условиях
, ,
где функции f(x) и g(x) заданы на всей числовой оси. Такая задача называется задачей с начальными условиями или задачей Коши.
Преобразуем волновое уравнение к каноническому виду, содержащему смешанную производную. Уравнение характеристик
распадается на два уравнения:
интегралами которых служат прямые
Введем новые переменные ξ=x – at, η=x + at и запишем волновое уравнение для переменных ξ и η.
, ,
,
,
и подставляя их в исходное уравнение, видим, что уравнение колебания струны в новых координатах будет
.
Интегрируя полученное равенство по η при фиксированном ξ, придем к равенству . Интегрируя это равенство по ξ при фиксированном η, получим
,
где φ и ψ являются функциями только переменных ξ и η соответственно. Следовательно, общим решением исходного уравнения является функция
. (8)
Найдем функции φ и ψ так, чтобы удовлетворялись начальные условия:
.
,
.
Интегрируя последнее равенство, получим:
,
где х0 и С – постоянные. Из системы уравнений
Таким образом, мы определили функции φ и ψ через заданные функции f и g, причем полученные равенства должны иметь место для любого значения аргумента. Подставляя в (8) найденные значения φ и ψ, будем иметь
.
Найденное решение называется формулой Даламбера решения задачи Коши для волнового уравнения
Пример. Решить уравнение при начальных условиях , .
Используя формулу Даламбера, сразу получаем
.
Решение волнового уравнения
методом разделения переменных
Метод разделения переменных применяется для решения многих задач математической физики. Пусть требуется найти решение волнового уравнения
, (9)
удовлетворяющее краевым условиям
u(x,0)=f(x), . (12),(13)
Частное решение уравнения (9), удовлетворяющее граничным условиям (10) и (11), ищут в виде произведения двух функций:
Подставляя функцию u(x,t) в уравнение (9) и преобразовывая его, получим
.
В левой части этого уравнения стоит функция, которая не зависит от x, а в правой – функция, не зависящая от t. Равенство возможно только в том случае, когда левая и правая части не зависят ни от x, ни от t, т. е. равны постоянному числу. Обозначим
, где λ>0. (14)
Из этих уравнений получаем два однородных дифференциальных уравнения второго порядка с постоянными коэффициентами
и . (15)
Общее решение этих уравнений
,
,
где A, B, C, D – произвольные постоянные.
Постоянные A и B подбирают так, чтобы выполнялись условия (10) и (11), из которых следует, что X(0)=X(l)=0, так как T(t)≠0 (в противном случае u(x,t)=0). Учитывая полученные равенства, находим
А=0 и .
Так как B≠0 (иначе, было бы X=0 и u=0, что противоречит условию), то должно выполняться равенство
,
.
Найденные значения λ называют собственными значениями для данной краевой задачи. Соответствующие им функции X(x) называются собственными функциями.
Заметим, что, если в равенстве (14) вместо – λ взять число λ (λ>0), то первое из уравнений (15) будет иметь решение в виде
.
Отличное от нуля решение в такой форме не может удовлетворять граничным условиям (10) и (11).
Зная , можем записать
.
Для каждого n получаем решение уравнения (9)
.
Так как исходное уравнение (9) линейное и однородное, то сумма решений также является решением, и потому функция
(16)
будет решением дифференциального уравнения (9), удовлетворяющим граничным условиям (10) и (11).
Найденное частное решение должно еще удовлетворять начальным условиям (12) и (13). Из условия (12) получим
.
Далее, дифференцируя члены ряда (16) по переменной t, из условия (13) будем иметь
.
Правые части двух последних равенств есть ряды Фурье для функций f(x) и φ(x), разложенных по синусам на интервале (0, l). Поэтому
. (17)
Итак, ряд (16), для которого коэффициенты Cn и Dn определяются по выписанным формулам, если он допускает двукратное почленное дифференцирование, представляет решение уравнения (9), удовлетворяющее граничным и начальным условиям.
Пример. Найти решение краевой задачи для волнового уравнения
, 0
Виды дифференциальных уравнений
Существует целый ряд задач, в которых установить прямую связь между величинами, применяемыми для описания процесса, не получается. Единственное, что можно сделать, это получить равенство, запись которого включает производные исследуемых функций, и решить его. Решение дифференциального уравнения позволяет установить непосредственную связь между величинами.
В этом разделе мы займемся разбором решений дифференциальных уравнений, неизвестная функция в которых является функцией одной переменной. Мы построили теоретическую часть таким образом, чтобы даже человек с нулевым представлением о дифференциальных уравнениях мог без труда получить необходимые знания и справиться с приведенными задачами.
Если какие-то термины окажутся для вас новыми, обратитесь к разделу «Определения и понятия теории дифференциальных уравнений». А тем временем перейдем к рассмотрению вопроса о видах дифференциальных уравнений.
Для каждого из видов дифференциальных уравнений применяется свой метод решения. В этом разделе мы рассмотрим все эти методы, приведем примеры с подробными разборами решения. После ознакомления с темой вам необходимо будет определять вид дифференциального уравнения и выбирать наиболее подходящий из методов решения поставленной задачи.
Возможно, прежде чем приступить к решению дифференциальных уравнений, вам придется освежить в памяти такие темы как «Методы интегрирования» и «Неопределенные интегралы».
Начнем ознакомление с темой мы с видов обыкновенных дифференциальных уравнений 1 -го порядка. Эти уравнения могут быть разрешены относительно производной. Затем перейдем в ОДУ 2 -го и высших порядков. Также мы уделим внимание системам дифференциальных уравнений.
Напомним, что y ‘ = d x d y , если y является функцией аргумента x .
Дифференциальные уравнения первого порядка
Простейшие дифференциальные уравнения первого порядка вида y ‘ = f ( x )
Начнем с примеров таких уравнений.
y ‘ = 0 , y ‘ = x + e x — 1 , y ‘ = 2 x x 2 — 7 3
Оптимальным для решения дифференциальных уравнений f ( x ) · y ‘ = g ( x ) является метод деления обеих частей на f ( x ) . Решение относительно производной позволяет нам прийти к уравнению вида y ‘ = g ( x ) f ( x ) . Оно является эквивалентом исходного уравнения при f ( x ) ≠ 0 .
Приведем примеры подобных дифференциальных уравнений:
e x · y ‘ = 2 x + 1 , ( x + 2 ) · y ‘ = 1
Мы можем получить ряд дополнительных решений в тех случаях, когда существуют значения аргумента х , при которых функции f ( x ) и g ( x ) одновременно обращаются в 0 . В качестве дополнительного решения в уравнениях f ( x ) · y ‘ = g ( x ) при заданных значениях аргумента может выступать любая функция, определенная для заданного значения х .
Наличие дополнительных решений возможно для дифференциальных уравнений x · y ‘ = sin x , ( x 2 — x ) · y ‘ = ln ( 2 x 2 — 1 )
Ознакомиться с теоретической частью и примерами решения задач таких уравнений вы можете в разделе «Простейшие дифференциальные уравнения 1 -го порядка».
Дифференциальные уравнения с разделяющимися переменными вида f 1 ( y ) · g 1 ( x ) d y = f 2 ( y ) · g 2 ( x ) d x или f 1 ( y ) · g 1 ( x ) · y ‘ = f 2 ( y ) · g 2 ( x )
Поговорим теперь об уравнениях с разделенными переменными, которые имеют вид f ( y ) d y = g ( x ) d x . Как следует из названия, к данному виду дифференциальных уравнений относятся выражения, которые содержат переменные х и у , разделенные знаком равенства. Переменные находятся в разных частях уравнения, по обе стороны от знака равенства.
Решить уравнения с разделенными переменными можно путем интегрирования обеих его частей: ∫ f ( y ) d y = ∫ f ( x ) d x
К числу дифференциальных уравнений с разделенными переменными можно отнести следующие из них:
y 2 3 d y = sin x d x , e y d y = ( x + sin 2 x ) d x
Для того, чтобы прийти от ДУ с разделяющимися переменными к ДУ с разделенными переменными, необходимо разделить обе части уравнения на произведение f 2 ( y ) ⋅ g 1 ( x ) . Так мы придем к уравнению f 1 ( y ) f 2 ( y ) d y = g 2 ( x ) g 1 ( x ) d x . Преобразование можно будет считать эквивалентным в том случае, если одновременно f 2 ( y ) ≠ 0 и g 1 ( x ) ≠ 0 . Если хоть одно из условий не будет соблюдаться, мы можем потерять часть решений.
В качестве примеров дифференциальных уравнений с разделяющимися переменными можно привести следующие из них: d y d x = y · ( x 2 + e x ) , ( y 2 + a r c cos y ) · sin x · y ‘ = cos x y .
К уравнениям с разделяющимися переменными мы можем прийти от ряда дифференциальных уравнений других видов путем замены переменных. Например, мы можем подставить в исходное уравнение z = a x + b y . Это позволит нам перейти к дифференциальному уравнению с разделяющимися переменными от дифференциального уравнения вида y ‘ = f ( a x + b y ) , a , b ∈ R .
Подставив z = 2 x + 3 y в уравнение y ‘ = 1 e 2 x + 3 y получаем d z d x = 3 + 2 e z e z .
Заменив z = x y или z = y x в выражениях y ‘ = f x y или y ‘ = f y x , мы переходим к уравнениям с разделяющимися переменными.
Если произвести замену z = y x в исходном уравнении y ‘ = y x · ln y x + 1 , получаем x · d z d x = z · ln z .
В ряде случаев прежде, чем производить замену, необходимо произвести преобразования исходного уравнения.
Предположим, что в условии задачи нам дано уравнение y ‘ = y 2 — x 2 2 x y . Нам необходимо привести его к виду y ‘ = f x y или y ‘ = f y x . Для этого нам нужно разделить числитель и знаменатель правой части исходного выражения на x 2 или y 2 .
Нам дано уравнение y ‘ = f a 1 x + b 1 y + c 1 a 2 x + b 2 y + c 2 , a 1 , b 1 , c 1 , a 2 , b 2 , c 2 ∈ R .
Для того, чтобы привести исходное уравнение к виду y ‘ = f x y или y ‘ = f y x , нам необходимо ввести новые переменные u = x — x 1 v = y — y 1 , где ( x 1 ; y 1 ) является решением системы уравнений a 1 x + b 1 y + c 1 = 0 a 2 x + b 2 y + c 2 = 0
Введение новых переменных u = x — 1 v = y — 2 в исходное уравнение y ‘ = 5 x — y — 3 3 x + 2 y — 7 позволяет нам получить уравнение вида d v d u = 5 u — v 3 u + 2 v .
Теперь выполним деление числителя и знаменателя правой части уравнения на u . Также примем, что z = u v . Получаем дифференциальное уравнение с разделяющимися переменными u · d z d u = 5 — 4 z — 2 z 2 3 + 2 z .
Подробный разбор теории и алгоритмов решения задач мы привели в разделе «Дифференциальные уравнения с разделяющимися переменными».
Линейные неоднородные дифференциальные уравнения первого порядка y ‘ + P ( x ) · y = Q ( x )
Приведем примеры таких уравнений.
К числу линейных неоднородных дифференциальных уравнений 1 -го порядка относятся:
y ‘ — 2 x y 1 + x 2 = 1 + x 2 ; y ‘ — x y = — ( 1 + x ) e — x
Для решения уравнений этого вида применяется метод вариации произвольной постоянной. Также мы можем представить искомую функцию у в виде произведения y ( x ) = u ( x ) v ( x ) . Алгоритмы применения обоих методов мы привели в разделе «Линейные неоднородные дифференциальные уравнения первого порядка».
Дифференциальное уравнение Бернулли y ‘ + P ( x ) y = Q ( x ) y a
Приведем примеры подобных уравнений.
К числу дифференциальных уравнений Бернулли можно отнести:
y ‘ + x y = ( 1 + x ) e — x y 2 3 ; y ‘ + y x 2 + 1 = a r c t g x x 2 + 1 · y 2
Для решения уравнений этого вида можно применить метод подстановки z = y 1 — a , которая выполняется для того, чтобы свести исходное уравнение к линейному дифференциальному уравнению 1 -го порядка. Также применим метод представления функции у в качестве y ( x ) = u ( x ) v ( x ) .
Алгоритм применения обоих методов приведен в разделе «Дифференциальное уравнение Бернулли». Там же можно найти подробный разбор решения примеров по теме.
Уравнения в полных дифференциалах P ( x , y ) d x + Q ( x , y ) d y = 0
Если для любых значений x и y выполняется ∂ P ( x , y ) ∂ y = ∂ Q ( x , y ) ∂ x , то этого условия необходимо и достаточно, чтобы выражение P ( x , y ) d x + Q ( x , y ) d y представляло собой полный дифференциал некоторой функции U ( x , y ) = 0 , то есть, d U ( x , y ) = P ( x , y ) d x + Q ( x , y ) d y . Таким образом, задача сводится к восстановлению функции U ( x , y ) = 0 по ее полному дифференциалу.
Выражение, расположенное в левой части записи уравнения ( x 2 — y 2 ) d x — 2 x y d y = 0 представляет собой полный дифференциал функции x 3 3 — x y 2 + C = 0
Для более подробного ознакомления с теорией и алгоритмами решения примеров можно обратиться к разделу «Уравнения в полных дифференциалах».
Дифференциальные уравнения второго порядка
Линейные однородные дифференциальные уравнения второго порядка с постоянными коэффициентами y ‘ ‘ + p y ‘ + q y = 0 , p , q ∈ R
Линейное однородное дифференциальное уравнение с постоянными коэффициентами обычно решается достаточно просто. Нам необходимо найти корни характеристического уравнения k 2 + p k + q = 0 . Здесь возможны три варианта в зависимости от различных p и q :
- действительные и различающиеся корни характеристического уравнения k 1 ≠ k 2 , k 1 , k 2 ∈ R ;
- действительные и совпадающие k 1 = k 2 = k , k ∈ R ;
- комплексно сопряженные k 1 = α + i · β , k 2 = α — i · β .
Значения корней характеристического уравнения определяет, как будет записано общее решение дифференциального уравнения. Возможные варианты:
- y = C 1 e k 1 x + C 2 e k 2 x ;
- y = C 1 e k x + C 2 x e k x ;
- y = e a · x · ( C 1 cos β x + C 2 sin β x ) .
Пример 13
Предположим, что у нас есть линейное однородное дифференциальное уравнение 2 -го порядка с постоянными коэффициентами y ‘ ‘ + 3 y ‘ = 0 . Найдем корни характеристического уравнения k 2 + 3 k = 0 . Это действительные и различные k 1 = — 3 и k 2 = 0 . Это значит, что общее решение исходного уравнения будет иметь вид:
y = C 1 e k 1 x + C 2 e k 2 x ⇔ y = C 1 e — 3 x + C 2 e 0 x ⇔ y = C 1 e — 3 x + C 2
Восполнить пробелы в теоретической части и посмотреть подробный разбор примеров по теме можно в статье «Линейные однородные дифференциальные уравнения 2 -го порядка с постоянными коэффициентами».
Линейные неоднородные дифференциальные уравнения второго порядка с постоянными коэффициентами y ‘ ‘ + p y ‘ + q y = f ( x ) , p , q ∈ R
Основным способом решение уравнений данного вида является нахождение суммы общего решения y 0 , которое соответствует линейному однородному дифференциальному уравнению y ‘ ‘ + p y ‘ + q y = 0 , и частного решения y
исходного уравнения. Получаем: y = y 0 + y
Способ нахождения y 0 мы рассмотрели в предыдущем пункте. Найти частное решение y
мы можем методом неопределенных коэффициентов при определенном виде функции f ( x ) , которая расположена в правой части записи исходного выражения. Также применим метод вариации произвольных постоянных.
К числу линейных неоднородных дифференциальных уравнений 2 -го порядка с постоянными коэффициентами относятся:
y ‘ ‘ — 2 y ‘ = ( x 2 + 1 ) e x ; y ‘ ‘ + 36 y = 24 sin ( 6 x ) — 12 cos ( 6 x ) + 36 e 6 x
Теоретические выкладки и подробный разбор примеров по теме можно найти в разделе «ЛНДУ 2 -го порядка с постоянными коэффициентами».
Линейные однородные дифференциальные уравнения (ЛОДУ) y ‘ ‘ + p ( x ) · y ‘ + q ( x ) · y = 0 и линейные неоднородные дифференциальные уравнения (ЛНДУ) второго порядка y ‘ ‘ + p ( x ) · y ‘ + q ( x ) · y = f ( x )
Линейные однородные и неоднородные дифференциальные уравнения и постоянными коэффициентами являются частными случаями дифференциальных уравнений этого вида.
На некотором отрезке [ a ; b ] общее решение линейного однородного дифференциального уравнения y ‘ ‘ + p ( x ) · y ‘ + q ( x ) · y = 0 представлено линейной комбинацией двух линейно независимых частных решений y 1 и y 2 этого уравнения, то есть, y = C 1 y 1 + C 2 y 2 .
Частные решения мы можем выбрать из систем независимых функций:
1 ) 1 , x , x 2 , . . . , x n 2 ) e k 1 x , e k 2 x , . . . , e k n x 3 ) e k 1 x , x · e k 1 x , . . . , x n 1 · e k 1 x , e k 2 x , x · e k 2 x , . . . , x n 2 · e k 2 x , . . . e k p x , x · e k p x , . . . , x n p · e k p x 4 ) 1 , c h x , s h x
Однако существуют примеру уравнений, для которых частные решения не могут быть представлены в таком виде.
Возьмем для примера линейное однородное дифференциальное уравнение x y ‘ ‘ — x y ‘ + y = 0 .
Общее решение линейного неоднородного дифференциального уравнения y ‘ ‘ + p ( x ) · y ‘ + q ( x ) · y = f ( x ) мы можем найти в виде суммы y = y 0 + y
, где y 0 — общее решение соответствующего ЛОДУ, а y
частное решение исходного дифференциального уравнения. Найти y 0 можно описанным выше способом. Определить y
нам поможет метод вариации произвольных постоянных.
Возьмем для примера линейное неоднородное дифференциальное уравнение x y ‘ ‘ — x y ‘ + y = x 2 + 1 .
Более подробно этот раздел освещен на странице «Линейные дифференциальные уравнения второго порядка».
Дифференциальные уравнения высших порядков
Дифференциальные уравнения, допускающие понижение порядка
Мы можем провести замену y ( k ) = p ( x ) для того, чтобы понизить порядок исходного дифференциального уравнения F ( x , y ( k ) , y ( k + 1 ) , . . . , y ( n ) ) = 0 , которое не содержит искомой функции и ее производных до k — 1 порядка.
В этом случае y ( k + 1 ) = p ‘ ( x ) , y ( k + 2 ) = p ‘ ‘ ( x ) , . . . , y ( n ) = p ( n — k ) ( x ) , и исходное дифференциальное уравнение сведется к F 1 ( x , p , p ‘ , . . . , p ( n — k ) ) = 0 . После нахождения его решения p ( x ) останется вернуться к замене y ( k ) = p ( x ) и определить неизвестную функцию y .
Дифференциальное уравнение y ‘ ‘ ‘ x ln ( x ) = y ‘ ‘ после замены y ‘ ‘ = p ( x ) станет уравнением с разделяющимися переменными y ‘ ‘ = p ( x ) , и его порядок с третьего понизится до первого.
В уравнении, которое не содержит аргумента х и имеет вид F ( y , y ‘ , y ‘ ‘ , . . . , y ( n ) ) = 0 , порядок может быть заменен на единицу следующим образом: необходимо провести замену d y d x = p ( y ) , где p ( y ( x ) ) будет сложной функцией. Применив правило дифференцирования, получаем:
d 2 y d x 2 = d p d y d y d x = d p d y p ( y ) d 3 y d x 3 = d d p d y p ( y ) d x = d 2 p d y 2 d y d x p ( y ) + d p d y d p d y d y d x = = d 2 p d y 2 p 2 ( y ) + d p d y 2 p ( y )
Полученный результаты подставляем в исходное выражение. При этом мы получим дифференциальное уравнение, порядок которого на единицу меньше, чем у исходного.
Рассмотрим решение уравнения 4 y 3 y ‘ ‘ = y 4 — 1 . Путем замены d y d x = p ( y ) приведем исходное выражение к уравнению с разделяющимися переменными 4 y 3 p d p d y = y 4 — 1 .
Более подробно решения задач по теме рассмотрены в разделе «Дифференциальные уравнения, допускающие понижение порядка».
Линейные однородные и неоднородные дифференциальные уравнения высших порядков с постоянными коэффициентами y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = 0 и y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = f ( x )
Решение уравнений данного вида предполагает выполнение следующих простых шагов:
- находим корни характеристического уравнения k n + f n — 1 · k n — 1 + . . . + f 1 · k + f 0 = 0 ;
- записываем общее решение ЛОДУ y 0 в стандартной форме, а общее решение ЛНДУ представляем суммой y = y 0 + y
— частное решение неоднородного дифференциального уравнения.
Нахождение корней характеристического уравнения подробно описано в разделе «Решение уравнений высших степеней». Для нахождения y
целесообразно использовать метод вариации произвольных постоянных.
Линейному неоднородному ДУ с постоянными коэффициентами y ( 4 ) + y ( 3 ) — 5 y ‘ ‘ + y ‘ — 6 y = x cos x + sin x соответствует линейное однородное ДУ y ( 4 ) + y ( 3 ) — 5 y ‘ ‘ + y ‘ — 6 y = 0 .
Более детальный разбор теории и примеров по теме вы можете найти на странице « Линейные однородные и неоднородные дифференциальные уравнения высших порядков с постоянными коэффициентами».
Линейные однородные и неоднородные дифференциальные уравнения высших порядков y ( n ) + f n — 1 ( x ) · y ( n — 1 ) + . . . + f 1 ( x ) · y ‘ + f 0 ( x ) · y = 0 и y ( n ) + f n — 1 ( x ) · y ( n — 1 ) + . . . + f 1 ( x ) · y ‘ + f 0 ( x ) · y = f ( x )
Найти решение ЛНДУ высших порядков можно благодаря сумме y = y 0 + y
, где y 0 — общее решение соответствующего ЛОДУ, а y
— частное решение неоднородного дифференциального уравнения.
y 0 представляет собой линейную комбинацию линейно независимых функций y 1 , y 2 , . . . , y n , каждая из которых является частным решением ЛОДУ, то есть, обращает равенство y ( n ) + f n — 1 ( x ) · y ( n — 1 ) + . . . + f 1 ( x ) · y ‘ + f 0 ( x ) · y = 0 в тождество. Частные решения y 1 , y 2 , . . . , y n обычно подбираются из известных систем линейно независимых функций. Подобрать их далеко не всегда просто и возможно, в этом и заключается основная проблема.
После того, как мы найдем общее решение ЛОДУ, найти частное решение соответствующего ЛНДУ можно благодаря методу вариации произвольных постоянных. Итак, y = y 0 + y
Получить более подробную информацию по теме можно в разделе «Дифференциальные уравнения высших порядков».
Системы дифференциальных уравнений вида d x d t = a 1 x + b 1 y + c 1 d y d t = a 2 x + b 2 y + c 2
Данная тема подробно разобрана на странице «Системы дифференциальных уравнений». Там же приведены примеры задач с подробных разбором.
http://pandia.ru/text/79/052/35879.php
http://zaochnik.com/spravochnik/matematika/differentsialnye-uravnenija/vidy-differentsialnyh-uravnenij/