НИЖЕГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ 
 

Кафедра: «Высшая математика»

Курсовая  работа

 

по  информатике 
 

Тема:

Вычисление определенного интеграла с заданной точностью методом Ньютона-Котеса (n=7) 

Вариант №  4 
 

      Выполнил:  студент группы 10-А3

Ераксин О.В.

                            

                                      Проверил:    Толкачёв И.Н.  
 

Нижний  Новгород

2011

Содержание

Стр.

    Введение…….……………………………………..…………………………………3

1. Постановка  задачи   ……………………………..…………………………………..6

1.1. Вычисление определенного интеграла с заданной точностью методом    Ньютона-Котеса ………………………………………………………………………..6

2. Методика решения  задачи ..……………………………………………………….10

2.1. Получение разрешающего уравнения ...…….………..…………………….......10

         2.2 Вычислительная схема метода Ньютона-Котеса ....………………….......11

3. Блок-схемы  алгоритма……………………………………………………………..12

         3.1. Блок-схема основной программы………………………………………...12

         3.2. Блок-схема процедуры……………………………………………….........14                          

         3.3. Блок-схема функции ……….........................................................................15

4. Программа  на языке Pascal………………………………………………………...16

4.1. Основная  программа..............................................................................................16

4.2. Процедура................................................................................................................19

5.Выводы....……………………………….…………………………………...……....206.Литература………………………………………………………………………......22 
 
 
 
 
 
 
 
 
 
 

Введение.

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

Методы  численного интегрирования классифицируются в зависимости от того, заданы ли значения аргумента через равные интервалы или нет. Краткий обзор некоторых из них.

Формулы Ньютона — Котеса получаются путем замены подынтегральной функции интерполяционным многочленом Лагранжа с разбиением отрезка интегрирования на n равных частей. Получающиеся формулы используют значения подынтегральной функции в узлах интерполяции и являются точными для всех многочленов некоторой степени, зависящей от числа узлов. Точность формул растет с увеличением степени интерполяционного многочлена. Заметим, что формулы прямоугольников, трапеций и Симпсона являются частными случаями формул Ньютона — Котеса.

Метод Гаусса не предполагает разбиения отрезка интегрирования на равные промежутки. Формулы численного интегрирования интерполяционного типа ищутся такими, чтобы они обладали наивысшим порядком точности при заданном числе узлов. Узлы и коэффициенты формул численного интегрирования находятся из условий обращения в нуль их остаточных членов для всех многочленов максимально высокой степени.

Формула Эрмита, являющаяся частным случаем формул Гаусса, использует многочлены Чебышева для вычисления интегралов вида

Получающаяся  формула характерна тем, что все  коэффициенты при  равны.

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

Формула Чебышева представляет интеграл в виде

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

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

При вычислении определенного интеграла, где - непрерывная на отрезке функция, иногда удается воспользоваться известной формулой Ньютона - Лейбница:

  (1)

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

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

    (2)

Подобный  подход удобен тем, что он приводит к алгоритмам, 
легко реализуемым на ЭВМ и позволяющим получать результат с достаточной точностью. При этом, естественно, предполагается, что отрезок разбит на n частей точками, наличие которых подразумевается при построении многочлена . 
 
 
 
 
 
 
 
 
 
 
 
 
 

Постановка  задачи.

Вычисление  определенного интеграла  с заданной точностью  методом Ньютона-Котеса.

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

,

 получим:

Таким образом:

(3)

где

(4)

По поводу полученных формул можно заметить,  что:

1)  коэффициенты  не зависят от функции , так как они составлены только с учетом узлов интерполяции;

2) если  - полином степени n, то тогда формула (2) точная, ибо в этом случае .

Как уже  отмечалось выше, применение формулы (2) предполагает построение на отрезке интегрирования системы узлов интерполяции (здесь ), которыми отрезок делится на n частей. Длина называется при этом шагом интегрирования. Естественно считать, 
что шаг h  постоянен, т. е. . В этом случае можно 
применить интерполяционную формулу Лагранжа для равноотстоящих узлов. Итак, с учетом

и

формула (5.21) для «весовых» коэффициентов

, принимает вид:   .  (5)

Перейдем  в этом интеграле всюду к переменной t. Из подстановки   получаем:   ,  т. е. .

При имеем ,  а при  будет . Тогда

, (6)

где

.    (7)

Числа (5) называют коэффициентами Котеса. Как видно, коэффициенты Котеса не зависят от функции , а только от числа n точек разбиения. С учетом формул (3) и (6), получаем следующий вид квадратурных формул Ньютона - Котеса:

  (8)

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

Последний член в формуле Ньютона — Котеса:

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

k— степень используемых многочленов, x*— некоторая точка в интервале , - значение k-й производной функции в точке x* (в этой точке принимает максимальное значение).

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

Коэффициенты  в формуле Ньютона - Котеса