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

книги из ГПНТБ / Курганов Р.А. Прогнозирование наклонного рассеивания радиоволн метеорными ионизациями

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

Орт ѵ = [п-и] лежит в части зеркальной плоскости, расположенной ниже плоскости горизонта для точки Р.

Положение вектора радианта s в зеркальной пло­ скости тс определяется углом ^ 0 между вектором радианта s и направлением орта и. Положение эле­

ментарной площадки

Ф = p-da-jip-dri определяется

углом у] между нормалью к ней t

и ортом

и.

 

Направление отсчета

щ,у\н$'

показано

на

рис.12,

где изображено сечение элемента объема зеркальной плоскостью тс.

d6 Линия горизонта \ точки Р

 

 

 

Рис. 12.

 

 

 

Изменим

порядок

интегрирования

выражения

(31), т.

е.

сначала

определим

„видимые" с

на­

правления

радианта размеры элемента

объема

dV,

т. е. площадь собирающей поверхности dV

для

метеоров

данного радианта, лежащего в

зеркальной

плоскости тс. Векторы

радиантов,

расположенные

под

углами 0 < -% < 180°, будут находиться над

горизонтом

точки Я, т. е. будут удовлетворять первому условию

интегрирования.

Границы

изменения

0 < у\0 < тс тож­

дественно

соответствуют границам изменения угла ß,

т. е. угла

между

вектором

радианта

и плоскостью

связи э 0 — - l < ß < ß 0 + JL, где ß0 — угол, соответствую­ щий положению вектора радианта в данной зеркальной плоскости, при котором зенитный угол радианта мини­ мален. Второе условие интегрирования, т. е. условие

60

интегрирования только по входящему в элемент объема метеорному потоку, определяет допустимые границы изменения у, при которых элементарная площадка da „освещается" со стороны рассматриваемого радианта снаружи. Эти границы будут т/0 — 90° < у < % + 90°.

„Видимые" размеры элементарной площадки dax=da-t-s или da± — delcosß', но ß' в данном случае будет являться функцией ч\. Так как углы г\ и ß — это углы между векторами, лежащими в одной зеркальной плоскости те, то df[ = dß и сф' = dß. После подстановки

в выражение

для

dN(x,

у,

z)

ß' = 7)0

-ц, üfß' = dß

и постоянных

для любого,

лежащего в

плоскости я,

вектора радианта

s пределов

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

ЧЧо — 7 < 7 ! <

+ у . получим

 

 

 

 

 

0 +

Т

 

 

 

 

 

 

 

 

rfßX

X f cos2

(ir)o -

^-1

Çn(s)-p(v)-dv,

 

но

 

 

 

 

 

ЧоJ+ я/2 COS2

(TJO -

i j ) fify = Y

,

 

т|0—я/2

 

 

 

 

Л Ѵ ( х ,

y, z) = it-ç>-dp(d<x)2

X

 

Х - Y

-

%)]-d?$n(s)-p(v).dv.

(33)

