СПОСОБ ЦИФРОВОЙ ГЕНЕРАЦИИ НЕПРЕРЫВНОГО СЛУЧАЙНОГО ПРОЦЕССА ПРОИЗВОЛЬНОЙ ДЛИНЫ

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

Генерация во временной области случайного сигнала x(t), имеющего нормальный закон распределения амплитуд, производится по известному алгоритму («Автоматическое управление вибрационными испытаниями», Библиотека по автоматике, вып. 579, М., Энергия, 1978 г.) путем формирования во временной области на основе заданного спектра амплитуд последовательности (временного блока a ₁) вещественных чисел a _1i (i=0, 1, …, N-1) в форме разложения Райса-Пирсона с шагом дискретизации Δt [сек] и преобразования полученной последовательности в непрерывный во времени случайный сигнал длиной T_Б=N·Δt=1/Δf [сек], обратно пропорциональной шагу квантования спектра амплитуд сигнала по частоте Δf.

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

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

Предлагается способ формирования непрерывного случайного процесса произвольной длины путем цифровой генерации во временной области на основе заданного спектра амплитуд сигнала последовательностей вещественных чисел a _ri (i=1, 2, …; i=0, 1, …, N-1) в форме разложения Райса-Пирсона с шагом дискретизации Δt [сек] (временных блоков a _r, длиной T_Б=NΔt); построения косинусоидального окна Ханна на длине T_HANN=(4·n-3) точек (n<<N, например, n=N/10):

;

формирования случайного процесса из последовательности следующих друг за другом временных блоков a _r (r=1, 2, …), состыкованным по последним a _r,N-1 и первым a _r+1.0 точкам блоков a _r и a _r+1 соответственно по следующему алгоритму: в конце блока a _r (r=1, 2, …) и в начале следующего блока a _r+1 выделяют последовательности точек, длиной n, к концу блока a _r добавляют (n-1) точек, симметричных последним n точкам блока a _r, умноженным на вторую половину окна Ханна, к началу блока a _r+1 добавляют (n-1) точек, симметричных первым n точкам блока a _r+1, умноженным на левую половину окна Ханна, полученные в блоках a _r и a _r+1 последовательности складывают для получения суммарной последовательности, производят замену последовательностей точек блоков так, что последние n точек блока a _r заменяют первыми n точками суммарной последовательности, первые n точек следующего блока a _r+1 заменяют последними n точками суммарной последовательности; создания для формирования на выходе из системы управления из отдельных следующих друг за другом блоков a _r непрерывного случайного сигнала цифрового механизма "ping-pong": в памяти компьютера организуют два буфера длиной N, в один из которых заносят сгенерированный блок данных a _r, который затем контроллером прямого доступа к памяти (DMA) без участия центрального процессора передается в память цифроаналогового преобразователя (ЦАП). За время выдачи первого блока a _r на выход ЦАП центральный процессор генерирует следующий блок данных a _r+1, осуществляет стыковку его начала с концом предыдущего a _r и переносит во второй буфер, который, в свою очередь, будет передан в ЦАП через DMA-канал. Затем буфера автоматически циклически меняются местами, и на выходе ЦАП генерируется непрерывный во времени случайный сигнал произвольной длины.

Отличительным признаком предлагаемого способа является:

построение косинусоидального окна Ханна на длине T_HANN=(4·n-3) точек (n<<N):

;

формирование случайного процесса из последовательности следующих друг за другом временных блоков a _r (r=1, 2, …), состыкованным по последним a _r,N-1 и первым a _r+1.0 точкам блоков a _r и a _r+1 соответственно по следующему алгоритму: в конце блока a _r (r=1, 2, …) и в начале следующего блока a _r+1 выделяют последовательности точек, длиной n, к концу блока a _r добавляют (n-1) точек, симметричных последним n точкам блока a _r, умноженным на вторую половину окна Ханна, к началу блока a _r+1 добавляют (n-1) точек, симметричных первым n точкам блока a _r+1, умноженным на левую половину окна Ханна, полученные в блоках a _r и a _r+1 последовательности складывают для получения суммарной последовательности, производят замену последовательностей точек блоков так, что последние n точек блока a _r заменяют первыми n точками суммарной последовательности, первые n точек следующего блока a _r+1 заменяют последними n точками суммарной последовательности; создание для формирования на выходе из системы управления из отдельных следующих друг за другом блоков a _r непрерывного случайного сигнала цифрового механизма "ping-pong": в памяти компьютера организуют два буфера длиной N, в один из которых заносят сгенерированный блок данных a _r, который затем контроллером прямого доступа к памяти (DMA) без участия центрального процессора передается в память цифроаналогового преобразователя (ЦАП). За время выдачи первого блока a _r на выход ЦАП центральный процессор генерирует следующий блок данных a _r+1, осуществляет стыковку его начала с концом предыдущего a _r и переносит во второй буфер, который, в свою очередь, будет передан в ЦАП через DMA-канал. Затем буферы автоматически циклически меняются местами, и на выходе ЦАП генерируется непрерывный во времени случайный сигнал произвольной длины.

Благодаря наличию указанного отличительного признака в совокупности с известными приобретается возможность формирования непрерывных случайных процессов x(t) заданной произвольной длины.

В результате поиска по источникам патентной и научно-технической информации решений, содержащих аналогичные признаки, не обнаружено.

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

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

Способ реализуется следующим образом. В конце блока a _r (для примера взято r=1) и в начале следующего a _r+1 выделяют последовательности длиной n точек (как правило, значение n является незначительным по отношению ко всей длине блока N точек, например, в данном примере n=N/10=10), на которых строится алгоритм стыковки фаз (Фиг.1) На длине T_HANN=(4·n-3) точек вычисляют косинусоидальное окно Ханна (Фиг.2):

.

К концу блока a _r,N-1 добавляют (n-1) точек, симметричных последним n точкам блока: a _1,N= a _1,N-2, …, a _1,N+n-2=a _1,N-n. Умножают полученную последовательность на вторую половину окна Ханна (Фиг.3). К концу последующего блока a _r+1 (точке a _2,0) добавляют (n-1) точек, симметричных первым n точкам блока: a _2,-1=a _2,1, …, a _2,1-n=a _2,n-1. Умножают полученную последовательность на первую половину окна Ханна (Фиг.4). Преобразованные последовательности (Фиг.3 и 4) складывают в суммарную последовательность (Фиг.5). Точками суммарной последовательности заменяют точки исходных блоков точек блоков a _r и a _r+1 так, что последние n точек блока a _r заменяют первыми n точками суммарной последовательности, а первые n точек следующего блока a _r+1 заменяют последними n точками суммарной последовательности (Фиг.6).