Краткие теоретические сведения. Действительная функция f(x) называется аналитической в точке e, если в некоторой окрестности ½x-e½<R этой точки функция разлагается в
Действительная функция f(x) называется аналитической в точке e, если в некоторой окрестности ½x-e½<R этой точки функция разлагается в степенной ряд (ряд Тейлора):
(1)
При e=0 получаем ряд Маклорена:
(2)
Разность (3)
называется остаточным членом и представляет собой ошибку при замене функции f(x) полиномом Тейлора.
Для ряда Маклорена
где 0<q<1. (4)
Таким образом, вычисление значения функции можно свести к вычислению суммы числового ряда
а1+а2+ . . . +an+ . . . . (5)
Известно, что числовой ряд называется сходящимся, если существует предел последовательности его частных сумм:
, (6)
где Sn= а1+а2+ . . . +an+ . . . .
Число S называется суммой ряда.
Из формулы (13) получаем S=Sn+Rn ,
где Rn - остаток ряда, причем R®0 при n®¥.
Для нахождения суммы S сходящегося ряда (5) с заданной точностью e нужно выбрать число слагаемых n столь большим, чтобы имело место неравенство
½Rn½<e.
Тогда частная сумма Sn приближенно может быть принята за точную сумму S ряда (5).
Приближенно n выбрать так, чтобы имело место неравенство
|Sn+1-Sn|<e или an<e.
Задача сводится к замене функции степенным рядом и нахождению суммы некоторого количества слагаемых при различных параметрах суммирования х . Каждое слагаемое суммы зависит от параметра х и номера n, определяющего место этого слагаемого в сумме.
Обычно формула общего члена суммы принадлежит одному из следующих трех типов:
а) ; ; ;
б) ; ; ;
в) ; ; .
В случае а) для вычисления члена суммы аn целесообразно использовать рекуррентные соотношения, т. е. выражать последующий член суммы через предыдущий: an+1=y(x, n)an. Это позволит существенно сократить объем вычислительной работы. Кроме того, вычисление члена суммы по общей формуле в ряде случаев невозможно (например из-за наличия n!).
В случае б) применение рекуррентных соотношений нецелесообразно. Вычисления будут наиболее эффективными, если каждый член суммы вычислять по общей формуле an=f(x, n).
В случае в) член суммы целесообразно представить в виде двух сомножителей, один из которых вычисляется по рекуррентному соотношению, а другой непосредственно an=f(x, n)*сn(x,n), где сn=cn-1y(x,n).
Постановка задачи
Для х изменяющегося от a до b с шагом (b-a)/k, где (k=10), вычислить функцию f(x), используя ее разложение в степенной ряд в двух случаях:
для заданного n;
для заданной точности e (e=0.0001).
Для сравнения найти точное значение функции.
Варианты
№ | функция | диапазон изменения аргумента | n |
№ | сумма |
Методические указания
Алгоритм решения задачи сводится к трем циклам, причем два из них вложены в третий. Внутренние циклы суммируют слагаемые при фиксированном параметре x, один (арифметический для заданного n), другой (итерационный для заданной точности e). При организации этих циклов следует обратить внимание на правильный выбор формулы для вычисления элемента ряда an и правильное присвоение начальных значений переменным цикла. Внешний цикл организует изменение параметра х.
Результаты расчетов представить в следующем виде:
Вычисление функции
X=...... SN=...... SE=..... Y=......
X=...... SN=...... SE=..... Y=......
..........
X=...... SN=...... SE=..... Y=......
Здесь X- значение параметра; SN- значение суммы для заданного n; SE- значение суммы для заданной точности; Y-точное значение функции.
Содержание отчета
1. Постановка задачи.
2. Вариант задания.
3. Математическая модель (формулы, по которым выполняются вычисления слагаемых ряда).
4. Листинг программы.
5. Полученные результаты.
6. Выводы.
Лабораторная работа № 4
"РАБОТА С ОДНОМЕРНЫМИ МАССИВАМИ"
Цель: Получение навыков обработки одномерных массивов.