Куницын, В. Е. Томография ионосферы / Куницын В. Е., Терещенко Е. Д. ; Акад. наук СССР, Кол. науч. центр, Поляр. геофиз. ин-т. – Москва : Наука, 1991. – 176 с.

координатеИСЗ- aQ (илиуглуместар). Еслиматрицапрямойзадачи LJM: FM "* непрерывна по углу спутника aQ , то вместо системы (3.1Б) можнополучить системупофазоразностнымилидоплѳровским даннымпутемдифференцирования (3.15) поуглуotg: А А - Dj ■ О'16) Здесь Dj з AIj/AaQ - доплѳровские данные, = ALjj/Aa 0 - конечно - разностное отношение (или производная) матрицы Lyj к приращениюугла. Доплѳровскиѳ данныеопределяются нетолькоизме­ нениемполнойфазы, связаннойсинтегральной электроннойконцент­ рациейполучу, нотакжеилокальнойэлектроннойконцентрацией Ng вточкеспутника. Поправканадоплѳровскиѳ данныеравнапроизведе­ нию Ns на составляющую скорости ИСЗ, направленную по лучу: А, reNBcos (a^ + р- aQ ). Этапоправкаможетбытьвнесенавитера­ ционныйалгоритмиполучающиесявпроцессеитерацийзначенияNg на границеh = hQ ионосферыбудутпостоянно "подправлять" измеряемые значения доплѳровской частоты. Величина изменения доплѳровской частоты, обусловленная локальной электронной концентрацией,может достигатьдолейгерца. Перейдемкрасчетуматрицы разностнойзадачи, которая, как уже говорилось, должна определяться по приращению матрицы L ^ , непрерывнойпоуглуспутника. Непрерывностьматрицы можетбыть обеспечена путем введения конечных треугольных элементов для представленияфункции F(h,T), т.ѳ. когдаискомаяфункция заменя­ ется кусочно-планарной аппроксимацией. Гладкая функция P(h,x) заменяется непрерывной многогранной поверхностью аппроксимации, производная по углу ИСЗ линейных интегралов по которой-уже непрерывная функция. Здесь мы ограничимся рассмотрением решения томографическойзадачистреугольнымиэлементами. Интерполяциябо­ леевысокогопорядкаможетбытьрассмотренапоаналогичной схеме, однаковыкладкистановятсяболеегромоздкими. Отметим, чтовычис­ лительныеметоды, основанныенапримененииконечных треугольныхи и граничных элементов, широко используются в различных областях наукиитехники [132, 133]. Треугольныеэлементыестественнополу­ чаются из сетки прямоугольников путем деления каждого из них диагональю пополам. Функция F( h , T) внутри каждого треугольного элементазаменяетсяналинейнуюаппроксимацию F(h,x) = а+ Ьт+ ch. Значения коэффициентов (а, Ь, с) в каждом конечном элементе определяютсяизсистемытрехлинейныхуравненийдлятрехграничных точек. Несложно сразу выписать выражения для коэффициентов в заданномконечномэлементе. Этивыражениянесколько отличаютсяу 79