Автономная система (математика)

редактировать

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

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

Автономные системы тесно связаны с динамическими системами. Любая автономная система может быть преобразована в динамическую систему, и, используя очень слабые предположения, динамическая система может быть преобразована в автономную систему.

Содержание

  • 1 Определение
  • 2 Свойства
  • 3 Пример
  • 4 Качественный анализ
  • 5 Методы решения
    • 5.1 Первый порядок
    • 5.2 Второй порядок
      • 5.2.1 Специальные case: x '' = f (x)
      • 5.2.2 Особый случай: x '' = x 'f (x)
    • 5.3 Высшие порядки
  • 6 См. также
  • 7 Ссылки

Определение

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

ddtx (t) = f (x (t)) {\ displaystyle { \ frac {d} {dt}} x (t) = f (x (t))}\ frac {d} {dt} x (t) = f (x (t))

где x принимает значения в n-мерном евклидовом пространстве ; t часто интерпретируется как время.

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

ddtx (t) = g (x (t), t) {\ displaystyle {\ frac {d} {dt}} x (t) = g (x (t), t)}\ frac {d} {dt} x (t) = g (x (t), t)

, в котором закон, регулирующий эволюцию системы, зависит не только от текущего состояния системы, но и от параметра t, который также часто интерпретируется как время; такие системы по определению не являются автономными.

Свойства

Пусть x 1 (t) {\ displaystyle x_ {1} (t)}x _ {1} (t) будет уникальным решением начального значения проблема для автономной системы

ddtx (t) = f (x (t)), x (0) = x 0 {\ displaystyle {\ frac {d} {dt}} x (t) = f (Икс (т)) \, \ mathrm {,} \ quad x (0) = x_ {0}}\ frac {d} {dt} x (t) = f (x (t)) \, \ mathrm {,} \ quad x (0) = x_0 .

Тогда x 2 (t) = x 1 (t - t 0) {\ displaystyle x_ {2} (t) = x_ {1} (t-t_ {0})}x_2 (t) = x_1 (t-t_0) решает

ddtx (t) = f (x (t)), x (t 0) = x 0 {\ displaystyle {\ frac {d} {dt}} x (t) = f (x (t)) \, \ mathrm {,} \ quad x (t_ {0}) = x_ {0}}\ frac {d} {dt} x (t) = f (x (t)) \, \ mathrm {,} \ quad x (t_0) = x_0 .

Действительно, обозначая s = t - t 0 {\ displaystyle s = t-t_ {0}}s = t-t_0 , мы имеем x 1 (s) = x 2 (t) {\ displaystyle x_ {1} (s) = x_ {2} (t)}x_1 (s) = x_2 (t) и ds = dt {\ displaystyle ds = dt}ds = dt , таким образом,

ddtx 2 ( t) = ddtx 1 (t - t 0) = ddsx 1 (s) = f (x 1 (s)) = f (x 2 (t)) {\ displaystyle {\ frac {d} {dt}} x_ { 2} (t) = {\ frac {d} {dt}} x_ {1} (t-t_ {0}) = {\ frac {d} {ds}} x_ {1} (s) = f (x_ {1} (s)) = f (x_ {2} (t))}\ frac {d} {dt} x_2 (t) = \ frac {d} {dt} x_1 (t-t_0) = \ frac {d} {ds} x_1 (s) = f (x_1 (s)) = f (x_2 (t)) .

Для начального условия проверкой является tr ivial,

x 2 (t 0) = x 1 (t 0 - t 0) = x 1 (0) = x 0 {\ displaystyle x_ {2} (t_ {0}) = x_ {1} (t_ {0} -t_ {0}) = x_ {1} (0) = x_ {0}}x_2 (t_0) = x_1 (t_0-t_0) = x_1 (0) = x_0 .

Пример

