Все выпуски

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

Представлена новая модификация парного фильтра Калмана, основанная на применении модифицированной взвешенной ортогонализации Грама–Шмидта и LD-разложении ковариационных матриц. Основными результатами данной работы являются: новая LD-модификация парного фильтра Калмана (теорема 1); новый алгоритм LD-ПФК для оценки состояния парных марковских моделей на основе прямой процедуры модифицированной взвешенной ортогонализации Грама–Шмидта и LD-разложения ковариационных матриц (алгоритм 2); результаты сравнительного анализа вычислительных свойств алгоритмов парной фильтрации. Полученные теоретические результаты дополняют теорию парной фильтрации в классе линейных дискретных парных марковских моделей с гауссовыми шумами.

Разработанный алгоритм реализован на языке MATLAB. Проведены серии численных экспериментов, результаты которых демонстрируют его работоспособность и численные преимущества перед другими существующими модификациями парного фильтра Калмана.

Представленные результаты могут быть в дальнейшем использованы для построения новых методов параметрической идентификации парных марковских моделей.

В работе представлен подробный и систематический анализ сходимости метода конечных элементов при моделировании распространения упругих волн в анизотропных средах, обладающих ортотропной симметрией тензора упругости. Пространственная дискретизация основана на применении метода конечных элементов с квадратичными базисными функциями $P_2$ на неструктурированных тетраэдральных сетках, что обеспечивает гибкость при работе со сложными геометриями. Для интегрирования по времени используется схема Ньюмарка с параметрами $\beta = 0,302$, $\gamma = 0,6$; при этом выбор именно этой схемы не оказывает влияния на пространственную сходимость, которая является предметом настоящего исследования.

Анализ сходимости выполнен на тестовой задаче с гладким локализованным начальным импульсом в полупространстве с граничными условиями свободной поверхности. Начальные условия подобраны таким образом, что решение остается гладким, а импульс не достигает границ расчетной области до момента сравнения решений. Это позволяет исключить влияние аппроксимации граничных условий на исследуемую сходимость и сосредоточиться исключительно на пространственной составляющей погрешности. Серия вычислительных расчетов проведена на шести последовательно сгущающихся неструктурированных сетках с числомтет раэдров от 39 728 до 372 904. Сравнение производилось с эталоннымрешением , полученным на существенно более мелкой сетке, что обеспечивает достаточную точность референсных данных. Для корректного вычисления нормы разности решения на всех сетках проецировались на единую грубую пространственную сетку.

Установлено, что порядок сходимости в норме $L_2$ составляет 2,314 при использовании медианной длины ребра тетраэдров в качестве характерного размера сетки и 2,878 при использовании обратного кубического корня из числа элементов. Полученные значения свидетельствуют о хорошей сходимости реализованной численной схемы и подтверждают корректность предложенного подхода.

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

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

Для автономных систем ОДУ рассмотрено простейшее подмножество двухстадийных схем Розенброка с комплексными коэффициентами, численная реализация которых требует одного LU-разложения и одного вычисления Якобиана за шаг интегрирования.

Проведено теоретическое исследование точности и устойчивости таких методов. Получены новые A-устойчивые методы 3-го порядка точности с различными свойствами и возможностью простой оценки главного терма локальной погрешности, что необходимо для автоматического выбора шага. Проведено тестирование новых методов.

Рассмотрен формализм квазиклассического приближения относительно малого коэффициента диффузии D, D→0, для многомерного уравнения Фоккера−Планка−Колмогорова с нелокальным и нелинейным вектором сноса в классе траекторно-сосредоточенных функций. Получена динамическая система Эйнштейна−Эренфеста типа (0, M), описывающая движение точки, в окрестности которой локализованы квазиклассические асимптотические решения. Построено семейство квазиклассических асимптотик с точностью O(D^(M+1)/2).

Настоящая работа посвящена периодическим и непериодическим полулокальным сглаживающим сплайнам или S-сплайнам класса C^p, состоящим из полиномов степени n.
Первые p + 1 коэффициентов каждого полинома задаются значениями предыдущего полинома и его p первых производных в точке склейки, остальные n − p коэффициентов при старших производных полинома определяются методом наименьших квадратов. Эти условия дополняются или начальными условиями (непериодический случай), или условием периодичности сплайн-функции на отрезке определения. В работе выписана система линейных уравнений, определяющих коэффициенты полиномов, составляющих сплайн. Матрица системы имеет блочный вид. Доказаны теоремы существования и единственности. Показано, что сходимость сплайнов к исходной функции зависит от величин собственных значений матрицы устойчивости. Приведены примеры устойчивых S-сплайнов.

Разработан алгоритм ускоренного моделирования КМОП СБИС (Сверх Больших Интегральных Схем с Комплементарной логикой на транзисторах Металл-Окисел-Проводник) под управлением точности. Алгоритм обеспечивает возможность проведения параллельного числительного эксперимента в много процессорной вычислительной среде. Ускорение расчета осуществляется за счет применения блочно-матричной и структурной (DCCC) декомпозиций. Особенность подхода состоит в выборе моментов и способов обмена параметрами и в применении многоскоростных методов интегрирования в процессе расчета подсистем. Благодаря этому имеется возможность оценивать и контролировать погрешность по требуемым характеристикам.

Данная работа посвящена сравнению двух критериев однородности: критерия χ², основанного на таблицах сопряженности признаков 2 × 2, и критерия однородности, основанного на асимптотических распределениях суммируемых квадратичных уклонений оценок функции распределения в модели зависимости «доза–эффект». Оценка мощности критериев производится при помощи компьютерного моделирования. Для построения функций эффективности используется метод ядерной оценки регрессии, основанный на оценке Надарая–Ватсона.

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