ВУЗ:
Составители:
Рубрика:
70
8.4. Исследовать зависимость относительной погрешности
численного дифференцирования от шага разностной сетки для
разностной производной второго порядка на трехточечном шаблоне ,
оценивая относительную погрешность как
()
%100⋅
′
′
−
′
=∆
теор
теорразн
y
yy
h :
а) функции xxxy 7)(
2
+= при
{
}
01,0;02,0;05,0;1,0;3,0;5,0
=
h в точке х=4;
б) функции xxxy 2sin)(
2
= при
{
}
01,0;02,0;05,0;1,0;2,0;4,0
=
h в точке х=1;
в) функции
2
2
7
)(
x
xx
xy
+
= при
{
}
01,0;02,0;05,0;1,0;3,0;5,0
=
h в точке х=2;
г) функции
xxxxy 34)(
23
−+=
при
{
}
025,0;05,0;1,0;25,0;5,0
=
h
в точке х=0.
Построить графики зависимостей
(
)
)( hfh
=
∆
.
8.5. На интервале [0;5] таблично задана функция )( xy :
x
i
0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0
y
i
0,1 0,7 0,9 1,2 1,5 1,6 1,7 1,9 2,0 2,1 2,2
Рассчитать и построить графики всех разностных производных
первого и второго порядков.
Сравнить результаты численного дифференцирования с графиками
теоретических производных
()
x
xy
2
1
=
′
и
()
3
4
1
x
xy −=
′′
.
8.2. Численное интегрирование
В большинстве научных и технических задач нахождение
определенного интеграла от элементарных функций , или от функций ,
заданных таблично , не может быть сведено к аналитическим выражениям .
В этих случаях применяется приближенное интегрирование .
Пусть надо вычислить определенный интеграл
∫
=
b
a
dxxfI )(
для
вещественной функции f(x), которая определена и ограничена на
замкнутом интервале [a,b]. Согласно геометрическому смыслу,
определенный интеграл I есть площадь, ограниченная кривой f(x), осью 0х
и прямыми х=а и х=b (рис.8.2).
Разбив интервал [a,b] на n
частей, найдем интегральную
сумму S:
()
∑
=
+
−=
n
i
iii
xxfS
1
1
)( ξ
,
которая при стремлении
интервалов
iii
xxx −=∆
+ 1
к нулю
для непрерывных функций
стремится к точному значению
интеграла:
IS
i
i
x
=
→∆ 0max
lim
.
x
x
0
=
a
x
1
x
2
x
i
x
i+1
ξ
i
x
n-1
x
n
=
b
y
f(x)
f(ξ
i
)
Рис. 8.2.
70 8.4. Исследовать зависимость относительной погрешности численного дифференцирования от шага разностной сетки для разностной производной второго порядка на трехточечном шаблоне, y′разн −yтеор ′ оценивая относительную погрешность как ∆(h ) = ⋅100% : ′ yтеор а) функции y ( x) = x 2 +7 x при h ={0,5; 0,3; 0,1; 0,05; 0,02; 0,01} в точке х=4; б) функции y ( x) =x 2 sin 2 x при h ={0,4; 0,2; 0,1; 0,05; 0,02; 0,01} в точке х=1; x 2 +7 x в) функции y ( x) = при h ={0,5; 0,3; 0,1; 0,05; 0,02; 0,01} в точке х=2; x2 г) функции y ( x) =x 3 +4 x 2 −3x при h ={0,5; 0,25; 0,1; 0,05; 0,025} в точке х=0. Построить графики зависимостей ∆(h ) = f (h) . 8.5. На интервале [0;5] таблично задана функция y(x) : xi 0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 yi 0,1 0,7 0,9 1,2 1,5 1,6 1,7 1,9 2,0 2,1 2,2 Рассчитать и построить графики всех разностных производных первого и второго порядков. Сравнить результаты численного дифференцирования с графиками 1 1 теоретических производных y ′(x ) = и y ′′(x ) =− . 2 x 4 x3 8.2. Численное интегрирование В большинстве научных и технических задач нахождение определенного интеграла от элементарных функций, или от функций, заданных таблично, не может быть сведено к аналитическим выражениям. В этих случаях применяется приближенное интегрирование. b Пусть надо вычислить определенный интеграл I = f ( x )dx для ∫ a вещественной функции f(x), которая определена и ограничена на замкнутом интервале [a,b]. Согласно геометрическому смыслу, определенный интеграл I есть площадь, ограниченная кривой f(x), осью 0х и прямыми х=а и х=b (рис.8.2). Разбив интервал [a,b] на n y частей, найдем интегральную f(ξi) n сумму S: S =∑ f (ξi )( xi +1 −xi ) , i =1 f(x) которая при стремлении интервалов i ∆ x = x i +1 −x i к нулю для непрерывных функций x0= a x1 x2 xi ξi xi+1 xn-1 xn= b x стремится к точному значению интеграла: lim S =I . Рис. 8.2. max ∆xi → 0 i
Страницы
- « первая
- ‹ предыдущая
- …
- 68
- 69
- 70
- 71
- 72
- …
- следующая ›
- последняя »