Асимптотическое интегрирование дифференциальных уравнений
Пусть имеем ряд (возможно и расходящийся)
 (53)
Обозначим через сумму первых членов ряда.
Будем говорить, что ряд (53) представляет собой асимптотическое разложение функции для достаточно больших , если выражение удовлетворяет условию
 (54)
( — любое фиксированное число), даже если ( — фиксировано). То обстоятельство, что данный ряд есть асимптотическое разложение функции (асимптотический степенной ряд), обозначается так:
Смысл асимптотического разложения состоит в том, что ряд (53) может служить источником приближенных формул
так что разность при будет бесконечно малой порядка выше , т.е.
Пример 11. Найти асимптотическое разложение функции
 (55)
Решение. Применяя раз интегрирование по частям, получаем
Обозначим и положим Имеем
так что ряд будет расходящимся для всех значений . Однако этот ряд может быть применен для вычисления значений при больших значениях . В самом деле, фиксируя некоторое значение , получим
отсюда, поскольку , будем иметь
 (56)
Для значений , достаточно больших, правая часть этого неравенства может быть сделана как угодной малой. Так, для будем иметь
поэтому значение функции может быть вычислено с большой точностью для больших значений , если взять сумму надлежащего числа членов ряда . Из оценки (56) следует, что
для всякого фиксированного , так что ряд дает асимптотическое разложение данной функции .
Если выполняется условие (54), то для коэффициентов ряда (53) из (54) получаем
 (57)
Отсюда следует, что если функция имеет асимптотическое разложение, то оно единственно.
Напротив, один и тот же ряд вида (53) может служить асимптотическим разложением для разных функций. Например, для функции
в силу (57) асимптотическим разложением является ряд (53), все коэффициенты , которого равны нулю. Очевидно, что этот же ряд является асимптотическим разложением и для функции . Говорят, что асимптотический ряд представляет не одну функцию, а класс асимптотически равных функций.
Операции над асимптотическими рядами
1) Если  (58) то  (59)
2) Если имеют место асимптотические разложения (58), то асимптотическое разложение функции может быть получено путем формального перемножения разложений (58).
3) Если функция имеет асимптотическое разложение
 (60)
начинающееся с члена , то имеет место асимптотическое разложение
 (61)
т.е. асимптотическое разложение (60) можно формально интегрировать почленно.
4) Формальное почленное дифференцирование асимптотического разложения, вообще говоря, недопустимо.
В самом деле, рассмотрим функцию ее асимптотическим разложением является ряд с коэффициентами , в то время как производная функции не имеет асимптотического разложения, поскольку даже не имеет предела при .
Однако если функция дифференцируема, а функция может быть разложена в асимптотический степенной ряд, то
Приложения к интегрированию дифференциальных уравнений
1. Рассмотрим дифференциальное уравнение
 (62) Ряд  (63)
расходящийся при всех значениях , формально удовлетворяет данному уравнению, в чем нетрудно убедиться непосредственной проверкой. Уравнению (62) удовлетворяет функция
причем интеграл в правой части находится при . Повторным интегрированием по частям находим где
При имеем
Следовательно, взяв первые членов ряда, мы совершим ошибку, меньшую (n+1)-го члена. Нетрудно видеть, что в данном случае
Поэтому построенный ряд является асимптотическим и может быть использован для вычисления интеграла и тем самым решения уравнения (62).
2. Если есть решение уравнения Бесселя , то, подставив вместо функцию , обнаружим, что будет удовлетворять уравнению
 (64)
Для больших это уравнение естественно попытаться заменить уравнением
 (65)
которое имеет решение .
Можно улучшить точность (для больших ) заменой постоянных и разложениями по отрицательным степеням :
Это означает, что решение уравнения (64) можно искать в виде
 (66)
Подставив выражение (66) в уравнение (64), получим
 (67)
Этот процесс можно продолжить и дальше. Существенно заметить, что эти выражения приводят к точному результату для (см. Бесселевы функции полуцелого индекса). Уравнение (65) является, как говорят, предельным уравнением для уравнения (64) (уравнение (65) получается из (64), если в коэффициенте при совершить предельный переход при ). Решение уравнения (65) для больших (особенно для ) достаточно хорошо определяет поведение решения исходного уравнения (64).
Примеры показывают, что асимптотическое поведение решений дифференциального уравнения не всегда можно вывести из поведения решений предельного уравнения.
Возьмем функцию  (68)
и сконструируем дифференциальное уравнение, для которого будет решением. Имеем
Попробуем, чтобы и были связаны соотношением
 (69)
Тогда заданная функция будет удовлетворять дифференциальному уравнению
![y''+ \left[\frac{\beta^2}{x^{2-2\beta}}+\frac{\alpha(\alpha-1)}{x^2}\right]y=0.](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAARsAAAA9BAMAAABvgdCNAAAALVBMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACttl6nAAAADnRSTlMAn7wBJoFZlEHSEOaRcLkE6kkAAAUiSURBVGje7ZrdaxtHEMBHt6dDPkugmKakhgbl4gR/NOHQCZPEIM7GKWncCtlUEONWXE1daFz8mqfC0ZpgQmqM0jy0oUU1FNoYgvFHcE1SjtAPSltjl9QPLTVNLIs0cfw3dFeSZcu53TvJK7cPngcx2rvx/W5mdmd2ZQCaaJoKnEXTtEpNvYpynDMNqlMUs0Jbz2ziAm/v1CmfG5XiXAP+osb44cjv6Lvm4Ygz/tHc/wmnB15wMpKp/tO54/RD0slonIrTwBdHXE4DKE5GJ5CtzzBDV5orzi9X58E342DTetRu9NzPeDgQ4onjMeA6dKQdbNpsH/beCsaR53jiXAW4Ltd7HVbdQbuygtToUbIi88T5FMTs13/1MMoR+cjl1rOljuBAHpUPjhjSPGPJjYfUu9tTJBizJH9Sd9J2ODGdH05gYOMbVuL47qp/6g0oi9UmfSLYvgOHpHF0kR+O0Lk8xKoQXUGIPheKPALwZ8CbOWTjnWicH85pgHCQcfNJHWqPGT6MEwuBlCXe6FeIGNXAQdj7QogxpZ6oEJ4C7B2UXIRTT0hT89UIEb0q3lnGawoDR8Y5LswAwt5ZMMG3nmcs9n7ccRoBsf4QwjgTq4DWAJI6BNbsZlbUcI8j3xWfsoJ1G88YxLDvS7cqq3F4CVMtRrrXXi4tEqNzVwBGy5joUjbwiIHja/qwo5ndWJ+42WRAr4oiyhFzsORehJMa1/O+MlZl/7wwycAJvF13yGGPgB+bW+o0pIL27DUVjrFxxOFtX2qCtYsMHMFlV+qJ0xfKeSaOrPRaW99u6bHCouvTbXBOU3Zc0s5XnKfiCHEmjieUqzAFeReWNrPAtMFppzwjvBPzAyrOErsbxOPJrSr0A57KDJxO1SWOREsN2WLjrJgwZOSjZsLlm5Bm4SC33nG5kxg/DpHSXF0w0WgunOFXpyG1bcqaZWz7KsVp6FNrxrAnRoqF5KQJeZyBmklxfW9xWuM34BaOXiSRE6OAk6tCwzFLym6GPpHoOZhI6EWcHy/ZC3mRLxOJvkTiTXLvJVfy/iaOr3M9lylILTaSBCefXiu65zEN5wBFrFKcA67kxWKwAhlYKO3mDm8GC1rU2mslwUJatXNHmIZPsKm4FaykCSu5S2gKYtYep3JbSCbpKhZTmfQLSV0kF6fgsLHHOLUhf7a0poVDmKONDA7e+zbthKMW90z4460/3ODIS4tUHOlOy+PS8TOp/mbSwWG4K1PAxpEPThWAiRIwB9zgeM6nqDjo1MTO3lK6kH/rmmlpe6sj2ZTQQHzCKii1QfQa9Bf35IywjBpDJg3nhjpEOyOMBbsMh+ZUmPRPbikWPA/bY0eRL/RRGg56IN6mdl6pRqdeWdaFsS3lVyXoLm370jTv/MQ4QNVctO4xs6iIF4fHXNGIM/TcKXtGluDIx+EsXq9eJ0pAh4yDKRqBc5STsN1ubH6r768Hj26QXZOm+nXDD/LfTkdyPZ8pIbleMrnjCGZLwyx8B4V1GysCeC0H0wb5Yduq8H2vzh1nGZJnn3o3NvLLBFFeub/kcAgmWr6sYHQoCvdgIQNazpiRwl4biHL/stO2QgNvxubQiYd3VPRg20QgioWcf3UJT7NrVuWpLGV3jFvOphejQdGoCs49f2bzhNq1eLNDcW9VcNA/E6uSVaal/+Om5u7qBOtIY6ouXaal/Mb5bqM6OBFVTJdtilRqFdglTgXVxUU3uNvTr32cfZx9HF44VfhZ//fEQqU4/80/PfwLDaee0EbB7LAAAAAASUVORK5CYII=) (70)
Предположим, что , например, . Тогда из условия (69) найдем и решение
![y(x)= \sqrt[\LARGE{4}]{x}\sin(\sqrt{x}+C)](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAMoAAAAaBAMAAAAEdd60AAAAMFBMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAlTPQ5AAAAD3RSTlMAgkHAm2gQ6CHQMVGwyAGUiHLCAAADtklEQVRIx+WW30sUURTHz+yvxt24tZlEvaxZpA9RiYmG4m6OWZS0K2b0lAMWIUS7KBkarGsF9tYPTAgqNvpFEKn98tElqHyI/oAgGLCQXnZxd21No+2c2bnrzv6qh+ilC7tz594753Pv95w5ZwD+szbxLyBfdv0tSyZtv4fyzI0UpnQWt8rCaNmzcr9eu/YEc9e2IMV8+xO2j6mBy3zGHOK9bpr+5NI/KEqJAz1nrqfvrf1ax9CfAzFvRkpPJ7VjqYeTAX5KmS+yq/N6BzL3wXDlTRjhq8HIj8WqcyjDAlLKMgYs57jxlvTQnYz5Kq7LIG4xDp/TJxxPb2pKyaa0IoWNZY6EtTXiIh8RMsSHddqGnai+uAC2tDxS2ra3PAti6KvrA/PavMHn570eJZdimifUN2DLeFZpAh6wURzY3nCiAWAu1x6e5asM96sFUnP9+9emGmVL7aXGABhDGbHWceDJjI7i/Ub/LmC7kXjXD3sYYm2t44MOBYwUtmKF2jS3bT8P7WCtP3V2OghCk9hrcyrDP1q7fbCK7wgFNR3fX39LR3HwPUgAW1f7WNwSBfgKkcMJBUxNpJJdbVrIVFSgFaM8d7UOt/8BmsGhsIQsLEGEU1DQdnA8i2dSmJP7aoMCVV6XbcGAlAuw0/MQ5fHncQBaMYPXhe+J8B2jDSkxEKIQ4XZQ0CpwK+TSSbt9t92+CSlD3MOoEHQFTYtiVL2jQBN+5qHMknbT6lPu5NgKhZ+lnfbep3ZXKEla361R3iklfkYUdpMCxUSRb81QrBJ/W1X7NG21vIjrKJY3kIpzcT4rxpxEKUtRWA1EXOT9e6/icB903qdsIwx5gNUDnD4+oODGjEF0JqdQPD5OYODgteNkFA7qKF34NhooAkrxtwccHgo2x8t5C1qby8pFJ6/4wIBj080xSqXGENvHHDKLKZ+jsBptHHVPwKwMLNG9ZA3pKCULAM/RCawRbx5JQwEoVaCyZWojntGbXUxMy2BEt4xIlRI6e1ZqextJzjuS0WRyl5kiZbxJFbStVpJkHcXiviaR49gCpYttA2qGsQQYKTQlZ1GsMYWssCB7SlqyHdhXGCaZcCrDGKMkKIiyRdb7BcTJbXSx4aoSP4VxsWzpLK8vVD/20rsVN2SJXJUlBk6vcT3CYxl8hTP/zhlXIQplfjZ0MVC0lFHmN9eoVeYIT3u5Nb4kVl7IgFrFukaL18sbWtHMqMhleUzFoGhFXrVYvNZ7/uiTgPUWnze7/sqXh1x8Ovx7C78AR18OJH/wmREAAAAASUVORK5CYII=) (71)
уравнения (70) при будет колеблющимся. С другой стороны,
поэтому предельным уравнением, соответствующим уравнению (70), будет
 (72) Его общее решение  (73) не содержит колеблющейся части.
Итак, асимптотическое поведение решения (71) дифференциального уравнения (70) нельзя "угадать" по поведению решения (73) предельного уравнения (72).
Приведем некоторые относящиеся сюда результаты. Пусть дано дифференциальное уравнение
 (74) где
 (75) так что
 (76)
Предельное уравнение в данном случае имеет вид
 (77)
и является уравнением с постоянными коэффициентами. Пусть — корни (которые мы для простоты предполагаем различными) характеристического уравнения
 (78)
Решения предельного уравнения — экспоненты .
Оказывается [16], что асимптотическое поведение решений уравнения (74) аналогично не поведению линейных комбинаций экспонент , а поведению линейных функций
 (79)
где показатели определяются формулами
 (80)
Функции (79) зависят не только от и , т.е. не только от предельных значений и при , но также и от коэффициентов , участвующих в правых частях равенств (75).
Теорема. Если характеристическое уравнение имеет различные корни и и если то уравнение
обладает линейно независимыми решениями и , которые можно представить асимптотическими рядами:
 (81)
Если корни характеристического уравнения совпадают, то может появиться логарифмический член. Решение можно представить асимптотическим рядом типа первого ряда (81), тогда как другое решение теперь представимо рядом вида
 (82)
Коэффициенты могут быть при этом найдены известным способом неопределенных коэффициентов путем подстановки выражений (81) или (82) в уравнение и приравнивания нулю коэффициентов при степенях . При этом формальное дифференцирование асимптотических разложений, законность которого априори неясна, приводит к правильным асимптотическим представлениям искомых функций.
Если заметили ошибку, опечатку или есть предложения, напишите в комментариях.
|