Квант

34

А. Егоров

Уравнения и предыды

  1. Простой пример

Для приближенного вычисления квадратного корня из положительного числа aa поступают следующим образом (см., например, книгу «Алгебра и начала анализа 9», с. 83).

В качестве первого приближения берет произвольное число a0>aa_0 > \sqrt{a}. Следующое приближение a1a_1 находится по формуле

a1=12(a0+aa0),a_1 = \frac{1}{2} \left( a_0 + \frac{a}{a_0} \right),

по a1a_1 находят a2a_2:

a2=12(a1+aa1)a_2 = \frac{1}{2} \left( a_1 + \frac{a}{a_1} \right)

и вообще

an=12(an−1+aan−1).a_n = \frac{1}{2} \left( a_{n-1} + \frac{a}{a_{n-1}} \right).

Послеовательность {an}\{a_n\} монotonно убивает и ограничена. По теореме Вейерштрасса она имеет предел. Пусть lim⁡n→∞an=s\lim_{n \to \infty} a_n = s. Тогда s>0s > 0 и s=12(s+as)s = \frac{1}{2} \left( s + \frac{a}{s} \right), то есть s2=as^2 = a, а так как s>0s > 0, то s=as = \sqrt{a}.

Если взять достаточно большое nn, то ana_n будет достаточно мало отличаться от a\sqrt{a}.

Уражнение

  1. Докажите, что послеовательность {an}\{a_n\} монotonно убивает.
  2. Оцените разность ∣an−a∣|a_n - \sqrt{a}|.

Проанализируем наш способ вычисления послеовательных приближений для a\sqrt{a}. Послеовательность {an}\{a_n\} задавалась рекуррентно: an=φ(an−1)a_n = \varphi(a_{n-1}), где φ(x)=12(x+ax)\varphi(x) = \frac{1}{2} \left( x + \frac{a}{x} \right). Установив сходимость последовательности {an}\{a_n\}, мы заключили, что ее предел ss удовлетворяет равенству s=12(s+as)s = \frac{1}{2} \left( s + \frac{a}{s} \right), то есть является корнем уравнения

φ(x)=x.\varphi(x) = x.

(1)

На рисунке 1 изображен график функции y=12(x+2x)y = \frac{1}{2} \left( x + \frac{2}{x} \right), x>0x > 0, и показаны три первых члена последовательности

{a0=2,an=12(an−1+2an−1).\begin{cases} a_0 = 2, \\ a_n = \frac{1}{2} \left( a_{n-1} + \frac{2}{a_{n-1}} \right). \end{cases}

Уражнение 3. Выразите ana_n через nn, если a0=1a_0 = 1 и a) an=a⋅an−1a_n = a \cdot a_{n-1}; 6) an=an−1aa_n = a^a_{n-1}; 7) an=a+an−1a_n = \sqrt{a + a_{n-1}}; 8) an=an−11+an−12a_n = \frac{a_{n-1}}{\sqrt{1 + a_{n-1}^2}}; 9) an=11+an−1a_n = \frac{1}{1 + a_{n-1}}.

  1. Уравнения и предыды

Пусть нам нужно решить уравнение (1), где φ\varphi — некоторая функция. Возьмем произвольное число a0∈D(φ)a_0 \in D(\varphi) и определим рекуррентно последовательность {an}:an=φ(an−1)\{a_n\}: a_n = \varphi(a_{n-1}). Предположим, что an∈D(φ)a_n \in D(\varphi) для всех n∈Nn \in \mathbb{N}.

Рис. 1.

35

Теорема 1. Если последовательность {an}\{a_n\} сходится к s и функция φ непрерывна при x=sx = s, то s является корнем уравнения (1).

Эта теорема и была фактически применена в п. 1. Выбор первого приближения (или «нулевого члена» последовательности) диктуется сообщениями удобства — a0a_0 естественно быть поближе к исходному корню.

Как будет видно из дальнейшего, уменне исследовать разрешимость уравнения (1) позволяет иногда решить вопрос о сходимости соответствующей последовательности {an}\{a_n\}.

Упражнения
4. Исследуйте сходимость последовательностей:

