4.3.4 Итерационные циклы

Среди циклов с неизвестным числом повторений большое место занимают итерационные. В этих циклах количество повторений заранее неизвестно. Выход из цикла происходит, если достигнута заданная точность результата или появилась стандартная ситуация (конец данных). Такие циклы используются при решении задач численными методами: приближенное вычисление функций, определенных интегралов, решение уравнений и т.д.

Основной особенностью итерационных циклов является то, что новое значение искомой величины (Y_след) вычисляется с использованием старого (Y_пред) по формуле: Y_след=f(Y_пред)

Обычно для вычисляемой величины Y задается некоторое начальное значение (приближение) Y_o и погрешность вычисления результата (Y) в виде: E=(0,5; 1)*10^-m, где m – целое число. Как правило, E=(0,1; 0,5, 0,01; 0,05, 0,001 и т.д.). Считается, что заданная точность достигнута, если найденные значения Y_пред и Y_след отличаются на величину, меньшую E, т.е. условие окончания вычислений можно записать так:

|Y_след – Y_пред|<E

Для программирования итерационных циклов в Паскале используются операторы While и Repeat.

Пример 4.11 Вычислить значение арифметического выражения:

Вычисление непрерывных радикалов производится в цикле, начиная от внутреннего радикала. В данной задаче начальное значение . Каждое следующее значение радикала будет вычисляться через предыдущее значение радикала по формуле , число изменяется от начального значения 5 до конечного значения 98 с шагом 3.

Var r,a:real;

begin

r:=sqrt(2); a:=5;

while a<=98 do

begin

r:=sqrt(a+r);

a:=a+3;

end;

writeln('R=',r);

end.

Приближенное вычисление функций

Пример 4.12 Вычислить значение функции SinX по приближенной формуле F = SinX = X - X³/3! + X⁵/5! - ...+(–1)^2n-1 X^2n-1/(2n-1)!, n=1, 2, 3 ..., с погрешностью E=10^-5. Функция F здесь представляет собой сумму членов степенного ряда: F = U₁ + U₂ + ... + U_n. Условие окончания вычисления функции: |F - F_след|< E

Заметим, что |F-F_след| = |U_след| – очередное слагаемое. Можно сказать, что вычисления прекращают, как только получают член ряда, меньший, чем заданная абсолютная погрешность.

Обратите внимание, что вычислять слагаемые ряда непосрественно по формуле (находить степень и факториал, делить их) крайне нерационально. В этих случаях целесообразно выражать последующий член через предыдущий. Покажем, как это делается.

Пусть речь идет о слагаемом с номером n. Тогда на основании общей формулы для члена ряда это слагаемое (предыдущее) без учета знака запишется так: .

Номер следующего слагаемого будет n+1, тогда само слагаемое запишется так: .

Разделив последнее выражение на предыдущее, получим выражение для вычисления следующего члена ряда по известному предыдущему

.

Обратите внимание, что в этой формуле n означает номер предыдущего члена. Если считать, что n имеет начальное значение 1 и изменяется с шагом 2, то знаменатель U_след может быть вычислен как произведение вида:

(n + 1) (n+2).

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

Вначале: n=1; U=X; F=0 {U – играет роль предыдущего и следующего члена}

Program NovSin;

Var

F, X, E, U : Real;

n : Integer;

Begin

Writeln('Введите аргумент Х и погрешность Е');

Readln(X,E);

n := 1;

U :=X ;

F:=0;

While Abs(U)>E Do

Begin

F:=F+U;

U:=-U*Sqr(X)/(2*N*(2*N+1));

N:=N+1;

End;

Writeln(X:8:6,' SinX=',F:10:6, ' Машинный SinX=',Sin(X):10:6);

End.

Количество повторений цикла заранее неизвестно. Оно зависит от требуемой точности и от значения аргумента. Если X=0, то цикл не выполняется ни разу. Количество повторений цикла в данном случае можно узнать по номеру N последнего вычисленного члена. В других случаях это можно сделать с помощью обычного счетчика:

k := 0;

while abs(U)>E do

begin

k := k+1;

F :=.......;

. . . . . .

end;

writeln('k= ',k);