Анализ периодических компонент интенсивности точечных процессов. Любушин А.А. - 12 стр.

                                              12
распределения значений R(ω ,τ | L) для заданной частоты (или периода) на искусственных
моментах времени и задать уровень значимости как квантиль этого распределения с
соответствующей доверительной вероятностью.
   5. Программная реализация метода.
   Программа PPPeriod (Point Processes Periodicity) написана на языке Compaq Visual
Fortran и представляет собой консольное приложение. Диалог с пользователем организован
по принципу последовательности «вопрос-ответ» на английском языке. Целью программы
является вычисление приращений максимума логарифмической функции правдоподобия для
модели интенсивности точечного процесса с заданным периодом по отношению к
максимуму логарифмической функции правдоподобия для модели чисто пуассоновского
процесса с постоянной интенсивностью (для нулевой гипотезы). Значения периода
сканируются в заданных пределах и те периоды, для которых эта разность имеет пик,
существенно возвышающийся над уровнем фоновых статистических флуктуаций оценки,
выделяют периодические компоненты потока событий.
   Входным файлом программе может служить любая символьная таблица чисел, в первой
колонке которой стоят монотонно неубывающие числа, интерпретируемые как значения
моментов времени. Прочие колонки входного файла, если таковые имеются, игнорируются.
   По выбору пользователя программа может работать в 2-х режимах:
   1) оценка по всей имеющейся выборке;
   2) оценка в скользящем временном окне заданной длины и с заданным взаимным
      смещением (в размерном времени).
   После запуска программы необходимо ответить на следующие вопросы:
   1. Имя входного файла. Файл должен находится в той же директории, что и загрузочный
      модуль программы. Если файла нет или он имеет неправильную структуру (например,
      в первой колонке входной таблицы значения убывают для каких-то строк), то
      программа завершает работу с соответствующим сообщением.
   2. Ввести значения Tmin и Tmax - минимальных и максимальных значений периодов
      сканируемого диапазона периодов в размерном времени, то есть в тех же единицах
      времени, которые присутствуют во входном файле.
   3. Ввести значение N p - числа пробных значений периодов, покрывающих диапазон

      [Tmin , Tmax ]   с   шагом,   равномерным    в   логарифмическом   масштабе,   то   есть
      рассматриваются периоды: