Генератор случайных функций

38639 О ОПИСАН ИЕ ИЗОБРЕТЕНИЯ К АВТОРСКОМУ СВИДЕТЕЛЬСТВУСоюз Советских Социалистических Республикф Зависимое от авт. свидетельства Ко -М, Кл. Ст 061 1/ Ло 1689552118-24) 1 цо Заявлено ЗО.Ч 111,1971 с присоединением зая Комитет по делам изобретений и аткрыти при Совете Министров СССРиоритет Опубликовано 14 Ч 1.1973. Бюллетень Мо 26Дата опубликования описания 8.Х.1973 ДК 681 142(088.8) о СЕОС ".",; "1.:":11 тс.с ." Гсо 1сГ 1 с 1,", . сС., с, СЭ,4 и с; Авторыизобретени алышев и Л. Т. Тота Городецкий,анрогский радиотехнический институт аявитель ЕНЕРАТОР СЛУЧАИНЫХ ФУНКЦИ Изобретение относится к области технической кибернетики и может быть использовано в различных системах связи, управления и контроля, а также в качестве специализированной приставки к ЭВМ, для генерирования случайных процессов с произвольными законами распределения и заданной функцией корреляции. Известны устройства, которые позволяют формировать случайные процессы с,нормальным и пуассоновскими законами распределения и заданными вероятностными характеристиками, содержащие генератор равномерно распределенных случайных чисел, цифровой функциональный преобразователь, цифровой фильтр, сумматор и цифровой, нелинейный функциональный преобразователь.Однако эти устройства не позволяют формировать случайные процессы с произвольными законами распределения и заданными вероятностными характеристиками.Цель предлагаемого изобретения заключается в создании цифрового генератора случайных функций, способного формировать случайные процессы с произвольными законами распределения и заданными вероятностными характеристиками (математическим ожиданием М, дисперсией о и функцией корреляции Я (т. Для реализации поставленной задачи впредлагаемое устройство дополнительно введены: устройство ввода, блок управления, вычислительное устройство, ассоциативное запо минающее устройство законов распределения,блоки памяти, причем выход цифрового фильтра через сумматор подключен к первой выходной клемме, а через цифровой нелинейный функциональный преобразователь - ко 10 второй выходной клемме, первый вход цифрового фильтра соединен с выходом генератора первичного нормального сигнала, второй вход - с первым выходом блока управления, второй выход которого соединен с вычисли тельным устройством, а третий выход черезассоциативное запоминающее устройство соединен со входами первого и второго блоков памяти, выходы первого блока памяти подключены к цифровому нелинейному функцио нальному преобразователю, выход второгоблока памяти через вычислительное устройство соединен с первым входом блока управления, второй вход которого соединен с устройством ввода.25 В основу метода генерирования случайныхпроцессов с заданными законами распределения и корреляционными функциями положен принцип нелинейного преобразования случайного процесса с предварительно сформпро ванной функцией корреляции.(3) а, /1 - г (6) 65 Блок-схема предлагаемого устройства .изображена на фиг. 1.Она содержит генератор 1 первичного нормального сигнала, генератор 2 равномерно распределенных чисел, цифровой функциональный преобразователь 3, цифровой фильтр 4, устройство ввода б, блок управления б, вычислительное устройство .7, ассоциативное запоминающее устройство 8 законов распределения, блоки памяти 9 и 10, сумматор 11, универсальный цифровой, нелинейный функциональный преобразователь 12.Математические предпосылки реализации предлагаемого универсального цифрового генератора случайных функций с заданными вероятностными характеристиками заключаются в следующем.Заданный закон распределения вероятностей , (У) выходного случайного процесса У(1) определяется кривой нелинейного преобразования У(х) нормального случайного процесса Х(1), осуществляемого преобразователем 12. Для получения необходимой функции корреляции выходного сигнала с законом распределения, отличным от нормального, требуется заданную функцию корреляции Р(т) пересчитать с выхода на вход преобразователя 12 с учетом найденного вида кривой преобразования, а затем полученную в результате пересчета корреляционную функцию Р(т) сформировать с помощью линейного фильтра из широкого нормального сигнала генератора 1. Полученный таким образом на выходе фильтра нормальный случайный процесс, пройдя преобразователь 12, будет иметь заданные закон распределения и корреляционную функцию.Алгоритм работы предлагаемого генератора включает в себя следующее:1. Определение вида и параметров кривой нелинейного преобразования случайного процесса Х с плотностью вероятности ,. (Х) в случайный процесс У с заданной плотностью распределения , (У) по следующему, справедливому для всех непрерывных функций соотношению:х,Д)1 х =(У) 1 у. (1) 2. Пересчет заданной функции корреляции Р(т) с выхода на вход блока 12 с учетом найденного вида кривой преобразования У(х),3, Расчет передаточной функции цифрового фильтра для формирования функции корреляции исходного нормального случайного процесса Х(1), которая равна:5 (ор)1 " (1")где 5(ж) и 5 (ю) - спектральные плотности входного и выходного случайного процесса,В соответствии с приведенным алгоритмом работы может быть реализован генераторслучайных функций, решающий поставленнуюзадачу при наличии формул перехода междунормальным законом распределения и любымдругим непрерывным законом,5 В целях пояснения решения поставленнойзадачи рассмотрим алгоритм работы предлагаемого генератора при формировании райсовских случайных процессов с заданнымифункциями корреляции.10 Уравнение кривой нелинейного преобразования определится из выражения:1 - Ф(х) = Я(И),15 где Ф(х) .= .е а 7,1) 2 рг- интеграл вероятностей;СО еегеггг)9(У, /)= ре е,(Г рдр- интеграл вероятностейРэлея-Райса.Используя известное разложение функциц Я(У, 1 р) в ряд, равенство (3) запишем в виде; 1 - Ф(х) - е Г 1+ Уф+ ( ) Уф4 64 30(4)После соответствующих преобразований равенства (4) получим выражение: З 5 У (х) = Сг 3/ - 1 п 1 - Э (х) ХкХВ Кг ( - 1 п 1 - Ф (х д, (5) которое является аналитической формулой нелинейного функционального преобразования нормальных случайных чисел Х в последовательность чисел (1= (1(х), распределенных по обобщенному закону Рэлея-Райса.4 Здесь С Кг и В - коэффициенты, соответствующие определенной плотности распреде.ления ,(У), которая определяется, в свою очередь, значением райсовского параметра КНа фиг. 2 показаны графики кривых преобразователя 12, полученные для преобразования нормального закона распределения Х(1) в райсовский для различных значений величин параметра 1 р.Пересчет заданной функции корреляцириР(т) с выхода преобразователя 12 осуществляется из выражения:Р=2 д 1 - г60 Х Кгде г,(т) - нормированная функция корреляции входного нормально го процесса Х(1);(7) с = ЬХУ(х)дх,СО2ХУ(х)дх 45 50 55 60 б 5 й,б,с,А - постоянные коэффициенты, определенные для фиксирован.ных значений райсовского параметра Р: ОЭ- х 2а =с/(х)ехр-- Йх,На фиг. 3 приведены кривые функциональной зависимости Рц(т) =рг(т) для фиксированных КАналогичным образом данная задача может быть решена практически для всех функций распределения вероятности, и.,еющих технические приложения.Предлагаемое устройс;во работает следующим образом,Генератор 1 первично:, н рмального случайного сигнала состоит:, енератора 2 равномерно распределенных ел и цифрового функционального преобр зтеля,З, который преобразует исходную,.в . черно распределенную последовательн,сть чисел в нормальный случайный сигнал У (аТ) с некоррелированными дискретными з.ачениями. Сигнал 17(пТ) поступает на цифр, вой фильтр, который формирует функцию корреляции г(пТ) этого сигнала, Цифровой филь;р 4 управляется командами, поступающими с блока управления б. Через устройство ввода б в буферную память блока управления записывается код, соответствующий определенному виду заданной функции корреляции Л(пТ) (где Т - интервал дискретности) и плотности вероятности (У) заданного закона распределения. Работа цифрового генератора происходит одновременно по двум каналам: каналу формирования заданной корреляционной функции Р,(пТ) и каналу выбора кривой не. линейного преобразования У(х) для формирования заданной плотности (У) закона распределения, В соответствии с этим код ( У) с блока б поступает на ассоциативное запоминающее усгройство 8 законов распределения, в памяти которого хранятся кодовые комбинации, соответствующие различным законам распределения. Запоминающее устройство осуществляет перебор всевозможных кодовых комбинаций и при совпадении одной из нихе с заданной выдает кодовую посылкублоки памяти 9 и 10. В блоке 9 хранятся значения коэффициентов С,К и В(п=0,1,2 Й), которые сигналом У, одновременно считываются на вход универсального цифрового нелинейного функционального преобразовате 5 10 15 20 25 30 35 40 6ля 12, реализующего формулу (5). В соответствии со значением выработанного с устройства д сигнала Гф и считанных из блока 9 коэффициентов СК, и Впо формуле (5), с помощью сдвиговых регистров и сумматоров, осуществляется перестройка передаточной функции так, чтобы реализовать необходимую для заданного вида плотности вероятности 1(У) кривую нелинейного преобразования У(х). В блоке 10 хранятся коды коэффициентов а Ь с и с 1 (7), которые по сигналу , считываются на вход вычислительного устройства 7, на другой вход которого с блока б засылается код, соответствующий заданному виду корреляционной функции Ку(пТ), Устройство 7 представляет собой цифровой арифметический вычислитель, реализующий выражение (6) и служащий для расчета значений нормированной функции корреляции г(пТ). С выхода устройства 7 рассчитанные ординаты г(пТ) засылаются в блок управления 6, с которого ординаты г,(пТ) поступают на вход фильтра 4. Фильтр 4, состоящий пз сдвиговых регистров и сумматора, формирует заданную функцию корреляции г(лТ) входного нормального сигнала 7(гТ). Нормальный сигнал Х(лТ) на выходе фильтра имеет сформированную функцию корреляции г(пТ) дискретных значений. Этот сигнал пропускается через преобразователь 12, на выходе которого образуется дискретный сигнал У(пТ), имеющий заданную плотность распределения и функцию корреляции К(пТ). Если необходимо иметь кривую нелинейного функционального преобразования в каком-либо квадранте, то ко входному нормальному процессу Х(пТ) с помощью сумматора можно добавить некоторое значение математиче. ского ожидания, например т.=о,2 о Зо, 4 о. Для получения дискретных нормальных случайных процессов с заданной корреляци. онной функцией организован Вых, 2,Предмет изобретенийГенератор случайных функций, содержа щий генератор первичного нормального сигнала, цифровой фильтр, сумматор и цифровой нелинейный функциональный преобразователь, отличающийся тем, что, с целью расширения функциональных возможностей устройства, он дополнительно содержит устройство ввода, блок управления, вычислительное устройство, блоки памяти, ассоциативное запоминающее устройство, при этом выход цифрового фильтра через сумматор подключен к первой выходной клемме, а через цифровой нелинейный функциональный преобразователь - ко второй выходной клемме, первый вход цифрового фильтра соединен с выходом генератора первичного нормального сигнала, второй вход - с первым выходом блока управления, второй выход которого соединен с вычислительным устройством, а третий выход через ассоциативное запоминающее устройст386390 го блока памяти через вычислительное устройство соединен с первым входом блока управления, второй вход которого соединен с устройством ввода. 1 П 7 до авитель Э чина Б. НанкинаБогдановапорожец и Н. А ак КорректорыЗаказ 2702/1 Тираж 647 Типография, пр апунова, 2 во соединен со входами первого и второго блоков памяти, выходы первого блока памяти подключены к цифровому нелинейному функциональному преобразователю, выход второ 0 ,0 О Изд.166ПодписноеНИИПИ Комитетазобретений и открьСовете Министровва, Ж, Раушская по деламтий приСССРнаб д, 4/5