Анализ высокоэнергетичных компонент космического излучения методами непараметрической статистики тема автореферата и диссертации по физике, 01.04.16 ВАК РФ

ВВЕДЕНИЕ

ГЛАВА I. ИДЕНТИФИКАЦИЯ АДРОНОВ КОСМИЧЕСКОГО ИЗЛУЧЕНИЯ

§ I. Байесовские решающие правила

§ 2. Краткий обзор методов непараметрического оценивания функции плотности вероятности

§ 3. Адаптивное КБС оценивание

§ 4. Автоматическая классификация адронов космического излучения

ГЛАВА П. АНАЛИЗ МНОГОМЕРШХ ДАНШХ ЭКСПЕРИМЕНТОВ

С РЕНТГЕН-ЭМУЛЬСИОННЫМИ КАМЕРАМИ

§ I. Меры совпадения и различия многомерных распределений

§ 2. Метода вычисления байесовского риска

§ 3. Определение доли семейств, инициированных различными первичными ядрами

§ 4. О возможности многомерного сравнения экспериментальных данных и теоретических моделей

ГЛАВА Ш. ОПРЕДЕЛЕНИЕ ЭНЕРГИИ МЮОНОВ КОСМИЧЕСКОГО ИЗЛУЧЕНИЯ

СПЕКТРОМЕТРИЧЕСКИМИ И КАЛОРИМЕТРИЧЕСКИМИ МЕТОДАМИ

§ I. Определение импульса мюона магнитными спектрометрами с железным сердечником

§ 2. Максимально-измеримый импульс магнитных спектрометров

§ 3. Коррекция котировочных констант магнитного спектрометра

§ 4. Сравнение различных методов оценивания энергии многослойными установками

§ 5. Подавление фона от ядерных взаимодействий

§ 6. Регрессионный метод определения энергии

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

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

Однако, как часто бывает, развитие новых вычислительных методов ставит много методологических вопросов. Среди них можно отметить вопрос об обоснованности статистических выводов о модели по результатам вычислительных экспериментов. В настоящее время ряд статистических процедур, используемых для анализа результатов машинного моделирования, весьма ограничен, поэтому проведение адекватного анализа часто оказывается настолько трудным, что иногда фактически игнорируется, и полученные оценки переносятся на физические явления без особых на то оснований /2/.

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

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

Этих опасностей старается избежать прикладная статистика -дисциплина, в которой мера адекватности выбранной математической модели и изучаемой действительности является решающим фактором, определяющим эффективность и действенность используемых затем методов статистической обработки /4/. В прикладной статистике стараются осуществить более реалистичную постановку задач. Это прежде всего касается отказа от априорного предположения об известной с точностью до неизвестного параметра функции распределения, то есть от параметрических методов. Непараметрическая статистика с самого начала предполагает, что функциональный вид распределений неизвестен, при этом стараются избегать ссылок на центральную предельную теорему, для обоснования нормальности ошибок, так как многочисленные исследования показали, что на практике реализуются несимметричные, многомодельные распределения, которые, в лучшем случае, лишь приближенно согласуются с нормальным /5/. Практическая работа с существенно ограниченными выборками поставила под сомнение актуальность для прикладных вопросов значительной части асимптотической теории математической статистики, как относящейся к пределам, которые фактически не достигаются. В /6/ даже отвергается любая статистическая техника, которую нельзя проверить с помощью машинного моделирования.

Однако, к сожалению, в физике космических лучей до сих пор широко применяются классические схемы обработки даншх, которые основываются на независимости и нормальности случайных ошибок, а поиск закономерностей ведется в предположении априорно постулируемой функциональной зависимости. Исторически эти метода впервые были применены к астрономическим и геодезическим наблюдениям, где в экспериментальных методах была достигнута большая точность и распределения ошибок носили гауссовский характер. Живучесть этих представлений связана и с тем, что абстрактная теория статистики использует нормальную модель как удобную идеализацию, позволяющук создать последовательную теорию. Однако, как показал еще К.Пирсо* /7/, миф о всеобщей нормальности является абсурдным, и вообще, очень редко .эмпирические распределения в экспериментальной физике можно заменить аналитическими. Как показано в /8/,использование сложных аналитических выражений может привести к тому, что функция будет хорошо описывать конкретные эмпирические данные, но -хуже искомую зависимость.

Характерной особенностью статистических решений в вычислительном эксперименте является то обстоятельство, что информация, необходимая для математической формулировки задачи обработки, заключена в наборах реализаций вероятностной модели.Эти наборы обычно называют обучающими выборками (ОВ). В ОВ отображаются всевозможные вариации наблюдаемых физических величин, порождаемые индивидуальными различиями и изменениями множества их свойств, допустимых в пределах данного состояния. В дальнейшем для генерации ОВ будет использован ряд машинных моделей. Это программа прохождения мюонов через Арагацский искровый калориметр /9/, алгоритм моделирования РПИ-детектора /10/, имитационная модель магнитного спектрометра /II/, теоретические банки данных сотрудничества "Памир" /12/ и программа моделирования отклика периодической системы детекторов ионизации /13/. Все модели осуществляют учет экспериментальной ситуации и на выходе генерируют величины, аналогичные наблюдаемым в эксперименте. Так как при получении ОВ контролируются входные параметры машинных моделей, то всегда известно, какой тип начальных данных (какая модель) реализовался в данном наборе. Использование ОВ с известной принадлежностью, как базы дальнейшего статистического анализа экспериментальных данных, носит название обучения с "учителем". Обучающие выборки отражают стохастический характер взаимодействий космических лучей, наличие многих каналов сильных и электромагнитных взаимодействий, а также случайный характер явлений, лежащих в основе функционирования регистрирующих приборов. Поэтому, было бы чрезмерным ожидать от распределений выходных характеристик машинных моделей "хорошего" поведения, укладавающегося в предположения классических процедур. К выбору априорных предположений, лежащих у начала анализа, следует подходить с известной осторожностью. Самым слабым предположением является предположение, что характерные черты изучаемого явления как то отражаются через структуру обучающей выборки. Согласно гипотезы компактности /14/ следует ожидать, что векторы признаков данного класса (модели) в некотором смысле будут ближе друг к другу, чем большинство векторов других классов. Локальное оценивание функции плотности вероятности в признаковом пространстве, как будет показано, приведет к содержательным и эффективным статистическим процедурам. Если же оценки плотности искать в предположении узкой структуры (например, предположить нормальность признака и оценить перше моменты его распределения), то можно прийти к совершенно неверным выводам.

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

Рассмотренное локальное оценивание является инструментом интерпретации данных. Допустим, существуют альтернативные модели, реализованные в обучающих выборках. Измеренное экспериментальное значение х" , как и реализации моделей в ОБ, являются точками N -мерного выборочного пространства. Первым шагом на пути решения, к какой из моделей эта точка наиболее близка, является интерпретация экспериментального значения. Если бы в выборочном пространстве был задан некий класс функций ^(х] , зависящий от параметра 6 , то лучшей интерпретацией была бы функция правдоподобия Щх/6), ОеЯ Заметим, что следует отличать метод максимального правдоподобия как метод интерпретации данных и как метод оценивания. В данном случае отнюдь не следует, что максимально правдоподобная оценка (оценка и , максимилизирующая ) является наилучшей оценкой параметра 6 /15/.

Так как в нашем случае нельзя ожидать, что ОБ можно описать функцией с не очень большим числом параметров, то единственной возможностью интерпретации экспериментального значения X являются локальные оценки плотности в точке X , с использованием альтернативных классов 0В-р(х/(0|) »в которых суммируется информация о взаимном расположении вектора X и классов обучающей выборки, , I = 4 , ; Ь - число классов ОВ.

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

Байесовском подходе /17/. Введение априорных плотностей, характерных для байесовского подхода, не встречает при обработке данных физики космических лучей серьезных затруднений. Например, при классификации адронов космического излучения в качестве априорной вероятности можно взять ЛГ/р отношение, измеренное в предыдущем наиболее достоверном эксперименте. Внесение априорной информации непосредственно в процедуру статистического решения позволяет учитывать, кроме того, и дополнительную информацию: энергию адро-на, измеренную ионизационным калориметром, и плотность электромагнитного сопровождения.

Многие задачи физики космических лучей имеют характер классификационных - требуется отнести зарегистрированное событие к одному из возможных состояний. При классификации большие объемы зашумленной, косвенной информации преобразуются в ценную информацию о отнесении события к конкретному типу. Задача определения первичного состава космического излучения при помощи анализа флуктуационных распределений мюонов и электронов в ШАЛ /18/ также относится к этому типу. И, наверно, было бы целесообразнее использовать обучающие выборки, полученные путем полного моделирования экспериментальной ситуации, непосредственно для классификации по ним экспериментальных данных, чем предварительно апрокси-мировав их аналитическими распределениями, решать итерационную г задачу минимизации X , для определения весовых коэффициентов. И вообще, чем меньше субъективного элемента вносится в анализ, тем достовернее будут вывода о природе физических явлений.

Классификационной является также задача выбора наиболее информативной системы признаков для возможного различия моделей сильного взаимодействия. Мерой различия классов ОВ может служить вероятность ошибочной классификации при использовании оптимального (байесовского) классификатора. Наилучшим будет считаться тот комплекс признаков, используя который, можно добиться наилучшей классификации, то есть минимума ошибок.

Другой тип статистических процедур связан с необходимостью приписать случайной величине не индекс класса, к которому она принадлежит (классификация), а величину, имеющую непрерывное распределение в некотором интервале. Если при классификации множество возможных решений было конечным, то теперь - оно бесконечно. К этому типу задач относится оценивание энергии мюона, распределенной согласно некоторому энергетическому спектру.

Следует опять-таки отметить, что и характеристика и метод решения задачи зависят от доступной априорной информации. Если энергетический спектр неизвестен, то задачу определения энергии мюона можно сформулировать как классификацию зарегистрированных событий по заданному числу энергетических интервалов. Обучающая выборка в этом случае состоит из N классов. Энергию каждого класса можно полагать равной средней в интервале или разыгрывать равномерно внутри каждого интервала (в работе /19/ было использовано 13 интервалов от 5 до 10000 ГэВ). Если же нас интересует не измерение энергетического спектра, а оценивание энергии каждого зарегистрированного мюона для исследования ядерных или электромагнитных взаимодействий, то тогда доступна априорная информация о энергетическом спектре и эту информацию следует использовать в статистической процедуре оценивания: обучающая выборка генерируется согласно известному энергетическому спектру.

Задача предсказания неизвестного значения переменной у , связанной с известной переменной у стохастической связью, решается методами восстановления регрессии /20/. Регрессией у на Г называется функция, задаваемая множеством упорядоченных пар х,7(х)). где у (xJ - математическое ожидание независимой переменной, при условии, что зависимая приняла определенное значение. Обычно предполагают форму регрессионной зависимости известной с точностью до конечного числа параметров, и оценивают эти параметры методами максимального правдоподобия или наименьших квадратов. Однако общее понятие регрессии шире, чем подгонка параметров известной кривой, и, следуя принципу минимальности используемой модели данных, мы не будем специфицировать вид функции регрессии, а будем ее восстанавливать локально по мере поступления экспериментальной информации.

В диссертационной работе проводится последовательное построение непараметрических методов решения задач обработки информации, получаемой с типичных экспериментальных установок в космических лучах. Нашей целью будет создание согласованной системы методов, позволяющих проводить многомерный анализ больших объемов информации, относящейся к различным характеристикам разнотипных частиц. Построение крупных комплексных установок (эксперимент АНИ) /21/ требует проведения большого цикла вычислительных экспериментов с машинными моделями прохождения излучения через атмосферу и экспериментальные установки. Только создание эффективной методики совместного анализа результатов вычислительных экспериментов с моделями и экспериментальных данных позволит делать выводы о характере сильного взаимодействия и составе первичного потока в области энергий вплоть до 10*® эВ.

Диссертация состоит из введения, трех глав, приложения и заключения. Всего 107 страниц текста, 14 рисунков и 7 таблиц. Список литературы содержит 94 ссылки.

Основные результаты диссертационной работы опубликованы в работах /19, 39, 43, 44, 45, 51, 52, 72, 73, 80, 81, 84, 85, 86, 88, 92/.

В заключение выражаю искреннюю признательность своим научный руководителям - профессорам Асатиани Т.Л. и Мамиджаняну Э.А. за всестороннюю помощь и поддержку в ходе выполнения работы. Благодарю чл.-корр. АН АрмССР Матиняна С.Г. за полезные и стимулирующие обсуждения. Считаю своим приятным долгом отметить ценное сотрудничество с коллективом ВЦ ЕрфИ и кандид. физ.-мат. наук Акопошм Н.З.

Постановке и решению задачи создания многомерных методов анализа данных РЭК экспериментов во многом способствовали многочисленные обсуждения с кандид. физ.-мат. наук Дунаевским A.M.

Вопросы создания статистических методов идентификации адро-нов космического излучения эффективно обсуждались с канд. физ.-мат. наук Авакяном В.В.

Новые методы обработки калориметрической информации и их применения к данным Арагацской установки обсуждались с канд. физ.-мат. наук Тер-Антоняном C.B.

Результаты расчетов параметров магнитного спектрометра эксперимента АНИ обсуждались с канд. физ.мат. наук Ивановым В.А. и Мнацаканяном Э.А.

Большой вклад в дело создания математического обеспечения внесли Галфаян С.Х. и Арутюнян С.Х.

Зазян М.З. и Чилингарян Н.С. активно участвовали в применении созданных непараметрических методов к обработке данных экспериментов "Памир" и "Пион".

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

Всем им автор выражает свою глубокую благодарность.

ЗАКЛЮЧЕНИЕ

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

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

1. Классификационный способ идентификации адронов космического излучения. Применение этого метода к обработке данных установки "Пион" позволило уменьшить ожидаемую ошибку идентификации до 13% и учесть априорную информацию об энергии частицы и плотности ливневого сопровождения.

2. Многомерный анализ данных рентген-эмульсионных экспериментов, включающий: а) поиск модельно-зависимых комбинаций признаков по результатам моделирования ядерно-электронных каскадов (одновременнс можно рассматривать наборы, включающие до 30 различных признаков) б) выбор теоретической модели, наиболее подходящей экспериментальным данным; в) определение массового состава первичного потока космического излучения в условиях неопределенной модели сильного взаимодействия. Анализ, проведенный в рамках сопоставления моделей

М А , Мб, Femi, , Femcxx ' показал' что Д°лю семейств, инициированных тяжелыми ядрами, можно будет определить с точносты не хуже 0,05; г) выделение событий определенного типа из первичной информации, например, бинокулярных событий.

3. Метод повышения максимально-измеренного импульса мюонных магнитных спектрометров путем коррекции юстировочных констант с помощью алгоритмов случайного поиска и применения непараметрических методов определения импульса мюона. Применение методов к данным Арагацского мюонного спектрометра (безфильмовый канал) позволило повысить максимально-измеримый импульс с I до 2.5 ТэВ/с,

4. Метод обработки данных калориметрических установок типа PAIR-METER и искровый калориметр, включающий в себя: а) классификационный способ измерения энергии мюонов.

Применение этого метода к данным Арагацского искрового калориметра позволило увеличить долю правильно классифицируемых мюонов больших энергий ( > I ТэВ) с 17 до 30%; б) регрессионный способ оценивания энергии, позволяющий учесть априорную информацию об энергетическом спектре мюонов; в) метод подавления фона от ядерных взаимодействий низко-энергетичных мюонов, позволивший довести уровень дискриминации с

ПГ1 ДО Ю"3;

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

Практическое преимущество новых методов продемонстрировано путем обработки данных машинных моделей и данных с экспериментальных установок в/с Арагац.

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

Результаты, вошедшие в диссертацию, докладывались на всесоюзных конференциях по космическим лучам (Якутск, 1977, Ереван, 1979), Международной конференции по математическим методам решения физических задач (Дубна, 1977), Международных конференциях по космическим лучам (Пловдив, 1977, Париж, 1981, Бангалор, 1983), на сессии ядерного отделения АН СССР (МИфИ, 1983), на всесоюзной конференции по математической теории распознавания образов (Звенигород, 1983), на всесоюзной школе-семинаре по программно-алгоритмическому прикладному статистическому анализу (Цахкадзор, 1983) и на рабочих совещаниях эксперимента АНИ (1980 - 1983 гг.).

1. Бунге M. Философия физики. M.: Прогресс, 1975. - 345 с.

2. Крейн М., Лемуан 0. Введение в регенеративный анализ моделей. М.: Наука, 1982. - 103с.

3. Бикел П., Доксам К. Математическая статистика. М.: Финансы и статистика, 1983. - 277с.

4. Айвазян С.А., Енюков И.С., Мешалкин Л.Д. Прикладная статистика. М.: Финансы и статистика, 1983. - 470с.

5. Стогов Г.В., Макшанов A.B., Мусаев A.A. Устойчивые методы обработки данных. Зарубежная радиоэлектроника, 1982, № 9, с.3-46.

6. Berkson J. Minimum CHI-SQUARE, not maximum Likelihood . Ann. Stat., 1980, v.8, p.457-487.

7. Pearson E.S. Studies 1n the history of probability and statistics. B1o-metrica, 1965, v.52, p.3-18.

8. Вапник В.H. Восстановление зависимостей по эмпирическим данным. М.: Наука, 1979, - 447с.

9. Асатиани Т.Л., Тер-Антонян C.B. Моделирование прохождения ju-мезонов через Арагацский искровой калориметр. Препринт ЕФИ-496 СII)—81, Ереван, 1981.

10. Авундкян А.Т., Казарян С.П., Оганесян А.Г. Алгоритмы моделирования РПИ-детекторов и некоторые его применения. Препринт ЕФИ-437(44)-80, Ереван, 1980.

11. Чилингарян A.A. Имитационная модель магнитного спектрометра. Тезисы докладов У конференции молодых ученых ЕФИ, Ереван,198I, с.64-65.

12. Dunaevsky A.M., Emelyanov Yu.A., Shorln B.F., Urysson A.V. The calculation of nuclear electromagnetic cascades. Preprint FIAN, No.149, 1980.

13. Kokoulin R.P. Monte-Carlo calculation for PAIR-METER. ICR-Report, 102-82-5, Tokyo, 1982.

14. Колмогоров A.H., Шомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1968. - 496с.

15. Efron В. Maximal likelihood and decision theory. Ann. Statist., 1982, v.10, p.340-357.

16. Хелстром К. Квантовая теория проверки гипотез и оценивание. М.: Мир, 1979. - 337с.

17. Райфа Г., Шлейфер Р. Прикладная теория статистических решений. М.: Статистика, 1977. - 306с.

18. Тер-Антонян С.В., Чилингарян А.А. Об одной возможности улучшения энергетического разрешения многослойных детекторов. ВАНТ, сер.Техника физ. эксп., Ереван: 1982, вып.3/12/,с.90-91.

19. Мостеллер Ф., Тьюки Дж. Анализ данных и регрессия. М,¡Финансы и статистика, - 1982, 235с.

20. Данилова Т.Е., Дунаевский A.M., Ерлыкин А.Д. и др. Проект эксперимента по исследованию взаимодействий адронов в области энергий Ю3-Ю5 ТэВ (эксперимент АНИ). Изв. АН Арм.ССР, сер.Физика, т.17, вып.3-4, с.129-232.

21. Асатиани Т.Л., Алчуджян С.В., Амроян К и др. Исследование горизонтального потока мюонов с помощью магнитного спектрометра и искрового калориметра. Изв.АН СССР, сер.физ., 1981, т.45, с.1265-1268.

22. Kltamura Т., Murakl Y., Okada A. et al. PAIR-METER, A new muon spectrometer. Proc. 14-th ICRC, Munich, 1975, v.6, p.2145-2149.

23. Алиханян А.И. Обзор работ по рентгеновскому излучению и детекторам на его основе. Труды международного симпозиума по переходному излучению частиц высоких энергий. Ереван, 1977, с.41-79.

24. Авакян В.В., Авундкян А.Т., Антонян К.Г. и др. Установка для исследования взаимодействий космических адронов с энергией, выше 500 ГэВ. Труды международного симпозиума по переходному излучению частиц высоких энергий. Ереван, 1977, с.277-286.

25. Леман Э. Проверка статистических гипотез. М.: Наука, 1979, 408 с.

26. Yuan L.C., Uto Н., Dell G.F., Alley R.W. Relat1v1st1c Particle separation by means of transition radiation utilizing multiwire proportional chambers. Phys.Lett., 1972, v.408, p.689-698.

27. Cherry M.L., Muller D., Prince Т.Д. The efficient identification of relativistlc particles by transition radiation. NIM, 1974, v.115, p.141-150.

28. Fukunaga К., Hostetier L.D. K-nearest-neighbor Bayes-r1sk estimation. -IEEE Trans. Inform. Theory, 1975, v.IT-21, p.285-293.

29. Rosenblatt M. Remarks on some nonparametrlc estimates of a density function. Ann. Math. Stat., 1956, v.42, p.1815-1842.

30. Фукунага К. Введение в статистическую теорию распознавания образов. М.: Наука, 1979. - 367с.

31. MahalanoMs P.C. On the generalized distance 1n statistics. Proc. of the National Institute of Sciences of India, 1963, v.12, p.49-55.

32. Cover ü.M., Hart P.E. Nearest neighbor pattern classification. IEEE Trans. Inform. Theory, 1967, v.IT-13, p.21-27.

33. Loftsgaarden D.O., Quesenberry C.P. A nonparametrlc estimate of a multivariate density function. Ann. Math. Stat., 1965, v.36, p.1049-1051.

34. Rabiner L.R., Levinson E.S., Rosenberg A.E., VHlpon J.G. Speaker-Independent recognition of Isolated words using clustering techniques. -IEEE Trans, on Acoust., Speech and Signal Process., 1979, v.ASSP-27, p.336-349.

35. Yunck T.R., Teteor F.B. Comparison of decision rules for automatic EEG classification. IEEE Trans. Pattern Analysis Machine Intell., 1980, v.PAMI-2, p.420-428.

36. Dublsson D., Lavlson R. Surveillance of a nuclear reactor by use of a pattern recognition methodology. IEEE Trans. Syst. Man and Cybernet., 1980, v. SMC 10, p.603-609.

37. Friedman. T.H. Data analysis techniques for high energy particle physics.-Proc. of the 1974 CERN school of computing, 1974, CERN 74-23, p.271-366.

38. Чшшнгарян А.А. Анализ методик интерпретации данных в применении к эксперименту АНИ. ВАНТ, сер.Техника физ. эксп., Харьков, 198I, вып. 2/8/, с.59-68.

39. Tapia R.A., Thompson J.R. Nonparametrlc probability density estimation. -The Johns Hopkins University Press, Baltimore and London, 1978, 176 p.

40. Chen Chi-Ru. Asymptotic mean square error of NN density estimation. -Chin. Ann. Math., 1981, v.2, p.425-430.

41. Патрик Э.А. Основы теории распознавания образов. М.: Советское Радио, 1Э80 - 407с.

42. Chllingaryan А.А. The development of statistical methods 1n cosmic ray physics. Proc. 18-th ICRC, Bangalor, 1983, v.5, p.524-526.

43. Чилингарян А.А. Идентификация частиц высоких энергий как задача распознавания образов. ВАНТ, серия Техника физ.эксперим., Ереван, 1983. вып.3/12/, с.61-68.

44. Chllingaryan А.А. On the efficiency of the high energy particle identification statistical methods. Preprint EPI-582(69)-82, Yerevan, 1982.

45. Авакян В.В., Авунджян А.Т., Багдасарян JI.C. и др. Исследование состава и спектров адронов космического излучения на высотах гор. ВАНТ, сер. Техника физ. эксперим., вып. 4/14/,с.25-44.

46. Авакян В.В., Авунджян А.Т., Гаряка А.П. и др. Определение сечения неупругого взаимодействия пионов и нуклонов с ядрами железа в интервале энергий 0.5 * 5.0 ТэВ. ВАНТ, сер. Техника физ. эксперим., Ереван, 1983, вып.4/14/, с.45-55.

47. Славатинский С.А. Ядерные взаимодействия в космических лучах. ВАНТ, серия Техника физ. эксперим., Ереван, 1982, вып.3/12/, с.36-48.

48. Славатинский С.А. Результаты исследований неупругих взаимодействий адронов в области энергий 10^-10® ГэВ. ВАНТ, серия Техника физ. эксперим., Ереван, 1983, вып.5/17/, с.3-16.

