Методика решения линейных обыкновенных дифференциальных уравнений
Уменьшение заказа это техника в математика для решения линейных обычный дифференциальные уравнения. Применяется, когда одно решение
известно и второй линейно независимый решение
желательно. Метод также применим к уравнениям n-го порядка. В этом случае анзац даст уравнение (n-1) -го порядка для
.
Линейные обыкновенные дифференциальные уравнения второго порядка
Пример
Рассмотрим общее однородное обыкновенное дифференциальное уравнение с линейными постоянными коэффициентами второго порядка. (ODE)
![ay '' (x) + через '(x) + cy (x) = 0, ;](https://wikimedia.org/api/rest_v1/media/math/render/svg/95fc4fa6a0bfed9411b1311a5b82cf4d993c1700)
куда
- действительные ненулевые коэффициенты. Два линейно независимых решения для этого ОДУ можно легко найти, используя характеристические уравнения за исключением случая, когда дискриминант,
, исчезает. В этом случае,
![{ displaystyle ay '' (x) + by '(x) + { frac {b ^ {2}} {4a}} y (x) = 0, ;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13984bf7fa751a3c7529fe92d1f7054f580a7abc)
из которых только одно решение,
![{ displaystyle y_ {1} (x) = e ^ {- { frac {b} {2a}} x},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4bc2de486bb26f707ef7a9a8c4af5ad5f6a8c2cf)
можно найти с помощью его характеристического уравнения.
Метод понижения порядка используется для получения второго линейно независимого решения этого дифференциального уравнения с использованием нашего одного известного решения. Чтобы найти второе решение, мы предполагаем
![у_ {2} (х) = v (х) у_ {1} (х) ;](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4a4dd9fed2215f00a85fb04e16e8b59164beaac)
куда
- неизвестная функция, которую предстоит определить. С
должен удовлетворять исходному ODE, мы подставляем его обратно, чтобы получить
![{ displaystyle a left (v''y_ {1} + 2v'y_ {1} '+ vy_ {1}' ' right) + b left (v'y_ {1} + vy_ {1}' вправо) + { frac {b ^ {2}} {4a}} vy_ {1} = 0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/286772aee7d28b3a0943e3efb2efeefa35db13b3)
Преобразуя это уравнение через производные от
мы получили
![{ displaystyle left (ay_ {1} right) v '' + left (2ay_ {1} '+ by_ {1} right) v' + left (ay_ {1} '' + by_ {1} '+ { frac {b ^ {2}} {4a}} y_ {1} right) v = 0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e0519ca16de6b33fd199df521823d15b39e6426)
Поскольку мы знаем, что
является решением исходной задачи, коэффициент при последнем члене равен нулю. Кроме того, подставив
в коэффициент второго члена дает (для этого коэффициента)
![2a left (- { frac {b} {2a}} e ^ {{- { frac {b} {2a}} x}} right) + be ^ {{- { frac {b} {2a }} x}} = left (-b + b right) e ^ {{- { frac {b} {2a}} x}} = 0.](https://wikimedia.org/api/rest_v1/media/math/render/svg/1bc18f61a859e4e32ed5b2cceea85eaa3dfa343c)
Следовательно, мы остаемся с
![ay_ {1} v '' = 0. ;](https://wikimedia.org/api/rest_v1/media/math/render/svg/bccc781f311e510f79840b9e20ebae182d13b5cc)
С
предполагается ненулевым и
является экспоненциальная функция (и, следовательно, всегда отличное от нуля), мы имеем
![v '' = 0. ;](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a6c8d2d1f242bdd0782bc4a72076aaf15f81131)
Его можно интегрировать дважды, чтобы получить
![v (x) = c_ {1} x + c_ {2} ;](https://wikimedia.org/api/rest_v1/media/math/render/svg/aca73549397e65d3b8d821f6b1af44bd9baf1f68)
куда
- константы интегрирования. Теперь мы можем записать наше второе решение как
![y_ {2} (x) = (c_ {1} x + c_ {2}) y_ {1} (x) = c_ {1} xy_ {1} (x) + c_ {2} y_ {1} (x ). ;](https://wikimedia.org/api/rest_v1/media/math/render/svg/6c808fc73c85c2e4aa3716fbb5a43ba2f7ffe46d)
Со второго срока в
является скалярным множителем первого решения (и, следовательно, линейно зависимым), мы можем отбросить этот член, получив окончательное решение
![y_ {2} (x) = xy_ {1} (x) = xe ^ {{- { frac {b} {2a}} x}}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/ddff71484668eb8e6f30c0129464b8276e1d3241)
Наконец, мы можем доказать, что второе решение
найденное этим методом, линейно не зависит от первого решения путем вычисления Вронскиан
![W (y_ {1}, y_ {2}) (x) = { begin {vmatrix} y_ {1} & xy_ {1} y_ {1} '& y_ {1} + xy_ {1}' end { vmatrix}} = y_ {1} (y_ {1} + xy_ {1} ') - xy_ {1} y_ {1}' = y_ {1} ^ {{2}} + xy_ {1} y_ {1} '-xy_ {1} y_ {1}' = y_ {1} ^ {{2}} = e ^ {{- { frac {b} {a}} x}} neq 0.](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8804e22e6b590ff02f5c8e8d8d712a964773e9b)
Таким образом
- второе искомое линейно независимое решение.
Общий метод
Учитывая общее неоднородное линейное дифференциальное уравнение
![у '' + p (t) y '+ q (t) y = r (t) ,](https://wikimedia.org/api/rest_v1/media/math/render/svg/787b0c619ee4e93b5d2107748c2992217b6369e3)
и единое решение
однородного уравнения [
], попробуем решить полное неоднородное уравнение в виде:
![y_ {2} = v (t) y_ {1} (t) ,](https://wikimedia.org/api/rest_v1/media/math/render/svg/fbb139fd9d5cb9a5dab4aa0a344e4127ab4bef59)
куда
- произвольная функция. Таким образом
![y_ {2} '= v' (t) y_ {1} (t) + v (t) y_ {1} '(t) ,](https://wikimedia.org/api/rest_v1/media/math/render/svg/1886c7be8de20b689cef0e272ed2c7cd5dac01e7)
и
![y_ {2} '' = v '' (t) y_ {1} (t) + 2v '(t) y_ {1}' (t) + v (t) y_ {1} '' (t). ,](https://wikimedia.org/api/rest_v1/media/math/render/svg/93bdd42d2adece4205be8a0b911ecd5c23d5113f)
Если они заменены на
,
, и
в дифференциальном уравнении, то
![y_ {1} (t) , v '' + (2y_ {1} '(t) + p (t) y_ {1} (t)) , v' + (y_ {1} '' (t) + p (t) y_ {1} '(t) + q (t) y_ {1} (t)) , v = r (t).](https://wikimedia.org/api/rest_v1/media/math/render/svg/c1544fe19da5a78be7f36856ce7fa32f9ede923d)
С
является решением исходного однородного дифференциального уравнения,
, поэтому мы можем сократить до
![y_ {1} (t) , v '' + (2y_ {1} '(t) + p (t) y_ {1} (t)) , v' = r (t)](https://wikimedia.org/api/rest_v1/media/math/render/svg/23a2a850ea943607959e260f40977a6f6a37252e)
которое является дифференциальным уравнением первого порядка для
(уменьшение порядка). Поделить на
, получение
.
В интегрирующий фактор является
.
Умножая дифференциальное уравнение на интегрирующий коэффициент
, уравнение для
можно свести к
.
После интегрирования последнего уравнения
найдено, содержащее одну постоянную интегрирования. Затем проинтегрируем
чтобы найти полное решение исходного неоднородного уравнения второго порядка, показывающее две постоянные интегрирования, как должно:
.
Смотрите также
Рекомендации