Уравнение y ′ = (2 - y) y {\ displaystyle y '= (2-y) y}y'=(2-y)yявляется автономным, поскольку независимая переменная, назовем ее x {\ displaystyle x}x , явно не появляется в уравнение. Чтобы построить поле наклона и изоклина для этого уравнения, можно использовать следующий код в GNU Octave / MATLAB

Ffun = @ ( X, Y) (2 - Y). * Y; % function f (x, y) = (2-y) y [X, Y] = meshgrid (0:.2: 6, - 1:.2: 3); % выберите размеры графика DY = Ffun (X, Y); DX = единицы (размер (DY)); % генерировать колчан значений графика (X, Y, DX, DY, 'k'); % построить поле направления черным цветом; контур (X, Y, DY, [0 1 2], 'g'); % добавить изоклины (0 1 2) в зеленый заголовок ('Поле наклона и изоклины для f (x, y) = (2-y) y')

Из графика видно, что функция ( 2 - y) y {\ displaystyle (2-y) y}(2-y) y x {\ displaystyle x}x -инвариантно, как и форма решения, т. Е. y (x) = y (x - x 0) {\ displaystyle y (x) = y (x-x_ {0})}y (x) = y (x-x_0) для любого сдвига x 0 {\ displaystyle x_ {0}}x_ {0} .

Символьное решение уравнения в MATLAB путем выполнения

y = dsolve ('Dy = (2-y) * y', 'x'); % решая уравнение символически

, получаем два равновесных решения, y = 0 {\ displaystyle y = 0}y=0и y = 2 {\ displaystyle y = 2}y = 2 , и третье решение с неизвестной константой C 3 {\ displaystyle C_ {3}}C_ {3} ,

y (3) = - 2 / (exp (C3 - 2 * x) - 1)

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

y1 = dsolve ('Dy = (2-y) * y', 'у (1) = 1', 'х'); % решить задачу начального значения символически y2 = dsolve ('Dy = (2-y) * y', 'y (2) = 1', 'x'); % для разных начальных условий y3 = dsolve ('Dy = (2-y) * y', 'y (3) = 1', 'x'); у4 = dsolve ('Dy = (2-y) * y', 'y (1) = 3', 'x'); y5 = dsolve ('Dy = (2-y) * y', 'y (2) = 3', 'x'); y6 = dsolve ('Dy = (2-y) * y', 'y (3) = 3', 'x'); ezplot (y1, [0 6]); ezplot (y2, [0 6]); % построить решения ezplot (y3, [0 6]); ezplot (y4, [0 6]); ezplot (y5, [0 6]); ezplot (y6, [0 6]); title ('Поле уклона, изоклины и решения для f (x, y) = (2-y) y') legend ('Поле уклона', 'Изоклины', 'Решения y_ {1..6}'); текст ([1 2 3], [1 1 1], strcat ('\ leftarrow', {'y_1', 'y_2', 'y_3'})); текст ([1 2 3], [3 3 3], strcat ('\ leftarrow', {'y_4', 'y_5', 'y_6'})); сетка включена;
Поле наклона с изоклинами и решениями

Качественный анализ

Автономные системы могут быть проанализированы качественно с помощью фазового пространства ; в случае одной переменной это фазовая линия .

Методы решения

Следующие методы применяются к одномерным автономным дифференциальным уравнениям. Любое одномерное уравнение порядка n {\ displaystyle n}n эквивалентно n {\ displaystyle n}n -мерной системе первого порядка (как описано in Обычное дифференциальное уравнение # Редукция к системе первого порядка ), но не обязательно наоборот.

Первый порядок

Автономное уравнение первого порядка

dxdt = f (x) {\ displaystyle {\ frac {dx} {dt}} = f (x)}\ frac {dx} {dt} = f (x)

является разделимым, поэтому его легко решить, преобразовав его в интегральную форму

t + C = ∫ dxf (x) {\ displaystyle t + C = \ int {\ frac {dx} {f (x)}}}t + C = \ int \ frac {dx} {f (x)}

Второй порядок

Автономное уравнение второго порядка

d 2 xdt 2 = f (x, x ′) {\ displaystyle {\ frac {d ^ { 2} x} {dt ^ {2}}} = f (x, x ')}\frac{d^2x}{dt^2} = f(x, x')

сложнее, но его можно решить, введя новую переменную

v = dxdt {\ displaystyle v = { \ frac {dx} {dt}}}v = {\ frac {dx} {dt}}

и выражая вторую производную от x {\ displaystyle x}x (через правило цепочки ) как

d 2 xdt 2 = dvdt = dxdtdvdx = vdvdx {\ displaystyle {\ frac {d ^ {2} x} {dt ^ {2}}} = {\ frac {dv} {dt}} = {\ frac {dx} {dt}} {\ frac {dv} {dx}} = v {\ frac {dv} {dx}}}\ frac {d ^ 2x} {dt ^ 2} = \ frac {dv} {dt} = \ frac {dx} {dt} \ frac {dv} {dx} = v \ frac {dv} {dx}

, так что исходное уравнение принимает вид

vdvdx = f (x, v) {\ displaystyle v {\ frac {dv} {dx}} = f (x, v)}v \ frac {dv} {dx} = f (x, v)

который является уравнением первого порядка uation не содержит ссылки на независимую переменную t {\ displaystyle t}tи, если решено, предоставляет v {\ displaystyle v}v как функцию от x {\ Displaystyle x}x . Затем, вспоминая определение v {\ displaystyle v}v :

dxdt = v (x) ⇒ t + C = ∫ dxv (x) {\ displaystyle {\ frac {dx} {dt}} = v ( x) \ quad \ Rightarrow \ quad t + C = \ int {\ frac {dx} {v (x)}}}\ frac {dx} {dt} = v (x) \ quad \ Rightarrow \ quad t + C = \ int \ frac {dx} {v (x)}

, которое является неявным решением.

Особый случай: x '' = f (x)

Особый случай, когда f {\ displaystyle f}f не зависит от x ′ {\ displaystyle x '}x'

d 2 xdt 2 = f (x) {\ displaystyle {\ frac {d ^ {2} x} {dt ^ {2}}} = f (x)}\ frac {d ^ 2 x} {dt ^ 2} = f (x)

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

Идея состоит в том, чтобы использовать тождество (исключая деление на ноль )

dxdt = (dtdx) - 1 {\ displaystyle {\ frac {dx} {dt}} = \ left ({\ frac {dt} {dx}} \ right) ^ {- 1}}\ frac { dx} {dt} = \ left (\ frac {dt} {dx} \ right) ^ {- 1}

который следует из правила цепочки . Отметим в стороне, что, инвертируя обе стороны автономной системы первого порядка, можно немедленно интегрировать по x {\ displaystyle x}x :

dxdt = f (x) ⇒ dtdx = 1 f (x) ⇒ t + С = ∫ dxf (x) {\ displaystyle {\ frac {dx} {dt}} = f (x) \ quad \ Rightarrow \ quad {\ frac {dt} {dx}} = {\ frac {1} { f (x)}} \ quad \ Rightarrow \ quad t + C = \ int {\ frac {dx} {f (x)}}}\ frac {dx} {dt} = f (x) \ quad \ Rightarrow \ quad \ frac {dt} {dx} = \ frac {1} {f (x)} \ quad \ Rightarrow \ quad t + C = \ int \ frac {dx} {f (x)}

, что является еще одним способом взглянуть на технику разделения переменных. Тогда возникает естественный вопрос: можем ли мы сделать что-то подобное с уравнениями более высокого порядка? Ответ положительный для уравнений второго порядка, но есть над чем поработать. Вторая производная должна быть выражена как производная по отношению к x {\ displaystyle x}x вместо t {\ displaystyle t}t:

d 2 xdt 2 = ddt (dxdt) = ddx (dxdt) dxdt = ddx ((dtdx) - 1) (dtdx) - 1 = = - (dtdx) - 2 d 2 tdx 2 (dtdx) - 1 = - (dtdx) - 3 d 2 tdx 2 = = ddx (1 2 (dtdx) - 2) {\ displaystyle {\ begin {align} {\ frac {d ^ {2} x} {dt ^ {2}}} = {\ frac {d} {dt}} \ left ({\ frac {dx} {dt}} \ right) = {\ frac {d} {dx}} \ left ({\ frac {dx} {dt}} \ right) {\ frac {dx} { dt}} = {\ frac {d} {dx}} \ left (\ left ({\ frac {dt} {dx}} \ right) ^ {- 1} \ right) \ left ({\ frac {dt} {dx}} \ right) ^ {- 1} = \\ = - \ left ({\ frac {dt} {dx}} \ right) ^ {- 2} {\ frac {d ^ {2} t} {dx ^ {2}}} \ left ({\ frac {dt} {dx}} \ right) ^ {- 1} = - \ left ({\ frac {dt} {dx}} \ right) ^ {- 3} {\ frac {d ^ {2} t} {dx ^ {2}}} = \\ = {\ frac {d} {dx}} \ left ({\ frac {1} {2}} \ left ({\ frac {dt} {dx}} \ right) ^ {- 2} \ right) \ end {align}}}\ begin {align} \ frac {d ^ 2 x} {dt ^ 2} = \ frac {d} {dt} \ left (\ frac {dx} {dt} \ right) = \ frac {d} {dx} \ left (\ frac {dx} {dt} \ right) \ frac {dx } {dt} = \ frac {d} {dx} \ left (\ left (\ frac {dt} {dx} \ right) ^ {- 1} \ right) \ left (\ frac {dt} {dx} \ справа) ^ {- 1} = \\ = - \ left (\ frac {dt} {dx} \ right) ^ {- 2} \ frac {d ^ 2 t} {dx ^ 2} \ left (\ frac {dt} {dx} \ right) ^ {- 1} = - \ left (\ frac {dt} {dx} \ right) ^ {- 3} \ frac {d ^ 2 t} {dx ^ 2} = \ \ = \ frac {d} {dx} \ left (\ frac {1} {2} \ left (\ frac {dt} {dx} \ right) ^ {- 2} \ right) \ end {align}

Чтобы еще раз подчеркнуть: достигнуто то, что вторая производная в t { \ displaystyle t}tесть пчела n, выраженное как производная в x {\ displaystyle x}x . Затем исходное уравнение второго порядка может быть окончательно интегрировано:

d 2 xdt 2 = f (x) {\ displaystyle {\ frac {d ^ {2} x} {dt ^ {2}}} = f (x) }\ frac {d ^ 2 x} {dt ^ 2} = f (x)
ddx (1 2 (dtdx) - 2) = f (x) {\ displaystyle {\ frac {d} {dx}} \ left ({\ frac {1} {2}} \ left ({\ гидроразрыва {dt} {dx}} \ right) ^ {- 2} \ right) = f (x)}\ frac {d} {dx} \ left (\ frac {1} {2} \ left (\ frac {dt} {dx} \ right) ^ {- 2} \ right) = f (x)
(dtdx) - 2 = 2 ∫ f (x) dx + C 1 {\ displaystyle \ left ( {\ frac {dt} {dx}} \ right) ^ {- 2} = 2 \ int f (x) dx + C_ {1}}\ left (\ frac {dt} {dx} \ right) ^ {- 2} = 2 \ int f (x) dx + C_1
dtdx = ± 1 2 ∫ f (x) dx + C 1 {\ displaystyle {\ frac {dt} {dx}} = \ pm {\ frac {1} {\ sqrt {2 \ int f (x) dx + C_ {1}}}}}\ frac {dt} {dx} = \ pm \ frac {1} {\ sqrt {2 \ int f (x) dx + C_1}}
t + C 2 Знак равно ± ∫ dx 2 ∫ е (Икс) dx + C 1 {\ displaystyle t + C_ {2} = \ pm \ int {\ frac {dx} {\ sqrt {2 \ int f (x) dx + C_ {1 }}}}}t + C_2 = \ pm \ int \ frac {dx} {\ sqrt {2 \ int f (x) dx + C_1}}

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

Особый случай: x '' = x 'f (x)

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

d 2 xdt 2 = (dxdt) nf (x) {\ displaystyle {\ frac {d ^ {2} x} {dt ^ {2}}} = \ left ({\ frac {dx} {dt}} \ right) ^ {n} f ( x)}\ frac {d ^ 2 x} {dt ^ 2} = \ left (\ frac {dx} {dt} \ right) ^ nf (x)

, где n {\ displaystyle n}n - некоторый параметр, не равный двум. Это будет работать, поскольку вторая производная может быть записана в форме, включающей степень x ′ {\ displaystyle x '}x'. Переписываем вторую производную, переставляем и выражаем левую часть как производную:

- (dtdx) - 3 d 2 tdx 2 = (dtdx) - nf (x) {\ displaystyle - \ left ({\ frac {dt } {dx}} \ right) ^ {- 3} {\ frac {d ^ {2} t} {dx ^ {2}}} = \ left ({\ frac {dt} {dx}} \ right) ^ {-n} f (x)}- \ left (\ frac {dt} {dx} \ right) ^ {- 3} \ frac {d ^ 2 t} {dx ^ 2} = \ left (\ frac {dt} {dx} \ right) ^ {- n} f (x)
- (dtdx) n - 3 d 2 tdx 2 = f (x) {\ displaystyle - \ left ({\ frac {dt} {dx}} \ right) ^ { n-3} {\ frac {d ^ {2} t} {dx ^ {2}}} = f (x)}- \ left (\ frac {dt} {dx} \ right) ^ {n - 3} \ frac {d ^ 2 t} {dx ^ 2} = f (x)
ddx (1 2 - n (dtdx) n - 2) = f (x) {\ displaystyle {\ frac {d} {dx}} \ left ({\ frac {1} {2-n}} \ left ({\ frac {dt} {dx}} \ right) ^ {n-2} \ справа) = е (х)}\ frac {d} {dx} \ left (\ frac {1} {2 - n} \ left (\ frac {dt} {dx} \ справа) ^ {n - 2} \ right) = f (x)
(dtdx) n - 2 = (2 - n) ∫ е (x) dx + C 1 {\ displaystyle \ left ({\ frac {dt} {dx}} \ справа) ^ {n-2} = (2-n) \ int f (x) dx + C_ {1}}\ left (\ frac {dt} {dx} \ right) ^ {n - 2} = (2 - n) \ int f (x) dx + C_1
t + C 2 = ∫ ((2 - n) ∫ f (x) dx + C 1) 1 n - 2 dx {\ displaystyle t + C_ {2} = \ int \ left ((2-n) \ int f (x) dx + C_ {1} \ right) ^ {\ frac {1} { n-2}} dx}t + C_2 = \ int \ left ((2 - n) \ int f (x) dx + C_1 \ right) ^ {\ frac {1} {n - 2}} dx

Справа будет +/-, если n {\ displaystyle n}n четное. Обработка должна быть другой, если n = 2 {\ displaystyle n = 2}n = 2 :

- (dtdx) - 1 d 2 tdx 2 = f (x) {\ displaystyle - \ left ({\ frac {dt} {dx}} \ right) ^ {- 1} {\ frac {d ^ {2} t} {dx ^ {2}}} = f (x)}- \ left (\ frac {dt} {dx} \ right) ^ {- 1} \ frac {d ^ 2 t} {dx ^ 2} = f (x)
- ddx (ln ⁡ (dtdx)) = е (Икс) {\ Displaystyle - {\ гидроразрыва {d} {dx}} \ влево (\ пер \ влево ({\ гидроразрыва {dt} {dx}} \ вправо) \ вправо) = е (х)}- \ frac {d } {dx} \ left (\ ln \ left (\ frac {dt} {dx} \ right) \ right) = f (x)
dtdx = C 1 e - ∫ е (x) dx {\ displaystyle {\ frac {dt} {dx}} = C_ {1} e ^ {- \ int f (x) dx}}\ frac {dt} {dx} = C_1 e ^ {- \ int f (x) dx}
t + C 2 знак равно C 1 ∫ е - ∫ е (x) dxdx {\ displaystyle t + C_ {2} = C_ {1} \ int e ^ {- \ int f (x) dx} dx}t + C_2 = C_1 \ int e ^ {- \ int f (x) dx} dx

высшие порядки

Аналогового метода решения автономных уравнений третьего и более высокого порядка не существует. Такие уравнения могут быть решены точно только в том случае, если они обладают каким-либо другим упрощающим свойством, например, линейностью или зависимостью правой части уравнения только от зависимой переменной (т. Е. Не от ее производных). Это не должно вызывать удивления, учитывая, что нелинейные автономные системы в трех измерениях могут производить поистине хаотическое поведение, такое как аттрактор Лоренца и аттрактор Ресслера.

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

См. Также

Ссылки

Последняя правка сделана 2021-06-12 19:28:03
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте