Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Алания М.В. Квазипериодические вариации космических лучей

.pdf
Скачиваний:
27
Добавлен:
23.10.2023
Размер:
5.11 Mб
Скачать

Следует заметить, что для осреднения по

формуле

(1.1.2) требуется точное знание реальных ошибок,

которыми

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

ведены в отдельных оригинальных работах последователей,

в следующем

параграфе мы предлагаем методику, дающую

возможность

достаточно правильно оценивать реальные

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

мим ошибкам at. Во всех

остальных

случаях вполне успеш­

но можно пользоваться

формулой

(1.1.1).

Станции

Моусон Т у л е

Резолыот-Бен Черчилль

Ипверкаргпл

Веллингтон

Оттава Дип-Ривер Чикаго Линкольн Вашингтон Сульфур

.Москва X псала Гёттинген

Херстмоисо Л и д с

Жесткость обрезания, Бв

Статистиче­ ская точ­ ность, %

Группа

 

I

 

0,57

0,38

 

0,00

0,50

 

0,00

0,50

I

0,11

0,54

 

1,81

0,48

 

1.71

0,37

 

0,96

0,52

 

0,87

0,31

 

1,54

0,48

 

1,99

0,52

 

1,98

0,28

11

0,94

0,23

 

2,03

0,38

 

1,17

0,44

 

2,38

0,41

 

2,30

0,31

 

1.710,56

 

Таблица 1

 

Жесткость обрезания, Бв

­

Группа

Станции

Статистиче скаяточ­ ность, %

 

 

 

Кланмакс

2,77

0,12

 

Беркли

4,01

0,45

 

Р и м

5,04

0,56

III

Цуг-Шпитце

3,33

0,15

Хер.чанюс

4,94

0,44

 

Ломннцкин Щит

3,15

0,30

 

Вепсенау

3,22

0,43

 

Буэнос-Айрес

10,70

0,51

 

Л а э

14,89

0,78

 

Рпо-де-Жанейро

11,47

0.63

 

У анкаио

14,18

0,17

IV

Макерере

14,90

0,28

 

Макапу

11,03

0,53

 

Норикура

9,13

0,20

 

§ 2. Определение ошибок данных наблюдений нейтронной и жесткой компоненты интенсивности космических лучей

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

11

п

V v r -

п

i = \

Прежде всего необходимо точное определение самих оши­ бок. ог данных наблюдений, Как показало наше исследова­

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

рактеризуются статистическими ошибками ст = 100/] «%, где п — среднее число частиц, зарегистрированных за единич­ ный интервал времени измерения без учета аппаратурного ВЛИЯНЛ1Я. Вероятность того, что за данный интервал време­ ни будет зарегистрировано п частиц (в предположении их независимости, это справедливо в первом приближении для нейтронных мониторов, счетчиковых телескопов и ионизаци­ онных камер), определяется распределением Пуассона. Од­ нако, для больших л, Пуассоновское распредление можно заменить распределением Гаусса.

Щп,п = р = ехр[ - (п - п)2/2п\.

(1.2.1)

Учитывая, что флюктуация интенсивности (в относи­ тельных единицах)

л — и

п

распределение по / будет иметь вид

Щ1' Т1) = y ^ f e x p [ - н Г -m -

(1.2,2)

Прежде всего найдем закон распределения эксперимен­ тальных значений по данным наблюдений различных стан­ ций. Мы проанализировали данные двухчасовых наблюдений нейтронной компоненты за период МГГ по ст. Клаймакс, Норпкура п ооре,д1Н0Ш1ые данные жесткой компоненты по ст. Токио и Норпкура.

12

В области малых изменений /, обусловленных в основ­ ном флуктуациями статистического характера, можно утвер­ ждать, что распределение по экспериментальным данным WaKc,i. должно быть близким к теоретическому распределе­ нию WTeor., которое определяется согласно выражению (1.2.2). На рис. 3 представлены кривые зависимости 1пЩп, й) Для различных отклонений I, согласно данным

Рис. 3. Логарифм теоретического и экспериментального частотного распре-

п —■п

деления разностей типа / = —— интенсивности космических лучей по дан-

п

ным нейтронной компоненты станций Клаймакс и Норикура и по осредненным данным жесткой компоненты станций Токио и Норикура.

Кривые, проведенные через черные точки, соответствуют экспериментально­ му частотному распределению, а кривые проведенные через кружки соответ­ ствуют теоретическому частотному распределению.

Обозначения: по оси ординат отложен InW(n, п), а по оси абсцисс валнчи-

п п

на разности типа I = ——— в процентах.

п

нейтронной компоненты станций Клаймакс, Норикура и осрвдненньш данным жесткой .компоненты станций Токио—Но­ рикура. На этом же рисунке также приведены теоретически

ожидаемые результаты распределения

In Ц7теор.

Анализ

показывает, что в основной области,

а

именно в [-10а, 10а]

ln^„cn.®lnU7ieop., так как практически отклонения /

больше,

чем ± 10а, чрезвычайно редки

(рис.

3). Если провести нор­

мировку для

экспериментального

распределения

Ц7ЭКСП.

внутри области

~ [—10а, 10а] изменения I и согласовать

теоретические распределение

№ТСоР.

с

экспериментальным

13

ИРэксп.

при

0, легко можно найти среднее

эффективное

значение

я, которому соответствует данное эксперименталь­

ное распределение W0|iC„. Действительно

 

 

 

 

W(I) = ( А / У Щ ехр[—пР/2],

 

 

 

 

j 2

 

 

 

 

In W (/) = In (A/V2riH)-

(1-2.3)

при / ~

0,

1п{А/У 2ш )= In W (о) н формула для

определения

среднего эффективного счета я, примет следующий вид:

(1.2.4)

п= (2[ln.W(o)~InW(l)})/P

Рис. 4. Изменение отношения In(lPai«:ii. /НПсор. )

п—п

в зависимости от величины разности типа 1= — —

 

 

 

 

 

 

 

 

п

по

данным

нейтронной

компоненты интенсив­

ности космических

лучей

станций 1) Клаймакс,

2)

Норикура и

3)

по

осреднениы.м

данным

жесткой

компоненты

станций Токио

п Но­

рикура.

По осп

ординат

отложено изменение

1п(Гэкс». /ЙУтеор.)

 

а по осп абсцисс—величина

 

разности

типа

I

 

п п

 

 

= — —— в процентах.

п

14

ria рис. 4. (кривые 1,2,3) показаны In (И7ТС0Г. /W3Kzn. ) соот­ ветственно для ст. Клаймакс, Норикура и Токио-Норикура. Лег­ ко заметить, что в области / Пт 1 % для ст. Клаймакс и Нори­ кура

Wmv.

Что же касается осреднениых значений жесткой компо­

ненты Токпо-Норпкура, это

различие выявляется

лишь при

/ > 1,4%. Можно полагать,

что при больших I

действуют-

факторы, не связанные со статистическими флюктуациями.

Таким образом из

рис.

4 видно,

что имеются

большие

с-тклонеиня, 'выходящие за

пределы

ошибок данных наб­

людений. Невидимому

они

вызваны

физическими

процес­

сами, которые в конечном счете обуславливают реальные ва­ риации космических лучей.

Согласно [16], статистическая ошибка двухчасовых изме­ рений для ст. Клаймакс составляет 0,12%, а для ст. Нори­ кура — 0,21%- Анализ же экспериментальных кривых час­ тотного распределения, приведенных на рис. 3 дает соответст­ венно значения 0,28% и 0,29% для ст. Клаймакс и Норнкура, т. е. существенно больше чем принимается в [163. Ана­ логичный вывод получается для осреднениых данных жест­ кой компоненты станций Токпо-Норпкура. Подобный анализ был проделан для данных наблюдений ряда приборов (ней­

тронных

мониторов, кубических

телескопов

и

ионизацион­

ных камер), периода МГГ расположенных

на

различных

пунктах

Земного

шара (см. табл. 2).

Найденные ошиб­

ки в 1,5—2

раза

больше

чем

статистические ошибки,

определяемые по формуле ст =

1 ]/д

, где я — среднее число им­

пульсов

за 2

часа

наблюдений.

Этот

вывод

находится в

согласии с работой [17], в которой а определялось по раз­ ностям показаний двух станций. Однако следует отметить, что метод, используемый в [17] для расчета ошибок данных наблюдения, относительно ограничен, так как не всегда уда­ ется подбирать необходимую для данного метода пару стан­ ций, близко расположенных друг от друга. Методом же, который был изложен здесь, можно определить реальные ошибки данных наблюдения для любой конкретной станции, независимо от другой.

15

 

 

 

 

 

 

Таблица 2

 

 

Статистичес­

Действитель­

Отношение

 

Станция

кая

ошибка

ная о шибка

а действ,

