40832

ВЫЧИСЛЕНИЕ ПРОИЗВОДНЫХ И ИНТЕГРАЛОВ

Лекция

Математика и математический анализ

При аппроксимации операторов численного дифференцирования и интегрирования наибольшее распространение ввиду своей простоты нашли интерполяционные формулы Ньютона. Формулы численного дифференцирования Формулы для расчета производной в точке x получаются следующим образом. Такие формулы называют простейшими формулами численного дифференцирования.3 получается три важные формулы второго порядка точности: 4.

Русский

2013-10-22

654 KB

24 чел.

ТЕМА 4. ВЫЧИСЛЕНИЕ ПРОИЗВОДНЫХ И ИНТЕГРАЛОВ

4.1. Как аппроксимируют операторы?

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

Аппроксимацию операторов обычно осуществляют следующим образом: функцию  аппроксимируют такой функцией , от которой оператор легко вычисляется, после чего полагают = [].

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

4.2. Формулы численного дифференцирования 

 Формулы для расчета производной  в точке x получаются следующим образом. Берется несколько близких к x узлов  (n m+1), называемых шаблоном (точка x может быть одним из узлов). Вычисляются значения  в узлах шаблона, строится интерполяционный многочлен Ньютона (3.6) и после взятия производной от этого многочлена  получается выражение приближенного значения производной (формула численного дифференцирования) через значения функции в узлах шаблона:.

При n=m+1 формула численного дифференцирования не зависит от положения точки x внутри шаблона, т.к. в этом случае m-я производная от полинома  m-й степени Pm(x) есть константа. Такие формулы называют простейшими формулами численного дифференцирования.

Анализ оценки погрешности вычисления производной

   (4.1)

показывает, что погрешность минимальна для значения x в центре шаблона и возрастает на краях. Поэтому узлы шаблона, если это возможно, выбираются симметрично относительно x. Заметим, что порядок погрешности при h0 равен n-m 1, и для повышения точности можно либо увеличивать n, либо уменьшать h (последнее более предпочтительно).

Приведем несколько важных формул для равномерного шаблона, получаемых с использованием (3.7):

  .   (4.2)

Простейшая формула (4.2) имеет второй порядок погрешности в центре интервала и первый по краям:

 .  (4.3)

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

        (4.4)

      (4.5)

       (4.6)

Для вычисления второй производной часто используют следующую простейшую формулу:

  ,  (4.7)

которая имеет второй порядок погрешности в центральной точке x2.

4.3. Формулы численного интегрирования 

Формулы для вычисления интеграла  получают следующим образом. Область интегрирования [a, b] разбивают на малые отрезки , в общем случае разной длины. Значение интеграла по всей области равно сумме интегралов на отрезках . Выбирают на каждом отрезке [xi, xi+1] 1–5 узлов и строят для каждого отрезка интерполяционный многочлен соответствующего порядка. Вычисляют интеграл от этого многочлена на отрезке. В результате получают выражение интеграла (формулу численного интегрирования) через значения подынтегральной функции в выбранной системе точек. Такие выражения называют квадратурными формулами.

 Приведем наиболее часто используемые квадратурные формулы для равных отрезков длиной

Формула средних

Формула средних получает-ся, если на каждом i-м отрезке взять один центральный узел  соот-ветствующий середине отрезка. Функция на каждом отрезке ап-проксимируется многочленом нулевой степени (константой) Заме-няя площадь криволинейной фигуры площадью прямо-угольника высотой  и основанием h, получим   (рис. 4.1):

  .   (4.8)

Погрешность формулы средних имеет второй порядок по h:

      (4.9)

Формула трапеций

Формула трапеций получается при аппроксимации функции f(x) на каждом отрезке [xi, xi+1] интерполяционным многочленом первого порядка, т.е. прямой, проходящей через точки , . Площадь криволинейной фигуры заменяется площадью трапеции с основаниями  и высотой h (рис.4.2):

. (4.10)

Погрешность формулы трапеций в два раза больше, чем погрешность формулы средних:

  .    (4.11)

Формула Симпсона

Формула Симпсона получается при аппроксимации функции f(x) на каждом отрезке [xi, xi+1] интерполяционным многочленом второго порядка (параболой) c узлами  . После интегрирования параболы получаем (рис.4.3)

   (4.12)

После приведения подобных членов формула (4.12) приобретает удобный для программирования вид:

Погрешность формулы Симпсона имеет четвертый порядок по h:

      (4.13)

Схема с автоматическим выбором шага  по заданной точности

Анализ формул (4.8),(4.10),(4.12) показывает, что точное значение интеграла находится между значениями  и, при этом имеет место соотношение

         (4.14)

Это соотношение часто используется для контроля погрешности вычислений. Расчет начинается с m = 2 и производится по двум методам, в результате получают . Если  – заданная погрешность), то шаг h уменьшают вдвое (m = m2) и расчет повторяют. Если точность достигается, то окончательное значение интеграла получается по формуле (4.14). При существенном уменьшении шага h начинают сказываться ошибки округления, поэтому шаг должен быть ограничен снизу некоторой величиной, зависящей от разрядной сетки ЭВМ (m  n).

   Схема алгоритма с автоматическим выбором шага представлена на рис. 4.4.


Формулы Гаусса

При построении предыдущих формул в качестве узлов интерполяционного многочлена выбирались середины и (или) концы интервала разбиения. При этом оказывается, что увеличение количества узлов не всегда приводит к уменьшению погрешности (сравни формулы средних и трапеций), т.е. за счет удачного расположения узлов можно значительно увеличить точность. Суть методов Гаусса порядка n состоит в таком расположении n узлов интерполяционного многочлена на отрезке [xi, xi+1], при котором достигается минимум погрешности квадратурной формулы. Детальный анализ показывает, что узлами, удовлетворяющими такому условию, являются нули ортогональнoго многочлена Лежандра степени n (см. подразд. 3.1). Так, для n = 1 один узел должен быть выбран в центре. Следовательно, метод средних является методом Гаусса порядка 1.

Для n=2 узлы должны быть выбраны следующим образом:

,

и соответствующая формула Гаусса 2-го порядка имеет вид

          .    (4.15)

Порядок погрешности этой формулы при h 0 равен p = 4, т.е. такой же, как у метода Симпсона, хотя используется только 2 узла!

Для n = 3 узлы выбираются следующим образом:

,

и соответствующая формула Гаусса 3-го порядка имеет вид

         .   (4.16)

Порядок погрешности этой формулы при h 0 равен p = 7, т.е. значительно выше, чем у формулы Симпсона, практически при одинаковых затратах на вычисления. Следует отметить, что формулы Гаусса, особенно широко применяются для вычисления несобственных интегралов специального вида, когда подынтегральная функция имеет достаточно высокие производные.

4.4. Варианты заданий

Для каждого варианта задан интервал [a, b], функция f(x) и указан метод вычисления интеграла. Вначале вычислить точные выражения для первой, второй производных  и для интеграла. Затем необходимо составить подпрограмму для вычисления первой и второй производных по формулам (4.4), (4.7) и подпрограмму вычисления интеграла указанным методом. Составить основную программу, которая вычисляет таблицу значений функции, ее

точных и приближенных производных  в точках  а также точное и приближенное значения интеграла. Расчет первой производной в крайних точках a и b выполнить по формулам (4.5) и (4.6) соответственно, вторую производную вычислять только во внутренних точках. Схема выполнения расчетов приведена на рис. 4.5.

Расчеты производной произвести для hp =0.2, 0.1 и 0.05. Расчеты интеграла произвести для m=10, 20 и 40.

При использовании алгоритма вычисления интеграла с автоматическим выбором шага по данной точности расчет произвести для =0.1, 0.01, 0.001 и получить зависимость m().

Проанализировать погрешность вычислений, для чего построить графики и вычислить погрешности производных и интеграла.

Таблица 4.1

N вар.

Функция f(x)

Интервал

Метод

Значение

a

b

интегрирования

1

–2

