Советов Б.Я., Яковлев С.А. Моделирование систем (3-е изд., 2001) (1186218), страница 28

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 28)

Запо­минается только накопленная сумма исходов и общее число ре114ализаций. Это немаловажное обстоятельство вообще является хара­ктерным при реализации имитационных моделей методом стати­стического моделирования на ЭВМ.4.2. ПСЕВДОСЛУЧАЙНЫЕ ПОСЛЕДОВАТЕЛЬНОСТИ И ПРОЦЕДУРЫИХ МАШИННОЙ ГЕНЕРАЦИИПри статистическом моделировании систем одним из основныхвопросов является учет стохастических воздействий. Количествослучайных чисел, используемых для получения статистически устой­чивой оценки характеристики процесса функционирования системыS при реализации моделирующего алгоритма на ЭВМ, колеблетсяв достаточно широких пределах в зависимости от класса объектамоделирования, вида оцениваемых характеристик, необходимой то­чности и достоверности результатов моделирования.

Для методастатистического моделирования на ЭВМ характерно, что большоечисло операций, а соответственно и большая доля машинноговремени расходуются на действия со случайными числами. Крометого, результаты статистического моделирования существенно за­висят от качества исходных (базовых) последовательностей случай­ных чисел. Поэтому наличие простых и экономичных способовформирования, последовательностей случайных чисел требуемогокачества во многом определяет возможность практического исполь­зования машинного моделирования систем [31, 37, 46].Рассмотрим возможности и особенности получения последова­тельностей случайных чисел при статистическом моделированиисистем на ЭВМ.

На практике используются три основных способагенерации случайных чисел: аппаратный (физический), табличный(файловый) и алгоритмический (программный).Аппаратный способ. При этом способе генерации случайные чис­ла вырабатываются специальной электронной приставкой — гене­ратором (датчиком) случайных чисел,— служащей в качестве одно­го из внешних устройств ЭВМ.

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