Аппаратура

 

о %

 

 

согласно

%

а статист.

 

 

 

 

(105)

 

 

 

 

 

 

 

 

Агуляр

Г.

 

0,15

0,36

2,4

NM

Беркли

г.

0,45

0,64*

1,4

ММ

Вашингтон

0,28

0,40

1,5

NM

Вейсенау

 

 

0,43

0,88

2.0

NM

Веллингтон

г.

 

0,37

0,80

2,2

NM

Клаймакс

г.

 

0,12

0,28

2,3

ММ

Норикура

г.

 

0,21

0,29

1.8

NM

Сульфур

г.

 

0,23

0,36

1,5

NM

Сульфур

г.

 

0,25

0,31

1,3

Т

Тбилиси

 

 

0.45*

0.64*

1.4

NM

Уанкайо

 

 

0.18

0,30

1,7

ММ

Цугспнтце

 

 

0,15

0,30

2,0

ММ

Цугспитце

 

 

0,20

0,36

1,8

Т

Чикаго

 

 

0,48

0,78

1,6

ММ

•Якутск

 

 

0,05

0,11

2,9

I

Обозначения: ММ— нейтронный монитор

Т

— кубический телескоп

 

/

— ионизационная

камера

 

*

— данные 1964

года

 

§ 3. Оценка ожидаемых величин амплитуд

первой и

■второй гармоник суточной вариации космических лучей

и сравнение с экспериментом

 

При исследовании звездно-суточных, солнечно-суточных,

.■полусуточных и других вариаций космических

лучей перио­

дического характера часто используется метод гармоническо­ го анализа.

Частотное распределение амплитуд и фаз первых и вто­ рых гармоник солнечно-суточных вариаций определялись по экспериментальным данным в ряде работ [18, 19, 20]. Ам­ плитуда и фаза гармоник имеют определенное распределе­

ние, обусловленное, по-существу, двумя факторами:

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

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

16

обусловленные первой и второй причинами. Количественное решение этой задачи особенно важно при исследовании су­ точных вариаций за небольшие интервалы времени, исчис­

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

Здесь мы будем рассматривать и воздействие статисти­ ческих и аппаратурных флуктуаций на частотное распреде­ ление амплитуды и фазы первой и второй гармоник при гар­ моническом анализе по 12 ординатам (21). Проведены так­

же оценки величин амплитуд первой и второй гармоник су­ точной вариации, ожидаемых за счет статистических флук­ туаций интенсивности космических лучей и ошибок аппара­

турного характера [15, 22]. Сраг.пение теоретических и экс­ периментальных результатов позволяет количественно найти

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

Представим дискретные значения изменений интенсив­ ности космических лучей в течение суток через I(t).

Тогда [23] f(t) можно выразить

 

со

 

 

2-М \

 

 

 

 

 

 

 

 

 

 

bhsin

)

 

k=\

 

 

Т

 

 

 

 

 

 

 

 

2Ш

\

(1.3.1)

о

I

О<sin

I т

+ Та j 1

где

 

k=\

 

 

 

 

 

 

 

О-

 

 

1"г

 

 

 

 

 

/ (x)dx ak =

1 ir"

 

 

п0 = ~ j

'| / (т) cos ki dx

 

0

 

2г.

6

 

 

 

 

1

 

 

 

 

bh ==

i*

 

 

k = \,

2. 3....

/ (т) sin kx dx

 

2xzkt

, a T = 1 сутки),

пли если дискретные значе­

(здесь x =

ния интенсивности xt (1, 2, З...2р)

задаются через равные проме­

жутки времени,

то

 

 

 

 

 

2. М. В. Алания, Л. X. Шаташвили

17

1

 

1

ki К

 

V ,.

ak =--

V

К,cos —

у i=

Р /=1

р

 

 

 

 

 

 

h =

. L V 1,./О1*1кЫ

 

 

Р i=i

р

 

 

Задача ставится так: три взаимно независимых 'измене­ ниях соответственно заданному статистическому распределе­ нию для заданного значения«(я—средняя скорость сче­ та интенсивности космических лучей за 2 часа) найти ожи­ даемое частотное распределение амплитуд гармоник rh н оценить их ожидаемые величины. Предположим, что имеют­ ся только случайные отклонения .V,- в скорости счета ог среднего значения я, распределенные по закону Гаусса. Пусть, '.алее W(rk, а) есть распределение амплитуд rh, для ко­ торою выполняется условие нормировки

— со

 

| W(rk, о) = 1,

 

 

 

 

(1.3.3)

 

-- X»

 

 

 

 

 

 

 

 

где ст =

1

 

 

 

 

 

 

 

 

, — .

 

 

 

 

 

 

 

 

 

V п

 

 

 

 

 

 

 

 

Распред.еление W (гк, а) можно представить через распределения

W(ah,a)

и W(bk,'т), при

условии,

что W(au,a) и

W(bh,a)

яв­

ляются

независимыми

друг

от друга

распределениями,

т. е.

 

 

Г 2

 

 

 

 

 

 

 

W(rlv а) =

J W (а*,а) W(\/7£^b<j)dr,

(1.3.4)

 

 

ri

 

 

 

 

 

 

 

где bk =

/ г \ -

a\.

 

 

 

 

 

 

 

Таким образом, найдя распределение W(ak, а), можно

определить W (гк, о).

Очевидно,

что

W(rh,o)

можно най­

ти путем интегрирования по

всем

возможным

изменениям

■И* *^’2) *

^2р"

 

 

 

 

 

 

 

 

Пользуясь соотношением

(1.8),

можно записать:

 

18

 

—|—СО

~\~'УО

 

W К , а) =

А dx3, ---- 4 x zp j

dxxexp[—xj(aft, хй, ■■• ,.\'3j0)/2a2jX

 

 

— oo

 

 

X

exp [ — (.v'i -f xl +

• • •

+ xlp)J2a-}.

(1.3.5)

Постоянная А определяется условием нормировки (1.3.3). Рассмотрим частный случай К = 1,Р = 6, (т. е. при гармони­ ческом анализе по 12 ординатам). Поскольку в этом случае

i=2

 

 

 

то

подставляя выражение х2

и

проводя

последовательные ин­

тегрирования, получим:

 

 

 

W (alt a) = А ехр[ — 3,1aj/a2].

(1-3.6)

Коэффициент 3,1 при

в

экспоненте выражения (1.3.6)

был уточнен Крымским Г. Ф. и др. в работе [24]. Согласно этой работы значения коэффициента 3,1 следует заменить на

3,0.

Учитывая выражение (1.3.6), (1.3.4) примет вид

Л>

1Р(Лст)= j W(a„ a)W(Vrl~ai , ajda,

(1.3.7)

= В exp [—Згу/cf-].

Из условий нормировки находится постоянная В,, хотя как

будет показано ниже, для нашей задачи нет необходимости знать ее значение.

Поскольку фаза ^ выражается через срх = arctg К

W(av a) и W(bv а)—независимые распределения при случайных хь легко понять, что при наших условиях фаза суточных гар­ моник за счет статистических флуктуаций не будет выделенной, т. е. любая фаза суток одинаково вероятна и распределение W{ср, а) должно носить изотропный характер.

Сравним теперь полученные результаты с эксперимен­ тальными данными. Для этого составим частотное распреде­ ление для первых гармоник на основе экспериментальных данных за весь период МГГ. Гармонический анализ прово-

19

лился с учетом влияния нециклических изменений для от­ дельных суток по данным нейтронной и жесткой компонент высокогорной станции Норнкура. Были проведены расчеты

1армонпческпх коэффициентов

и Ь1 н найдены серин

зна­

чений а'[! и

(/С =

1,2 • • • 529

— номер суток,

где у =

1 со­

ответствует

двухчасовым данным

нейтронной

компоненты,

а / = 2 — жесткой

компоненты).

Затем были

вычислены

Л' 11 Фг7 !I построены их частотные распределения (рис. 5).

Рис. 5. Логарифм экспериментального и теоретического частотного рас­ пределения амплитуды первой гармоники суточной вариации а) по данным нейтронной и б) по данным жесткой компоненты станции Норнкура.

Обозначения: крестики со стрелками соответствуют экспериментально­ му. а сплошная кривая — теоретическому частотному распределению

На рис. 5 (а, б) крестики изображают частотное рас­ пределение амплитуды первой гармоники солнечно-суточной вариации по данным а) нейтронной н б) жесткой компонен­ ты космических лучей согласно наблюдениям станции Норн­

кура. Сплошная кривая — ожидаемое частотное распреде­ ление по формуле (1.3.7), причем нормировка этой кривой проводилась путем сшивания с экспериментальными данны­ ми при гх 0. Такая нормировка теоретической кривой поз­

20

Соседние файлы в папке книги из ГПНТБ