21446

Обыкновенные дифференциальные уравнения

Лекция

Математика и математический анализ

Функция называется решением (или интегралом) д.у., если она раз непрерывно дифференцируема на некотором интервале и при удовлетворяет уравнению. Процесс нахождения решения д.у. называется его интегрированием...

Русский

2013-08-02

438.5 KB

3 чел.

Лекция 1.

Обыкновенные дифференциальные уравнения.

Общие понятия. Примеры.

Обыкновенными дифференциальными уравнениями (о.д.у.) называются уравнения вида:

                                                      (1)

где -известная функция,  независимая переменная,  неизвестная функция.

Порядком дифференциального уравнения (д.у.) называется наивысший порядок производной функции , входящей в уравнение.

Функция   называется решением (или интегралом) д.у. (1), если она  раз непрерывно дифференцируема на некотором интервале  и при  удовлетворяет уравнению. Процесс нахождения решения д.у. называется его интегрированием .

Пример 1. Пусть  непрерывная на интервале  функция ,  её первообразная, тогда

                                                                                (2)

т.е. для отыскания первообразной мы имеем о.д.у. I-ого порядка. Его решения:

          ,        

где   произвольная постоянная.

Д.у. (2) имеет бесчисленное множество решений и это справедливо для всех о.д.у. Чтобы выделить единственное решение уравнения (2), достаточно задать значение первообразной   в какой-либо точке, например, . Тогда решение единственно и равно

             

Основные элементарные функции являются решениями о.д.у.

Пример 2.

Уравнение:

                                                                     (3)   

имеет решения  что легко проверить. Функция  удовлетворяет условиям  а  удовлетворяет условиям

К о.д.у. приводят многие задачи естествознания.

Пример 3.

Движение материальной точки массы  описывается II законом Ньютона . Если точка движется по оси  и - её положение в момент времени , то функция  удовлетворяет о.д.у. II-ого порядка:

                                                          (4)

Чтобы найти единственное решение уравнения (4), необходимо  задать начальные данные

                                                               (5)

  Если задача об отыскании всех решений о.д.у. сводится к вычислению конечного числа интегралов и производных от известных функций и к алгебраическим операциям, то говорят, что уравнение интегрируется в квадратурах. Класс таких уравнений узок, Поэтому для исследования д.у. как правило применяются приближенные и численные методы.

         Пример 4.

Модель типа хищник-жертва. Рассмотрим динамику популяции 2-х видов, взаимодействующих между собой по типу хищник-жертва. При этом предполагается, что жертва может найти достаточно пищи, но и при каждой встрече с хищником последний убивает жертву.

Пусть  и - соответственно количество жертв и хищников в момент времени . Пусть также коэффициент рождаемости жертв  и коэффициент естественной смертности (т.е. без учета уничтожения хищником)  являются константами, причем . Таким образом, в отсутствие хищников популяция жертв будет расти  со скоростью . Кроме того, предположим, что число случаев, когда хищник убивает жертву зависит от вероятности их встречи и, следовательно, пропорционально произведению .

Таким образом, популяция жертв удовлетворяет следующему д.у.:

где   , а

Относительно хищников предположим, что их число при отсутствии жертв по естественным причинам убывает со скоростью . В результате  встреч с жертвами число хищников увеличивается. Таким образом, имеем уравнение:

 

где   , а

Итак, модель хищник-жертва сводится к системе 2-х нелинейных о.д.у.

                                                   (a)

0, 0, 0, 0. Это уравнения Лотки-Вольтерра. Система должна быть дополнена начальными условиями

 .                                                               (b)

Явное аналитическое решение задачи (a), (b) неизвестно. Рассмотрим метод решения этой задачи, основанный на теории возмущений.

Первый шаг состоит в нахождении стационарных состояний (, ) или точек равновесия. Эти точки находятся из условий:

 

откуда

.                                                  (с)

Кроме того, точкой равновесия будет так называемое тривиальное решение или точка покоя

     .

Разложим правые части уравнений (а) в ряды Тейлора в окрестности точки (,). С учетом (с) получим

 .

Таким образом, в окрестности этой точки исходные уравнения (а) можно аппроксимировать линейными

                                                                      (d)

Положим

      

Тогда система (d) примет вид:

                                                                           (e)

Продифференцируем по t, например, I-е уравнение  и подставим в него второе

                  .

Легко проверить, что общим решением этого уравнения является

                      ,                                       (f)

