Математическое моделирование систем и процессов. Тюмиков Д.К. - 20 стр.

UptoLike

Составители: 

Рубрика: 

0 0.5 1 1.5 2 2.5 3 3.5 4
0
2
4
6
8
10
ax()
bx()
cx()
dx()
x
двумя координатами и определяется количество точек m, находящихся под
графиком функции и количество точек N всего сгенерированных по заданному
квадрату. Доказывается, что отношение m/N при достаточно большой выборке
стремится к значению определенного интеграла. Метод применяется при
сложных функциях, для которых трудно или невозможно найти первообразную.
Используемые функции:
Runif(n,a,b) – вектор N случайных чисел, имеющих равномерное распределение
(b, a) – границы интервала, a<b.
Протокол
Пример№1. Введем функции для ограничения области вычисления:
1:)(
2:)(
9:)(
6:)(
4..01.0,0:
2
=
=
=
=
=
xxd
xxc
xxb
xxa
x
Построение графиков заданных функций приведено на рис 8.1. Область
вычисления находится в пределах:
20
5.21
y
x
Подынтегральное выражение:
2
)(:),( yxyxf +=
Рис.8.1. Графики функций a(x),b(x), c(x), d(x)
Число повторений N:= 10000
Определим векторы случайных значений на отрезках:
)16,0,(:
)2,0,(:
)5.2,1,(:
NrunifZ
NrunifY
NrunifX
=
=
=
Функция подсчета:
=
=
=
1
0
),,(*)
1
(*)016(*)02(*)15.2(:
)),((*)1,)(,2(*))((*))((*))((:),,(
N
n
nnn
zyxk
N
I
yxfzyxdxifyxcyxbyxazyxk