3

средних

5.983

2

0

3

трапеций

6.699

3

1

8

Симпсона

8.896

4

4

7

автомат

6.118

5

5

8

Гаусса 2

6.067

6

3

6

Гаусса 3

3.367

7

1

4

средних

0.100

8

0

4

трапеций

0.153

9

–8

2

Симпсона

713.3

10

–2

5

автомат

69.42

11

–5

3

Гаусса 2

167.6

12

–1

4

Гаусса 3

22.09

13

1

7

средних

3.533

14

–2

5

автомат

154.73

15

–4

2

Симпсона

20.375

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

1. Как производится аппроксимация оператора общего вида?

2. Получите формулы (4.2), (4.3) и (4.7).

3. Дайте геометрическую интерпретацию методов: средних, трапеций, Симпсона. Какой порядок погрешности они имеют?

4. В чем суть метода Гаусса?

PAGE  39


 

А также другие работы, которые могут Вас заинтересовать

35350. Организация видеосвязи. Видеоконференцсвязь 1.56 MB
  Видеоконференцсвязь применяется для личного и делового общения, проведения совещаний, конференций и других важных мероприятий, в которых важно присутствие тех субъектов общения, которые по тем или иным причинам не могут находиться в месте проведения мероприятия лично.
35351. Организация аудио связи по локальной сети между рабочими станциями в дуплексном режиме 584.5 KB
  Аудиоконференцсвязь — область информационной технологии, обеспечивающая одновременно двустороннюю передачу, обработку, преобразование и представление аудиоинформации на расстояние в режиме реального времени с помощью аппаратно-программных средств вычислительной техники.
35352. Манипуляционные роботы 778 KB
  История механики богата примерами, которые свидетельствуют о постоянном стремлении человека создать механизмы и устройства, подобные живым существам. Это стремление обусловлено многими причинами...
35353. Создание серверной и клиентской частей системы проведения аудиоконференций с возможностью работы в локальной или глобальной сети 890.5 KB
  Серверная часть должна обеспечивать эффективную работу с клиентами в условиях одновременного подключения нескольких пользователей, а так же обладать дружественным интерфейсом и отвечать следующим основным требованиям...
35354. Управління теками, файлами і ярликами 310 KB
  9 створити в теці Петренко ярлик антивірусної програми для перевірки дискети і запустити її на виконання 10 додати в меню Пуск пункт для запуску антивірусної програми; 1 перетянуть ярлык в главное меню; 2 Правой кнопкой мыши вызвать контекстное меню выбрать там закрепить в меню Пуск. Як створити теку на робочому столі Як створити теку в теці Як скопіювати зміст теки розміщеної на диску D: на робочий стіл' Як виділити групу суміжних об'єктів Як виділити групу несуміжних об'єктів Як виділити всі об'єкти...
35355. Управління папками, файлами і ярликами 109.5 KB
  Індивідуальне завдання Для того щоб освоїти прийоми роботи з теками і файлами необхідно виконати наступне: 1 відкрити вікно папки диска D: і створити в ній скажімо папку Petrenko букви латинські; открываем иконку с названием мой компьютер в окне находим иконку диска С: и отрываем ее в окне диска С: создаем новую папку даем ей имя латинскими буквами Petrenko 2 перейменувати папку Petrenko в папку Петренко букви кирилиці; Левой клавишей мыши выделяем папку под названием Petrenko и выбираем из предоставленного списка...
35358. Тема: Використовування програми Провідник Мета: придбати уміння і навик роботи з програмою Провідник. 291.5 KB
  1 відкрити вікно теки диска D: і створити в ній скажімо теку Petrenko букви латинські; в левой части окна выбираем диск Д в строке меню выбираем пункт файл создать папку 2 перейменувати теку Petrenko в теку Петренко букви кирилиці; в правой части окна выделяем папку Petrenko в строке меню выбираем пункт файлпереименовать 3 з теки диска D: скопіювати в теку Петренко якийнебудь файл або теку; в правой части окна выбираем необходимый файл выделяем его в строке меню...