49. Ройнишвили Н.Н. Анализ экспериментальных данных получаемых методом рентгеноэмульсионных камер на основе моделирования прохождения ЯЭК в атмосфере. ВАНТ, сер.Техника физ.эксперим., вып. 3/12/, Ереван, 1982, с.49-60.

50. Чилингарян А.А. О новой методике обработки данных рентгеноэмульсионных камер. ВАНТ, сер. Техника физ. экспер., Ереван, 1983, вып. 1/13/, с.81-84.

51. Chillngaryan A.A., Galfayan S.Kh., Zazyan M.Z., Dunaevsky A.M. Selection of model-dependent features 1n photon-hadron famines, registered by roentgen-emulslon chambers. Proc. 18-th ICRC, Bangalor, v.5, p.487-490.

52. Дунаевский A.M. Емельянов 10.А., Урысон А.В., Шорин Б.Ф. Препринт ФИАН, № 18, 1980.

53. Никольский С.И. Спектр первичного космического излучения по данным о широких атмосферных ливнях. Изв.АН СССР, сер.Физика, 1971, т.35, с.2117-2122.

54. Дунаевский A.M., Емельянов Ю.А., Пашков С.В. и др. Расчеты ядерно-электромагнитных каскадов высокой энергии. Труды ФИАН, 1984, т.154, с.142-217.

55. Dunaevsky A.M., Pashkov S.A., Slavatinsky S.A. Some properties of hadron-nuclear Interaction at 1014 1016 eV. - Proc. Intern. Semln. on Cosmic Ray and Particle Phys., Tokyo, 1984.

56. Асатиани Т.Л., Иванов В.А., Мнацаканян Э.А., Тер-Антонян С.В. О возможной программе исследования мюонов в широких атмосферр хных ливнях с энергиями 10-10 ТэВ. ВАНТ, сер. Техника физ. эксперим., 1982, выи. 3/12/, с.79-89.

57. Кутин Г.И. Методы ранжировки комплексов признаков. Зарубежная радиоэлектроника, 1980, № 4-.

58. Джонсон Н., Лион Ф. Статистика и планирование экспериментов в технике и науке. М.:Мир, 1980. - 609с.

59. Иберла К. Факторный анализ. М.: Статистика, 1980. - 395с.

60. Раудис Ш. Влияние объема выборки на точность выбора модели в задаче распознавания образов. Стат. проблемы управления, вып.50, Вильнюс, 198I. с.9-29.

61. Toussaint 6.Т. Bibliography of m1sclassif1cat1on. IEEE Trans, on Inform. Theory, 1974, v. IT-20, p.472-479.

62. Котюков В., Буторин А. Экспериментальное сравнение различных оценок ошибки классификации. Сб. Стат. проблемы управления, Вильнюс, 198I, вып.50, с.78-88.

63. Alpgard К., Ansorge R.E., Asman В. et al. First results on complete events from pp collisions at the CM energy of 540 GeV. Phys.Lett., 1981,1. V.107 B, p.310-314.

64. Galsser Т.К., Stanev T. Cosmic ray composition and particle physics1014 1016 eV. - Workshop on Cosmic Ray Interaction High Energy Results,1.-Paz, 1982, p.281-304.

65. Jodh 6.B., Ellsworth R.W., Stanev T., Galsser T.K. Cosmic rays 1n the atmosphere and characteristics of high energy Interactions. Nucl.Phys., 1981, v. B183, p.12-28.

66. Балаян Г.JI., Оганесян А.Г., Ходжамирян А.Ю. Прямые мюоны от распадов тяжелых кварков в области энергий космических лучей. ВАНТ, серия Техника физ.эксперим., Ереван, 1983, вып.5/17/, с.42-48.

67. Asatiani T.L., Ter-Antonyan S.V. The Investigation of hadroproductlon by cosmic ray muons on lead nuclei with the spark calorimeter method. Proc. 18-th ICRC, Bangalor, 1983, v. 7, p.87-90.

68. Matsuno S., Kajlno F., Kawashlma Y. et al. Cosmic ray muon spectrum up to 20 TeV at 89° zenith angle. Phys.Rev. D, 1984, v.29, p.1-23.70

