Движение которых определяется дифференциальными уравнениями

Интегрирование дифференциальных уравнений прямолинейного движения материальной точки

Краткое изложение результатов

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

Сила, зависящая от времени

Если на материальную точку действует сила, зависящая от времени , то дифференциальное уравнение прямолинейного движения вдоль оси Ox имеет вид:
.
Вводим ускорение и интегрируем это уравнение.
.
Здесь и далее A и B – произвольные точки на оси Ox . Заменим . Получаем закон изменения скорости от времени:
.
Интегрируя уравнение , получаем закон движения точки :
;
.

Сила, зависящая от скорости

Пусть на точку действует сила, зависящая от скорости . Составляем дифференциальное уравнение движения и интегрируем его:
.
Последнее уравнение дает в неявном виде зависимость . Решаем его. После чего интегрируем уравнение , как описано выше.

Есть второй способ интегрирования уравнения движения в случае зависимости силы от скорости. Для этого переходим от переменных x и t к переменным и x . Считаем, что скорость является функцией от координаты x :
;
.
Последнее уравнение дает в неявном виде зависимость . Далее интегрируем уравнение :
.
Это уравнение дает в неявном виде закон движения точки .

Сила, зависящая от перемещения

Пусть на точку действует сила, зависящая от перемещения . Составляем уравнение движения, переходим от переменных x и t к переменным и x , и интегрируем полученное дифференциальное уравнение:
;
;
.
Это уравнение представляет собой закон сохранения механической энергии для прямолинейного движения. Из него находим зависимость скорости от перемещения . После чего интегрируем уравнение , как это описано выше.

Дифференциальное уравнение движения точки

Рассмотрим прямолинейное движение материальной точки, находящейся под действием постоянных и переменных сил. Направим ось Ox системы координат вдоль линии движения точки. Пусть на нее действуют n сил, проекции которых на ось Ox мы обозначим как . Положение точки, при прямолинейном движении, однозначно определяется ее координатой x . Нам нужно определить закон движения точки , то есть закон изменения ее координаты со временем.

Уравнение движения точки определяется вторым законом Ньютона, который в случае прямолинейного движения имеет вид:
(1) .

Вместо того, чтобы в каждом уравнении выписывать все n сил, введем их равнодействующую, проекция которой, на ось x равна сумме проекций всех сил на эту ось:
.
Тогда задача сведется к движению материальной точки под действием одной силы . При этом уравнение движения примет наиболее простой вид:
(2) .
В дальнейшем, проекцию равнодействующей мы будем называть просто силой, действующей на точку.

Сила может быть как постоянной, так и зависеть от времени t , координаты x и от скорости . К сожалению, если зависит от всех перечисленных факторов, то не всегда возможно решить уравнение (2) аналитически. Поэтому мы рассмотрим те случаи, когда возможно получить аналитическое решение этого уравнения. Заметим, что если сила является постоянной, то уравнение (2) можно решать любыми, приводимыми ниже, способами.

Почему мы обозначаем в виде проекции силы на ось x , хотя рассматриваем только движение вдоль одной оси? – Потому что под обозначением силы R в виде одной буквы, часто подразумевается ее абсолютная величина: . Она имеет неотрицательные значения: . А когда мы пишем силу как проекцию , то подразумеваем, что эта величина может быть как положительной (если сила направлена вдоль оси x ), так и отрицательной (когда она направлена противоположно оси x ). В теоретической механике, в подобных случаях, иногда также говорят, что есть алгебраическое значение силы. Это относится не только к силе, ни и ко всем другим, рассматриваемым далее, векторным величинам.

Интегрирование уравнения движения с силой, зависящей от времени

Вначале рассмотрим случай, когда задан закон изменения силы со временем: . Перепишем уравнение (2), явно указав эту зависимость:
(t1) .

В этом уравнении время t является независимой переменной; координата x – зависимой переменной; – это вторая производная координаты по времени: . Масса m – это постоянная, то есть заданное число. С математической точки зрения, уравнение (t1) есть дифференциальное уравнение второго порядка, не содержащую зависимую переменную x в явном виде.

Решение такого уравнения выполняется с помощью подстановки
.
Тогда
.
Подставляя в (t1), мы получаем дифференциальное уравнение первого порядка:
(t2) .
Выполняя подстановку, мы ввели новую переменную , равную производной координаты x по времени t . Эта производная является проекцией скорости точки на ось Ox . Таким образом, процесс решения разбивается на две части. Сначала мы, решаем уравнение (t2), и находим закон изменения скорости со временем: . Затем, используя уравнение , находим закон изменения координаты .

Упростим уравнение (t2), разделив его на массу m :
(t3) ,
где – ускорение точки. Поскольку зависимость силы от времени известна, то и зависимость ускорения от времени также известна.

Уравнение (t3) является дифференциальным уравнением первого порядка с разделяющимися переменными. Разделяем переменные и интегрируем.
;
(t4) ;
(t5) .
Здесь – постоянная интегрирования. Чтобы ее определить, нужно знать значение скорости в какой-либо момент времени . Если мы сможем выразить интеграл через известные функции, то подставив в (t5) значения времени и значение скорости в этот момент, мы сможем определить постоянную .

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

Прямолинейное движение точки M под действием силы R_x.

Пусть нам известно, что в момент времени , точка M находилась в положении A, имела координату и скорость . Рассмотрим произвольный момент времени . Пусть в этот момент времени точка M находится в положении B, с координатой и скоростью . Величины и нам пока не известны. Наша задача их найти.

Перепишем (t4) явно указав, что есть функция от t :
(t6) .
Интегрируем (t6) от момента времени до :
.
Слева – интеграл от полного дифференциала. Поэтому он интегрируется элементарно:
.
Здесь мы учли, что . В результате получаем:
;
.

Этот результат можно получить и несколько иначе, если в интеграле сразу перейти к переменной . Тогда пределы интегрирования станут и . В результате получим тоже самое:
.

Итак, мы нашли значение скорости в произвольный момент времени :
(t7) .
Заменим обозначение момента времени на t . В результате получим закон изменения скорости со временем t :
(t8) .

