3.4. Рекуррентные уравнения

Рассмотрим обобщение последовательностей Фибоначчи. Формула

un + r  = c1 un + r – 1 + c2 un + r – 2 + …+ cr un

называется однородным  линейным  рекуррентным уравнением  порядка r. Ее решением является последовательность {un}, однозначно  определенная начальными значениями u0, u1, u2, …∙, ur –1. Решение такого уравнения называется возвратной, или рекуррентной  последовательностью порядка r.

Пример 1

Геометрическая прогрессия является решением уравнения un+1 = qun. Ее члены описываются формулой un = u0qn. Отсюда геометрическая прогрессия является возвратной последовательностью порядка 1.

Пример 2

Арифметическая прогрессия un = u0 + nd удовлетворяет соотношению:

un+1 - un = un+2 - un+1.

Получаем однородное рекуррентное уравнение:

un+2 = 2un+1 - un.

Начальные данные задаются значениями u0 и  u1 = u0 + d. Отсюда, арифметическая прогрессия является возвратной последовательностью порядка 2.

Пример 3

Произвольная периодическая последовательность является возвратной последовательностью порядка p, удовлетворяющей рекуррентному  соотношению:

un + p=un.

Здесь p – период последовательности.

Для заданного рекуррентного уравнения

un + r  = c1 un + r – 1 + c2 un + r – 2 + …+ cr un

найдем производящую функцию  возвратной последовательности {un}.

Обозначим:

K(x) = 1- c1x - c2x2 -… - crxr.

Теорема 1

Произведение u(x)K(x) = D(x) является многочленом степени, меньшей r.

Доказательство

Вычислим коэффициент ряда D(x) при xn+ r.  Он при  n ≥ 0 будет равен:

un + r - (c1 un + r – 1 + c2 un + r – 2 + …+ cr un) = 0.

Отсюда D(x) – многочлен степени, меньшей чем n.

Замечание

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

Пример 4

Применим доказанную теорему к  решению рекуррентного уравнения

un+2  = 5 un+1  - 6 un

при начальных условиях u= u1 = 1.

Здесь K(x) = 1 - 5x + 6x2.

Вычислим:

D(x) = K(x)u(x) =  (1 - 5x+6x2)(u0 + u1x + … ) = (u- 5u0 x + u1x +…).

По теореме 1 коэффициенты при   x2,   x3…   равны  нулю,   а u= u1 = 1. Следова

тельно,

D(x) = (1  -5 x +  x) = 1 - 4x.

Получаем:

.

Следующий шаг – разложение знаменателя K(x) в произведение

(1 - a1x) (1- a2x).

В данном случае это можно сделать с помощью формулы Виеты. Поскольку имеют место равенства:

 a1 + a2 = 5;                 a1 a2 = 6,

то числа  a1  и  a2  будут корнями квадратного уравнения   a2 - 5a +6 = 0. Решая это квадратное уравнение, получаем:

  .

Приходим к формуле

.

Теперь найдем разложение в сумму простых дробей методом неопределенных коэффициентов:

 .

Получим систему линейных уравнений:

  

имеющую решение A = 2, B = -1. Отсюда

.

Это приводит к ответу

un = 2n+1-3n.

В общем случае числа aI   в разложении 

K(x) = (1 - a1x) (1- a2x) … (1 - arx)

являются корнями уравнения:

F(a)=a r- c1a r -1  -- cr-1a - cr = 0,

ибо K(x) = . Это уравнение  F(a) = 0 называется характеристическим.

Если все корни уравнения  F(a)=0  действительны и различны, то получаем:

,

откуда

.

Это позволяет составить систему линейных уравнений для нахождения Ai с помощью известных значений u0, u1, ∙, ur-1.

Пример 5

Рассмотрим рекуррентное уравнение:

при заданных  u0 = 2,    u1=1.

Составим   характеристическое  уравнение:

F(a) = a2 + a - 2 = 0.

Его корни a 1 = 2, a 2 = -1 не равны между собой и являются вещественными числами. Поэтому решение рекуррентного уравнения можно искать в виде:

 .

Значения un  известны при n = 0, 1, откуда получаем систему линейных уравнений для нахождения коэффициентов A1, A2:

Решаем эту систему уравнений:

получаем A1 = A2 = 1. Приходим к следующему ответу:

.

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

,

где k – кратность корня .

Пример 6

Решим рекуррентное уравнение:

,      u0 = 1, u1 = 6.

С этой целью составим характеристическое уравнение:

a2- 6a +9=0.

Оно имеет кратные корни a 1 = a 2 =3. Поэтому решение рекуррентного уравнения нужно искать в виде

un  A13n + A2 n 3n

Подставляя известные значения u0  = 1 и u1  = 6, приходим к системе уравнений

Она имеет решение A1 = A2 = 1. Получаем ответ:

un = 3n(1 + n).

Упражнения

Свойства производящих функций

1. Сколько делителей, включая само число и 1, имеет число 720?

2. Доказать, что при n > 0 имеют место соотношения:

1) ;

2) .

3. Найти производящую функцию последовательности an = 10ncos(2n).

4. Найти производящую функцию последовательности an = 10nsin(2n).

5. Найти производящую функцию последовательности an = 10nn.

6. Найти производящую функцию последовательности an = 10nn.

7. Найти производящую функцию последовательности an = 10nn(n+1).

8. Найти производящую функцию последовательности an = 10nn2.

9. Найти производящую функцию последовательности an = 10n/(n+1).

10. Найти производящую функцию последовательности an = 10n/(n+1)(n+2).

11. Доказать рекуррентное соотношение для производящих функций  последовательностей nk:

.

12. Найти производящие функции последовательностей an = nk,  при  k = 1, 2, 3, 4.

Решение рекуррентных уравнений

13. un+2 = un+1 + 2un,  u0 = 1,  u1 = 1.

14. un+2 = -un+1 + 2un,  u0 = 0,  u1 = 1.

15. un+2 = -3un+1 – 2un,  u0 = 1,  u1 = 1.

16. un+2 = un+1 + 6un,  u0 = a,  u1 = b.

17. un+2 – 4un+1 – 5u= 0,  u0 = 8,  u1 = 10. 

Ответ: un = 7+3n

18. un+3 – 3un+2 +un+1– 3u= 0,  u0 = 1,  u1 = 3, u2 = 8.     

Ответ: u2n= (9×32n+(-1)n)/10,   u2n+1= (9×32n+1+3×(-1)n)/10, n≥0.

19. un+3 + un+2 – un+1– u= 0,  u0 = 1,  u1 = 2, u2 = 3.        

Ответ: un=(-1)n(n-1)+2

20.  un+2 – 4un+1 + 4u= 0,  u0 = a,  u1 = b.             

21.  un+2 = 2un+1 – 2un,  u0 = 1,  u1 = 2.

Указание: Производящая функция последовательности un:

u(x)= .

Ответ:

22. un+3 = 3un+2 – 3un+1+ un,  u0 = 1,  u1 = 3, u2 = 3;           

Указание: Искать решение в виде: un=A+ Bn + Cn2.

23. un+3 =  – 3un+2 – 3un+1– un,  u0 = 0,  u1 = 1, u2 = –2.

24. un+2 = un+1 – un,  u0 = 1,  u1 = 9.                                      Ответ:   

25. un+2 = 2un+1 – 4un,  u0 = 1,  u1 = 2.                                   Ответ: