Вычислить определенный интеграл

Лабораторная работа № 8

Вычисление производной функции

Вычисление определенных интегралов

Цель работы:

выработать навык использования табличного процессора MS Excel для численного (приближенного) вычисления производных функций (дифференцирование) и определенных интегралов (интегрирование) заданных в табличном виде.

Задание 1

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

 

Методика выполнения задания.

Способы вычисления производной функции, заданной в табличном виде:

 

способ правой конечной разности         


 

 


способ левой конечной разностью               


способ центральных разностей     

 

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

 

 

Пример

Вычислить производную функции     на интервале [0; 5].

Аналитически производная выражается формулой .

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

 

 

Порядок решения примера:

1. В ячейку А2 вводим 0 (начальное значение интервала) и с шагом 0,5 заполняем диапазон до значения 5 (конечная ячейка А12).

2. Вводим формулы:

в ячейку В2 =1/(1+ A 2) (функция 1/(1+x)) и копируем до ячейки В12;

в ячейку С2 (Правые разности) =(В3-В2)/( A 3- A 2) и копируем до ячейки С11;

в ячейку D 3 (Левые разности) =(В3-В2)/( A 3- A 2) и копируем до ячейки D 12;

в ячейку Е3 (Центральные разности) =(В4-В2)/( A 4- A 2) и копируем до ячейки E 11;

в ячейку F 2 (Точное значение) =−1/(1+ A 2)^2 (аналитически вычисленная производная -1/(1+x)^2) и копируем до ячейки F 12.

 


На рисунке приведено расположение формул в ячейках таблицы.

 

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

 

 

Вычислить производные

 

В вариантных примерах шаг табулирования функции внутри интервала выбрать по формуле:

(максимальное значение − минимальное значение) / 10.

 

 

Вариант Функция Производная Интервал
1 [10; 110]
2 [1; 20]
3 [0; 10]
4 [0; 3]
5 [4; 12]
6 [0; 0,5]
7 [0; p]
8 [0; p]
9 [1; 8]
10 [2; 6]
11 [1; 4]
12 [0; 10]
13 [0; 2]
14 [1; 5]
15 [10; 80]
16 [1; 5]
17 [0; 3]
18 [0; 1]
19 [2; 12]
20 [1; 2]

 

 

Задание 2

Для примера своего варианта, используя MS Excel, вычислить определенный интеграл тремя способами: правостороннее приближение, левосторонние приближение и по формуле трапеций. Провести анализ вычислений для всех способов. Сделать выводы.

Методика выполнения задания.

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


   Правостороннее приближение.                 Левостороннее приближение.

Приближение трапеций

 

Если все интервалы имеют одинаковую длину h , формулы имеют вид  для право­стороннего приближения и   для левостороннего приближения. Однако обычно для вычисления интегралов используется формула трапеций. Она несколько слож­нее по сравнению с первыми двумя, но зато позволяет получать, как правило, более точные результаты. Идея состоит в том, что соседние узловые точки функции соеди­няются прямой, поэтому вся площадь под графиком функции состоит как бы из трапе­ций. Соответственно, сама площадь равна сумме площадей этих трапеций. Формула для интервалов равной длины  имеет вид .

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

Далее рассмотрим, как описанные методы могут использоваться для вычисления интегралов на практике.

 

 

Пример.

Вычислить интеграл  .

Этот интеграл может быть вычислен аналитически. В частности, очевидно, что  

Именно это значение и попытаемся получить в результате интегрирования с использованием численных методов.

 

1. Оформляем таблицу, как показано на рисунке.

2. В ячейку А4 вводим число 0 как начальное значение для переменной интегрирования. После этого в ячейку А5 вводится формула =А4+ПИ()/10, согласно которой значе­ние последующего узла получается прибавлением к предыдущему узлу десятой части от длины диапазона интегрирования (команда ПИ()/10) и копируем ее до А14.

3. В ячейку В4 вводим заданную формулу = SIN ( A 4) и копируем ее до ячейки В14.

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


В ячейку С4 вводим формулу =(А5 - А4)*В5 (правосторонние приближение), в ячейку D 4 вводим практически такую же формулу =(А5 - А4)*В4 (левосторонние приближение), и для метода трапеций в ячейку Е4 вводится формула =(А5 - А4)*(В5+В4)/2.

 

Заполнять данные для элементарных площадей следует вплоть до предпоследней строки диапазона (т.е. ячейки С14, D 14, Е14 заполнять не следует — для этих ячеек введенные ранее формулы опреде­ления элементарных площадей не работают).

Осталось только вычислить сумму значений в столбцах С, D и Е. Для этого выде­ляем ячейку С15 и вводим туда формулу =СУММ( C 4: C 13).

В ячейку D 15 вводится формула =СУММ( D 4: D 13), а в ячейку Е15 — формула =СУММ( E 4: E 13). Окончательный результат показан на рисунке.

Результат вычислений для всех методов одинаковый.

Абсолютная точность: |2−1,9835| = 0,0165.

Относительная точность: 0,0165 / 1,9835 = 0,0083.

 


                                                                                                                

Вычислить определенный интеграл

Вариант Интеграл Вариант Интеграл
1 11
2 12
3 13
4 14
5 15
6 16
7 17
8 18
9 19
10 20

Контрольные вопросы

1. Какие формулы используются для вычисления производной функции?

2. Какими способами вычисляется определенный интеграл?

 

 


Дата добавления: 2019-11-16; просмотров: 90; Мы поможем в написании вашей работы!

Поделиться с друзьями:




Мы поможем в написании ваших работ!