Программная реализация итерационных методов

3.5 Программная реализация итерационных методов

Рисунок 12. Решение системы уравнений методом простых итераций

Рисунок 13. Решение уравнения методом Зейделя

Раздел 4. Сравнительный анализ методов численного дифференцирования и интегрирования

4.1 Методы численного дифференцирования

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

Предположим, что в окрестности точки xi функция F(x)дифференцируема достаточное число раз. Исходя из определения производной:

используем для её вычисления две приближенные формулы:

 (1)

(2)

Формулы (1) и (2) называют правыми и левыми разностными производными.

Для оценки погрешностей формул численного дифференцирования используется формула Тейлора:

откуда можно вычислить:

 (3)

Выражение (3) имеет погрешность порядка (x-xi), следовательно, формулы правых и левых разностных производных имеют погрешность одного порядка с h , где

h=xi-xi-1

Такая точность достаточно невысока, поэтому применяется так называемая центрально-симметричная форма производной, погрешность которой одного порядка с h2

 (4)

Хотя очевидно, что формула (4) используется для внутренних точек отрезка.

Для примера возьмём ряд точек:

Вычислим производную функции f(x)=sin(x) в одной из них двумя способами.

Очевидно, что h=

По центрально-симметричной формуле:

По формуле левой разностной производной:

Табличное значение =cos()=0.8660 ,т.е. значение производной, полученное по центрально-симметричной формуле ближе к истинному.

4.2 Методы численного интегрирования

В прикладных исследованиях часто возникает необходимость вычисления значения определённого интеграла

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

Суть приближенного вычисления заключается в двух операциях:

1. в выборе конечного числа вместо n

2. в выборе точки  в соответствующем отрезке.

В зависимости от выбора  мы получаем различные формулы для вычисления интеграла:

Формулы левых и правых прямоугольников (5),(6)

 (5)

 (6)

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

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

где

m=n/2

h=b-a/n

b, a- концы рассматриваемого отрезка.

Для сравнения результатов вычисления вышеизложенными формулами численного интегрирования вычислим 3-мя способами следующий интеграл, разделив отрезок [0, ] на 6 равных отрезков:

 h=

По формуле левых прямоугольников:

По формуле трапеции:

По формуле Симпсона:

А результат полученный аналитически равен

=1

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

Раздел 5. Численные методы решения обыкновенных дифференциальных уравнений

Инженеру часто приходится иметь дело с техническими системами и технологическими процессами, характеристики которых непрерывным образом меняются во времени. Соответствующие явления как правило подчиняются физическим законам, которые формулируются в виде дифференциальных уравнений. Одной из основных математических задач, которые приходится решать для таких уравнений, является задача Коши (начальная задача). Чаще всего к ней приходят тогда, когда известно начальное состояние физической величины системы в некоторый момент времени t0 (x0,y0) и требуется предсказать её поведение в момент времени t>t0 ( x>x0). В курсе математического обеспечения САПР, мы рассматривали методы решения задачи Коши с помощью решения обыкновенного дифференциального уравнения первого порядка .

Напомним, что решением обыкновенного дифференциального уравнения (ОДУ) первого порядка является функция y , которая при подстановке в уравнение , превращает его в тождество.

Обыкновенными дифференциальными уравнениями называются такие уравнения, которые содержат одну или несколько производных от искомой функции y=y(x). Их можно записать в виде

,

где х – независимая переменная.

Наивысший порядок n входящей в уравнение производной называется порядком дифференциального уравнения.

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

Графические методы используют геометрические построения.

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

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

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

Похожие работы

Численные методы линейной алгебры

Скачать

18618

0

16

... уравнений (2) сводится к последовательному решению двух следующих систем уравнений с треугольными матрицами коэффициентов L Y = B; (6) U X = Y (7) линейный алгебраический уравнение численный где Y =  - вектор вспомогательных переменных. Такой подход позволяет многократно решать системы линейных ...

Сравнительный анализ методов оптимизации

Скачать

42464

5

31

... 4 - график унимодальной, но не выпуклой функции Таким образом, кроме перечисленных свойств, выпуклые функции обладают также и всеми свойствами унимодальных функций. 2. Прямые методы безусловной оптимизации Для решения задачи минимизации функции f (х) на отрезке [а; b] на практике, как правило, применяют приближенные методы. Они позволяют найти решение этой задачи с необходимой точностью ...

Численное интегрирование функций

Скачать

7913

0

3

... – остаточный член, характеризующий погрешность формулы. Заметим, что формулы вида (2) называют квадратурными формулами. Геометрический смысл численного интегрирования состоит в вычислении площади криволинейной трапеции, ограниченной графиком функции f(х), осью абсцисс и двумя прямыми х = а и х = b. Приближенное вычисление площади приводит к отбрасыванию в квадратурных формулах остаточного члена ...

Автоматизированный анализ проектирования на микроуровне

Скачать

47503

2

14

... задачи, а именно: 1. Создана расчетная схема анализа на основании сравнительного анализа численных методов, а также программных и технических средств их осуществления; 2. Создан выбор метода автоматизированного анализа объекта проектирования; 3. Спланирован и проведен эксперимент, анализируя результаты которого, приходим к выводу, что данная модель может использоваться с параметрами: r = 5 R = ...