где  и  - const, определяемые из начальных условий (b). Из (e) получаем:

                       .                  (g)

Исключая из (f) и (g) время, имеем:

                                       

Возвращаясь к исходным обозначениям

                          ,                                                   (h)   

где мы обозначили.

Соотношение (h) – это уравнение эллипсов (в зависимости от с2, т.е. от начальных условий) с центром в точке ).

Стрелки указывают направление, соответствующее возрастанию времени. Из (d) видно, что при , растет  (т.к. >0) и наоборот.

Видно, что изменение популяций носит циклический характер: через определенное время популяция возвращается к первоначальному уровню.

Уравнения первого порядка, разрешенные относительно производной.

      Такого рода уравнение имеет вид:

                                                                                     (6)

В дальнейшем будет доказано, что при некоторых ограничениях, налагаемых на функцию  уравнение (6) имеет единственное решение, удовлетворяющее условию:

                                      .

Его общее решение, содержащее все без исключения решения, зависит от одной произвольной постоянной.

Дифференциальное уравнение  устанавливает зависимость между координатами точки и угловым коэффициентом касательной  к графику решения в этой точке. Зная x и y, можно вычислить . Следовательно, дифференциальное уравнение рассматриваемого вида определяет поле направлений и задача интегрирования д.у. заключается в том, чтобы найти кривые, называемые интегральными кривыми1, направление касательных, к которым в каждой точке совпадает с направлением поля. Пример 5.

                          

Для построения поля направлений найдем геометрическое место точек, в которых касательные к искомым интегральным кривым сохраняют постоянное значение. Такие линии называются изоклинами. Уравнение изоклин получим, считая, что

                                ,

где k- постоянная. Таким образом,

                        ,

т.е. в данном случае изоклинами являются окружности с центром в начале координат, причем угловой коэффициент касательной к искомым интегральным кривым равен радиусу этих окружностей. Для построения поля направлений дадим постоянной k некоторые определенные значения. После этого можно приближенно провести искомые интегральные кривые.

      Пример 6.

                          

Изоклинами являются гиперболы или причем при k=1 гипербола распадается на пару прямых x=0 и y=0. При k=0 получаем изоклину ; эта гипербола разбивает плоскость на

части, в каждой из которых y сохраняет постоянный знак. Интегральные кривые  пересекая гиперболу

переходят из области возрастания функции y(x) в область её убывания или наоборот, и, следовательно, на ветвях этой гиперболы расположены точки максимума и минимума интегральных кривых.

Определим теперь знаки второй производной в различных областях плоскости:

    или        

Кривая    или

                          (7)

разбивает плоскость на две части, в одной из которых y<0, и, следовательно, интегральные кривые выпуклы вверх, а в другой y>0, и значит, инте

гральные кривые выпуклы вниз. При переходе через кривую (7)

интегральные кривые переходят от выпуклости к вогнутости, и, следовательно, на этой кривой расположены точки перегиба интегральных кривых.

В результате проведенного исследования известны области возрастания и убывания интегральных кривых, известно расположение точек максимума и минимума, известны области выпуклости и вогнутости и расположение точек перегиба, известна изоклина k=1. Этих сведений вполне достаточно для того, чтобы сделать набросок расположения интегральных кривых.

      Рассмотрим некоторые уравнения I порядка, которые интегрируются в квадратурах.

  1.  Уравнение с разделяющимися переменными.

Это уравнения вида:

                         

                                                                                     (13)

будем считать, что функции ,   непрерывны при , а, кроме того ,  т.к. если при , , то из (13) следует, что  есть решение.

Запишем уравнение (13) в виде:

                   

и проинтегрируем

                                                                     (14)

это и есть общее решение уравнения (13).

Пример 7.

                   .

Здесь

          ,

так что , где С - произвольная постоянная. Есть еще решение y0.

Решением задачи Коши

при  

здесь будет:

 .

  1.  Уравнения, приводящиеся к уравнениям с разделяющимися переменными.

а) Рассмотрим уравнение

 ,

где  .

Перейдем к новым переменным , в результате получим

                    т.е.       

или

       

т.е. переменные разделились. Интегрируя, получим:

                   

б) Рассмотрим теперь однородное д.у. I порядка

     

однородная функция нулевого порядка, т.е.  

                 

Положим , т.е. , тогда

 

  .   

в) Уравнение вида

 