Интеграл справа записан не вполне корректно, хотя так часто пишут. Рассмотрим пример определенного интеграла . Он зависит от пределов интегрирования a и b , но не зависит от переменной интегрирования t . Можно сказать, что переменная t принимает заданные значения из отрезка , которые применяются только для вычисления интеграла. Поэтому для переменной интегрирования t можно использовать любое обозначение. Например, можно использовать переменную . Тогда .

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

Теперь найдем закон изменения координаты x от времени. Интегрируем уравнение
.
Разделяем переменные:
(t10) .
Здесь мы также можем выполнить интегрирование от A до B, но мы продемонстрируем другой способ, как получить результат в удобном виде, применяя неопределенный интеграл. Поскольку неопределенный интеграл определен с точностью до постоянной, то запишем его с нижним пределом интегрирования . Интегрируем (t10):
(t11) .
Найдем значение постоянной интегрирования . Для этого подставим сюда :
.
Далее учитываем, что значение координаты точки в момент времени нам известно: . Также учитываем, что интеграл в правой части имеет равные пределы интегрирования и поэтому равен нулю. В результате получаем:
.
Отсюда находим значение постоянной интегрирования: . В результате получаем закон движения точки:
(t12) .

Итак, мы нашли, что если на точку действует сила , то для определения ее закона движения, нужно сначала определить закон изменения скорости со временем:
(t7) .
А затем определить закон движения:
(t12) .
При этом мы полагаем, что нам известны скорость и координата в некоторый момент времени . Если бы мы проводили интегрирование через неопределенные интегралы в общем виде, то и были бы постоянными интегрирования и .

Постоянная сила

Разберем случай, когда действующая на точку сила имеет постоянное значение: . В этом случае ускорение также постоянно: . Интегрируем, используя таблицу неопределенных интегралов. Из (t7) находим закон изменения скорости со временем:
;
(t14) .
Мы видим, что скорость линейно изменяется со временем.

Подставляем в (t12) и находим закон движения точки:

;
(t15) .

Если в начальный момент времени , скорость точки была , а координата , то . Из (t14) и (t15) получаем:
;
.

Равномерное движение

Если проекция силы на ось Ox равна нулю: , то ускорение также равно нулю: . В этом случае из (t14) находим, что скорость точки постоянна:
.
Из (t15) находим, что координата линейно меняется со временем:
.

Если в начальный момент времени , скорость точки была , а координата , то ;
;
.

Интегрирование уравнения движения с силой, зависящей от скорости

Разберем случай прямолинейного движения материальной точки, когда действующая сила зависит от скорости . Такие задачи встречаются при движении в жидкой или газообразной среде, когда на точку помимо постоянных сил, действует сила трения, зависящая от скорости. В этом случае, уравнение движения имеет вид:
(v1) .
Разделим обе части уравнения на массу m :
(v2) ,
где – ускорение точки. Теперь нам известна зависимость ускорения точки от ее скорости. Уравнение (v2) не содержит в явном виде как зависимую переменную x , так и независимую переменную t . Поэтому его можно решать двумя способами.

Решение уравнения, определяя v_x(t)

Применим к уравнению (v2) метод решения дифференциального уравнения, не содержащего зависимую переменную в явном виде. Для этого, как и в предыдущем случае, делаем подстановку
.
Тогда
.
Подставляя в (v2), мы получаем дифференциальное уравнение первого порядка:
(v3) .

Пусть, как и в предыдущем случае, в момент времени , точка находилась в положении A, имела координату и скорость . И пусть в произвольный момент времени , точка находится в положении B с координатой и скоростью . Нам нужно найти величины и .

Разделяем переменные.
;
.
Перепишем это уравнение, указав, что скорость является функцией от времени:
.
Интегрируем по времени от до :
.
В левой части сделаем замену переменной. От переменной t перейдем к переменной . При этом изменим пределы интегрирования учитывая, что при ; и при :
(v4) .

Заменим обозначения переменных , и переменной интегрирования . Подставим в (v4):
(v5) .
Это уравнение, в неявном виде, дает закон изменения скорости от времени t . Вычислив интеграл, и выполнив преобразование, мы можем выразить через t : .

Далее, по формуле (t12) ⇑ определяем закон движения материальной точки:
(t12) .

Решение уравнения, определяя v_x(x)

Выпишем уравнение (v2) еще раз.
(v2) .
Для применения этого метода, в качестве независимой переменной возьмем координату x , а в качестве зависимой – скорость . То есть считаем, что скорость является функцией от координаты: .

Выразим через переменные x и вторую производную координаты по времени:
.
Подставим в (v2) и разделяем переменные:
;
.
В левой части в явном виде запишем как функцию от x , и интегрируем по x от положения A до B:
;
.
В интеграле слева переходим от переменной x к :
(v6) .

Переобозначим переменные:
(v7) .
Это уравнение дает в неявном виде зависимость скорости от координаты:
.
Подставив сюда , получим для x дифференциальное уравнение первого порядка. Решаем его методом разделения переменных:
;
.
Интегрируем от положения A до B:
;
.
Заменим переменные:
(v8) .

Уравнение (v8) дает в неявном виде закон движения материальной точки .

Интегрирование уравнения движения с силой, зависящей от перемещения

Наконец рассмотрим случай прямолинейного движения материальной точки, когда действующая сила зависит от перемещения x . Такие задачи встречаются при движении в потенциальных полях – в гравитационных или электрических. Сюда также относится движение груза, прикрепленного к упругой пружине.

Выписываем уравнение движения для этого случая:
(x1) .
Это дифференциальное уравнение второго порядка. Оно не содержит независимую переменную t в явном виде. Также как и в предыдущем случае, применяем метод решения дифференциального уравнения, не содержащего независимую переменную в явном виде.

Перейдем к новым переменным. В качестве независимой переменной возьмем координату x , а в качестве зависимой – скорость . Считаем, что скорость является функцией от координаты: .

Выразим вторую производную координаты по времени через переменные x и :
;
Подставим в (x1) и разделяем переменные:
(x2) ;
.
Интегрируем по x от A до B:
(x3) .
Вычисляем интеграл, используя таблицу неопределенных интегралов:
;
.
Подставляем в (x3):
(x4) . Нетрудно видеть, что слева стоит изменение кинетической энергии материальной точки. Справа – работа, которую совершает сила при перемещении материальной точки из A в B. Само уравнение (x4) представляет собой теорему об изменении кинетической энергии точки для прямолинейного движения.

Вернемся снова к уравнению (x2).
(x2) .
Его можно проинтегрировать и другим способом.

Для этого представим правую часть в виде производной по координате:
,
где – координата произвольной заранее выбранной точки C .
Левую часть также представим в виде производной по координате:
.
Тогда (x2) можно записать в виде:
.

Поскольку производная по x от выражения в скобках равна нулю, то само выражение является постоянной, не зависящей от x величиной:
.
Такая форма записи, когда некоторая функция от переменных приравнивается постоянной, называется интегралом дифференциального уравнения. Перепишем его в следующем виде:
(x5) .
Здесь – кинетическая энергия точки; – потенциальная энергия, отсчитываемая от, произвольным образом выбранной, точки C ; E – постоянная интегрирования, которая в данном случае имеет определенный физический смысл – это полная механическая энергия материальной точки. Поэтому мы ее обозначили привычной для этого случая буквой E . Само уравнение (x5) представляет собой закон сохранения механической энергии. С математической точки зрения, энергия E является интегралом дифференциального уравнения, или, как говорят в механике, интегралом движения точки. То есть величиной, сохраняющей при движении постоянное значение.

Выше мы пришли к выводу, что постоянная интегрирования E не зависит от координаты x , но ничего не сказали о ее зависимости от времени. Однако, для одномерного движения, со временем может изменяться только одна координата x . Поскольку постоянная E от нее не зависит, то она не зависит также и от времени t . Поэтому полная механическая энергия сохраняет постоянное значение и в различные моменты времени.

Нетрудно видеть, что формулировки (x4) ⇑ и (x5) ⇑ эквивалентны. Для доказательства, приравняем механическую энергию точки для двух положений A и B:
;

.
Здесь мы разбили интеграл от до на два интеграла – от до ; и от до . Интегралы от до сократились.

Найдем зависимость скорости точки от координаты. При этом мы считаем, что скорость точки в положении A нам известна. Рассмотрим два положения: A и B. Из (x4) ⇑ имеем:
,
где – работа, которая производит сила при перемещении точки из A в B. Наконец, заменим на x , и на . В результате получим искомую зависимость:
(x6) ,
где – работа, которая производит сила при перемещении материальной точки из A в точку с координатой x . Скорость определена с точностью до знака (плюс или минус). Знак нужно выбирать из начальных условий и исследования движения. Если в точке , то при достаточно малых значениях . Далее точка может остановиться и начать движение в обратную сторону. Тогда нужно выбрать знак минус, чтобы скорость стала отрицательной.

Теперь, зная зависимость , находим закон движения материальной точки. Для этого интегрируем уравнение:
;
;
;
.
Это уравнение дает в неявном виде зависимость координаты x от времени t .

Приложение к движению в пространстве

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

Пусть нам известно, что в момент времени , материальная точка находилась в точке A, и имела скорость . Выберем трехмерную систему координат Oxyz , и распишем эти начальные условия по компонентам:
При ;
При ;
При .

Сила в пространстве, зависящая от времени

Пусть на материальную точку действует сила, зависящая от времени: . Составим уравнения ее движения:
.

Выпишем уравнение для координаты x с начальными условиями:
; при .
Здесь все необходимые величины известны, и они не зависят от значений других координат. Мы можем найти закон изменения координаты x со временем, применяя интегрирование уравнения движения с силой, зависящей от времени ⇑ для прямолинейного движения.

Выпишем уравнение для координаты y с начальными условиями:
; при .
Здесь также известны все необходимые величины, и они не зависят от значений других координат. Мы также можем найти закон изменения координаты y со временем, применяя интегрирование, как для прямолинейного движения.

Точно также мы можем найти закон изменения координаты z со временем. В этом случае говорят, что переменные разделились. Уравнения движения, составленные для каждой из координат, вместе с начальными условиями, не зависят от значений других координат. Поэтому каждое такое уравнение можно проинтегрировать отдельно. В результате мы получим закон движения материальной точки в трехмерном случае: .

Силы, приводящие к разделению переменных

Пусть теперь на точку действуют три взаимно перпендикулярные силы. И пусть одна из них зависит только от времени; вторая – от проекции скорости на направление силы; третья – от проекции радиус-вектора на направление силы.

Выберем систему координат Oxyz , оси которой направим вдоль направлений действующих сил. Тогда в этой системе координат отличными от нуля будут только три проекции сил: . Составляем уравнения движения:
;
;
.
Мы видим, что и в этом случае переменные разделились. Каждое из этих уравнений зависит только от одной переменной. И мы можем решить его, применяя изложенные выше методы. Все это применимо и к случаю, когда любая из этих сил является постоянной.

И, разумеется, тут могут быть различные вариации, приводящие к разделению переменных. Например, если зависящая от времени сила лежит в плоскости xy , а перпендикулярная ей сила зависит только от координаты z . В этом случае переменные также разделяются.

Автор: Олег Одинцов . Опубликовано: 12-10-2020

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

Содержание:

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

При решении многих задач математики, техники, экономики и других отраслей науки бывает трудно установить закон, связывающий искомые и известные переменные величины. Но удается установить связь между производными или дифференциалами этих переменных, которая выражается уравнениями или системами уравнений. Такие уравнения называют дифференциальными уравнениями. Термин «дифференциальное уравнение» введен в 1676 году В. Лейбницом.

Мы рассмотрим только уравнения с функциями одной переменной и обычными производными, которые называют обычными дифференциальными уравнениями.

Основные понятия о дифференциальных уравнениях

Определение. Дифференциальным уравнением называется уравнение, которое связывает независимую переменную x, искомую функцию y = f (x) и еепроизводные или дифференциалы разных порядков, то есть уравнение
(7.1)

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

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

Так, уравнение y’ – 2 xy 2 + 5 = 0 является дифференциальным уравнением первого порядка, а уравнения y» + 2 y’ – y – sin x = 0 — дифференциальным уравнением второго порядка.

Определение. Решением дифференциального уравнения (7.1) называется такая функция y = φ (x), которая при подстановке в уравнение (7.1) превращает его в тождество.

Например, для дифференциального уравнения
y’- 2 x = 0 (7.2)
решением является функция y = x 2 . Найдем производную y’= 2x и подставим в уравнение, получим: 2x – 2x = 0, 0 ≡ 0.

