Програмування: рекурентні послідовності та співвідношення (реферат)

Реферат з інформатики

Програмування: рекурентні послідовності та співвідношення

1.Рекурсії.

Для обчислення степеня в алгоритмі накопичування добутку (див. П. 3.3)
змінна p приймала значення 1, a, a2, a3, … , an. У цій послідовності
перший член 1, а кожний наступний дорівнює попередньому, помноженому на
a. Позначивши члени послідовності через p0, p1, p2, … pn, маємо
рівність: pi=pi-1*a при i=1,2,…,n. Така рівність, що виражає член
послідовності через попередні (один або кілька), називається рекурентним
співвідношенням.

«Рекурентний» означає «зворотний». Справді, елемент послідовності тут
визначається через попередні, і для його обчислення треба повернутися до
них. Усім добре відомі рекурентні співвідношення вигляду an=an-1+d або
bn=bn-1* q – їм задовольняють члени відповідно арифметичних або
геометричних прогресій. Конкретна ж прогресія, тобто послідовність
чисел, задається першим членом a1 і різницею d (або знаменником q).
Власне, послідовність степенів у прикладі p0, p1, p2, … – геометрична
прогресія: вона визначається першим членом p0=1 і рекурентним
співвідношенням pi=pi-1*a при будь-якому i>0. Послідовність, члени якої
задовольняють деяке рекурентне співвідношення, також називається
рекурентною.

Приклад . Розглянемо послідовність {f} чисел 1, 1, 2, 3, 5, 8, 13, …
, у якій f1=f2=1, а наступні члени задаються рекурентним співвідношенням

fn=fn-2+fn-1, n>2.

Вона називається послідовністю чисел Фібоначчі – за прізвиськом Леонардо
Пізанського, який першим її описав. За першими двома її членами можна
обчислити третій. Для обчислення четвертого перший член уже не
потрібний, тому що f4=f2+f3. Для обчислення п’ятого достатньо пам’ятати
лише третій і четвертий тощо. Обчислюючи члени послідовності один за
одним, ми дістанемося будь-якого, почавши з перших двох. При цьому
щоразу ми використовуємо лише два останніх значення і, обчисливши
наступне, «забуваємо» перше з двох використаних.

Нехай дано номер n, n>2, і треба обчислити fn. Опишемо ці обчислення. З
попередніх міркувань випливає, що потрібні дві змінні для двох сусідніх
членів і третя для наступного (назвемо їх fa, fb і fc), а також змінна
m для зберігання номера останнього з обчислених членів.

Спочатку fa=1, fb=1, m=2, (*)

потім обчислимо fc:=fa+fb і збільшимо m на 1. Якщо значення fb і fc
зробити відповідно значеннями fa і fb (fa:=fb, fb:=fc), то обчислення
четвертого члена можна задати таким самим оператором fc:=fa+fb. Отже,
поки m1.

Умови продовження обчислень можуть бути різними, наприклад, >d або >d
для деякого d>0.

Розглянемо друге з них. Оскільки в ньому вказано два сусідніх члени,
потрібні дві змінні для їх збереження,

причому обидві повинні мати різні значення вже перед першою перевіркою
умови продовження. Після того, як

вона виявляється істинною, для обчислення нового члена передостанній
член уже не потрібний, тому що

рекурентне співвідношення має порядок 1. Тому в тілі циклу треба
спочатку вказати переприсвоювання, а потім

обчислення нового члена. Номера членів послідовності нас не цікавили,
тому алгоритм має вигляд:

X:=(a+1)/2; Y:=0.5*(X+a/X);

while abs(X-Y)>d do

begin

X:=Y; Y:=0.5*(X+a/X);

end;

{ abs(X-Y)<=d, значення Y вважається шуканим, адже |Y-a| ||. Тому,

якщо додати всі члени від першого до останнього з таких an, що |an|>d за
деякого d>0, то одержана сума

відрізняється від sinx не більш, ніж на d.

Отже, треба обчислити sn=, де n невідомо, а відомо лише, що |an|>d,
|an+1|? d. Очевидно,

sn=sn-1+an за будь-якого n>1, а s1=a1=x. Ці рівності виражають
залежність значення суми від попередньої суми і

відповідного доданка, тобто послідовність значень сум рекурентна.
Помітимо, що при d<|x| доданок a1 не треба додавати до суми, і результатом повинна бути "сума без доданків". Тому до послідовності сум додамо s0=0; тепер sn=sn-1+an для n>0.

Знайдемо рекурентне співвідношення для послідовності доданків ,
виразивши an через

an-1. Для цього у виразі для an побудуємо вираз, яким задається an-1:

=

= .

Отже, при n>1, a1=x. Запишемо одержані рекурентні співвідношення в
систему:

Побудуємо за нею алгоритм обчислення. Оскільки порядок обох
співвідношень 1, достатньо двох змінних, S і A,

для збереження членів послідовностей. Спочатку A:=x; S:=0. Далі перед
кожним обчисленням S:=S+A треба

спочатку перевірити, що A>d. Після додавання A до S обчислюється новий
доданок (значення A), і все

повторюється. Таким чином, цикл складений діями в такому порядку:

перевірка умови A>d,

додавання S:=S+A,

обчислення нового значення A.

Нехай змінна I зберігає номер останнього обчисленого доданка; спочатку
I=1. Оскільки при обчисленні нового

доданка використовується його номер, то цей номер треба попередньо
збільшити. Тепер алгоритм очевидний:

S:=0; A:=x; I:=1;

while A>d do

begin

S:=S+A; I:=I+1;

A:=A*(-x*x)/((2*I-2)*(2*I-1));

end

{A<=d, і воно не додано до S; значення S – шукане} Оформлення цього алгоритму у вигляді функції з параметром x і розумно підібраним значенням d залишається вправою.? Задачі 4.4.* Написати функцію обчислення кількості десяткових цифр натурального числа. 4.5.* Один із варіантів алгоритму Евкліда обчислення найбільшого спільного дільника чисел a і b (НСД(a,b)) грунтується на обчисленні рекурентної послідовності {p}, де p1=max{a,b}, p2=min{a,b}, pn=pn-2 mod pn-1 при n>2. Шуканим є останнє ненульове

значення послідовності. Уточнити алгоритм Евкліда у вигляді функції.

4.6. Послідовність {xn}, задана співвідношеннями

x1=(a+m-1)/2,

xi=( (m-1)xi-1 + a/x)/m при i > 1,

сходиться до a1/m. Запрограмувати обчислення a1/m при довільному
додатному дійсному a з точністю e , тобто за потрібне число

приймається перше xn таке, що | xn-xn-1 |

Добавить комментарий

Ваш e-mail не будет опубликован. Обязательные поля помечены *