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

оптимальныхметодовприменительнокконкретнымсхемамлучевойра- диотомографии. Отметимодноважноеобстоятельство, связанноесвыборомразме­ ровсеткиреконструкции. Ограничениеснизунаразмерсетки, обус­ ловленноевыбраннымприближениемлинейнойтомографии, приводилось выше (3.22). Крометого, самивычислительныепроцедурынакладывают ограничения на размеры сетки. С одной стороны, с ростом числа элементов разбиения заданной области растет точность вычисления интегралов (3.21) иэлементовматрицы А (3.16) или L. Сдругой стороны, сувеличениемразмерности системырастет числообуслов­ ленности матрицы, что увеличивает погрешность реконструкции, вследствие неточностей в данных и коэффициентах матрицы [139, 1401. Следовательно, должно существовать оптимальное число элементовразбиения, приводящеекнаименьшейпогрешности реконст­ рукции. Нахождениеминимума ошибкиреконструкцииисоответствую­ щего размера сетки возможно после ряда расчетов томографической реконструкции. Моделирование спутникового радиотомографичѳского эксперимента было осуществлено в серии численных экспериментов на различных сетках (1441. Результаты моделирования показали возможность реконструкции двумерной структуры ионосферы и небольшим числом приемников. Вдальнейшем, имеяввидуприложениякэкспериментам, будем вести речь о моделировании томографической реконструкции тремяприемниками. Достаточно очевидным следствием моделирования является выводо зависимостиконечного результата от начального приближения. 0 том, насколько точно "попадание" с начальным приближением, можно судить по характеру и скорости изменения распределения (в заданной метрике) с итерациями. Результаты моделирования показали, что очень хорошо восстанавливаются высококонтрастные, т.ѳ. резко отличающиеся от окружающего фона, детали, даже если их размер составляет один интервал дискретизации, причем в начальном приближении эти детали могут полностьюотсутствовать. Реконструкцияобщегофонапроходитменее точно и более медленно, нередко после большого числа итераций возникают характерные артефакты в виде "волн" или "биений", избавиться от которых можно путем мѳжитѳрационного сглаживания. Цривосстановлениираспределенияэлектроннойконцентрациисущест­ венноиспользованиесвойстваположительностиреконструируемойфун­ кции. Посколькуситерациямискоростиизменениязначенийконцент­ рации в областях с высокой и низкойконцентрациямиотличаются, следует использовать в этих областях различные константы ре­ лаксации. Рассмотрим результаты сравнения реконструкции фазовым и фазоразностнымметодами. Заметим, чтоиспользованныевмодѳлирова- 85