130
Вестник СамГУ. 2015. № 6(128)
УДК 517.923
И.В. Филимонова, Т.С. Хачлаев1
ОБ АСИМПТОТИЧЕСКИХ СВОЙСТВАХ ОПРЕДЕЛЕННЫХ НА ПОЛУОСИ РЕШЕНИЙ ОДНОГО ПОЛУЛИНЕЙНОГО ОДУ2
В статье рассматриваются решения обыкновенного дифференциального полулинейного уравнения, коэффициенты которого зависят от нескольких вещественных параметров. Если коэффициент выбрать так, что уравнение не будет содержать производной первого порядка от неизвестной функции, то это будет случай уравнения Эмдена — Фаулера. Асимптотическое поведение решений уравнения Эмдена — Фаулера при неограниченно больших значениях переменного описано в книге Ричарда Беллмана. Рассматриваемые в статье уравнения, содержащие первую производную от неизвестной функции, встречаются в некоторых задачах для эллиптических уравнений с частными производными в неограниченных областях. От того, с каким знаком первая производная входит в уравнение, существенно зависит описание решений. Частично результат этой статьи может быть получен из работ И.Т. Кигурадзе. Для описания асимптотического поведения решений нелинейного уравнения используются леммы о поведении решений линейных уравнений с достаточно сильно (слабо) растущим потенциалом.
Ключевые слова: обыкновенные дифференциальные уравнения, нелинейные уравнения, полулинейные уравнения, уравнения Эмдена — Фауле-ра, асимптотика решений, положительные решения, существование решений, принцип максимума.
В работе рассматриваются решения уравнения
где а > 1, к,р — произвольные постоянные. Изучается асимптотическое поведение при х ^ решений этого уравнения, определенных при х > хо > 0. Знак
параметра к, и величина параметра р существенным образом влияют на поведение решений.
Так как коэффициент перед нелинейным членом в уравнении (1) отрицательный, то решения не могут иметь положительных максимумов и отрицательных минимумов, значит, при больших х решение уравнения (1) есть монотонная функция. Поэтому можно ограничиться рассмотрением положительных решений, и нелинейный член писать коротко хрха.
х© Филимонова И.В., Хачлаев Т.С., 2015
Филимонова Ирина Владимировна ([email protected]), кафедра дифференциальных уравнений, МГУ имени М.В. Ломоносова, 119991, Российская Федерация, г. Москва, Ленинские горы, 1.
Хачлаев Тимур Султанович ([email protected]), кафедра высшей математики, МИРЭА, 119454, Российская Федерация, г. Москва, пр. Вернадского, 78.
2Статья подготовлена по докладу конференции "СамДиф-2015".
х" + кг' - хр\х\а-1х = 0,
(1)
Случай к = 0 полностью рассмотрен в книге [4], в главе 7 посвященной асимптотическим свойствам решений уравнения Эмдена — Фаулера. Для случая к = 0 результат нашей работы сформулирован в теореме 1.
Отметим, что при к < 0 уравнение (1) заменой в = екх преобразуется в: х88 = = (1п в)рв-2ха. Решениям, определенным при больших х, будут отвечать решения, определенные при больших в. Решения уравнения х88 = а(в)ха рассматривались в работе [6], но только один наш результат для случая к < 0, р < — 1 может быть получен как следствие теорем работы [6], а именно это теорема 4 [6] при а(в) = (1п в)Рв-2.
Некоторые свойства решений уравнения (1) при к = 0 были приведены в работе [1]. В случае р = 0, к < 0 в работе [2] получено асимптотическое разложение решения в ряд. В этих работах рассмотрение уравнения (1) носило вспомогательный характер для изучения решений полулинейного дивергентого эллиптического уравнения в цилиндрических областях [1] и решений недивергентного эллиптического уравнения во внешней к компакту области [2].
Асимптотическое поведение решений уравнений с частными производными не всегда соответствует поведению решений предполагаемого обыкновенного уравнения [3]. По совету проф. И.В. Асташовой мы изучаем свойства решений уравнения (1) и надеемся это поможет и в некоторых других задачах из уравнений с частными производными.
Нетривиальное решение х(х) уравнения (1), определенное при больших х, не может очень быстро расти.
Лемма 1. Существует М(р, а), что любое решение х(х) уравнения (1), определенное при х > хо для всех х > х\, оценивается как
\х\ < Мх-^т .
Приведем схему доказательства леммы 1. Сначала доказывается, что уравнение
х'' + кг' — ха =0 (2)
имеет решение хо (х) с двумя асимптотами. Пусть это х = 0, х = I. Введем обозначения Мо, ¡о: Мо = штхе(о,г) хо(х) = х(1о). Тогда для решения х(х), уравнения (2), определенного на отрезке [0,1], будет выполнено х(1о) ^ Мо. Отсюда выводится,
что если на отрезке [0,1] определено решение у(х) уравнения у'' + ку' — Б(х)уа =
_ р
= 0, где Б(х) > Во, то выполнено у (¡о) < МоБо . Оценка решения х(х) уравнения (1), определенного при х > хо, в точке х такой, что х > хо + ¡о, получится, если рассмотреть х(х) как решение линейного уравнения на отрезке [х — ¡о,х + + I — ¡о].
Сформулируем основной результат работы:
Теорема 1. Нетривиальное решение х(х) уравнения (1), определенное при больших х, имеет при х ^ то вид:
если к < 0, р> —1, то х(х) = ±(к—)—х-(1+р)/(1-а)(1 + о(1)),
если к < 0, р = —1, то х(х) = ±( —)— 1п-1/(а-1) х(1 + о(1)),
если к > 0 , р > —1, то х(х) = С(е-кх + 0(е-ах)), где С = 0, а € (к, ка) можно
выбрать произвольно,
если к > 0, р < —1, то есть три вида решений :
х(х) = С(1 + 0(хр+1)), где С = 0,
х(х) = С(е-кх + 0(е-ах)), где С = 0, а € (к,ка) можно выбрать произвольно, х(х) = ±( )1/(-1)х(Р+1)/(1-)(1 + о(1)),
если k < 0 , p < -1, то z(x) = C(1 + O(xp+1)), где C = 0. Более того, такого вида решения действительно существуют.
Асимптотические формулы для решений, стремящихся к ненулевой константе и экспоненциально убывающих, получаются из формул для представления решения линейного уравнения с асимптотически мало отличающимися от постоянных коэффициентами.
Доказательство формул в случае степенной асимптотики требует более тонких методов.
Лемма 2. Пусть Q(x) — ж при x — ж. Тогда для z(x) решения уравнения
z" + kz' - Q^z = 0 (3)
x
такого, что при x — ж выполнено z(x) — ж, для любых в = const > 0, C = = const > 0 найдется xi, что при x > x\ выполняется неравенство z(x) > Cxe.
Лемма 3. Пусть Q(x) — 0 при x — ж, Q(x) > 0, k > 0. Тогда для z(x) решения уравнения (3), для любых в = const > 0, C = const > 0 найдется xi, что при x > xi выполняется неравенство z(x) < Cxe.
Замечание. Доказать отсутствие неограниченного решения в условиях леммы 3 невозможно, так как функция z(x) = ln x имеет решение уравнения (3) с функцией Q(x) = (k — x-1) ln-1 x.
В данной статье обоснуем асимптотические формулы в наиболее сложном случае k > 0, p < —1.
< Так как z(x) монотонна при x > xo, то есть три возможности —^^о z(x) = = ж, li^^x—^^о z(x) = C > 0, li^^x—^^о z(x) = 0.
В двух последних случаях функция z(x) ограничена при x > xo, поэтому
z'' + kz' — a(x)z = 0, где a(x) = xpza-i = O(xp) при x — ж, так как p < —1, то
z(x) = C1(1 + O(xp+1)) + C2e-kx(1 + O(xp+1)) при x — ж,
т. е. либо z — C1 =0, либо z ~ C2e-kx при x — ж.
Пусть теперь limx—TO z(x) = ж. Покажем, что уравнение (1) имеет верхнее и нижнее решение вида Cxa, а = (p + 1)/(1 — v). Действительно подставляя это выражение в левую часть уравнения (1), найдем
1 а( а — 1)
Cx«-1( ka — C -1)
Поэтому, если Cа 1 > ka, то при x > max(0, Са— ) эта функция верхнее
€°-1-ка>
решение, если Са-1 < ка при х > шах(0, кО^-с"--1) эта функция нижнее решение.
Покажем теперь, что если г(х) ^ то при х ^ то, то г(х)/ха = К(х) имеет конечный ненулевой предел.
Допустим, что предела не существует, тогда найдутся числа С > С2, что К(х) принимает значения большие С и меньшие С2 на любом интервале (т, то). Выберем из интервала (С\,С2) число С так, чтобы Са-1 = ка. Тогда как отмечено выше, функция Сха на интервале (х1, то) будет верхним (или нижним) решением уравнения (1), график которого пересекается с графиком рассматриваемого решения г(х), бесконечное число раз на интервале (х1, то), что противоречит определению верхнего (нижнего) решения.
Поэтому отношение R(x) имеет предел. Допустим, что это бесконечный предел. Отметим, что xpz7-1 = R7-1x-1. Поэтому по лемме 2 z(x) растет быстрее хв при любом в, что противоречит оценке z < леммы 1.
Допустим теперь, что предел R(x) при x ^ ж равен нулю. Тогда xpz7-1 = = R7-1x-1. Поэтому по лемме 3 z(x) растет медленнее xe при любом в > 0. Выбрав в < а, найдем, что xpz7-1 = O(x-1-(a-e)(7-1)) при x ^ ж. И тогда получим, что
z(x) = С1(1 + о(1)) + C2e-kx(1 + о(1)) при x ^ ж,
что противоречит неограниченности функции z(x).
Таким образом, доказано, что R(x) = z(x)/xa ^ C = 0 при x ^ ж. Ясно, что С7-1 = ka. □
Существование экспоненциально убывающих решений проводится с помощью теоремы Лерэ — Шаудера о существовании неподвижной точки отображения выпуклого множества в свое предкомпактное подмножество. Построение решений, стремящихся к ненулевой константе, проводится с помощью итерационного процесса, подобно работе [5].
Литература
[1] Хачлаев Т.С. Асимптотическое поведение решений полулинейного эллиптического уравнения с растущим коэффициентом в цилиндрической области // Успехи математических наук. 2004. Т. 59. Вып. 2. С. 185-186.
[2] Сурначев М.Д. Асимптотическое поведение положительных решений уравнений Эм-дена — Фаулера // Проблемы математического анализа. 2011. № 1. С. 129-176.
[3] Филимонова И.В. О поведении решений полулинейного параболического или эллиптического уравнения, удовлетворяющих нелинейному краевому условию, в цилиндрической области // Труды семинара им. И.Г. Петровского. 2007. Вып. 26. С. 369-390.
[4] Беллман Р. Теория устойчивости решений дифференциальных уравнений. М.: Иностр. лит., 1954.
[5] Atkinson F.V. On second-order non-linear oscillations // Pacific Journal of Mathematics. 1955. V. 5(1). P. 643-647.
[6] Кигурадзе И.Т. Асимптотические свойства решений одного нелинейного дифференциального уравнения типа Эмдена — Фаулера // Изв. АН СССР. Сер.: Математика. 1965. Т. 29. Вып. 5. С. 965-986.
References
[1] Khachlaev T.S. Asymptotic behavior of solutions of a semilinear elliptic equation with increasing coefficient in a cylindrical domain. Uspekhi matematicheskikh nauk [Russian Mathematical Surveys], 2004, Vol. 59. Issue 2, pp. 185-186 [in Russian].
[2] Surnachev M.D. Asymptotic behavior of positive solutions to Emden-Fowler type equations. [Problems of Mathematical Analysis], 2011, no. 1, pp. 129-176 [in Russian].
[3] Filimonova I.V. On the behavior of solutions of a semilinear parabolic or elliptic equation satisfying a nonlinear boundary condition in a cylindrical domain. Trudy seminara im. I.G. Petrovskogo [Proceedings of the seminar named after I.G. Petrovsky], 2007, Issue 26, pp. 369-390 [in Russian].
[4] Bellman R. Stability theory of solutions of differential equations. M., IL, 1954 [in Russian].
[5] Atkinson F.V. On second-order non-linear oscillations. Pacific Journal of Mathematics, 1955, Vol. 5(1), pp. 643-647 [in English].
[6] Kiguradze I.T. Asymptotic properties of solutions of a nonlinear differential equation of Emden-Fowler type. Izv. AN SSSR. Seriia: Matematika [Proceedings of the Academy of Sciences of the USSR. Series: Mathematics], 1965, Vol. 29, Issue 5, pp. 965-986 [in Russian].
I.V. Filimonova, T.S. Khachlaev3
ON ASYMPTOTIC PROPERTIES OF SOLUTIONS, DEFINED ON THE HALF OF AXIS OF ONE SEMILINEAR ODE
The paper deals with the solutions of ordinary differential semi-linear equation, the coefficients of which depend on several real parameters. If the coefficient is chosen so that the equation does not contain the first-order derivative of the unknown function, it will be the case of Emden — Fowler equation. Asymptotic behavior of Emden — Fowler equation solutions at infinity is described in the book of Richard Bellman. The equations with the first-order derivative, considered in this work, erase in some problems for elliptic partial differential equations in unbounded domains. The sign of the coefficient in first-order derivative term essentially influences on the description of solutions. Partly the result of this paper can be obtained from the works of I.T. Kiguradze. In present work we use lemmas about the behavior of solutions of the linear equations with a strongly (weakly) increasing potential.
Key words: ordinary differential equations, nonlinear equations, semilinear equations, Emden — Fowler equation, asymptotic behavior of solutions, positive solutions, existence of solutions, maximum principle.
Статья поступила в редакцию 4/V7/2015. The article received 4/VJ/2015.
3 Filimonova Irina Vladimirovna ([email protected]), Department of Differential Equations, Lomonosov Moscow State University, 1, Leninskie Gory, Moscow, 119991, Russian Federation.
Khachlaev Timur Sultanovich ([email protected]), Department of Mathematics, Moscow State Institute of Radio Engineering, Electronics and Automation, 78, Vernadskogo Street, Moscow, 119454, Russian Federation.