Технические средства и программное обеспечение систем автоматизации научных исследований в геофизике / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. – Апатиты : [б. и.], 1982.

Для вычисления корреляционной функции г ( I ) выражение ( 1 ) мож но использовать в следующем виде: N -H I I ( х ;+е~ < х> Х У г < У » ■/ N-IH ( 2 ) N - ltl у Х (х;+Г <Х»2 E(yi-<y>f где Е - величина переменного сдвига одного исследуемого ряда от носительно другого ( о < /л/ 4 N —l ) . Однако расчеты по формуле ( 2 ) занимают слишком много ма шинного времени. Поэтому для сокращения последнего корреляцион ные функции вычисляются с помошью суммарных статистик: Алгоритм, определяемый ( 3 ) , реализован нами на ЭВМ типа ДОС АСВТ М -4 0 3 0 . Машинной обработке подвергаются одновремен но до десяти рядов исходных данных по СДВ - и геофизическим изме рениям. Каждая обрабатываемая выборка представляет собой матри цу В ( М , М ) , где N - длина строки матрицы ( N m a x = 5 0 0 ), М - дойна столбца матрицы ( М т а х = 1 0 ). Программы вычислений по алгоритму ( 3 ) составлены на языке ФО РТ РАН -1У . Анализ и конт роль правильного ввода входного числового массива осуществляется программно, отдельно от машинного счета. После сформирования мат рицы В ( М , М ) в основной программе дальнейшие вычисления выпол няются многократным обращением к сл&дуюшим подпрограммам: а) формирование корреляционных пар; б ) сдвиги строк; в) расчет ко эффициентов парной линейной корреляции; г ) нахождение минимально го и максимального (по абсолютной величине) значений коэффициен тов корреляции - r m; n ( О , г т а х ( 0 , причем £max = N/2 ; д) блок печати. Выходная информация представляет собой автокорреляционные и корреляционные функции для всех заданных временных рядов. В про- цессе обработки вспомогательными операциями можно изменять фор мат и содержание выходной информации, выводя на печать в том чис ле ряд полезных промежуточных результатов (суммарные статистики, выборочные средние и дисперсии и т .д .). Все указанные программы находятся в фазовой библиотеке ЭВМ, откуда их можно вызывать при обработке. Данные наблюдений считы- N-W N-ltl N -ltl ( 3 ) 42