Следует заметить, что y = x 2 не единственное решение уравнения. Это уравнение имеет бесконечное множество решений, которые можно записать так: y = x 2 + C.

Дифференциальные уравнения первого порядка

Определение. Дифференциальным уравнением первого порядка называется уравнение, которое связывает независимую переменную x, искомую функцию y = f (x) и ее первую производную:
F (x, y, y’) = 0. (7.3)

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

Если уравнение (7.2) решить относительно у’, то оно будет иметь вид:
y’= f (x, y) или . (7.4)

Простые примеры показывают, что дифференциальное уравнение может иметь бесконечное множество решений. Это мы видим на примере уравнения (7.2). Легко убедиться также, что дифференциальное уравнение имеет решениями функции y = Cx, а дифференциальное уравнение — функции где C — произвольное число.

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

Определение. Общим решением дифференциального уравнения (7.3) называется функция
у = φ (х, С), (7.5)
которая зависит от одной произвольной постоянной и удовлетворяет дифференциальное уравнение при произвольном значении C.

Если функция (7.5) выражается неявно, то есть в виде
Ф (х, у, С) = 0, (7.6)
то (7.6) называется общим интегралом дифференциального уравнения.

Определение. Частным решением дифференциального уравнения (7.3) называется такое решение, которое получается из общего решения (7.5) при некотором конкретном значении постоянной C.

Ф (х, у, С₀) называется частным интегралом дифференциального уравнения.

На практике при решении конкретных задач часто приходится находить не все решения, а решение, которое удовлетворяет определенным начальным условиям. Одной из таких задач является задача Коши, которая для дифференциального уравнения первого порядка формулируется так: среди всех решений дифференциального уравнения (7.3) найти такое решение y, которое при заданном значении независимой переменной x = x₀ равна заданному значению y₀ , то есть y (x₀) = y₀ или (7.7)

Условие (7.7) называется начальным условием решения.

Покажем на примере, как найти частное решение дифференциального уравнения, когда известно общее решение и задано начальное условие.

Мы видим, что дифференциальное уравнение имеет общее решение y = Cx. Зададим начальное условие . Подставим эти значения в общее решение, получим 6 = 2С, откуда С = 3. Следовательно, функция y = 3x удовлетворяет и дифференциальное уравнение, и начальное условие.

Ответ на вопрос о том, при каких условиях уравнение (7.4) имеет
решение, дает теорема Коши.

ТЕОРЕМА (о существовании и единственности решения). Если функция f (x, y) и ее частная производная определены и непрерывные в области G, которая содержит точку M0 (x0; y0) , то существует единственное решение y = φ (x) уравнения (7.4), которое удовлетворяет начальному условию: y (x₀) = y₀.

Теорема Коши дает достаточные условия существования единого решения дифференциального уравнения (7.4). Заметим, что в условии теоремы не требуется существования частной производной .

График произвольного частного решения дифференциального уравнения называется интегральной кривой. Общему решению отвечает семья кривых. Так мы проверили, что уравнение имеет общее решение y = Cx, то ему соответствует семья прямых,
которые проходят через начало координат (рис. 1).

Уравнение имеет общее решение, ему соответствует семья равносторонних гипербол (рис. 2).

Если задано начальное условие то это означает, что задана точка M₀ (x₀;y₀), через которую должна проходить интегральная кривая, отвечающая искомому частному решению. Таким образом, отыскание частного решения дифференциального уравнения по заданному начальному условию геометрически означает, что из семьи
интегральных кривых мы выбираем проходящую через точку M₀ (x₀; y₀).

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

Общего метода решения дифференциальных уравнений первого порядка не существует. Рассмотрим некоторые методы решения отдельных типов дифференциальных уравнений.

Дифференциальные уравнения с разделенными переменными

Определение. Уравнение вида
f₁ (y) dy = f₂ (x) dx, (7.8)
где f₁ (y) и f₂ (x) — заданные функции, называется дифференциальным уравнением с разделенными переменными.

В этом уравнении каждая из переменных находится только в той части уравнения, где находится ее дифференциал. Уравнение dy = f (x) dx является частным случаем уравнения (7.8). Чтобы решить уравнение (7.8), надо проинтегрировать обе его части:
.

Понятно, что произвольную постоянную С можно записывать в любой части равенства.

Пример 1. Решить дифференциальное уравнение:
, удовлетворяющее начальному условию

Решение. Проинтегрируем левую и правую части уравнения, причем для удобства потенцирования, произвольную постоянную запишем в виде ln |C| получим:


— это общее решение дифференциального уравнения.
Подставляя в общее решение начальное условие, найдем С: 2 = С.
Итак,
является частным решением данного уравнения.

Дифференциальные уравнения с разделяющимися переменными

Определение. Уравнение вида
f1 (x) f2 (y) + g1 (x) g2 (y) = 0 (7.9)
называется дифференциальным уравнением с разделяющимися переменными.

В этом уравнении переменные еще не разделены, но, поделив обе части уравнения на произведение f₂ (y) g₁ (x), получим уравнение с разделенными переменными:

Интегрируя это уравнение, запишем
.

Получили общий интеграл данного уравнения.

Пример 2. Решить дифференциальное уравнение
x (y + 1) dx – (x 2 + 1) ydy = 0.

Решение. Поделим обе части этого уравнения на (y + 1) (x 2 + 1), после чего получим
.

Интегрируя, получим

— общий интеграл дифференциального уравнения.

Пример 3. Найти частное решение дифференциального уравнения (1 + x 2 ) dy + ydx = 0, удовлетворяющее начальному условию y (0) = 1.

Решение. Отделим переменные, поделив уравнение на y ⋅ (1 + x 2 ), и проинтегрируем данное уравнение:

Получили общий интеграл дифференциального уравнения.

Используя начальное условие, найдем произвольную постоянную С:
ln 1 + arctg 0 = C, откуда C = 0.

Найденную постоянную подставим в общий интеграл и отыщем частное решение:
откуда

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

Определение. Функция двух переменных f (x, y) называется однородной n- го измерения, если выполняется условие

Например, f (x, y) = x 2 + y 2 , f (tx, ty) = t 2 f (x 2 + y 2 ) — однородная функция второго измерения.

Определение. Дифференциальное уравнение
y ‘= f (x, y) (7.10)
называется однородным, если функция f (x, y) однородная нулевого измерения.

Покажем, что это уравнение можно свести к уравнению с разделенными переменными.
Рассмотрим функцию f (tx, ty). Сделаем замену будем иметь:

Тогда уравнение (7.10) запишется в виде (7.11)
В общем случае переменные в однородном уравнение не разделяются сразу. Но, если ввести вспомогательную неизвестную функцию u = u (x) по формуле
или y = xu, (7.12)
то мы сможем превратить однородное уравнение в уравнение с разделенными переменными.

Из формулы (7.12) найдем y’ = u + xu’ и уравнение примет вид: u + xu’ = φ (u),
то есть , откуда .

После интегрирования получим
Отсюда находим выражение для функции u, возвращаемся к переменной y = xu и получим решение однородного уравнения.

Чаще всего не удается найти функцию u явно выраженной, тогда, после интегрирования, в левую часть следует подставить вместо u.
В результате получим решение уравнения в неявном виде.

Пример 1. Найти решение однородного уравнения

Решение. Заменой y = xu сведем заданное уравнение к уравнению
или .

Отделяя переменные, найдем
откуда или , то есть
.
Возвращаясь к переменной y, получим общее решение: .

Линейные дифференциальные уравнения

Определение. Линейным дифференциальным уравнением первого порядка называется уравнение, которое содержит искомую функцию и ее производную в первой степени без их произведения:
y’ + P (x) y = Q (x). (7.13)

Здесь P (x), Q (x) — известные функции независимой переменной x. Например, y’ + 2 xy = x 2 .

Если Q (x) = 0, то уравнение (7.13) называется линейным однородным и является уравнением с разделяющимися переменными.

Если Q (x) ≠ 0, то уравнение (7.13) называется линейным неоднородным, которое можно решить несколькими способами.

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

Решение дифференциального уравнения (7.13) ищем в виде y = u (x) v (x) или y = uv, (7.14)
где u (x), v (x) — неизвестные функции. Одну из этих функций можно взять произвольную, а другая определяется из уравнения (7.13).

Из равенства y = uv найдем производную y’:
y’= u’ ⋅ v + u⋅ v’.

Подставим y и y’ в уравнение (7.13):
u’v + uv’ + P (x) ⋅ u⋅ v = Q (x) или u’v + u (v’ + P (x) ⋅ v) = Q (x).

Выберем функцию v такой, чтобы v’ + P (x) v = 0. (7.15)
Тогда для отыскания функции u получим уравнение:
u’v = Q (x). (7.16)

Сначала найдем v из уравнения (7.15).
Отделяя переменные, имеем , откуда

Под неопределенным интегралом здесь будем понимать какую-то одну первообразную от функции P (x), то есть v будет определенной функцией от x.

Зная v, находим u из уравнения (7.16):

откуда

Здесь мы уже берем для u все первообразные.

Найденные функции u и v подставляем в (7.14) и получаем общее решение линейного дифференциального уравнения:
(7.17)

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

Пример 1. Решить дифференциальное уравнение .
Решение. Решение ищем в виде y = uv, тогда y’= u’ ⋅ v + u⋅ v’.
Подставим y и y’ в уравнение: или
. (7.18)

Выражение, стоящее в скобках, приравниваем к нулю, имеем
или

Отделим переменные, домножив обе части уравнения на , тогда .
После интегрирования, получим ln |v| = ln |x| (здесь ограничимся одной первообразной), откуда v = x.
Подставим v = x в уравнение (7.18):

Общее решение запишется:
y = x (x + C) = x 2 + Cx.

Пример 2. Найти частное решение дифференциального уравнения который удовлетворяет начальному условию y (0) = 0.

Решение. Заданное уравнение — это линейное неоднородное уравнение первого порядка, решение которого ищем в виде y = u⋅v.
Тогда

Подставим v в уравнение и найдем u:

Общее решение дифференциального уравнения будет:

Подставляем начальные условия в найденное решение и находим С:

Из общего решения получаем частное решение
.

Дифференциальное уравнение Бернулли

Определение. Уравнения вида
(или )
называется дифференциальным уравнением Бернулли.

Данное уравнение отличается от уравнения (7.13) только множителем y» (или x») в правой части. Для того, чтобы права часть данного уравнения была такой, как в (7.13), разделим его левую и праву часть на y»:

Сделаем замену:
Домножим левую и правую части полученного уравнения на (n + 1) и, используя замену, получим:

Мы получили линейное дифференциальное уравнение относительно новой переменной

Пример 1. Найти общее решение дифференциального уравнения xy’ + y = y 2 ln x.

Решение. .
Сделаем замену Тогда

Данное уравнение решим, сделав замену z = u (x) ⋅ v (x).

Выбираем функцию v (x) так, чтобы выражение в скобках равнялось нулю, и эта функция была бы частным решением уравнения

Тогда .

Проинтегрировав правую часть этого уравнения по частям, получим , а при y -1 = z = uv, имеем

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

Обыкновенным дифференциальным уравнением называется любое соотношение, связывающее независимую переменную искомую функцию и производные искомой функции до некоторого порядка включительно.

Обыкновенное дифференциальное уравнение может быть приведено к виду

Здесь — известная функция, заданная в некоторой области

Число т. е. наивысший из порядков производных, входящих в (1), называется порядком уравнения.

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

По этой ссылке вы найдёте полный курс лекций по высшей математике:

Основные понятия и определения

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

В этой главе мы будем рассматривать уравнение, разрешенное относительно производной:

Наряду с этим уравнением мы всегда будем рассматривать перевернутое уравнение

используя последнее в окрестности тех точек, в которых обращается в бесконечность.

Во многих случаях оказывается целесообразным «место уравнении (2) и (2′) рассматривать одно равносильное им дифференциальное уравнение

Обе переменные и входят в это уравнение уже равноправно, и любую из них мы можем принять за независимую переменную.

Умножая обе части уравнения (3) на некоторую функцию получаем более симметричное уравнение:

где Обратно, всякое уравнение вида (4) можно переписать в виде уравнений (2) или (2′), разрешая его относительно или так что уравнение (4) равносильно следующим двум уравнениям:

Иногда уравнение записывают *з так называемой симметрической форме:

Возможно вам будут полезны данные страницы:

Решение уравнения. Предположим, что правая часть уравнения (2), определена на некотором подмножестве вещественной плоскости Функцию определенную в интервале мы будем называть решением уравнения (2) в этом интервале*, если:

1. Существует производная для всех значений из интервала (Отсюда следует, что решение представляет собою функцию, непрерывную ею всей области определения).
2. Функция обращает уравнение (2) в тождество:

справедливое для всех значений из интервала Это означает, что при любом из интервала точка принадлежит множеству и

Так как наряду с уравнением (2) рассматривается перевернутое уравнение (2′), то и решения этого перевернутого уравнения естественно присоединять к решениям уравнения (2).

В этом смысле в дальнейшем мы будем для краткости называть решения уравнения (2′) решениями уравнения (2).

Примеры с решением

Пример 1.

является решением уравнения

в интервале ибо она определена и дифференцируема в эгои интервале, и, подставляя се в уравнение (9), получаем тождество:

справедливое при всех значениях

Пример 2.

Функция есть решение равнения в интервале

Пример 3.

является решением уравнения

в интервале

Иногда функцию обращающую уравнение (2) в тождество (7), т. е. решение уравнения (2), называют интегралом этого уравнения. Мы будем употреблять термин интеграл только в смысле п. 16.

Системы обыкновенных дифференциальных уравнений

При решении многих задач нужно найти функции y₁ = y₁ (x), y₂ = y₂ (x), . y_n = y_n (x), которые удовлетворяют системе дифференциальных уравнений, содержащих независимую переменную x , искомые y₁ , y₂ , . y_n и их производные.

Пример. Пусть материальная точка массы m имеет криволинейную траекторию движения в пространстве. Определить положение точки в любой момент времени t, когда на нее действует сила .

Положение точки в любой момент времени t определяется ее координатами x, y, z; следовательно, x, y, z являются функциями от t. Проекциями вектора скорости точки на оси координат будут производные x’ , y’ , z’.
Положим, что сила, а соответственно и ее проекции F_x, F_y, F_z зависят от времени t, от положения x, y, z точки и от скорости движения точки, то есть от . Искомыми неизвестными функциями в этой задаче будут три функции x = x (t), y = y (t), z = z (t). Эти
функции определяются из уравнений динамики:

Мы получили систему трех дифференциальных уравнений второго порядка. В случае движения, когда траектория является плоской кривой, лежит, например, в плоскости Оxy, получим систему двух уравнений для определения неизвестных функций x (t) и y (t):

Рассмотрим простейшие системы дифференциальных уравнений.

Системы дифференциальных уравнений первого порядка

Система n уравнений первого порядка с n неизвестными функциями имеет вид:
(7.38)

где x — независимая переменная, y₁, y₂, . y_n — неизвестные функции.

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

Решением системы называется совокупность функций y₁, y₂, . y_n, которые превращают каждое уравнение системы в тождество относительно x.

Задача Коши для системы (7.38) состоит в нахождении функций y₁, y₂, . y_n , удовлетворяющих систему (7.38) и заданные начальные условия:
(7.39)

Интегрирование системы (7.38) делают следующим образом. Дифференцируем по x первое уравнение системы (7.38):

Заменим производные
их выражениями f₁, f₂, . f_n из уравнений системы (7.38), получим уравнение

Дифференцируем полученное уравнение и, подставив в это равенство значения производных из системы (7.38), найдем

Продолжая дальше таким образом, получим

В результате получаем следующую систему уравнений:
(7.40)

Из первых (n-1) уравнений определим y₂, y₃, . y_n:
(7.41)

и подставим их значения в последнее уравнение системы (7.40) для определения y₁:

Продифференцируем это выражение (n-1) раз, определим
как функции от x, C₁, C₂, . C_n. Подставим эти функции в (7.41), найдем
(7.43)

Для того, чтобы полученное решение удовлетворяло заданным начальным условиям, остается только найти значение произвольных постоянных из уравнений (7.42) и (7.43) так, как мы это делали для одного дифференциального уравнения.

Пример 1. Проинтегрировать систему

когда заданы начальные условия
Решение. Дифференцируем по x первое уравнение, имеем:
. Подставляем сюда значение и из системы, получим

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

Общим решением этого уравнения является
(*)
и тогда (**)

Подберем постоянные С₁ и С₂ так, чтобы выполнялись начальные условия. На основании (*) и (**) имеем:
1 = С₁ – 9; 0 = С₂ – 2С₁ + 14, откуда С₁ = 10, С₂ = 6.
Таким образом, решением системы, которое удовлетворяет заданным начальным условиям, будет:

Системы линейных дифференциальных уравнений с постоянными коэффициентами

Система дифференциальных уравнений:
(7.44)
где коэффициенты a_ij — постоянные числа, t — независимая переменная, x₁ (t), . x_n (t) —
неизвестные функции, называется системой линейных дифференциальных уравнений с постоянными коэффициентами.

Эту систему можно решать путем сведения к одному уравнению n-го порядка, как это было показано выше. Но эту систему можно решить и другим способом. Покажем, как это делается.

Будем искать решение системы (7.44) в виде:
(7.45)

Надо определить постоянные α₁, α₂, . α_n и k так, чтобы функции (7.45) удовлетворяли систему (7.44). Подставим функции (7.45) в систему (7.44):

Сократим на e kt и преобразуем систему, сведя ее к такой системе:
(7.46)

Это система линейных алгебраических уравнений относительно α₁, α₂, . α_n. Составим определитель системы:

Мы получим нетривиальные (ненулевые) решения (7.45) только при таких k, при которых определитель превратится в ноль. Получаем уравнение n-го порядка для определения k:

Это уравнение называется характеристическим уравнением для системы (7.44).

Рассмотрим отдельные случаи на примерах:

1) Корни характеристического уравнения действительны и различны. Решение системы записывается в виде:

Пример 2. Найти общее решение системы уравнений:

Решение. Составим характеристическое уравнение:
или k 2 – 5k + 4 = 0, корни которого k₁ = 1, k₂ = 4.

Решение системы ищем в виде

Составим систему (7.46) для корня k₁ и найдем и :
или

Откуда Положив получим
Итак, мы получили решение системы:

Далее составляем систему (7.46) для k = 4:

Откуда
Получим второй решение системы:
Общее решение системы будет:

