Главная – Журналы – Сибирский журнал вычислительной математики 2024 номер 3

С целью решения проблемы, связанной с тем, что матрица коэффициентов многомерной модели c ошибками в переменных (ММОП) содержит постоянные столбцы, модель ММОП расширена до частичной многомерной модели с ошибками в переменных (Ч-ММОП) и предложен новый алгоритм модели Ч-ММОП, основанный на принципе частичной модели с ошибками в переменных (Ч-МОП) и непрямой корректировке. Алгоритм прост и легко реализуем. Для проверки используется пример преобразования координат, а результаты сравниваются с существующим алгоритмом модели ММОП; они показывают эффективность предлагаемого алгоритма. Наконец, алгоритм Ч-ММОП применяется к многоточечной серой модели (МСМ(1,N)) мониторинга оседания грунта. Результаты показывают, что модель Ч-ММОП, предлагаемая в данной статье, лучше учитывает влияние ошибок точек мониторинга, а результаты хорошо соответствуют реальной ситуации.

Юнитоидом называется квадратная матрица, которая может быть приведена к диагональному виду посредством конгруэнтного преобразования. Среди различных диагональных форм юнитоида A имеется лишь одна (с точностью до порядка, принятого для главной диагонали), все ненулевые диагональные элементы которой суть числа с модулем 1. Она называется канонической формой матрицы A относительно конгруэнций, а аргументы ее ненулевых диагональных элементов называются каноническими углами этой матрицы. Если A не вырождена, то ее канонические углы тесно связаны с собственными значениями матрицы A^-∗A, называемой коквадратом матрицы A. Хотя определение юнитоида напоминает понятие диагонализуемой матрицы в теории подобий, кажущаяся аналогия между этими двумя матричными классами обманчива. Мы показываем, что жорданова клетка J_n(1), которая в теории подобий рассматривается как антипод диагонализуемости, является юнитоидом. Более того, ее коквадрат C_n(1) имеет n различных унимодулярных собственных значений. Мы погружаем матрицу J_n(1) в семейство жордановых клеток J_n(λ) с параметром λ, меняющимся в диапазоне (0,2]. В некоторой точке, расположенной левее единицы, J_n(λ) перестает быть юнитоидной матрицей. Мы подробно обсуждаем этот момент в попытке понять, как может произойти подобная трансформация. Обсуждаются и аналогичные моменты, соответствующие меньшим значениям λ. Указаны некоторые примечательные факты, связанные с собственными значениями коквадратов и числами обусловленности этих значений.

Изучение распространения в пространстве и времени парниковых газов, а также оценка потоков с поверхности Земли этих газов с помощью системы усвоения данных представляет собой актуальную задачу мониторинга состояния окружающей среды. Одним из подходов к оценке потоков парниковых газов является подход, основанный на предположении, что потоки постоянны в заданной подобласти и на заданном временном интервале (порядка недели). Это обусловлено как необходимостью эффективной реализации алгоритма, так и свойствами используемых в таких задачах данных наблюдений. Современные задачи оценки потоков парниковых газов с поверхности Земли имеют большую размерность, поэтому обычно рассматривается вариант, в котором оцениваемой переменной являются потоки, а модель переноса и диффузии входит в оператор наблюдения. При этом возникает проблема использования больших окон усвоения, в пределах которых оцениваются значения потоков на нескольких временных интервалах. В работе рассматривается алгоритм оценки потоков по данным наблюдений из заданного временного интервала. Алгоритм является вариантом алгоритма ансамблевого сглаживания, широко применяемого в таких задачах. Показано, что при использовании окна усвоения, в котором происходит оценка значений потоков для нескольких временных интервалов, алгоритм может становиться неустойчивым, при этом нарушается условие наблюдаемости.

В статье рассматривается новый класс нелинейных интегро-дифференциальных уравнений Вольтерра второго рода с ядром типа свертки. Используя теорему Шаудера о неподвижной точке, мы получаем некоторые условия, достаточные для существования и единственности решений. Кроме того, для получения приближенного решения предлагаемого уравнения Вольтерра используется метод Нистрема. Приведены численные примеры для подтверждения полученных результатов.

Дробные операторы переменного порядка могут использоваться в различных физических и биологических приложениях, где скорость изменения представляющей интерес величины может зависеть от пространства и/или времени. В данной статье мы предлагаем явную конечно-разностную аппроксимацию для пространственно-временного дробного волнового уравнения Рисса-Капуто переменного порядка с начальными и граничными условиями в конечной области. Предлагаемая схема является условно устойчивой и имеет глобальную ошибку усечения O(τ²+h²). Также представлен численный эксперимент для проверки эффективности предлагаемой схемы.

В работе исследуется технология расчета разностных задач с внутренними граничными условиями баланса потоков, построенными с помощью односторонних многоточечных разностных аналогов первых производных произвольного порядка точности. Предлагаемая технология одинаково подходит для любых типов решаемых дифференциальных уравнений и допускает однотипную реализацию при любых порядках точности. Она, в отличие от аппроксимаций, опирающихся на продолженную систему уравнений, не приводит к осложнениям при расщеплении многомерных задач на одномерные. Сформулированы достаточные условия разрешимости и устойчивости реализации алгоритмов методом прогонки для граничных условий произвольного порядка точности. Доказательство основано на приведении многоточечных граничных условий к виду, не нарушающему трехдиагональную структуру матриц, и установлении условий диагонального преобладания в преобразованных строках матрицы, соответствующих внешним и внутренним граничным условиям.

В работе рассматривается имплементация алгоритма построения адаптивной вычислительной сетки в численное решение прямой одномерной задачи магнитотеллурического зондирования (задачи Тихонова-Каньяра) методом локальных интегральных уравнений, предложенным авторами ранее. Конструирование адаптивной вычислительной сетки основано на геометрических принципах, рассматривающих оптимизацию кусочно-постоянного интерполянта аппроксимируемой функции электрической проводимости среды. Проведены численные эксперименты для исследования и иллюстрации эффективности комбинированного метода. Апробация осуществлялась на модели Като -Кикучи с известным точным решением.