Блок-схема итерационного алгоритма метода расслоенной выборки.

.Если для области Gk выполним Nk опытов, получим оценку математического ожидания искомого показателя для данной области:

. (11)

Результирующая оценка должна рассматриваться как случайная дискретная величина, значения которой наблюдаются с вероятностями pk . Тогда результирующая оценка определяется усреднением:

. (12)

.Определим дисперсию оценки (9), имея в виду, что все N1+ N2+ N3+…+ N10 слагаемые - независимые случайные величины:

. (13)

Дисперсия случайной величины может быть оценена следующим образом:

(14)

3. Введя в рассмотрение доли от общего количества опытов, соответствующие областям Gk,

,

на основе (10) получим соотношение для определения количества опытов, необходимого для получения результата с погрешностью не выше :

(15)

При удачном разбиении области G и удачном выборе соотношения количества опытов для отдельных областей Gk дисперсия оценки (13) может быть существенно снижена. Оптимальные значения должны быть пропорциональны произведениям

Провели начальную серию опытов N=200. После проведения данной серии опытов были получены следующие результаты:

· Оценка математического ожидания для каждой из 10 областей на основании (11):

· Результирующая оценка математического ожидания по (12):

· Дисперсия для каждой из 10 областей по (14):

· Дисперсия оценки математического ожидания по (13):

· Требуемое количество опытов, рассчитанное по (15):

опытов.

Алгоритм повторялся до тех пор, пока не выполнилось условие . Данное условие выполнилось после третьей итерации алгоритма.

После второй итерации получили:

· N = 1982 опытов.

· Оценка математического ожидания для каждой из 10 областей:

· Результирующая оценка математического ожидания:

· Дисперсия для каждой из 10 областей:

· Дисперсия оценки математического ожидания:

(

· Требуемое количество опытов:

опытов.

После третьей итерации алгоритма:

· N = 2191 опытов.

· Оценка математического ожидания для каждой из 10 областей:

· Результирующая оценка математического ожидания:

· Дисперсия для каждой из 10 областей:

· Дисперсия оценки математического ожидания:

· Требуемое количество опытов:

опытов.

Дифференциальное уравнение (1) решается численным интегрированием методом Эйлера первого порядка [4] с шагом 0.001. Программа, реализующая данный метод снижения трудоемкости, написана на языке Delphi 7 [5].

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