рефераты конспекты курсовые дипломные лекции шпоры

Реферат Курсовая Конспект

Численное интегрирование

Численное интегрирование - раздел Математика, ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА Определенным Интегралом Функции F(X), Взятом В Интервале От A До B, Называ...

Определенным интегралом функции f(x), взятом в интервале от a до b, называется предел, к которому стремится интегральная сумма при стремлении всех промежутков ∆xi к нулю:

.

При приближенном вычислении определенного интеграла шаг интегрирования h=∆x выбирается конечным: , где Ii - элемент интегральной суммы. Заменяя подынтегральную функцию на каждом шаге отрезками линий нулевого, первого и второго порядков, получаем приближенные формулы для вычисления интеграла методами прямоугольников, трапеций и Симпсона соответственно.

1. Приближенные формулы для вычисления интеграла методом прямоугольников (рис 12.1).

Рис. 12.1 - вычисление интеграла методом прямоугольников

Правило прямоугольников (n=0). Заменяем график функции F(x) горизонтальной линией (линий нулевого порядка) и вычисляем значение элемента интегральной суммы как площадь прямоугольника

, где h - шаг интегрирования, у0 - значение функции в точке х=х0

у(х0)=у0


2. Приближенные формулы для вычисления интеграла методом трапеций (рис. 12.2).

Рис. 12.2. - вычисление интеграла методом трапеций

Правило трапеций (n=1). Заменяем график функции F(x) прямой, проходящей через две точки (х00) и (х0+h,у1), и вычисляем значение элемента интегральной суммы как площадь трапеции

3. Приближенные формулы для вычисления интеграла методом Симпсона (рис. 12.3).

Рис. 12.3. - вычисление интеграла методом Симпсона

Правило Симпсона (n=2). Заменяем график функции F(x) квадратичной параболой, проходящей через три точки с координатами (х00), (х0+h,у1), (х0+2h,у2). Расчетную формулу для вычисления элемента интегральной суммы получим, используя интерполяционный многочлен Лагранжа, в виде: y(x)=y0A0(x)+y1A1(x)+y2A2(x), где:

При x0=0; x1=h; x2=2h, получим:

При интегрировании на отрезке [a,b] расчетные формулы для методов прямоугольника, трапеций и Симпсона имеют вид:

,

где h - шаг по x, fa, fi, fb - значения функции при x равном a, xi, b соответственно.

Для метода прямоугольников приведены две расчетные формулы, так как площадь прямоугольника на каждом шаге интегрирования может определяться по левой или правой стороне. Суть метода прямоугольников для отрезка [a, b] проиллюстрирована на рисунке, при этом площадь под кривой f(x) (вспомните геометрический смысл определенного интеграла) заменена суммой площадей заштрихованных прямоугольников (рис. 12.4).

Рис. 12.4 - суть метода прямоугольников для отрезка [a, b]

 

Задание для самостоятельной работы

1. Вычислить значение определенного интеграла аналитически и численно четырьмя методами для пяти значений N, где N – число разбиений интервала интегрирования N=10; 20; 50; 100; 1000. Результаты расчета вывести на экран и распечатать в виде таблицы. Представление результатов расчета:

N Аналит. Значение Метод прямоуг. 1 Метод прямоуг. 2 Метод трапеций Метод Симпсона
         
         
         
         
         

 


 

2. Построить графики функций I=F(N).

Варианты интегралов приведены в таблице:

Вар. Вид интеграла Вар. Вид интеграла

 


– Конец работы –

Эта тема принадлежит разделу:

ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА

C А СИНЮТИН... ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА... УДК A...

Если Вам нужно дополнительный материал на эту тему, или Вы не нашли то, что искали, рекомендуем воспользоваться поиском по нашей базе работ: Численное интегрирование

Что будем делать с полученным материалом:

Если этот материал оказался полезным ля Вас, Вы можете сохранить его на свою страничку в социальных сетях:

Все темы данного раздела:

Таганрог 2007
УДК 681.326.3 – 181.48(075.8)   Синютин С.А. Вычислительная математика. Учебное пособие по дисциплине «Вычислительная математика» – Таганрог: Изд-во ТТИ ЮФУ, 2007. – 123

Ввод матриц
Вы можете вводить матрицы в MATLAB несколькими способами: • вводить полный список элементов; • загружать матрицы из внешних файлов; • генерировать матрицы, используя встр

Индексы
Элемент в строке i и столбце j матрицы А обозначается A(i,j). Например, А(4,2) - это число в четвертой строке и втором столбце. Таким образом, можно вычислить сумму элементов в четвер­том столбце м

Функции
MATLAB предоставляет большое количество элементарных математических функций, таких как abs, sqrt, exp, sin. Вычисление квадратного корня или логарифма отрицательного числа не является ошибко

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

Создание графика
Функция plot имеет различные формы, связанные с входными параметрами, например plot(y) создает кусочно-линейный график зависимости элементов у от их индексов. Если вы задаете два вект

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

Печать графики
Опция Printв меню Fileи команда print печатают графику MATLAB. Меню Printвызывает диалоговое окно, которое позволяет выбирать общие стандар

Команда help
  Команда help – это самый основной способ определения синтаксиса и поведения отдельных функций. Информация отображается прямо в командном окне. Например, help magic

Сценарии и функции
MATLAB - это мощный язык программирования, также как и интерактивная вычислительная среда. Файлы, которые содержат код на языке MATLAB, называются М-файлами. Вы создаете М-файлы, используя текстово

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

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

Теория погрешностей и машинная арифметика.
  Пусть - точное значение, - приближенное значение некоторой величины

