книги из ГПНТБ / Курганов Р.А. Прогнозирование наклонного рассеивания радиоволн метеорными ионизациями
.pdfОрт ѵ = [п-и] лежит в части зеркальной плоскости, расположенной ниже плоскости горизонта для точки Р.
Положение вектора радианта 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, $)*=$a-£a-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'
теряет смысл. Форма закона распределения регистри-