назад Оглавление вперед


[Старт] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] [167] [168] [169] [170] [171] [172] [173] [174] [175] [176] [177] [178] [179] [180] [181] [182] [183] [184] [185] [186] [187] [188] [189] [190] [191] [192] [193] [194] [195] [196] [197] [198] [199] [200] [201] [202] [203] [204] [205] [206] [207] [208] [209] [210] [211] [212] [213] [214] [215] [216] [217] [218] [219] [220] [221] [222] [223] [224] [225] [226] [227] [228] [229] [230] [231] [ 232 ] [233] [234] [235] [236] [237] [238] [239] [240] [241] [242] [243] [244] [245] [246] [247] [248] [249] [250] [251] [252] [253] [254] [255] [256] [257] [258] [259] [260] [261] [262] [263] [264] [265] [266] [267] [268] [269] [270] [271] [272] [273] [274] [275] [276] [277] [278] [279] [280] [281] [282] [283] [284] [285] [286] [287] [288] [289] [290] [291] [292] [293]


232

Следующие действия демонстрируют применение процедуры метода отбора с использованием последовательности случайных чисел из табл. 18.1.

Этап 1. Использование числа R = 0,0589 приводит к случайному значению х = 0,0589, соответствующему плотности h(x).

Этап 2. Выбираем из табл. 18.1 следующее число R = 0,6733.

ЭтапЗ. Так как 7(0,0589)/g(0,0589) = 0,3326/1,5= 0,2217 меньше R = 0,6733, мы отбрасываем значение* = 0,0589.

Для получения второго значения повторяем действия.

Этап 1. Использование числа R = 0,4799 (из первого столбца табл. 18.1) приводит к случайному значению х = 0,4799, соответствующему плотности h(x).

Этап 2. Выбираем из табл. 18.1 следующее число R = 0,9486.

Этап 3. Так как /0,4799)/g(0,4799) = 0,9984 больше R = 0,9486, мы принимаем значение х = 0,4799 как соответствующее бета-распределению.

Из этого примера видно, что итерации метода отбора должны повторяться с новыми значениями случайных чисел, равномерно распределенных на интервале [0, 1], до тех пор, пока не будут выполнены условия этапа 3.

Эффективность метода отбора повышается в результате выбора такой мажорирующей функции g(x), которая "облегала" бы функцию f(x) как можно более плотно, порождая в то же время приемлемую с аналитической точки зрения аппроксимирующую функцию h(x). Например, метод будет более эффективным, если прямоугольную мажорирующую функцию g(x) на рис. 18.6 заменить ступенчатой (см. упражнение 18.3.4.2). С увеличением числа ступеней функция g(x) становится все более прилегающей к функции f(x), и, следовательно, возрастает вероятность принятия полученного случайного значения как искомого. Однако получение тесно прилегающей мажорирующей функции влечет за собой дополнительные вычисления, которые могут стать чрезмерными, что, в свою очередь, может свести на нет преимущества высокой вероятности принятия полученных значений.

УПРАЖНЕНИЯ 18.3.4

1. В примере 18.3.6 продолжите выполнение процедуры метода отбора до получения следующего приемлемого случайного значения.

2. Рассмотрим плотность вероятности бета-распределения из примера 18.3.6. Определите двухступенчатую "пирамидальную" мажорирующую функцию g(x) с двумя равными скачками величины 0,75. Получите случайное значение, соответствующее бета-распределению, с использованием новой мажорирующей функции и тех же случайных чисел из интервала [0, 1] (табл. 18.1), которые использовались в примере 18.3.6. Вывод, который можно здесь сделать, сводится к тому, что использование более точной мажорирующей функции повышает вероятность принятия полученного значения как искомого. Заметьте, однако, что объем вычислений, связанных с использованием новой функции, увеличился.



3. Определите функции g(x) и h(x) для применения метода отбора к плотности распределения следующего вида:

. sin(x) + cos(x) „ ж

w 2 2

Используйте случайные числа из первого столбца табл. 18.1 для получения двух значений, соответствующих плотности распределения f(x). (Совет. Для удобства используйте прямоугольную функцию g(x) над областью определения функции f(x).)

4. Время между приходом клиентов в парикмахерскую описывается следующим распределением:

/(/)=у, 12 < г < 20.

Время стрижки является случайной величиной, плотность вероятности которой

/,(/) = £, 18<г<22.

Константы kx и кг выбираются из условия, что функции ft(x) и f2(x) являются плотностями вероятностей. Используйте метод отбора (и случайные числа из табл. 18.1), чтобы определить время ухода первого клиента из парикмахерской и время прихода второго клиента. Предположите, что первый клиент приходит в момент времени Т = 0.

18.4. ГЕНЕРИРОВАНИЕ СЛУЧАЙНЫХ ЧИСЕЛ

Равномерно распределенные случайные числа из интервала [0, 1] играют ключевую роль в получении выборок из любого вероятностного распределения. Истинные случайные числа из интервала [0, 1] можно генерировать лишь с помощью электронных приборов. Так как имитационные модели реализуются на компьютере, использование электронных приборов для генерации случайных чисел слишком бы замедлило процедуру имитационного моделирования. Кроме того, электронные приборы активизируются случайным образом. Следовательно, невозможно по желанию воспроизвести одну и ту же последовательность случайных чисел. Этот факт чрезвычайно важен, так как для отладки, проверки и утверждения имитационной модели часто требуется дублирование одной и той же последовательности случайных чисел.

В имитационном моделировании единственным подходящим методом генерации случайных чисел из интервала [0, 1] является метод, основанный на арифметических операциях. Такие числа не являются истинно случайными, так как они могут быть определены заранее, поэтому их называют псевдослучайными.

Наиболее часто используется мультипликативный метод сравнений, который генерирует случайные числа из интервала [0, 1] с использованием арифметических операций. В соответствии с этим методом псевдослучайное число Rn при заданных значениях параметров и0, b,c,nm можно вычислить по следующей формуле:

ип = {bun x +c)mod(/w), « = 1,2,..., „и.

*.=-. « = 1.2, .»



Начальное значение параметра и0 обычно называют начальным числом генератора случайных чисел.

Варианты мультипликативного метода сравнений, которые генерируют случайные числа с улучшенными статистическими характеристиками, описаны в книге [2].

Пример 18.4.1

Используя мультипликативный метод сравнений, сгенерируем три случайных числа при следующих начальных данных: 6 = 9, с = 5, и0=11ияз=12.

ы, =(9xll + 5)modl2 = 8, R, = = 0,6667,

и2 =(9x8 + 5)modl2 = 5, Л2 =- = 0,4167,

M3=(9x5 + 5)modl2 = 4, R,= - = 0,3333.

Шаблон Excel chl8RandomNumberGenerator.xls реализует описанный мультипликативный метод сравнений.6 На рис. 18.7 показаны случайные числа, сгенерированные в соответствии с начальными данными этого примера. Отметим, что сгенерированная последовательность случайных чисел состоит только из четырех различных чисел, затем эти числа повторяются. Конкретный выбор параметров и0, Ь, с, и т является решающим фактором, определяющим статистические качества генератора случайных чисел, а также длины его цикла (по окончании цикла генерируемая последовательность начинает повторять себя, как на рис. 18.7). Использование параметров, выбранных "наобум", не дает хорошего генератора случайных чисел. Надежные генераторы, наряду с достаточно большой длиной цикла генерируемых случайных чисел, должны пройти соответствующие статистические проверки, чтобы гарантировать равномерное распределение на интервале [0, 1] полученной последовательности. На это условие нужно обращать внимание при использовании непроверенного программного обеспечения в качестве генератора случайных чисел.

Multiplicative Congruential Method

Input data

u0 =

How many numbers?

Output results

Preii to <j»nei.

Generated random numbers

0.66667

0 41667

0.16667

0 91667

0.66667

16 17

0 41667

016667

0.91667

0.66667

20 ,

0.41667

Рис. 18.7. Сгенерированные в Excel случайные числа

Excel обладает несколькими встроенными средствами генерирования случайных чисел. Так, функция СЛЧИС генерирует равномерно распределенные на интервале [0, 1] числа, средство Генерация случайных чисел из надстройки Пакет анализа позволяет генерировать случайные числа, имеющие различные вероятностные распределения, в том числе равномерное, нормальное, Пуассона. - Прим. ред.

[Старт] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] [167] [168] [169] [170] [171] [172] [173] [174] [175] [176] [177] [178] [179] [180] [181] [182] [183] [184] [185] [186] [187] [188] [189] [190] [191] [192] [193] [194] [195] [196] [197] [198] [199] [200] [201] [202] [203] [204] [205] [206] [207] [208] [209] [210] [211] [212] [213] [214] [215] [216] [217] [218] [219] [220] [221] [222] [223] [224] [225] [226] [227] [228] [229] [230] [231] [ 232 ] [233] [234] [235] [236] [237] [238] [239] [240] [241] [242] [243] [244] [245] [246] [247] [248] [249] [250] [251] [252] [253] [254] [255] [256] [257] [258] [259] [260] [261] [262] [263] [264] [265] [266] [267] [268] [269] [270] [271] [272] [273] [274] [275] [276] [277] [278] [279] [280] [281] [282] [283] [284] [285] [286] [287] [288] [289] [290] [291] [292] [293]