Математическое моделирование. Рубрика в журнале - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование
Статья научная
Рассматривается задача маршрутизации перемещений с ограничениями и функциями стоимости, допускающими зависмость от списка заданий. Предполагается, что начальное условие процесса с дискретным временем может выбираться в пределах метрического пространства, удовлетворяющего условию полной ограниченности. По постановке задачи предполагается посещение конечной системы мегаполисов (непустых конечных множеств) с выполнением тех или иных работ, стоимости которых зависят всякий раз от пункта прибытия и пункта отправления. Стоимости перемещений и выполняемых работ агрегируются аддитивно. Для решения используется вариант широко понимаемого динамического программирования, обеспечивающий нахождение ε-оптимального решения при любом значении ε>0.
Бесплатно
Ослабление сферических ударных волн в гетерогенных средах
Статья научная
В данной работе приводится анализ численного исследования взаимодействия сферических ударных волн (УВ) с механическими и тепловыми неоднородностями, расположенными на плоской поверхности. Результаты расчетов показали, что наличие на плоской поверхности гетерогенного слоя приводит к ослаблению падающей УВ внутри этого слоя в зависимости от объемной доли конденсированной фазы. Достаточно большая объемная доля твердой фазы в гетерогенном слое приводит к деформации головного фронта УВ и сильное ее торможение, в результате чего отражения от стенки практически не происходит и ножка Маха не образуется, однако при a2>=5*10-2 - она образуется над поверхностью гетерогенного слоя. Увеличение температуры газа в слое приводит к возрастанию скорости (УВ), в то время как амплитуда ее снижается. Приведенные в работе расчеты показали, что увеличение температуры газа в слое способствует увеличению момента силы, действующей на решетку из твердых частиц. При высоких температурах газа в гетерогенном слое, несмотря на силовое и тепловое взаимодействие газа и частиц, происходит образование предвестника.
Бесплатно
Особенности математического моделирования процесса гидродинамического исследования нефтяных пластов
Статья научная
Оценка запасов нефтяного месторождения проводится специалистами геологических служб на основании гидродинамических исследований пластов. Существующие методики направлены на определение гидропроводности нефтеносного пласта и продуктивности скважин по данным их кратковременной эксплуатации. Для обработки результатов измерений используют различные методы, в частности, те, которые основаны на численном решении прямых и обратных задач фильтрации. При решении задачи нахождения коэффициента гидропроводности численными методами необходимо учитывать особенности задач подземной гидромеханики. Эти особенности нужно учитывать при составлении математической модели рассматриваемого процесса и при разработке алгоритмов ее численного решения. Ряд условий позволяет сформулировать задачу определения коэффициента гидропроводности как обратную нелинейную задачу гидродинамики. Существенно важным для решения данной задачи является доказательство единственности ее решения. В данной работе формулируются условия для обратной задачи фильтрации со смешанными граничными условиями, гарантирующие единственность ее решения.
Бесплатно
Статья научная
Рассматривается задача определения систематических ошибок нескольких РЛС по реальным данным измерений движущихся объектов (самолетов). В случае, когда модель пространственной зависимости систематических ошибок не известна полностью, их нахождение сводится к некорректно поставленной задаче оценивания. В работе предлагается подход, позволяющий в этих условиях получить разумную оценку. Его основу составляет локальная аппроксимация неизвестных систематических ошибок, рассматриваемых как функции геометрического положения. Пространство положений разбивается на систему достаточно малых областей. Внутри каждой области ищется вектор, локально приближающий сдвиг от систематических ошибок. Из-за некорректности задачи может быть определено только множество неопределенности, которое содержит все векторы сдвига, которые могли бы дать одинаковые измерения. Набор построенных множеств неопределенности можно рассматривать как многозначную функцию геметрического положения. Далее производится выборка однозначной функции систематических ошибок из многозначной на основе критерия, минимизация которого позволяет выделить наиболее плавно изменяющуюся функцию. Алгоритм опробован на реальных данных траекторного наблюдения.
Бесплатно
Оценка погрешности численного метода решения одной обратной задачи
Статья научная
Рассмотрен линейный дифференциальный оператор и система краевых условий, задаваемая линенйыми в пространстве n раз непрерывно дифференцируемых функций линейно-независимыми функционалами. Функция Грина для краевой задачи, определенной этим оператором и упомянутыми функционалами, строится как решение интегрального уравнения Фредгольма II рода, параметры которого определяются функцией Грина вспомогательной задачи. Предложенный метод обращения дает возможность эффективно решить как прямую (т.е. задачу нахождения решения), так и обратную (т.е. задачу нахождения правой части уравнения по экспериментально полученному решению) задачи. Обсуждены особенности численной реализации метода и возможности оценки точности полученных решений.
Бесплатно
Статья научная
Настоящая работа посвящена изучению одного класса систем дифференциальных уравнений нейтрального типа. Указаны области притяжения нулевого решения и установлены оценки экспоненциального убывания решений на бесконечности. В частности, из этих оценок вытекает асимптотическая устойчивость нулевого решения рассматриваемых систем. Результаты получены с использованием модифицированного функционала Ляпунова - Красовского.
Бесплатно
Статья научная
Построенные в ВЦ РАН динамические модели региональной экономики содержат большое число неизвестных параметров, идентифицировать которые можно с помощью минимизации нелинейных критериев близости расчетных и статистических данных. В работе для этого применяется параллельный индексный метод глобальной оптимизации, разработанный в ННГУ. Метод использует редукцию размерности на основе кривых Пеано и информационно-статистический подход, дополненный различными модификациями.
Бесплатно
Параметризация алгоритмов идентификации электрического диполя
Статья научная
Рассматривается задача идентификации параметров положения произвольно-ориентированного электрического диполя над плоскостью с бесконечной проводимостью по его электромагнитному полю, индуцируемому в точке наблюдения. Данная задача входит в комплекс математических моделей практически важной проблемы прогнозирования развития грозовых очагов. Для решения задачи, из-за ее плохой обусловленности, предлагается строить множество оценок с использованием параметризованного семейства алгоритмов, а окончательное решение принимать по результатам статистического анализа. Среди возможных методов определения параметров положения произвольно ориентированного электрического диполя, в работе рассмотрена параметризация прямого метода и экстремального параметрического метода. Однако, для повышения статистической значимости результирующей оценки параметров следует использовать по возможности большее число алгоритмов. Приведенные в работе результаты вычислительного эксперимента подтверждает эффективность подхода, однако большое разнообразие результатов измерения не позволяет строить достаточно точные гарантированные результирующие оценки параметров положения. Для повышения качества результирующих оценок, предлагается фильтрация построенного множества оценок, основанная на применении фильтрующих функционалов и их комбинации. Для построения фильтрующих функционалов использованы метод антитрасс и метод проектирования.
Бесплатно
Построение интегральной модели на примере динамики ветроэнергетической установки
Статья научная
Данная работа посвящена применению интегро-степенных рядов Вольтерра к описанию нелинейных динамических систем типа вход - выход. Универсальность данного математического аппарата позволяет создавать программное обеспечение для проведения экспериментов на компьютере. Работа продолжает исследования в области идентификации ядер Вольтерра, начатые в Институте систем энергетики им. Л.А. Мелентьева СО РАН. В первой части статьи рассматривается новый алгоритм идентификации полиномов Вольтерра второй степени для систем, допускающих активный эксперимент с помощью тестовых наборов возмущающих воздействий. Во второй части работы приведены результаты численных расчетов для эталонной динамической системы. Эталоном послужила модель ветроэнергетической установки с горизонтальной осью вращения. Построены квадратичные полиномы Вольтерра, описывающие нелинейную динамику угловой скорости вращения элементов ветроустановки от угла наклона лопастей и скорости ветра. Практическая идентификация ядер Вольтерра производилась относительно некоторого выбранного стационарного состояния моделируемой системы.
Бесплатно
Статья научная
Рассматривается алгоритм построения квазиструктурированных сеток, которые состоят из равномерных прямоугольных подсеток и строятся в два этапа. На первом из них расчетная область покрывается равномерной прямоугольной макросеткой, а на втором в каждом макроэлементе задается своя прямоугольная равномерная подсетка. Существенным является то, что подсетки могут быть несогласованными. За счет регулировки плотности узлов подсеток достигается адаптация квазиструктурированной сетки к неоднородностям внутри области. Для адаптации сетки к внешней границе подсетки подвергаются локальной модификации, состоящей в сдвиге приграничных узлов на границу. Излагается алгоритм локальной модификации для построения качественной квазиструктурированной сетки, который не нарушает структурированности подсеток. Предлагаемые квазиструктурированные сетки выгодно отличаются от структурированных сеток тем, что не требуют введения лишних узлов, которые необходимы лишь для поддержки структурированности, и от неструктурированных сеток тем, что не требуют хранения большого объема информации. Решение краевых задач на квазиструктурированных сетках ищется методом декомпозиции расчетной области на подобласти, сопрягаемые без наложения. Данный метод легко распараллеливается и поэтому может быть применим для проведения расчетов на многопроцессорных суперЭВМ.
Бесплатно
Построение наблюдения для задачи оптимального динамического измерения по искаженным данным
Статья научная
Теория оптимальных динамических измерений базируется на минимизации разности значений виртуального наблюдения, т.е. полученного с помощью расчетной модели, и экспериментальных данных, которые обычно искажены некоторыми помехами. В статье приведено описание математической модели оптимального динамического измерения при наличии помех разного вида. Кроме того, в статье предлагается алгоритм построения значений наблюдения по значениям, полученным в ходе эксперимента, которые предполагаются искаженными некоторыми случайными воздействиями. Предполагается, что на экспериментальные данные воздействует "белый шум", который понимается как производная Нельсона - Гликлиха от винеровского процесса. Для построения значений наблюдения используется априорная информация о форме функции, описывающей значения наблюдения. Сама процедура построения наблюдения состоит из двух этапов. На первом этапе формулируется критерий определения положения экстремальной точки сигнала с использованием статистики специального вида. А на втором этапе описывается процедура построения значений сигнала на основе информации о положении точки экстремума и форме выпуклости сигнала.
Бесплатно
Представление суммы Минковского для двух полиэдров системой линейных неравенств
Статья научная
Любой выпуклый полиэдр представим как множество решений некоторой системы линейных неравенств. Алгебраическая сумма по Минковскому выпуклых полиэдров X, Y С R n также является выпуклым полиэдром, и, следовательно, также представим как множество решений некоторой системы линейных неравенств. В статье предложен полиномиальный алгоритм решения указанной задачи, основанный на формировании ряда избыточных ограничений в представлении слагаемых и их трансляции в результирующее представление. Предложен эффективный способ использования параллельных и распределенных вычислений для реализации алгоритма.
Бесплатно
Применение дробно-рациональных интерполяций для решения краевых задач с особенностями
Статья научная
Статья посвящена разработке, реализации и тестированию нового метода решения сингулярно-возмущенных краевых задач для нелинейных уравнений с частными производными второго порядка в прямоугольной области. Для приближения решения в методе использованы прямые (тензорные) произведения дробно-рациональных функций, полученных из интерполяционных полиномов с узлами Чебышева, записанных в барицентрической форме, с помощью специальной замены переменной. Замена делается с целью адаптации положения узлов интерполяции к особенностям искомой функции и приводит к их сгущению в окрестности больших градиентов решения. Для аппроксимации нелинейных уравнений используется сочетание итерационного метода установления и метода коллокаций, что позволяет свести задачу на каждой итерации к решению матричного уравнения Сильвестра. Такой подход приводит к существенному снижению времени вычислений. Высокая эффективность метода продемонстрирована на примере тестовой краевой задачи в квадрате, решение которой имеет пик в центре области, обусловленный наличием у неизвестной функции полюса в комплексной плоскости.
Бесплатно
Применение нечеткой логики для создания системы управления устойчивостью двухколесного мотоцикла
Статья научная
Рассмотрена возможность применения методов управления на основе нечеткой логики для создания системы управления устойчивостью двухколесного транспортного средства на примере мотоцикла. Реализован программный комплекс эмулирующий физическое поведение мотоцикла, управляющий рулем и отображающий графически процесс управления в реальном времени.
Бесплатно
Разрешимость краевой задачи для вырождающихся уравнений соболевского типа
Статья научная
Целью работы является доказательство существования и единственности регулярных решений первой краевой задачи для систем уравнений соболевского типа с эллиптико-параболическими операторами с пространственным вырождением. А.И. Кожановым были рассмотрены начально-краевые задачи для уравнений соболевского типа с эллиптико-параболическими операторами второго порядка, действующими по пространственным переменным. В его работах были доказаны существование решения при выполнении условий «характеристической выпуклости» границы области относительно пространственных операторов. Техника, используемая в настоящей работе, будет близка к технике работ вышеуказанного автора. Для исследования вырождающихся систем уравнений соболевского типа используется также сочетание метода регуляризации и метода априорных оценок. С помощью метода регуляризации строится семейство приближенных решений вырождающихся уравнений. Анализ интегральных неравенств, при получении априорных оценок, основан на интегрировании по частям, применении неравенств Коши - Буняковского и Гельдера и неравенства Юнга. Также применяются свойства весовых соболевских пространств.
Бесплатно
Разрешимость нестационарной задачи теории фильтрации
Статья научная
Рассмотрена одна задача для класса неклассических уравнений математической теории волн. Отличительной особенностью этой задачи является зависимость от времени функциональных коэффициентов эллиптического оператора в правой части уравнения. Методом ее исследования является редукция к задаче Коши для нестационарного уравнения соболевского типа. Уравнения соболевского типа с зависящим от времени оператором в данной постановке рассматриваются впервые. Введено в рассмотрение понятие относительно спектрально ограниченной оператор-функции. Условия, гарантирующие выполнение этого свойства задачи, позволяют также выделить подпространство начальных значений, для которых существует единственное решение задачи Коши. Это подпространство мы назвали обобщенным фазовым пространством решений для нестационарного уравнения соболевского типа. Решение такой задачи для уравнений соболевского типа, а также и в исходной постановке, получено с помощью рекурсивной формулы.
Бесплатно
Расчет динамики баллистической модели ракет
Статья научная
Для обеспечения безопасности испытаний баллистических моделей ракет в гидродинамических бассейнах нашли широкое применение гидравлические улавливающие устройства в виде заполненной водой трубы с глухим днищем. Для ликвидации явления гидроудара в торце трубы предусматривается воздушный колокол. Разработанные математическая модель и метод расчета динамики баллистической модели в гидравлическом улавливающем устройстве позволяют выбрать геометрические параметры улавливающего устройства и проводить торможение модели в заранее заданном расчетном режиме. Уравнение для продольного движения модели получено из уравнения Лагранжа. По предлагаемому методу создана программа на ЭВМ и выполнены примеры расчетов. Проведено сравнение расчетов с экспериментальными данными, полученными при испытаниях в гидробассейне. Наблюдается удовлетворительное соответствие расчетных и экспериментальных данных, что служит подтверждением достоверности и надежности разработанной математической модели тормозного гидродинамического устройства. Разработанная математическая модель позволяет при заданном числе Эйлера и для заданной массы модели выбирать необходимые для торможения основные параметры улавливающего тормозного устройства. Предложенный метод расчета может быть использован для определения геометрических параметров тормозного и улавливающего устройства при проведении испытаний баллистической модели в гидробассейне.
Бесплатно
Реконструкция распределенных управлений в гиперболических системах динамическим методом
Статья научная
Рассматривается задача о восстановлении априори неизвестных распределенных управлений в гиперболических динамических системах по результатам приближенных измерений состояний (скоростей) наблюдаемого движения системы. Задача решается в динамическом варианте, когда для определения текущего приближения неизвестного управления разрешено использовать только поступившие в данный момент приближенные измерения, реконструкция управления должна осуществляться в динамике (по ходу процесса, по ходу движения системы). Рассматриваемая задача некорректна. Для решения задачи предлагается воспользоваться методом динамической регуляризации, разработанным Ю.С. Осиповым и его школой. Построены новые модификации динамических регуляризующих алгоритмов решения задачи, которые в отличие от традиционных алгоритмов позволяют получить усиленную сходимость регуляризованных приближений, в частности получить кусочно-равномерную сходимость. Выполнена конечномерная аппроксимация задачи. Приводятся результаты численного моделирования, которые позволяют судить о способности модифицированных алгоритмов восстанавливать тонкую структуру искомых управлений.
Бесплатно
Решение обратных задач достижения сверхразрешения при использовании нейронных сетей
Статья научная
Рассматривается актуальная проблема получения приближенных численных решений обратных задач в виде интегральных уравнений Фредгольма первого рода для систем радио- и гидролокации и дистанционного зондирования. Полученные решения дают возможность существенно повысить точность измерений, а также довести угловую разрешающую способность до значений, превышающих критерий Рэлея. Это позволяет: получать детализированные радиоизображения различных объектов и зондируемых областей; определять количество отдельных малоразмерных объектов в составе сложных целей, которые раздельно не фиксировались без представляемой обработки сигналов; получать координаты таких малоразмерных объектов с высокой точностью; повысить вероятности получения правильных решений задач распознавания и идентификации объектов. Метод применим для современных многоэлементных измерительных систем. Он основан на экстраполяции сигналов, принимаемых всеми элементами, за пределы самой системы. Решена задача создания необходимой для этого нейронной сети и ее обучения. В итоге, синтезируется новая виртуальная измерительная система значительно большего размера, что позволяет резко повысить угловое разрешение и тем самым повысить качество приближенных решений рассматриваемых обратных задач. На примерах демонстрируется эффективность метода, оценивается адекватность и устойчивость получаемых решений. Исследуется степень превышения виртуальной угломерной системой критерия Рэлея в зависимости от отношения сигнал/шум.
Бесплатно
Статья научная
В работах авторов были найдены линейные формулы, позволяющие находить приближенные собственные значения дискретных полуограниченных операторов. Используя их можно находить собственные значения дискретных операторов с любым порядковым номером. При этом снимаются многие вычислительные проблемы, возникающие в классических методах связанные с порядковым номером вычисляемых собственных значений и вопросов корректности производимых операций при их нахождении. Сравнение полученных результатов вычислительных экспериментов показали, что собственные значения, найденные по линейным формулам и методом Галеркина, хорошо согласуются. Причем, по мере увеличения порядкового номера собственных значений отличия уменьшаются. Используя линейные формулы, позволяющие вычислять собственные значений дискретных полуограниченных операторов, в статье изложен метод решения обратных спектральных задачах для операторов Штурма - Лиувилля, заданных на последовательных геометрических графах с конечным числом звеньев. Алгоритм апробирован на последовательном двухреберном графе. Результаты многочисленных экспериментов показали хорошую точность и высокую вычислительную эффективность разработанного метода.
Бесплатно