Локализация корней
  ПРИМЕР 1. Локализация корней (рис. 4.1). % Локализовать корни уравнения f(x)=0, где f(x)= x^3 - cos(x) + 1 % Введём функцию f(x) f = inline('x.^3 - cos(x)

Метод бисекции
Пусть[a,b] v отрезок локализации. Предположим, что функция f(x) непрерывна на [a,b] и на концах принимает значения разных знаков .

Метод Ньютона (метод касательных)
  Расчетная формула метода Ньютона имеет вид: . Геометрически метод Ньютона означает, что следующее приближение к кор

Метод простой итерации (метод последовательных повторений)
Для применения метода простой итерации следует исходное уравнение преобразовать к виду, удобному для итерации

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

Интервал неопределенности корня.
Если функция непрерывна, то найдется такая малая окрестность корня

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

Методика решения алгебраического уравнения
Мы остановимся здесь подробнее на методике решения алгебраического уравнения, т.е. уравнения вида: , левую часть которого будем обозначать т

Нормы векторов и матриц
Обозначим через - точное решение системы, а через - приближенное решение системы. Д

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

Метод Гаусса
Рассмотрим метод Гаусса (схему единственного деления) решения системы уравнений. Прямой ход состоит из m-1 шагов исключения. 1. Шаг. Исключим неизвестное

LU разложение матрицы
Представим матрицу A в виде произведения нижней треугольной матрицы L и верхней треугольной U.   Введем в рассмотрение матрицы

Метод Гаусса с выбором главного элемента по столбцу
На k-ом шаге прямого хода в качестве ведущего элемента выбирают максимальный по модулю коэффициент при неизвестной

Метод Холецкого
Если матрица системы является симметричной и положительно определенной, то для решения системы применяют метод Холецкого (метод квадратных корней). В основе метода лежит алгоритм специального LU

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

Метод Якоби
Самый простой способ приведения системы к виду удобному для итерации состоит в следующем: из первого уравнения системы выразим неизвестное x1, из второго уравнения системы выраз

Метод Зейделя
Метод можно рассматривать как модификацию метода Якоби. Основная идея состоит в том, что при вычислении очередного (n+1)-го приближения к неизвестному xi при i >1 используют

Метод простой итерации.
Если A - симметричная и положительно определенная матрица, то систему уравнений часто приводят к эквивалентному виду: x = x -

Постановка задачи приближения функции по методу наименьших квадратов
Пусть функция y=f(x) задана таблицей своих значений: , i=0,1,-n. Требуется найти многочлен фиксированной степени m

Определение параметров эмпирической зависимости
Часто из физических соображений следует, что зависимость между величинами хорошо описывается моделью вида

Многочлены Бернштейна
Предположим, что функция задана в отрезке [0,1] в точках , при некотором фиксирован

Постановка задачи интерполяции функций
Пусть функция y = f(x) задана таблицей своих значений: , i=0,1,...n. Требуется найти многочлен степени n

Оценка погрешности интерполяции
Если функция n+1 раз на отрезке [a,b] , содержащем узлы интерполяции , i=0,1,...n, то для погрешности интерполяц

Глобальная и кусочно-полиномиальная интерполяция
Пусть функция f(x) интерполируется на отрезке [a, b]. Метод решения этой задачи с помощью единого многочлена для всего отрезка

Интерполяция сплайнами
Пусть отрезок [a,b] разбит точками на n частичных отрезков . Сплайном степени m называется функция

Первая производная. Двухточечные методы
Для двухточечных методов при вычислении производных используется значение функции в двух точках. Приращение аргумента задается тремя способами, откладывая Δx = h вправо, влево и в обе стороны

Вычисление производных второго порядка
Вторая производная вычисляется как первая производная от первой производной. Для следующей пятиточечной схемы (рис. 11.2):  

Вычисление производных третьего порядка
Производные третьего порядка вычисляются как первая производная от производной второго порядка. Для рассмотренной пятиточечной схемы расчетная формула имеет вид:

Постановка задачи Коши для дифференциального уравнения первого порядка
Пусть дано дифференциальное уравнение первого порядка . Требуется найти функцию , у

Численное решение задачи Коши методом Эйлера
Численное решение задачи Коши состоит в построении таблицы приближенных значений в точках

Оценка погрешности метода Эйлера
Локальной погрешностью метода называется величина . Найдем величину локальной погрешности метода Эйлера:

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

Решение систем дифференциальных уравнений методом Эйлера.
Пусть требуется решить нормальную систему дифференциальных уравнений: с начальными условиями:

Численное интегрирование систем дифференциальных уравнений
Системой дифференциальных уравнений называется система вида где x - независимый аргумент; yi

Аппроксимация функции по Фурье
Пусть функция задана в интервале . В этом случае (при наличии у нее соответствующих

Преобразование Фурье
Так называется действие, с помощью которого по заданной в интервале функции строитс

Быстрое преобразование Фурье.
В предыдущем пункте было описано дискретное преобразование Фурье - сопоставление набору значений функции набора коэффициентов

Уравнение Лапласа (эллиптическое уравнение)
Это уравнение описывает некоторое установившееся стационарное состояние в пространстве (x,y) в некоторой области G (16.1):

Уравнение теплопроводности (параболическое уравнение)
Типичная задача, описывающая теплоперенос вдоль одномерного стержня единичной длинны при заданном начальном распределении температуры от

Хотите получать на электронную почту самые свежие новости?
Education Insider Sample
Подпишитесь на Нашу рассылку
Наша политика приватности обеспечивает 100% безопасность и анонимность Ваших E-Mail
Реклама
Соответствующий теме материал
  • Похожее
  • Популярное
  • Облако тегов
  • Здесь
  • Временно
  • Пусто
Теги