Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
шпора информатика.doc
Скачиваний:
13
Добавлен:
24.09.2019
Размер:
689.66 Кб
Скачать

31. Вычисл-е определ-го интеграла по одной из фор-л. Алг-м реализ-ии на эвм выбранной формулы.

Форм-лы для приближ-го вычисл-я интегралов примен-ся очень часто. Дело в том, что для большого числа элементар ф-ий первообразные уже не выражаются через элементарные функции, в результате чего нельзя вычислить определённый интеграл с помощью формулы Ньютона – Лейбница:

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

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

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

достаточно мала, то интегральная сумма должна уже мало отличаться от величины интеграла. Таким образом, мы получаем приближённое равенство (форм.1):

которое и является формулой прямоугольников. Здесь через y0, y1,…, yn обозначены значения функции y = f(x) в точках деления x0, x1,…, xn .

Аналогичная формула прямоугольников получится и в том случае, если брать для интегральной суммы значение функции не в левых, а в правых концах участков разбиения. Тогда формула примет вид (форм. 2):

Для функции, монотонной на отрезке интегрирования, всякая интегральная сумма, а значит, и определённый интеграл заключены между приближёнными значениями, указанными в правых частях формул 1 и 2. Геометрическую иллюстрацию этого факта можно видеть на рисунке:

Слагаемые, входящие в различные суммы, показаны пунктиров и штриховкой. Благодаря этому представление о погрешностях формулы прямоугольников можно получить, рассматривая разность результатов, полученных по формулам 1 и 2.

Если функция имеет на отрезке интегрирования конечное число экстремумов, то отрезки интегрирования можно разбить на участки монотонности и получить, таким образом, оценку погрешности формулы прямоугольников.

Алг-тм:

1) вв. a, b, и n, f(x);

2) ; I:=0; - знач-е интеграла

3) пока a<=x<=b повторять (формула левых прямоуголь-ков)

I:=I+h*f(x);

x:=x+h;конец цикла

4) вывести зн-е I.

32,Приближенные методы решения обыкновенных дифференциальных уравнений первого порядка: постановка задачи, один из методов решения (Эйлера, Рунге-Кутта). Алгоритм для реализации на ЭВМ выбранного метода.

Простейшее диф Ур-ние: (1)

Решением Ур (1) назыв такая фор-ла , (x=x0, y=y0)- начальные условия.

Задача нахождения ф-и есть задача Коши

Однако не смотря на это методы нахождения найти ее не всегда удается. Поэтому возникает необходимость в знакомстве с методами кот позвол найти прбл-ое решение Ур-я(1). Это прибл реш-е м быть выражено в разл формах и в зависим от этого часто условно прибл методы м выразить 3 вариантами:

1)аналитич-ие методы кот позвол представить приближен реш-е диф урав-я в виде некот аналитич выражения

2)графич прибл методы, кот-е позвол выразить реш-е дифф урав-я в виде графика.

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

Метод Рунге-Кутта:

Много общего с методом Эйлера.

, где x=x0,y=y0. [a,b];h

yi,(i=0,1,2..n)- точное значение в т. xi[0,b]

То нахождение приближ-ий закл-ся

, где

, (i=0,1,2,3…n-1)

При построении таблицы знач реш-я данного диф Ур-я методом Рунге-Кутта достаточно вносить вычисления значений:

( , , , ) Точность реш-я диф Ур-я методом Р-К опред-ся по фор-ле: ,

где y(xn)-точное знач-е реш-я y(x) в т. xn

y*n – это приближенное значение реш-я в т. Xn и вычеслен в шаге h/2

yn- приближенное значение реш вычислен при шаге h

В некоторых задачах задается точность  с кот нужно получить реш-е

Если  задания точность решения то n- число точек разб-ния [a,b] выбирается так, чтобы шаг h=(b-a)/h удовлетворял условию h4<

Алгоритм

Алг<r-k>

1 шаг) вводим а,b,h,y

2 шаг) опишем функцию f(x,y) {ввод f(x)}

3 шаг) x:=a; y:=y0;

4 шаг) выведем «х0 = …» «у0 = …»

5 шаг) вычислим последовательность уi

i:=1; Организуем цикл пока 1<=i<=n повтор

k1 = h*f(x,y);

k2(h) = h*f(x+h/2,y+k1/2);

k3(h) = h*f(x+h/2,y+k2/2);

k4(h) = h*f(x+h,y+k3)

d:=(k1+2k2+2k3+k4)/6; y:=y+d; x:=x+h;

вывод хi и yi

end;

6 шаг) конец программы

6