Ковариация

.

Достаточным условием неотрицательности ковариации, дающим выигрыш в оценке разности средних, является одинаковая упорядоченность условных средних , относительно векторного аргумента , т.е. выполнение неравенства

(7.1)

для любых значений векторных аргументов , .

Действительно, учитывая, что, находим

Так как f()>0 для всех v, то при выполнении (7.1)имеем В12³0, что и требовалось доказать.

Когда в качестве результатов моделирования выступают вероятности событий A1, А2 для вариантов S1 и S2 системы, то условные значения

,

,

где условная вероятность, i=1, 2.

Тогда достаточное условие неотрицательности ковариации запишется в виде

³0, (7.2)

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

Одинаково упорядоченными являются монотонно возрастающие или монотонно убывающие функции m1(v) и m2(v) скалярного аргумента v, а также одинаковые функции m1(v)=m2(v) независимо от их монотонности. Пример одинаково упорядоченных возрастающих (а) и убывающих (б) функций m(v) показан на рис. 7.3.

 

 

а б

m m

m1 m2

 

 

m1

 

m2

 

 

n n

Рис.7.3. Пример одинаково упорядоченных функций

 

Если положительные функции mj(), , одинаково упорядочены, то произведение любой комбинации этих функций mk()ms()…mm() одинаково упорядочено с произведением любой комбинации ml()mq()…mp(). Это же можно сказать и об условных вероятностях P(Aj/), .

 

Пример 7.4. Пусть методом статистического моделирования на ЭВМ необходимо сравнить результаты моделирования двух вариантов S1 и S2 системы, составленных из одинаковых блоков B1B4 (структура системы показана на рис. 7.4) и сравниваемых по критерию надежности с учетом случайных изменений внешней температуры. События A1 и A2 соответствуют безотказной работе вариантов S1 и S2 системы в течение заданного времени Т. Вероятность безотказной работы Bi при заданной температуре v можно определять как

,

где li(v) - интенсивности отказов, являющиеся возрастающими функциями температуры.

 

S1 S2

 
 

 

 


Рис. 7.4. Структуры сравниваемых вариантов систем S1 и S2

 

Таким образом, функции P(Bi/v) являются одинаково упорядоченными убывающими функциями. Можно показать, что функции

также одинаково упорядочены и убывают с ростом температуры v. Поэтому, используя при машинном эксперименте с вариантами S1 и S2 системы одни и те же реализации v случайной температуры v, получим в результате моделирования большую точность сравнения вероятностей Р(А1) и Р(А2), чем при раздельном моделировании S1 и S2 системы с использованием независимых реализаций v.

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

Когда независимые компоненты в воздействиях внешней среды Е отсутствуют, т.е. 1=2=, условные средние m1=M[q1/], m2=M[q2/] преобразуются в детерминированные зависимости критериев от случайных воздействий q1=f1(), q2=f2().

При этом условия одинаковой упорядоченности становятся еще более жесткими.

Так, например, условия (7.2) выполняются лишь тогда, когда для всех значений исключено одно из состояний: A12 или 1А2. Другими словами, положительная корреляция B12 и связанные с ней преимущества гарантируются лишь тогда, когда вариант системы S1 равномерно лучше (хуже) варианта S2. В принятых в п. 6.3 обозначениях это соответствует pC = 0 или pD = 0.

Состояния C = A12 или D =1А2 вариантов систем S1 и S2 возможны лишь при наличии двух неисправных блоков Вi, , состояние A = A1A2 возможно при отсутствии неисправностей или при одной неисправности, а состояние В =12 при трех или четырех неисправностях. Обозначив через ij ситуацию с неисправностями блоков Bi и Bj, находим соответствие между состояниями и убеждаемся в отсутствии состояния D.

 

Следует помнить, что условия одинаковой упорядочен-
ности (7.1) и (7.2) являются достаточными, но не необходимыми и достаточными условиями неотрицательности корреляции. Поэтому, обнаружив в конкретной схеме проведения имитационного эксперимента нарушениеэтих условий при некоторых реализациях входных воздействий , следует более детально рассмотреть процедуру сравнения средних значений или вероятностей. Например, при сравнении вероятностей, задаваясь значениями Dp = p1 p2, pA и pD, необходимо рассчитать значения
р2 = рAD, р1 = p2+Dp, рC = pD+Dр и вычислить коэффициенты корреляции и «выигрыша» соответственно:

;

,

где Nн, и Nз объемы выборки, необходимые для получения заданной точности оценки Dр при использовании независимых и зависимых реализаций.

Таким образом, использование зависимых испытаний дает возможность значительно сократить затраты машинного времени на моделирование. Рассмотренная методика сравнения характеристик вариантов при синтезе системы с учетом их корреляции является формальной. Однако основа для получения с помощью этой методики практических преимуществ неформальная операция выбора такой схемы имитации, при которой искусственно создавалась бы требуемая корреляция.

Оценка результатов моделирования системы. Рассмотрим возможность оценки при обработке результатов моделирования абсолютных значений характеристик процесса функционирования системы S. Пусть исследование одного из вариантов системы, например S2, выполнено аналитическим методом и определено среднее значение m2 критерия q2. Тогда оценка = m2 d среднего значения m1 имеет дисперсию

D[] = D[] = (D[]+D[])/gm = (1+a)D[]/gm,

где gm коэффициент выигрыша, получаемого при оценке разности средних значений d = m2 m1 за счет зависимости испытаний; a = D[]/D[]. Оценка точнее, если (1+a)/gm<1.

Однако затраты машинного времени для получения оценки , которые обозначим как t12, превышают при заданном N затраты машинного времени t1, необходимого для автономной оценки . Поэтому при заданной точности оценки среднего оценка дает выигрыш по затратам машинного времени на имитацию только в том случае, если (1+a)t12/(gmt1)<1.

Для нормально распределенных критериев q1 и q2 оценка дисперсии =D2+D. Выигрыш в затратах машинного времени на имитационное моделирование по сравнению с автономной оценкой будет лишь при условии (1+a)t12/(gDt1)<1, где gD коэффициент выигрыша, получаемого при оценке разности дисперсии Dза счет зависимых испытаний.

Рассмотренные методы сравнения вариантов S1 и S2 моделируемой системы можно использовать в алгоритмах оптимизации на этапе проектирования системы S, т.е. при ее синтезе, по результатам имитационного эксперимента с ее машинной моделью Мм.

При синтезе системы S на основе проведения машинных экспериментов с моделью Мм возникает задача анализа чувствительности модели к вариациям ее параметров. Под анализом чувствительности машинной модели Мм понимают проверку устойчивости результатов моделирования, т.е. характеристик процесса функционирования системы S, полученных при проведении имитационного эксперимента, по отношению к возможным отклонениям параметров машинной модели от истинных их значений .

Анализ чувствительности позволяет сравнивать методические погрешности, полученные при построении машинной модели Мм, с неточностями задания исходных данных, что особенно важно при практической реализации для целей синтеза системы S.

Малым отклонениям будут соответствовать изменения характеристик , которые в практических расчетах можно оценить величиной

,

где ; r0 остаточный член второго порядка малости относительно вариации, который используется для проверки точности решения.

Частная производная определяется в точках, соответствующих номинальным значениям параметров . Если , где
оптимальные параметры системы по показателю , то =0 и необходимо проводить оценку с использованием второй производной . Таким образом, частные производные , количественно характеризуют чувствительность машинной модели Мм к изменениям ее параметров.

Большие отклонения характеристик при малых вариациях свидетельствуют о неустойчивости модели Мм по отношению к этим вариациям. Для получения оценок показателя удобно рассматривать зависимые реализации внешних воздействий при различных и проводить соответствующую обработку результатов машинного эксперимента с моделью Мм.

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