- •Лабораторная работа № 1 Процедурные типы данных
- •Теоретическая часть
- •Процедурные типы
- •Процедурные переменные
- •Параметры процедурного типа
- •Практическая часть
- •Способ Ньютона (метод касательных)
- •Приближенное решение алгебраических и трансцендентных уравнений
- •Cпособ половинного деления
- •Cпособ пропорциональных частей (метод хорд)
- •Способ итераций
- •Приближенное интегрирование функций
- •Формула прямоугольников
- •Формула трапеций
- •Формула Симпсона
- •Задания Написать программу, реализующую итерационные вычислительные процессы, с использованием процедурных типов данных.
Способ итераций
Способ итераций заключается в следующем. Уравнение f( x ) переписывается в виде x = g( x ). После этого выбирают начальное приближение x1, подставляют его в правую часть уравнения и полученное значение x2 = g( x1 ) принимают за второе приближение для корня. Вообще, если найдено приближение xn, то cледующее приближение xn+1 определяется по формуле xn+1 = g( xn ). Если с заданной степенью точности выполняется равенство xn+1 = xn , то xn и принимают за искомое значение корня.
Приближенное интегрирование функций
Задана предельная абсолютная погрешность . Вычислить интеграл функции f( x ) непрерывной на отрезке [a, b] с точностью . В данной задаче вычисление с точностью означает следующее. Отрезок интегрирования [a, b] разбивается на n равных частей, а функцию f( x ) на рассматриваемом отрезке заменяют интерполирующей или аппроксимирующей функцией g( x ) простого типа и строится сумма Sni, которая является приближенным значением интеграла. Если | Sni+1 - Sni | < , Sni считается значением интеграла с точностью . (Здесь ni < ni+1 ( i = 1, 2, ... )).
Формула прямоугольников
Формула прямоугольников для приближенного вычисления определенного интеграла имеет вид
S = h*(y0 + y1 + ... + yn-1),
где h = ( b-a )/n,
xj = a + j*h ( j = 0, 1, 2, ..., n ),
yj = f( xj ).
Формула трапеций
Формула трапеций для приближенного вычисления определенного интеграла имеет вид
S = h/2*( y0 + 2*( y1 + y2 + ... + yn-1 ) + yn ),
где h = ( b-a )/n,
xj = a+ j*h ( j = 0, 1, 2, ..., n ),
yj = f( xj ).
Формула Симпсона
Формула Симпсона для приближенного вычисления определенного интеграла имеет вид
S = h/3*(( y0 + y2*k ) + 4*( y1 + y3 + ... + y2*k-1 ) + 2*( y2 + y4 + ... + y2*k-2 )),
где n - четно ( 2*k ),
h = ( b-a )/n,
xj = a+ j*h ( j = 0, 1, 2, ..., n ),
yj = f( xj ).
Задания Написать программу, реализующую итерационные вычислительные процессы, с использованием процедурных типов данных.
Варианты 1, 8, 15, 22. По способу половинного деления решить уравнения:
x2 = sin(x).
x3 = sin(x).
x = arcsin((x+1)/4).
Варианты 2, 9, 16, 23. По способу хорд решить уравнения:
x3 - 5*x + 1 = 0.
x3 - 9* x2 + 20*x - 1 = 0.
x3 - 3* x2 - 3*x + 10 = 0.
Варианты 3, 10, 17, 24. По способу касательных решить уравнения:
x5 + 5*x + 1 = 0.
sin( x ) + x = 1.
x2 - 10*lg( x ) - 3 = 0
Варианты 4, 11, 18, 25. По способу итераций решить уравнения:
x2 = ln( x+1 ).
ln( x ) = 4 - x2.
lg( x ) = 0.1*x.
Варианты 5, 12, 19, 26. Вычислить интегралы по формуле прямоугольников:
f(x) = 1/( x(1/3) + 1 ); [0, 1.2].
f(x) = ( 4 + x2 )(1/2); [0, 3].
f(x) = 1/sqrt( 1-0.25*sin( x ) ); [0, pi/4].
Варианты 6, 13, 20, 27. Вычислить интегралы по формуле трапеций:
f(x) = exp( -x )*cos( pi*x/4 ); [0, 2].
f(x) = tg( x )(1/2) ; [0, pi/6];
f(x) = (1 + cos( x2 )(1/2) ); [0, pi].
Вариант 7, 14, 21, 28. Вычислить интегралы по формуле Симпсона:
f(x) = (x2-5) *(10-x); [0, 10].
f(x) = exp(2*x) * sin(2*x2 + 1); [0, 1].
f(x) = arccos(exp(-(3*x)(1/3) ); [0.2, 0.56].