69. Allkofer O.C., Carstensen K., Bella 6. et al. Muon spectra from Dels up to 7 TeV. Proc. 17-th ICRC, Paris, 1981, v.10, p.321-324.

70. Kokoulin R.P., Petrukhln A.A. Comparative Analysis of possibilities of various techniques of muon spectrum measurements. Proc. 18-th ICRC, Paris, 1981, v.5, p.17-20.

71. Асатиани Т.Jl., Прохоров В.Н., Чилингарян A.A. Программное обеспечение безшильмового телескопа Арагацкого мюонного спектрометра. Препринт ЕФИ-239(32)-77, 13с., Ереван, 1977.

72. Chilingaryan A.A., Ter-Antonyan S.V. On the possibility of Improving the multi-layered detector resolution. Preprint EPI-519(6)-82, Yerevan, 1982.

73. Асатиани Т.JI., Иванов В.А, Мнацаканян Э.А. Магнитный спектрометр установки АНИ. ВАНТ, сер. Техника физ.эксперим., Харьков, 198I, вып.2/8/.

74. Rossi В. High energy particles. Prentice-Hall, Englewood Cliffs, New Jersey, 1952.

75. La May L., Masek G.E., Maung T. et al. A TeV energy cosmic ray spectrometer. NIM, 1972, v.99, p.547-552.

76. E1ch1ngler H., Regler M. Review of Track-fitting methods 1n counter experiments. Preprint CERN, Mo.81-06, Geneva, 1981, p.89.

77. Хьюбер П. Робастность в статистике. М.:Мир, 1984. - 303с.

78. Hodges J.L., Lehmann E.L. Estimates of location based on rank tests. -Ann. Math. Stat., 1963, v.34, p.598-617.

79. Асатиани Т.JI., Иванов В.А., Мнацаканян Э.А. и др. Использование магнитного спектрометра АНИ для исследования потока мюонов космического излучения. ВАНТ, серия Техника физ.эксперт*. , Ереван, 1983, вып. 2/13/, с.85-89.

80. Прохоров В.Н., Чилингарян А.А. Юстировка Арагацкого магнитного спектрометра. Препринт ЕФИ~24-1(34-)-77, Ереван, 1977.

81. Говорун Н.И., Иваненко И.М., Чвыров А.С. Определение параметров бесфильмовых искровых камер. Н.С. ОИЯИ Р5-5397, Дубна, 1970.

82. Asatlani T.L., Alchudjian S.V., Gazarian К.A. et al. Momentum spectrum and charge ratio of cosmic ray muons at zenith angle 84°. Proc, 15-th ICRC, Plovdiv, 1977, v.11, p.362-366.

83. Асатиани T.JI., Алчуджян С.В., Газарян К.А. и др. Исследование характеристик высокоэнергичных мюонов космического излучения. Изв.АН СССР, сер. физ., 1980, с.590-593.

84. Asatiani T.L., Abrahamian A.V., Alchudjian S.V. et al. The Investigation of muon groups at large zenith angles with the Aragats spark calorimeter. -Proc. 18-th ICRC, Bangalor, 1983, v.7, p.60-63.

85. Chilingaryan A.A. Programmed methods to Increase maximum detectable momentum of magnetic spectrometers. Proc. 17-th ICRC, Paris, 1979, v.10, p.364-367.

86. Барнавели Т.Т., Халдеева И.В., Эристави Н.А. Методика определения средней энергии проникающих частиц в многослойныхдетекторах. Препринт ИФ АН ГССР 7, 1982.

87. Devroge L.P. The uniform convergence of nearest neighbor regression and their application in optimization. IEEE Trans, on Inform. Theory, 1978, v. IT-24, p.142-151.

88. Асатиани Т.Л., Тер-Антонян С.В. Учет спектра при определении энергии мюонов многослойными установками. Препринт ЕФИ-653(43)-83, Ереван, 1983.

89. Чилингарян А.А. Развитие статистических методов в физике космических лучей. Материалы всесоюзной школы-семинара.

90. Программно-алгоритмическое обеспечение прикладного многомерного статистического анализа." Москва, 1983, с.367-370.

91. Сборник научных программ на фортране. М., Статистика, 1974-г 535с.

92. Гласс Р., Нуазо Р. Сопровождение программного обеспечения. М.: Мир, 1983. 156с.