2) Корни характеристического уравнения различны, но среди них есть комплексные:

k₁ = α + iβ, k₂ = α – iβ. Этим корням будут отвечать решения:

(7.47)

(7.48)

Можно доказать также, что истинные и мнимые части комплексного решения также будут решениями. Таким образом, получим два частных решения:
(7.49)
где — действительные числа, которые определяются через .

Соответствующие комбинации функций (7.49) войдут в общий решение системы.

Пример 3. Найти общее решение системы

Решение. Составляем характеристическое уравнение:
или k 2 + 12k + 37 = 0, корни которого k₁ = –6 + i, k₂ = –6 – i .

Подставляем поочередно k₁, k₂ в систему (7.46), найдем

Запишем уравнение (7.47) и (7.48) для наших данных

Перепишем эти решения в таком виде:

За частные решения можно взять отдельно действительные и отдельно мнимые части:

Общим решением системы будет

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

Классическая механика: о диффурах «на пальцах»

Введение

В этой статье я продолжаю тему цифровой обработки сигналов. В ней я постараюсь простым языком рассказать о концепции игровой механики (физики) с использованием подхода на основе дифференциальных уравнений. В будущем я собираюсь оценить действительно ли реализация такого подхода приведет к резкому увеличению вычислительной нагрузки. В рамках этой статьи не получится — слишком большой объем. В этой я собираюсь описать назначение коэффициентов, входящих в математическую модель динамического объекта, описать их физический смысл, т.е. их влияние на поведение динамического объекта.
Начнем пожалуй…

Физический смысл

В своей статье про Фильтр Калмана я описал что такое уравнения вида «вход-выход», передаточная функция и операторная форма записи дифференциальных уравнений (см. раздел «Основные понятия» в [2]). Примером части математической модели динамического объекта в операторной форме записи может служить следующее уравнение:
(1)
Это распространенная упрощенная модель динамических систем. Опережая напишу, что модель тела в движке «Farseer» использует урезанный аналог представленной выше модели (динамического звена второго порядка). Ниже представлено описание принятых в ней обозначений.

• а0, а1, а2 — коэффициенты инерции, демпфирования и жесткости, соответственно.
• b0, b1 — коэффициенты входного воздействия.
• s — оператор Лапласса (d/dt).
• α(s), β(s) — выходная и входная переменные, как функции оператора Лапласса.

Представленное уравнение описывает динамическую систему типа «один вход — один выход» (SISO). Ее можно использовать для описания динамики объекта по одной из его степеней свободы. Как Вам, возможно, известно, у свободного тела есть шесть степеней свободы — три поступательные (линейное движение вдоль трех осей системы координат (СК) ) и три вращательные (повороты вокруг осей СК). Таким образом, полная модель физического тела будет описываться шестью такими уравнениями (или четырьмя для 2D случая). Вы сразу можете сказать, что уже это свидетельствует о слишком высокой сложности такого подхода. Но на самом деле в Farseer, к примеру, класс тела (Body) содержит и линейные координаты тела (по сути это пара α(s) по OX и OY), и линейные скорости (пара s*α(s) по OX и OY) и параметры ориентации и угловой скорости. Эти параметры обсчитываются раздельно для каждой из осей, т.е. количество уравнений такое же — два уравнения по оси OX (линейное и угловое движение) и два для оси OY. Разница лишь в форме уравнений.
Алгоритм в движке Farseer — приближенный и упрощенный, но позволяет работать с варьируемыми квантами времени. В функцию обсчета параметров движения (Island.Solve(ref TimeStep step, ref Vector2 gravity)) передается время, прошедшее после последнего обсчета параметров. Это позволяет при недостаточной производительности компьютера держать скорость течения игрового времени примерно постоянным в ущерб плавности и реалистичности движения игровых объектов.
При построении дискретной модели на основе диференциальных уравнений мы четко завязываемся на фиксированный квант времени. Уравнения интегрируются для изначально заданной частоты дискретезации, и если с момента последнего обсчета по какой-то причине прошло времени больше, чем заданный квант (в англоязычной литературе его называют «time sample»), то мы либо должны произвести обсчет несколько раз, либо получим замедление движения объекта. Последнее я как раз и наблюдал в игре «Command Cortex» на слабой машине. Движения акторов были плавными но медленными (акторы, управляемые человеком, получают преимущество). Таким образом, нельзя говорить об исключительном преимуществе одного из этих подходов.
Теперь о том, за что отвечают коэффициенты представленного выше уравнения. Это уравнение описывает движение физического тела относительно положения равновесия при α(s) = 0. Это еще одна из причин кажущегося неудобства применения такой модели в игровой механике. При отсутствии приложенных внешних сил данная модель рано или поздно возвратит (при условии устойчивости модели) тело в положение равновесия. Представьте игровой мир, наполненный шариками, которые все время стремятся в начало координат (например, в левый верхний угол экрана). К такому поведению приводит наличие коэффициента жесткости (см. а2 выше). Представьте, что тело соединено с началом координат пружиной. Пока на тело действуют силы, пружина растянута, но стоит убрать внешнее воздействие и тело устремится к нулю. Тела в движке Farseer таким поведением не обладают. Если мы зададим коэффициент a2 равным нулю, то и в данном случае тела не будут стремиться к началу координат (см. выше я писал, что модели в Farseer по сути урезанные варианты этой модели). Ну и зачем этот коэффициент тогда нужен, спросите Вы. Если раскрыть скобки в левой части уравнения (1) и вместо слагаемого
a2*α(s)
напишем
a2*(α(s) — α0)
то через α0 мы получим возможность задать положение, к которому игровой объект будет стремиться. Величина коэффициента а2 отвечает за то, насколько быстро тело переместится в заданное положение равновесия (чем больше значение, тем выше жесткость пружины). Как такое реализовывается в Farseer я пока не выяснил, но думаю придется создать дополнительный источник воздействия.
Теперь коэффициент a1. Это коэффициент демпфирования. Чем больше значение этого коэффициента, тем быстрее гасится скорость (линейная или угловая). Аналогия из жизни — вязкие жидкости, такие как масло, мед, эпоксидная смола. Эти жидкости очень вязкие (имеют большое значение коэффициента демпфирования). Чем выше скорость движения тела в них, тем выше сопротивление этому движения. Если медленно двигать в них ложку, например, то преодолеть сопротивление большого труда не составит, а вот если ударить с размаху, то удар будет жестким.
Величина коэффициента а0 характеризует инерцию объекта. При описании линейного движения в качестве коэффициента а0 используется масса. Чем выше его значение, тем медленнее тело набирает скорость при приложении к нему внешних сил.
Теперь о коэффициентах в правой части уравнения (1). Тут нужно заметить, что данная модель расширенная на случай, когда входное воздействие определяется не только самим значением внешней силы, но и ее изменением. Для описания динамики игровых объектов это, возможно, будет излишним. Однако в промышленных системах управление встречаются и такие модели. Каков же их физический смысл? Коэффициент b1 это по сути коэффициент передачи внешней силы вовнутрь объекта. Обычно он равен еденице, т.е. сила передается как есть.
Коэффициент b0 интересен. Он играет роль форсирующего коэффициента. Представьте очень инерционный объект, к которому прикладывают силу, плавно нарастающую со временем. Если скорость нарастания и конечная величина силы будут малы, то объект очень медленно набирать скорость. Но если силу сделать большой, то после достижения внешней силой заданного значения объект не остановится в каком-то положении, а будет колебаться под действием инерции. Форсирование — это воздействие, пропорциональное скорости нарастания внешней силы. Если мы выберем его большим, то даже при малой скорости нарастания внешней силы наш объект будет достаточно быстро набирать скорость, а когда внешняя сила достигнет заданного значения, форсирование отключится. Вот такой вот хитрый этот «b0».

Динамика в картинках

Чтобы наглядно показать влияние коэффициентов дифф. уравнения на поведение динамического объекта решил построить графики переходного процесса при ступенчатом (step response) и импульсном (impulse response) входных воздействиях. Всего представлено 6 групп графиков (по одной группе для каждого коэффициента). Графики построены в пакете Octave (v. 3.4) с установленным пакетом «Signal Processing».
Итак, в качестве исходной возьмем модель вида:
=========================================
>>> w = tf([1 1],[1 1 1])

Transfer function «w» from input «u1» to output…

y1: (s + 1)/(s^2 + s + 1)

Continuous-time model.
=========================================


Код «w = tf([1 1],[1 1 1])» в символьном виде имеет вид:
>>> w = tf([b0 b1],[a0 a1 a2])
На скриншотах внизу-справа — примерное время стабилизации (коридором стабильности считаем ± 5% от заданной величины).

Попробуем поиграться с коэффициентом жесткости a2.
>>> w1 = 0.1*tf([1 1],[1 1 0.1])
y1: (s + 1)/(s^2 + s + 0.1)

>>> w2 = 10*tf([1 1],[1 1 10])
y1: (s + 1)/(s^2 + s + 10)
Примечание: пришлось подшаманить с коэффициентами усиления, чтобы результирующий коэффициент усиления был равен единице.


Что видно на графиках? Слева-направо представлены графики для w, w1 и w2, соответственно. Графики w1 более плавные и медленнее достигают установившегося значения. Графики w2 имеют более колебательный характер, но быстрее достигают установившегося значения. Вывод: жестче пружина — больше колебаний, но короче переходный процесс.

Попробуем поиграться с демпфированием (а1).
>>> w1 = tf([1 1],[1 0.25 1])
y1: (s + 1)/(s^2 + 0.25s + 1)

>>> w2 = tf([1 1],[1 2 1])
y1: (s + 1)/(s^2 + 2s + 1)


Сразу вывод: больше вязкость — быстрее затухают колебания.

Попробуем поиграться с инерцией (а0).
>>> w1 = tf([1 1],[0.1 1 1])
y1: (s + 1)/(0.1s^2 + s + 1)

>>> w2 = tf([1 1],[2 1 1])
y1: (s + 1)/(2s^2 + s + 1)


Вывод: меньше масса чугуняки — меньше болтанки и короче переходный процесс.

Перейдем к правой части и поиграемся с b1.
>>> w1 = 10*tf([1 0.1],[1 1 1])
y1: (10 s + 1)/(s^2 + s + 1)

>>> w2 = 0.25*tf([1 4],[1 1 1])
y1: (0.25 s + 1)/(s^2 + s + 1)


Вроде бы разница еле заметна, если смотреть на графики Step Response. Но на графиках Impulse Response хорошо виден эффект этого коэффициента. Если он равен единице, то график импульсного переходного процесса начинается с единицы (на самом деле он выходит из нуля, но не суть важно — второе значение в графике еденица). График w1 «начинается» со значения 10 (обратная величина от 0.1), а график w2 — начинается со значения 0.25 (обратное к 4). Таким образом, коэффициент b1 можно «обозвать» коэффициентом эффективности управления (входного воздействия).

И напоследок вкусненькое — игры с коэффициентом b0. Это хитрый коэффициент, потому и сравнение будет не таким как было выше. Чтобы показать его эффект придется варьровать несколько коэффициентов.
>>> w1 = tf([6 1],[1 1 1])
y1: (6 s + 1)/(s^2 + s + 1)

>>> w2 = tf([6 1],[1 3 1])
y1: (6 s + 1)/(s^2 + 3 s + 1)


Чем отличаются друг от друга w1 и w2? У w2 в три раза больше коэффициент демпфирования. В результате получаем интересные выводы. Графики w1 и w2 раньше пересекают уровень установившегося значения чем дефолтный график. Однако график w1 сохраняет форму дефолтного с его колебательностью, а график w2 за счет увеличенного демпфирования более сглаженный. Таким образом, играясь с форсированием и демпфированием мы можем заставить даже чугунный утюг порхать по рингу как бабочка без колебаний туда-сюда.

На правах PS

В данной статье я рассматривал лишь положительные значения коэффициентов. Их положительность — необходимое условие устойчивости мат. модели. Однако можно попробовать поиграться и с отрицательными значениями. Неустойчивой системой также можно управлять. Вспомните о самолетах пятого поколения (например, наш Беркут). Обратная стреловидность крыла — это неустойчивый планер, но зато высокая маневренность. Автоматика способна скорректировать эту неустойчивость и при этом, когда нужно, закладывать крутые виражи.
Если получится, состряпаю игрушку, с которой можно будет наглядно увидеть все эти эффекты.