Поделив dN(x, y, z) на величину элементарного объема dV — np2dp(daf, где p = 62-/?0, и преобразовав

61

f1 ~72 0 - 5 " ) ] к виду [l ~ Й 'где a

=

T '

ловина длины

трассы,

получим

 

 

 

ß0 + n,'2

 

 

 

Д(х, 3-, z) =

- J [\-^-d^n(s)-p(v).dv.

(34)

 

 

Po—^.'2

г/

 

 

Численное интегрирование этого выражения по dv, ufß позволяет получить искомую объемную плотность точек, отражающих сигналы, превышающие пороговый уровень на время Г ; > 0 .

Для расчета количества отражений с длитель­ ностью больше некоторой Г(. > 0, расчета средней длительности отражений и коэффициента заполнения необходимо определить объемную плотность точек, отражения от которых превышают пороговый уровень на время, большее некоторого Tt.

Образовавшийся в элементе объема след разру­ шается за счет одновременного действия амбиполяр­ ной диффузии и прилипаний и если отраженный им сигнал имеет в данном случае длительность превыше­ ния порогового уровня Tt, то его длительность при разрушении следа только за счет амбиполярной диф­

фузии

была бы Tgl

7^-ехр (~~^ .

где тп

— постоян­

ная

времени прилипаний,

но

Tgl

= -zg'ln — или

Г ; = гг -1п—— для

недоуплотненных

следов,

если

 

 

аотѵ

 

 

 

 

 

 

а ; < а с

и Tgi = 1,13- 10-14г-<Х;,

если

а ; > а с ,

т. е. след

в точке

отражения

переуплотненный.

Отсюда

вели­

чина

электронной

плотности в

точке

отражения

следа,

создающего

сигнал

длительности

Т0

будет

а / = а о » '

ехр (^--е4^

для

а, < а.

и

а(. = 1ДЗХ

 

 

iL*

 

 

 

 

 

 

X 1СГ14- — п для а ^ а . . По величине а, может быть

определено по формуле (28) давление в_точке макси­ мальной ионизации следа pml, затем nL{s) по форму-

62

ле (30) и интегрированием выражения (34) интересую­ щая нас объемная плотность точек, отражающих сигналы, превышающие пороговый уровень на время большее Tt.

Высота над поверхностью земли, соответствующая

максимальному

 

значению АДл, у,

z), т.

е. условию

dt±i (х, у, z)

А

есть характеристическая

высота отра-

———-—- = 0,

 

dz

 

 

 

 

жений с длительностью больше Ті

для объема метеор­

ной зоны

ионосферы с площадью основания Дл;-Ду = 1.

Объемная

плотность

точек,

отражающих

сигналы

с длительностью от

Г,-

до

Т1+ъ

\1+1{х,

у, z) =

— \(х,у,

z) kl+](x,

у,

z).

Количество

отражений

с длительностью от Т, до Ті+и зарегистрированных во всем диапазоне высот с объема с площадью осно­ вания dx-dy = 1,

W/. і+\ (•*. У) =» S [Д/ (•*, У, z) - Д/ + 1 (Л:, у, *)]• Дг.

z

Общее количество отражений, зарегистрированных на трассе с длительностью от Ті до Т1+1,

X у

Двумерная плотность распределения числа отражений

с длительностью

от Т{

до Т1+х

в опорной

плоскости

Рі, 1+1 (х, У) = Ni'1+1

yS>

Общее

количество отраже-

ний, зарегистрированных

во всем диапазоне высот

с объема

с площадью основания Дл>Ду =

1,

 

N(x,

y) =

^Nhi+l(x,

у).

 

 

 

 

 

І

 

 

Часовое

число отражений на трассе

 

N = Z%N(x,y)-àx-Ay.

 

X

у

 

 

Двумерная

плотность

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

численности

отражений

в опорной

плоскости pN(x,

у)—

ЛГ(^)>

Плотность

вероятности распределения

длительности

63

 

 

 

 

 

 

 

 

1

) =

' ' 1 + 1

. Сум-

марная

длительность

отражений,

зарегистрированных

со всего диапазона

высот с объема с площадью

осно­

вания

 

Д х . Д у = 1 ,

Т(х,

y)

=

^Nu+l(x,

у)'1^1^

.

 

 

 

 

 

 

 

 

 

і

 

 

 

Общая

длительность

отражений,

зарегистрированных

на трассе

Т = £ £ Г ( х ,

_у)-Дл>Ду. Коэффициент

запол-

 

 

 

т

X у

 

 

 

 

 

 

 

 

нения

•>) =

 

 

 

 

плотность распределения

——. Двумерная

 

 

 

3600

 

 

 

 

 

4

 

 

 

коэффициента

заполнения

 

в

опорной

плоскости

Рп іхі

У

)

Т

 

Средняя

длительность отраже-

~т

НИИ т =

.

N

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

Получаемое при прогнозе часовое число метеорных радиоотражений N есть среднее для данного часа значение. Вследствие случайности величины плотности потока метеорных частиц, падающего с элемента площади небесной сферы Ь, величины их скоростей ѵ и момента падения t регистрируемое за данный час количество метеоров N является сложной функцией 3-х случайных величин b, v, t. Форма и параметры закона распределения N при среднем значении, равном прогнозируемому N, также должны быть определены в результате прогноза.

Если считать, что путем принятия соответствую­ щих мер аппаратурная погрешность измерений может быть уменьшена до любой достаточно малой величины, то регистрируемое за отрезок времени àT на трассе число метеорных радиоотражений

^дг= S п(х> У> z' Р).

X, V, г, ß

64

где ti{x,

у , z, ß) —число

зарегистрированных за

время — AT в элементе объема

dV'= dx-dy-dz

метеор­

ной зоны

ионосферы отражений от метеоров,

падаю­

щих с элемента небесной сферы, центр которого имеет

радиант,

расположенный

под углом

ß в

зеркальной

для точки X, у , z

плоскости,

 

 

 

 

 

 

§b(x,

 

 

 

 

п(х,_у,

z,

ß) =

jb V

у, z,

Щ-а{х,

у, z,

ß) X

Хт~и5(х,

y, z,

v)-p(b)-p(v)

-db-dv-dx-dy-dz-d$-bT.

 

 

 

 

 

 

 

(35)

После интегрирования по значениям

b, ѵ,

характери­

зующим данный элементарный поток метеоров, получим

п (х,

у,

z,

ß) = b (х,

у,

z,

ß) (х, у,

z, ß) X

 

 

Xm-h5(x,

 

у, z,

$)-dx-dy-dz-d$-bT.

 

 

 

(36)

Тогда n(x, y, z, ß) можно представить как

 

 

(38)

 

P),,v„P- b V

 

 

 

|.( px

 

 

 

Хп(х,:

 

 

)=n

 

 

t

 

 

 

 

 

или

 

у,

z, §)p{b)-p{v)-db-dvW

 

 

 

 

(37)

 

JJP(*)-|*(^-n\x, y,

 

 

 

 

 

 

 

n{x,y,z,

 

z,

ß)-p (b) -p(v)

 

-db-dv,

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß > ) = l ,

 

£(v)*=\.

 

 

 

 

 

 

В этих

выражениях

 

 

Л

Л _ . _

Л

 

 

 

под­

под b, m

, п(х, у , ß, z)

разумеваются

средние

значения за интервал

времени,

по которому

производится

усреднение,

т. е.,

напри­

мер, за полмесяца

при

прогнозе.

Если

за

отрезок

времени Д7\ т. е. время, соответствующее

одной

выборке,

было

зарегистрировано

п

 

метеоров,

то это

значит, что величины

b, m 1 , 0

определены как сред­

ние по п реализациям,

т. е.

 

 

 

 

 

 

 

 

 

 

 

У, г,

=

b

(—У'^^У, ß,

z)

 

 

 

 

 

 

 

 

\_m I

 

 

 

 

 

 

 

 

B-395.-5

65

ИЛИ

^

n(x, y, z, $)*=$aa-n(x,

y, z, ß).

(39)

Так как всегда п< оо, то в общем случае ß„¥=ß,

лл

ц.я =7^= (г, т. е. существует дисперсия регистрируемой численности п(х, у, z), обусловленная дисперсией экспериментально измеренных за время ДГ средних

лл

значений ß„ и ц„. Очевидно, что

п-Ь'2

п (яг- 1 -5 )2

Db— дисперсия плотности потока метеорных частиц, падающих с пределов данной единичной площади небесной сферы. Величины Db исследуется в настоя­ щее время Ю. А. Пупышевым по дисперсии числа метеорных радиоотражений, зарегистрированных р/л станцией в разных азимутах [100].

D{m-15)1,5ч =

^

^ >

Т

-Dv,

(40)

 

 

дѵ

JA A

 

 

где Dv — дисперсия

 

 

V,

В

скоростей

метео­

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

ров данного элементарного

потока.

первых измерений

Как установлено

на

основании

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

Так как входящие в выражение (39) величины

лл л

ß, [А, п(х, у, z, ß) — независимы, то

Dn = f-^-Dn + n2-D^n + n2-Di>.n.

(41).

л

Дисперсия Dn обусловлена случайностью падения ме­ теоров во времени, т. е. является Пуассоновской диспер-

Лл

сией Dn = n. Дисперсия общего числа метеоров, за­ регистрированных за время ДГ на трассе

/Ѵ= £

п,(х,

у, z, ß) или УѴ = 2

 

\п4п-Ч

і (х,у,

г, р)

I

'

1

 

 

 

 

66

равна

i

или, T. к.

л л

 

D ^ = £ (£Ц + «2-Dß„. +

).

(42)

 

i

 

 

 

 

 

Сумма

первых членов £

= X я / равна JV—среднему

за единицу времени

числу

отражений на трассе,

т. е.

 

 

 

 

 

Л

 

прогнозируемому

значению часового числа 2J Dnt

= N.

 

X

л

л

i

 

 

 

Тогда

DN = N+

ZJ

пЦО^-^

D\xn),

где ^ — прогнози­

руемое

значение

численности для

элемента объема

dx-dy-dz — dV и направления ß.

Подставляя в это выражение значения Db и Dv, соответствующие регистрируемому элементарному потоку, можно вычислить величину DN, а зная /V, DN и закон распределения — спрогнозировать параметры закона распределения регистрируемого на трассе числа метеорных радиоотражений за время Д7\

Однако исследование законов распределения дву­ мерной плотности падающего метеорного потока и распределения скоростей метеоров элементарных потоков является чрезвычайно сложной и длительной

процедурой,

требующей

набора,

обработки и анализа

огромного

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

материала,

и не пред­

ставляется

выполнимым

в ближайшее время.

Поэтому

рассмотрим

метод

прогноза

параметров

законов

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

 

регистрируемой численности

по имеющимся результатам

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

непосред­

ственно

на метеорных

радиотрассах. Входящие в (42)

 

„ X

Л

Л

 

 

 

величины 2Jti\-Dpn

и \п\-и^п

могут быть

преобразо-

 

 

I

 

i

 

 

 

 

 

5*

67

ваны к виду

S t\ Л — о л

n).D}nrN'.D^

і

и

здесь

— относительная дисперсия регистрируемого числа отражений на трассе, обусловленная неравновесными

л л

вкладами

относительной дисперсии

Dß„

и

ДА„ эле-

ментарных

метеорных потоков.

 

і

і

 

 

 

Следовательно

 

 

 

DN=N

+ yV2-Z)-|L + N2-Dp~ .

(43)

Если дисперсия

УѴ определяется

по

„к"

выборкам,

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

чество реализации /V, то по

определению

дисперсии,

как математического

ожидания квадрата

отклонения

величины от

ее

среднего

значения

 

 

 

 

A

k

 

л

 

 

 

A*ff = E >"0Ч)--Аѵг ,

 

 

где m (Л^) — частость

появления

реализаций

в вы-

 

 

 

 

 

k

 

 

борке. Так как

D^N

=

то D^

= ^^-^'D}x

или

і

68

£)(J.jv — Dp

При

достаточном количестве (боль­

ше 10) выборок

=

(-4-^ и £Ѵдр = ^

. Аналогично

=

. В данном случае £>ß_ и

представляют

собой относительную дисперсию средних по N реа-

Лл

лизациям значений Ь- и mh? по отношению к средним

ÏЛ - 1 5

значениям

опт

для

интервала

суток, из

которого

взяты nku

выборок. Дисперсия

среднего

по

N реали-

 

 

 

Л

 

 

 

 

 

 

зациям значения Ь- относительно

среднего

значения

b за данный день -=- на

порядок

меньше

величины

Dp-

и ее

можно

не

учитывать.

 

час из „k"

Пусть „k" выборок

взяты

за

данный

дней.

Так

как

по

имеющимся

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

данным не сделано вывода о наличии месячного и се­ зонного изменения закона распределения геоцентри­ ческих скоростей метеоров [101], то можно считать, что, в первом приближении, этот закон одинаков для

разных суток.

Тогда, как и

в

случае

выборок в пре­

делах одного

дня величина

Dp- будет

определяться

как Dp-=y=~,

т. е. как дисперсия выборочных сред­

них по /V реализациям,

взятых

из одной

генеральной

совокупности.

 

 

 

 

 

 

С другой стороны,

определенно

установлено [99],

что среднее значение плотности падающего на землю

потока метеорных

частиц

имеет заметный

месячный

и сезонный ход.

Поэтому

нельзя считать,

что „ku

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

л

совокупности. Это значит, что величина Dß_ в данном

случае будет представлять собой относительную дис-

л

Персию средних за выборку значений Ь~ относительно среднего значения за „k" дней и формула Z)ß_ = -^-

09'

теряет смысл. Форма закона распределения регистри-

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