a) an=a+an−1,a0=0a_n = \sqrt{a + a_{n-1}}, a_0 = 0;
b) an=11+an−1,a0=1a_n = \frac{1}{1 + a_{n-1}}, a_0 = 1.
5. Укажите способ нахождения положительного корня уравнения
x3−x−1=0.x^3 - x - 1 = 0. (1)
6. (С пособ приближенного в вычисления a3\sqrt[3]{a}.) Последовательность {an}\{a_n\} строится по следующему правилу:
{a0=a>0,an=13(2an−1+aan−1).\begin{cases} a_0 = a > 0, \\ a_n = \frac{1}{3}(2a_{n-1} + \frac{a}{a_{n-1}}). \end{cases}
Докажите, что lim⁡n→∞an=a3\lim_{n \to \infty} a_n = \sqrt[3]{a}.
7. Пусть a>0,b>0,p>0,q>0a > 0, b > 0, p > 0, q > 0. Расмотрим последовательность
ab,bpa+qb,pa+qbpb+q(pa+qb),…,anbn,…\frac{a}{b}, \frac{b}{pa+qb}, \frac{pa+qb}{pb+q(pa+qb)}, \dots, \frac{a_n}{b_n}, \dots
где an+1=bn,bn+1=pan+qbna_{n+1} = b_n, b_{n+1} = pa_n + qb_n. Существует ли у этой последовательности предел?
(Эту задачу предложил наш читатель В. Фищук.)

  1. Главный пример

Остальная часть статьи посвящена вопросу о сходимости последовательности
{α0=1,αn=aαn−1,\begin{cases} \alpha_0 = 1, \\ \alpha_n = a^{\alpha n-1}, \end{cases} (2)
то есть 1, a,aa,aaa,aaaa,…a, a^a, a^{a^a}, a^{a^{a^a}}, \dots. Эта сходимость зависит, конечно, от параметра aa.

Оказывается, для больших aa выяснить существование предела сравнительно легко, а для маленьких —

гораздо труднее. Приведем несколько примеров.

При a=2a = 2 получаем последовательность
1; 2; 4; 16; 65536; …
(она, очень/no, расходится).

При a=1,5a = 1,5 — последовательность
1; 1,5; 2,106; 2,349; 2,592; …
α11≈10,99;α12≈86,163;α13>1016\alpha_{11} \approx 10,99; \alpha_{12} \approx 86,163; \alpha_{13} > 10^{16}
(она тоже расходится, но здесь с самого начала это далеко не так очень/no).

При a=2=1,4142a = \sqrt{2} = 1,4142 …— последовательность
1; 2\sqrt{2}; 1,632; 1,761; …
α8≈1,966;α25≈1,9999\alpha_8 \approx 1,966; \alpha_{25} \approx 1,9999 (тут, как мы увидим, предел существует и равен 2).

При a=1a = 1
1, 1, 1, …
(предел, очень/no, существует и равен 1).

При a=0,5a = 0,5
1; 0,5; 0,707; 0,613; 0,654; …
α12≈0,641206;α13≈0,641177\alpha_{12} \approx 0,641206; \alpha_{13} \approx 0,641177.

При a=0,01a = 0,01
1; 0,01; 0,955; 0,0123; 0,945; 0,013; …
α15≈0,013;α16≈0,941\alpha_{15} \approx 0,013; \alpha_{16} \approx 0,941 (поведение двух последних последовательностей мы прокомментируем в конце статьи).

Из теоремы 1 мы знаем, что если lim⁡n→∞an=s\lim_{n \to \infty} a_n = s, то s является корнем

36

Уравнение
ax=xa^x = x
(функция x→axx \rightarrow a^x непрерывна во всех точках числовой прямой). Следовательно, если уравнение (3) при некотором aa не имеет корней, то соответствующая последовательность (2) расходит.

На рисунке 2 изображен график функции y=axy = a^x и первые члены последовательности (2) — при a>1a > 1 и предположении, что уравнение (3) имеет корни; на рисунке 3 — то же при 0<a<10 < a < 1.

  1. a>1a > 1
    При a>1a > 1 последовательность (2) монотно возрастает (проверьте!).

Теорема 2. При a>1a > 1 последовательность (2) сходится моеда и molъко моеда, когда уравнение (3) имеет хотя бы один корень. В этом случае ее предель равен меньшему из корней уравнения (3).

Доказательство. Мы только что отмечали, что если последовательность (2) сходится, то ее предель является корнем уравнения (3). Пусть теперь уравнение (3) имеет корни. Обозначим через s0s_0 наименьший корень. Покажем, что αn<s0\alpha_n < s_0 при всех n∈Nn \in \mathbb{N}. Легко видеть, что s0>1=α0s_0 > 1 = \alpha_0. Если же s0>αns_0 > \alpha_n, то s0=as0>aαn=an+1s_0 = a^{s_0} > a^{\alpha_n} = a_{n+1}. Из теоремы Вейерштрасса следует, что последовательность (2) сходится, причем, очевидно, к s0s_0. Теорема доказана.

  1. 0<a<10 < a < 1
    В этом случае последовательность (2) ограничена: 0<αn<10 < \alpha_n < 1. К сожалению, она не монотна. В самом деле, поскольку 0<a<10 < a < 1, функция y=axy = a^x убивает, и α2=aa<1=α0\alpha_2 = a^a < 1 = \alpha_0, следовательно, если уравнение (3) при некотором aa не имеет корней, то есть α3>α1\alpha_3 > \alpha_1. Далее, α4=aα3<aα1=α2\alpha_4 = a^{\alpha_3} < a^{\alpha_1} = \alpha_2, то есть 1=α0>α2>α41 = \alpha_0 > \alpha_2 > \alpha_4. Пользовать метод математической индукции, легко доказать, что подпоследовательность α1,α3,α5,…,α2n−1,…\alpha_1, \alpha_3, \alpha_5, \dots, \alpha_{2n-1}, \dots нашей последовательности монотно возрастает, а подпоследовательность α0,α2,α4,…,α2n,…\alpha_0, \alpha_2, \alpha_4, \dots, \alpha_{2n}, \dots монотно убивает. При этом всегда α2n>α2n−1\alpha_{2n} > \alpha_{2n-1}. По теореме Вейерштрасса каждая из этих подпоследовательностей имеет предел. Пусть lim⁡n→∞α2n−1=s1\lim_{n \to \infty} \alpha_{2n-1} = s_1 и lim⁡n→∞α2n=s2\lim_{n \to \infty} \alpha_{2n} = s_2. Ясно, что s2≥s1s_2 \geq s_1 и подпоследовательность (2) будет сходящейся в том и только в том случае, когда s1=s2s_1 = s_2 (доказательство!).

Что можно сказать про числа s1s_1 и s2s_2? Поскольку α2n+2=aaα2n\alpha_{2n+2} = a^{a^{\alpha_{2n}}}, α2n+1=aaα2n−1\alpha_{2n+1} = a^{a^{\alpha_{2n-1}}} и функция x→aaxx \rightarrow a^{a^x} непрерывна во всех точках числовой прямой, числа s1s_1 и s2s_2 являются корнями уравнения (3).

Теорема 3. Числа s1s_1 и s2s_2 являются соотношением, меньшим и большим корнями уравнения (4).

Доказательство. Пусть ss — произвольный корень уравнения (4). Для доказательства теоремы достаточно установить, что α2n−1<s<α2n\alpha_{2n-1} < s < \alpha_{2n} для всех n∈Nn \in \mathbb{N}.

Прежде всего заметим, что s<1s < 1. В самом деле, если s≥1s \geq 1, то, ввиду as>0a^s > 0, aas<1a^{a^s} < 1, то есть aas≠sa^{a^s} \neq s. Кроме того, очень удобно, s>0s > 0. Из 0<s<10 < s < 1 следует a0>as>a1a^0 > a^s > a^1, то есть α0>as>α1\alpha_0 > a^s > \alpha_1. Значит, aα0<aas<aα1a^{\alpha_0} < a^{a^s} < a^{\alpha_1}, то есть α1<s<α2\alpha_1 < s < \alpha_2. Индукционный шаг проделывается аналогично.

Следствие 1. Если уравнение (4) имеет более одного корня, то последовательность (2) расходится.

Следствие 2. Если уравнение (4) имеет более одного корня, то последовательность (2) сходится, и предель

37

Рис. 4.

этой последовательности есть корень уравнения (3).

Легко видеть, что всякий корень уравнения (3) удовлетворяет также и уравнению (4). Но при 0<a<10 < a < 1 уравнение (3), очевидно, имеет корень (нарисуйте графики!). Значит, при таких aa и уравнение (4) обязаτεльно имеет хотя бы один корень.

Упражнение 8: Рассмотрим последовательность α0=x>0,αn=aαn−1\alpha_0 = x > 0, \alpha_n = a^{\alpha n-1}. При каких xx эта последовательность сходится? (Рассмотите отдельно случай a>1a > 1 и 0<a<10 < a < 1).

Итак, вопрос о сходимости последовательности (2) совпается к вопросу о числе корней уравнения (3) при a>1a > 1 и уравнения (4) — при 0<a<10 < a < 1.

  1. Исследование уравнения (3) при a>1a > 1

Уравнение (3) равносильно уравнению

ln⁡xx=ln⁡a.\frac{\ln x}{x} = \ln a. (рис. 4). Эта функция определяна при x>0x > 0, достигнет максимума, равного 1/e1/e, при x=ex = e, убивает при x>ex > e и возвращает при 0<x<e0 < x < e.

Упражнение 9. Проведите подобное исследование функций y=ln⁡xxy = \frac{\ln x}{x}.

Следовательно, при ln⁡a>1/e\ln a > 1/e, то есть a>e1/ea > e^{1/e}, уравнение (5) корней не имеет.

Рис. 5.

При 0<ln⁡a<1/e0 < \ln a < 1/e, то есть 1<a<e1/e1 < a < e^{1/e}, уравнение (5) имеет два корня.

При a=e1/ea = e^{1/e} оно имеет один корень.

Заметим для дальнейшего, что при ln⁡a≤0\ln a \leq 0, то есть 0<a≤10 < a \leq 1, уравнение (5) имеет единственный корень.

  1. Исследование уравнения (4) при 0<a<10 < a < 1

Мы уже знаем, что все корни уравнения (4) заключены в интервале [0;1][0; 1]. Уравнение (4) равносильно уравнения ax=log⁡axa^x = \log_a x, x=log⁡alog⁡axx = \log_a \log_a x или x=1ln⁡aln⁡ln⁡xln⁡ax = \frac{1}{\ln a} \ln \frac{\ln x}{\ln a}. Обозначим через φ\varphi функцию

φ(x)=1ln⁡aln⁡ln⁡xln⁡a−x.\varphi(x) = \frac{1}{\ln a} \ln \frac{\ln x}{\ln a} - x.

Мы хотим выяснить, сколько корней в зависимости от параметра aa имеет уравнение φ(x)=0\varphi(x) = 0, равносильно уравнению (4). Чтобы это сделать, нужно посчитать, сколько существует точек пересечения графика функций φ\varphi с осью абсцисс.

Прежде всего заметим, что к этим точкам пересечения заведomo принадлежит точка s0s_0, где s0s_0 — единственный корень уравнения (3), (5),— см. заменяния в конце пп. 6 и 5.

Исследуем функцию φ\varphi.

Так как ln⁡a<0\ln a < 0, функция φ\varphi определяна при 0<x<10 < x < 1. При x→0x \rightarrow 0 (справа) φ(x)→−∞\varphi(x) \rightarrow -\infty; при x→1x \rightarrow 1 (слева) φ(x)→+∞\varphi(x) \rightarrow +\infty.

38

Найдем экстремумы функции ff; для этого исследуем ее производуно:

φ′(x)=1xln⁡xln⁡a−1.\varphi'(x) = \frac{1}{x \ln x \ln a} - 1.

Приравняв ее к нулю, мы получим уравнение

xln⁡x=1ln⁡a.x \ln x = \frac{1}{\ln a}.

(6)

На рисунке 5 изображен график функции y=xln⁡xy = x \ln x.

Упражнение 10. Исследуйте функцию y=xln⁡xy = x \ln x и построите ее график.

Мы видим, что при x=1/ex = 1/e функция y=xln⁡xy = x \ln x достигает минимума, равного −1/e-1/e. Поэтому, если 1/ln⁡a<−1/e1/\ln a < -1/e, то есть ln⁡a>−e\ln a > -e или a>e−ea > e^{-e}, то уравнение (6) корней не имеет. Так что при e−e<a<1e^{-e} < a < 1 производаная функция ff строго больше нуля, то есть при таких значениях aa функция ff строго возрастает. Значит, уравнение f(x)=0f(x) = 0 при e−e<a<1e^{-e} < a < 1 имеет единственный корень s0s_0 (примерный график функции ff для e−e<a<1e^{-e} < a < 1 изображен на рисунке 6).

Если a=e−ea = e^{-e}, то φ′(x)=0\varphi'(x) = 0 при x=1/ex = 1/e и φ′(x)>0\varphi'(x) > 0 при всех остальных x∈[0;1/e]∪[1/e;1]x \in [0; 1/e] \cup [1/e; 1].

В этом случае ff тоже строго возрастает во всей области определения, и уравнение f(x)=0f(x) = 0 лю-прежнему имеет единственный корень x=1/ex = 1/e (график функции ff для a=e−ea = e^{-e} изображен на рисунке 7).

Итак, если e−e≤a<1e^{-e} \leq a < 1, уравнение (4) имеетровно один корень.

Если же 0<a<e−e0 < a < e^{-e}, то уравнение (6) имеет два корня: x1x_1 и x2x_2 (рис. 5). Причем при x∈[0;x1]x \in [0; x_1] x1x_1 [∪] x2x_2; 1 [производаная φ′(x)\varphi'(x) положительная, а при x∈[x1;x2]x \in [x_1; x_2] — отрицательна, так что x1x_1—точка максимума, а x2x_2—точка минимума функции ff; значит, φ\varphi при x∈[0;x1]x \in [0; x_1] [возрастает и при x∈[x1;x2]x \in [x_1; x_2] 1 [возрастает, а при x∈[x1;x2]x \in [x_1; x_2] убывает (разберитесь в этом!).

Покажем теперь, что s0∈[x1;x2]s_0 \in [x_1; x_2]. Для этого достаточно доказать, что φ′(s0)<0\varphi'(s_0) < 0. Имеем

φ′(s0)=1s0ln⁡s0ln⁡a−1=1(ln⁡s0)2−1\varphi'(s_0) = \frac{1}{s_0 \ln s_0 \ln a} - 1 = \frac{1}{(\ln s_0)^2} - 1

(так как ln⁡s0s0=ln⁡a\frac{\ln s_0}{s_0} = \ln a). Докажем, что ln⁡s0<−1\ln s_0 < -1 (тогда φ′(s0)<0\varphi'(s_0) < 0). Для этого обратимся вновь к функции y=ln⁡xxy = \frac{\ln x}{x} (рис. 4). Поскольку она при 0<x<10 < x < 1 возрастает, а ln⁡a<−e\ln a < -e, корень s0s_0 уравнения (5) лежит левее корня уравнения ln⁡xx=−e\frac{\ln x}{x} = -e, то есть левее точки e−1e^{-1}. Значит, s0<1/es_0 < 1/e и ln⁡s0<−1\ln s_0 < -1.

Итак, про функцию ff нам теперь известно, что

  1. φ(x)→−∞\varphi(x) \rightarrow -\infty при x→0x \rightarrow 0;
  2. φ(x)→+∞\varphi(x) \rightarrow +\infty при x→1x \rightarrow 1;
  3. φ\varphi возрастает на промежутке ∣0;x1∣|0; x_1|, в точке x1x_1 имеет максимум, убывает на промежутке ∣x1;x2∣|x_1; x_2|, в точке x2x_2 имеет минимум и снова возрастает на промежутке ∣x2;1∣|x_2; 1|;