Радиотомография глобальных ионосферных структур / В. Е. Куницын, Е. Д. Терещенко, Е. С. Андреева [и др.] ; Акад. наук СССР, КНЦ, ПГИ. – Препр. ПГТ 90-10-78. – Апатиты : [б. и.], 1990.

линейных интегралов по которой уже непрерывная функция. Можно приме нять и интерполяцию более высокого порядка, когда функция F представ ляется рядом по степеням ^ и Т более чем первого порядка. Значения функции F в каждом конечном элементе определяются ее значениями е ря де точек дискретной сетки. Линейный интеграл (12) представляет собой сумму интегралов по всем конечным элементам, которые пересекает / -ый луч I j = Z j X ( k ) F ( h . , T ) c f k (16) где у (h ) = ( R+k ) [ R z 4 lh 2 fi + 2 R h ^ Подставляя в (16) соответствующую интерполяционную формулу для F(h,T) можно проводить интегрирование. При этом линейный интеграл (16) по конечным элементам будет состоять из ряда слагаемых (трех в случае треугольных конечных элементов при кусочно-планарной аппроксимации). Каждое слагаемое содержит сомножителем значение{/*Ь граничных точках конечного элемента. После интегрирования по всем конечным элементам, связанным с J -лучом, суммируя коэффициенты при FM , получим элемент матрицы Zj/kj Выполняя численное интегрирование по всем лучам, получим матрицу L jM. Матрица Z j M связана с набором {о(„} положений ИСЗ и соответствую щей серией лучей. Можно вычислить матрицу и для другого набора близких положений ИСЗ с заданным приращением {aC,*AU0} . После чего определяется матрица для фазоразностной томографической задачи/?JM = = ( L' jm - L JM ) / A do Сведение томографической задачи к линейным интегралам (12) и да лее к интегралам виде (14 ), (15) является следствием линеаризации за дачи. Линеаризации в том смысле, что возмущение среды считается доста точно малым и дает лишь поправку.в подынтегральные функции линейных интегралов (2 ), но траектории лучей при этом не изменяются, т.е. ос таются прямыми. Представляется необходимым оценить, насколько такое приближение допустимо в экспериментах по ионосферной радиотомографик. Опенку влияния искривления траектории получим из рассмотрения модель ной задачи распространения в параболическом по k слое ЛЬ (h)=nm * с полутолщиной ^ и максимумом на высоте Из реше ния лучевых уравнений /20/ несложно получить величину смещения траек тории А X и отклонения траектории A k от прямолинейной с углом 6„ к вертикали к _ л - ( J « a z л Т ~ }C O iP во ( J- ' ' iCOtSGo I / •' ( I ? )