4.8, а. Здесь ИШ — источник шума; КС — клю­чевая схема; ФИ — формирователь импульсов; ПС — пересчетнаясхема. При усилении шумов на выходе ИШ получается напряжение115Ыш(0, которое является слу­чайным процессом, пока­занным на временной диа­грамме рис. 4.8, б. Причемотрезок шумовой реализа­ции и, (0, сформированныйна интервале времени (О, Т)г tс помощью КС, содержитслучайное число выбросов.Сравнениенапряжения1 I'll I i 4UI III IЫх(0 с пороговым Un позво­ляет сформировать на вы­ходе ФИ серию импульсов«ф(0- Тогда на выходе ПС4 *i+;tможет быть получена по­Рис.

4.8. Аппаратный способ получения слу­следовательность., случай­чайных чиселных чисел x,(t). Например,если провести масштабирование и принять длину интервала (0, 7)за единицу, то значения интервалов времени At,=ti+l — t) междусоседними импульсами u$(t) будут случайными числами дс(е(0, 1).Возможны и другие схемные решения аппаратных генераторов слу­чайных чисел [29, 37].

Однако аппаратный способ получения случай­ных чисел не позволяет гарантировать качество последовательностинепосредственно во время моделирования системы S на ЭВМ,а также повторно получать при моделировании одинаковые после­довательности чисел.Табличный способ. Если случайные числа, оформленные в видетаблицы, помещать во внешнюю или оперативную память ЭВМ,предварительно сформировав из них соответствующий файл (мас­сив чисел), то такой способ будет называться табличным.

Однакоэтот способ получения случайных чисел при моделировании системна ЭВМ обычно рационально использовать при сравнительно небо­льшом объеме таблицы и соответственно файла чисел, когда дляхранения можно применять оперативную память. Хранение файлаво внешней памяти при частном обращении в процессе статистичес­кого моделирования не рационально, так как вызывает увеличениезатрат машинного времени при моделировании системы S из-занеобходимости обращения к внешнему накопителю. Возможныпромежуточные способы организации файла, когда он переписыва­ется в оперативную память периодически по частям. Это уменьшаетвремя на обращение к внешней памяти, но сокращает объем опера­тивной памяти, который можно использовать для моделированияпроцесса функционирования системы S.Алгоритмический способ.

Способ получения последовательностейслучайных чисел основан на формировании случайных чисел в ЭВМс помощью специальных алгоритмов и реализующих их программ.v/sA^V^A/116Каждое случайное число вычисляется с помощью соответствующейпрограммы по мере возникновения потребностей при моделирова­нии системы на ЭВМ.Таблица 4.1НедостаткаДостоинстваТребуется периодическая про­Запас чисел не ограниченРасходуется мало операций веркаНельзя воспроизводить после­вычислительной машиныНе занимается место в памяти довательностиИспользуетсяспециальноемашиныустройствоНеобходимы меры по обеспе­чению стабильностиТабличныйЗапас чисел ограниченТребуется однократная про­Занимает много места в опе­веркаМожно воспроизводить по­ ративной памяти или необходи­мо время на обращение к вне­следовательностишней памятиАлгоритми­Требуется однократная про­Запас чисел последовательно­ческийверкаМожно многократно воспрои­ сти ограничен ее периодомзводить последовательности чи­Существенные затраты ма­селшинного времениЗанимает мало места в памятимашиныНе используются внешниеустройстваСпособАппаратныйДостоинства и недостатки трех перечисленных способов получе­ния случайных чисел для сравнения представлены в табл.

4.1. Изэтой таблицы видно, что алгоритмический способ получения слу­чайных чисел наиболее рационален на практике при моделированиисистем на универсальных ЭВМ.Генерация базовой последовательности. При моделировании си­стем на ЭВМ программная имитация случайных воздействий лю­бой сложности сводится к генерированию некоторых стандартных(базовых) процессов и к их последующему функциональному преоб­разованию. Вообще говоря, в качестве базового может быть принятлюбой удобный в случае моделирования конкретной системы ^про­цесс (например, пуассоновский поток при моделировании Q-схемы).Однако при дискретном моделировании базовым процессом явля­ется последовательность чисел {xi}=x0, xlt ..., xN, представляющихсобой реализации независимых, равномерно распределенных на ин­тервале (0, 1) случайных величин {Q = £0, £u ..., £N или — в стати­стических терминах — повторную выборку из равномерно распре­деленной на (0, 1) генеральной совокупности значений величины« [29, 37, 46].Непрерывная случайная величина £ имеет равномерное рас­пределение в интервале (а, Ь), если ее функции плотности117(рис.

4.9, а) и распределения (рис. 4.9, б) соответственно примутвидО,{х<я,(х-а)/(Ь-а),а^х4:Ь,1,х>6.Определим числовые характеристики случайной величины {,принимающей значения х,— математическое ожидание, дисперсиюи среднее квадратическое отклонение соответственно:ььM\Z\4 xf(x) dx=\xdxl(b-a) = (a+b)l2;D\Z\=$(x-M\£\)2f(x)dx=(b-a)2ll2;<rt=+y/D\e\ = Q>-a)l<fr/3).При моделировании систем на ЭВМ приходится иметь дело сослучайными числами интервала (0, 1), когда границы интервалаа = 0 и Ь=\.

Поэтому рассмотрим частный случай равномерногораспределения, когда функция плотности и функция распределениясоответственно имеют вид{0, х<0,fix) \o, x<0, x>l;х, 0 < х < 1 ,1, х > 1 .Такое распределениеимеет математическое ожидание М[€] = 1/21и дисперсию D[£\ = I12Это распределение требуется получить на ЭВМ. Но получить егона цифровой ЭВМ невозможно, так как машина оперирует с празрядными числами. По­этому на ЭВМ вместо непре­рывной совокупности равно­мерных случайных чисел ин­тервала (0, 1) используютдискретную последователь­ность 2" случайных чисел то­го же интервала.

Закон рас­пределения такой дискрет­ной последовательности наРис. 4.9. Равномерное распределение слузывают квазиравномернымчайной величины118распределением.Случайная величина £, имеющая квазиравномерное распределе­ние в интервале (0, 1), принимает значения JC,=Z/(2"—1) с вероят­ностями pi= 1/2", 1=0, 2я—1.Математическое ожидание и дисперсия квазиравномерной слу­чайной величины соответственно имеют видг 1~'11'2'-1 .Щ&= У - = - - = — - — У i=L J^*->•)« 1 t*= 0 2"-12"~ 0 2"|_2"-1(2"-1)2",toЛ"I V } " А-"2j_ 1 /(2"-1)2"(2" +1 -1)~2"\6(2"-1)2(2"-1)2"1(2"-l)2"222" ,Г 0 \(2"-1) 22"-l+V(2"-l)2"l\_l2"+l2(2"-1)у \ 2 ? - \Таким образом, математическое ожидание квазиравномернойслучайной величины совпадает с математическим ожиданием рав­номерной случайной последовательности интервала (0, 1), а диспер­сия отличается только множителем (2"+1)/(2" —1), который придостаточно больших п близок к единице.На ЭВМ невозможно получить идеальную последовательностьслучайных чисел хотя бы потому, что на ней можно оперироватьтолько с конечным множеством чисел.

Характеристики

Список файлов книги