будет однородным, если

 ,

т.е. ,однородные функции одного порядка, (), и будет уравнением с разделяющимися переменными, если

                       .

1) Пусть  –это решение уравнения (6) на отрезке [a,b] оси х. График этой функции   называется интегральной кривой уравнения (6).


 

А также другие работы, которые могут Вас заинтересовать

42995. Разработка привода и натяжной станции подземного ленточного конвейера 5.59 MB
  Современное массовое и крупносерийное производство продукции разнообразных отраслей промышленности выполняется поточным методом с широким использованием автоматических линий. Поточный метод производства и работа автоматической линии основаны на конвейерной передаче изделий от одной технологической операции к другой. Следовательно конвейеры являются составной и неотъемлемой частью современного технологического процесса – они устанавливают и регулируют темп производства, обеспечивают его ритмичность, способствуют повышению производительности труда и увеличению выпуска продукции. Конвейеры являются основными средствами комплексной механизации и автоматизации транспортных и погрузочно-разгрузочных работ и поточных технологических операций.
42996. Расчет подстанции (п/ст) «Симахинская» 1.85 MB
  Питание данной подстанции осуществляется воздушной линией электропередач 110 кВ от подстанции Таежная. Описание существующей схемы электрических соединений подстанции Схема электрических соединений подстанции рис. Главными признаками определяющими тип подстанции являются её местоположение назначение и роль в энергосистеме число и мощность установленных трансформаторов их тип и высшее напряжение. Все подстанции можно разбить на три основные категории [78]: – по упрощенным схемам как правило без выключателей на стороне высокого...
42998. Обеспечение аварийным источником электроснабжения потребителей объекта по адресам: г. Санкт-Петербург, Лиговский пр.37, ул. Восстания д.1, ул. Восстания д.6 2.69 MB
  Проектом предусматривается установка стационарного дизельэлектрического агрегата мощностью 350 кВА и мобильного дизельэлектрического агрегата мощностью 150 кВА для аварийного электроснабжения потребителей объекта. Дизельэлектрические агрегаты размещаются по адресу: г. Запас дизельного топлива для каждого дизельэлектрического агрегата расположен во встроенном топливном баке. Емкости топливных баков обеспечивают время работы не превышающее 8 часов работы дизельэлектрических агрегатов на номинальном режиме.
42999. Розрахунок деталі вал-шестерня 2.02 MB
  При проектуванні дискових шеверов необхідно прагнути до вибору максимальних зовнішніх діаметрів шевера. Якщо пройняти до уваги що вказані параметри гвинтового зачіпляє роблять вплив на розміри дискового шевера то останні повинні бути різними у міру переточування зубів шевера. Це зумовлює одну з особливостей розрахунку дискового шевера – необхідність визначення розмірів шевера при різних ступенях його сточенності. Це умова – забезпечення повної обробки активної частини профілю колеса – є основним при розрахунку шевера.
43001. Круглый фасонный резец с радиальной подачей для обработки деталей 1.1 MB
  Форма режущего лезвия определяется формой обрабатываемого изделия формой профиля обрабатываемой поверхности. 57 квалитеты; Идентичность формы причем точность фасонного резца на одиндва класса выше детали; Высокая производительность обработки за счет экономии времени – обрабатываются одновременно все участки фасонного профиля детали; Большой срок службы – за счет увеличения количества переточек по передней поверхности; Применение фасонных резцов не требует высокой квалификации рабочего. При изучении исходных данных следует...
43002. Динамический и силовой анализ механизма 99.5 KB
  Динамический анализ механизма включает в себя определение движущего момента такого, чтобы звенья механизма двигались с заданными скоростями при заданных нагрузках и массах действующих на механизм. Формула для определения движущего момента
43003. Расчет заготовки коробчатой формы 364.5 KB
  Основной критерий оборудования − это номинальное усилие пресса. 1 операция Вырубка заготовки При вырубке круглой заготовки необходимо учитывать следующие усилия: Pвыр− усилие затрачиваемое на вырубку заготовки пуансоном; Рпрот− усилие затрачиваемое на проталкивание заготовки пуансоном через отверстие в матрице; Ртр− усилие затрачиваемое на трение отхода о пуансон. Усилие операции: Pоп = Pвыр Рпрот Ртр. Усилие вырубки заготовки: k=1113 – коэффициент учитывающий притупление режущих кромок неравномерность зазора...