Об оценке устойчивости и погрешности разностной схемы в задаче оптимального управления для линейного уравнения Шредингера
Рассмотрение задачи оптимального управления для линейного уравнения Шредингера с критерием качества Лионса. Проведение оценки устойчивости и погрешности разностной схемы. Исследование скорости сходимости разностных аппроксимаций по функционалу.
Рубрика | Программирование, компьютеры и кибернетика |
Вид | статья |
Язык | русский |
Дата добавления | 26.04.2019 |
Размер файла | 347,6 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Н. М. Махмудов
Размещено на http://www.allbest.ru//
68
Размещено на http://www.allbest.ru//
Об оценке устойчивости и погрешности разностной схемы в задаче оптимального управления для линейного уравнения Шредингера
Н. М. Махмудов
Рассматривается задача оптимального управления для линейного уравнения Шредингера с критерием качества Лионса. К этой задаче применяется разностный метод. Устанавливается оценка устойчивости и погрешности разностной схемы, исследуется вопрос об оценке скорости сходимости разностных аппроксимаций по функционалу. Подобные вопросы освещены в работах [1-4] и др.
Постановка и дискретизация задачи© Н. М. Махмудов, 2010
оптимальный управление уравнение шредингер
Рассмотрим задачу о минимизации функционала
(1)
на множестве
при условиях
,(2)
, (3)
, (4)
,(5)
где - заданные числа, - ограниченная измеримая функция, удовлетворяющая условиям
,(6)
а функции , удовлетворяют условиям
, (7)
, (8)
- заданные числа.
При каждом заданном под решением редуцированной задачи (2)-(5) понимаются функции
,
,
удовлетворяющие соотношениям (2)-(5) для почти всех , .
Исходя из результатов работы [2] можем утверждать, что редуцированная задача (2)-(5) имеет единственное решение и верны оценки
, (9)
(10)
для , где и - постоянные величины, которые не зависят от , но зависят от остальных исходных данных задачи. Из условий (6)-(9) следует, что функции , , , являются гладкими, поэтому при с помощью результатов работы [2] можем определить следующие оценки:
,
(11)
(12)
для .
Произведя дискретизацию при каждом натуральном , рассмотрим задачу о минимизации функции
(13)
на множестве
при условиях
(14)
, (15)
, (16)
, (17)
где сеточные функции , определены следующими формулами:
, (18)
, , (19)
.(20)
Эти функции определены на сетке
, ,,
.
Обозначим
С помощью сумматорных тождеств нетрудно доказать справедливость утверждения.
Теорема 1. Для решения разностной схемы (14)-(17) при верна оценка
,(21)
для любого , где - постоянная величина, не зависящая от и .
Оценка погрешности разностной схемы
Оценим погрешность аппроксимации. С этой целью рассмотрим следующие усреднения решения редуцированной задачи (2)-(5) при :
,
.(22)
Определим оператор на множестве по формуле
. (23)
Обозначим . Ясно, что будет решением следующей системы:
(24)
, (25)
, (26)
, (27)
Где
. (28)
Теорема 2. Пусть выполнено условие согласования , где , - постоянные, не зависящие от и . Тогда верны оценки
,
, (29)
для , где
. (30)
Доказательство. Аналогично оценке (21) легко доказать справедливость оценки
, (31)
.
Оценим правую часть неравенства (31). Используя формулу для , эту оценку можно представить в виде
,(32)
Где
, (33)
,
, (34)
,
. (35)
Используя формулы (34),(35), а также (20),(21) для слагаемых ,, можно получить справедливость неравенства
, (36)
. (37)
Оценим слагаемое . Предварительно оценим значения при , . С помощью формулы (33) находим
.(38)
Используя формулу для , нетрудно получить соотношение
,(39)
Оценим значения
.
С помощью формулы для получим соотношение
(40)
Пусть справедливо равенство
(41)
для фиксированного . Тогда согласно формуле (40) становится справедливым соотношение
. (42)
Обозначим через , а - через и сделаем замену переменных в виде . Тогда функционал примет вид
, (43)
где . Очевидно, что есть линейный функционал относительно . Кроме того, этот функционал ограничен в пространстве . Действительно, в силу оценки (11) для фиксированного функция принадлежит пространству . Поэтому будет принадлежать пространству для фиксированного . Используя (43), получим неравенство
, (44)
где - постоянная, не зависящая от и . Кроме того, функционал обращается в нуль на многочленах второй степени. Действительно, если возьмем и подставим его в формулу для , то получим
(45)
Таким образом, выполнены все условия леммы Брэмбла-Гильберта [1, с.29]. Тогда в силу этой леммы и оценки (44) получим неравенство
. (46)
При выполнении обратной замены переменных и возврате к прежним обозначениям согласно (46) получим
,
. (47)
Рассмотрим соотношения для , при . Используя формулу, определяющую величины , получим
.(48)
Оценим значения и . Легко видеть справедливость неравенств
(49)
.(50)
Очевидно неравенство
(51)
.
С учетом (51) и (50) получим
(52)
Используя лемму Брэмбла-Гилберта, аналогично неравенству (47) для получим
,
. (53)
В силу справедливости соотношений (49), (52) и (53) и (48) получим
,
. (54)
Аналогично оценке слагаемого оценим величины :
, . (55)
Аналогично оценкам (54),(55) и в силу граничных условий для можно записать соотношения
, . (56)
, . (57)
Используя (39), (47) для при , получим
+
. (58)
По аналогичной схеме можно получить следующую оценку:
.(59)
В силу справедливости неравенств (36), (37), (52), (57), оценок (11), (12) и условия согласования, используя (32), получим соотношение
,
. (60)
Из последних неравенств и оценки (31) находим
,
. (61)
Выбирая , получим формулу, доказывающую теорему.
Оценка скорости сходимости разностных аппроксимаций по функционалу
Оценим разность исходного функционала и дискретной функции.
Теорема 3. Пусть выполнены условия теоремы 2. Тогда для любого и имеет место оценка
, , (62)
Доказательство. Используя формулу (1) и (13), получим
,(63)
Согласно виду формулы, определяющей , справедливо неравенство
. (64)
В силу утверждения теоремы 2 получаем
. (65)
Рассмотрим разность . С учетом соотношения (22) находим
.(66)
С учетом (66) получим
(67)
Согласно (11) и (67)
. (68)
Суммируя (65) и (66), определяем
. (69)
Аналогично получим следующее соотношение:
. (70)
Неравенства (63) и (70) доказывают теорему.
Сформулируем две вспомогательные леммы.
Лемма 1. Пусть выполнены условия теоремы 3. Пусть, кроме того, оператор определяется формулой (23). Тогда и имеет место оценка
, . (71)
Доказательство этой леммы проводится с использованием утверждения теоремы 3.
Пусть оператор определяется формулой , (72)
где
(73)
Лемма 2. Пусть выполнены условия теоремы 3. Пусть, кроме того, оператор определяется формулой (72). Тогда и имеет место оценка
,
. (74)
Доказательство. Нетрудно установить, что .
Поэтому в записи теоремы 3, выбирая вместо и проводя аналогичное доказательство, получим оценку
,
. (75)
Легко видеть, что
Отсюда получим
.
Это неравенство и соотношение (75) доказывают лемму.
Сформулируем и докажем теорему о скорости сходимости разностных аппроксимаций по функционалу.
Теорема 4. Пусть выполнены условия леммы 1 и 2. Пусть, кроме того, и являются решениями задач (1)-(5) и (13)-(17) соответственно, т. е.
.
Тогда последовательность разностных задач (13)-(17) аппроксимирует задачу (1)-(5), т.е.
(76)
и справедлива оценка о скорости сходимости:
, .(77)
Доказательство этой теоремы проводится с использованием леммы 1 и 2.
Список литературы
Потапов М.М., Разгулин А.В., Шамеева Т.Ю. Аппроксимация и регуляризация задачи оптимального управления типа Шредингера // Вестн. Моск. ун-та. Сер.15. Вычисл. матем. и киберн. 1987. №1. С.8-13.
Ягубов Г.Я. Оптимальное управление коэффициентом квазилинейного уравнения Шредингера: дис. доктор. наук. Киев, 1994. 318 с.
Ягубов Г.Я. Разностный метод решения задачи оптимального управления коэффициентом квазилинейного уравнения Шредингера с интегральным критерием качества по границе области // Проблемы матем. модел. и опт. управления. Баку, 2001. С.37-48.
Ягубов Г.Я. Сходимость разностного метода решения задачи оптимального управления для нелинейного уравнения Шредингера с интегральным критерием качества // Вестн. Сумгаит. гос. ун-та. 2001. №1. С.37-42.
Размещено на Allbest.ru
...Подобные документы
Понятие разностных схем, сеточная функция, пространство и нормы. Аппроксимация дифференциальных операторов. Корректность разностной схемы и сходимость. Одномерное уравнение переноса с переменными и постоянными коэффициентами. Схема бегущего счета.
дипломная работа [388,3 K], добавлен 11.11.2009Математическое описание алгоритмов схемы и операций для уравнения Лапласа. Изучение разностной схемы "крест" для нахождения численного решения эллиптического уравнения, задача Дирихле. Использование указателей в среде Matlab для решений методом Гаусса.
дипломная работа [859,3 K], добавлен 23.10.2014Решение конечно-разностной задачи Дирихле для уравнения Лапласа в прямоугольной области. Погрешность замены дифференциального уравнения разностным. Использование схемы узлов при получении сеточных уравнений. Сущность метода Зайделя. Листинг программы.
курсовая работа [348,5 K], добавлен 26.04.2011Определение оптимального плана выпуска продукции частного предприятия по изготовлению мебели с применением метода линейного программирования (симплекс-метод). Построение схемы движения информации в подсистеме оптимального плана выпуска продукции.
лабораторная работа [301,5 K], добавлен 08.06.2009Теоретическая основа линейного программирования. Задачи линейного программирования, методы решения. Анализ оптимального решения. Решение одноиндексной задачи линейного программирования. Постановка задачи и ввод данных. Построение модели и этапы решения.
курсовая работа [132,0 K], добавлен 09.12.2008Алгоритм решения задач линейного программирования симплекс-методом. Построение математической модели задачи линейного программирования. Решение задачи линейного программирования в Excel. Нахождение прибыли и оптимального плана выпуска продукции.
курсовая работа [1,1 M], добавлен 21.03.2012Методы определения оптимального плана производства (приобретения) продукции с учетом ограниченного обеспечения ресурсами различного вида. Технология поиска оптимального решения задач линейного программирования (ЗЛП) с помощью итоговой симплекс-таблицы.
лабораторная работа [42,8 K], добавлен 11.03.2011Тестирование модуля отыскания корня уравнения методом половинного деления. Схема алгоритма тестирующей программы. Численное интегрирование по методу Симпсона с оценкой погрешности по правилу Рунге. Проверка условий сходимости методов с помощью MathCAD.
курсовая работа [1,1 M], добавлен 04.02.2011Анализ метода линейного программирования для решения оптимизационных управленческих задач. Графический метод решения задачи линейного программирования. Проверка оптимального решения в среде MS Excel с использованием программной надстройки "Поиск решения".
курсовая работа [2,2 M], добавлен 29.05.2015Решение задачи линейного программирования симплекс-методом: постановка задачи, построение экономико-математической модели. Решение транспортной задачи методом потенциалов: построение исходного опорного плана, определение его оптимального значения.
контрольная работа [118,5 K], добавлен 11.04.2012Разработка структурной схемы и обобщенного алгоритма работы прибора. Оценка максимальной погрешности линейного датчика давления и нормирующего усилителя. Разработка элементов принципиальной электрической схемы микропроцессорной системы сбора данных.
курсовая работа [2,1 M], добавлен 08.02.2015Анализ решения задачи линейного программирования. Симплексный метод с использованием симплекс-таблиц. Моделирование и решение задач ЛП на ЭВМ. Экономическая интерпретация оптимального решения задачи. Математическая формулировка транспортной задачи.
контрольная работа [196,1 K], добавлен 15.01.2009Разработка цифрового нерекурсивного и рекурсивного фильтров с заданными параметрами. Проектирование фильтра в программе Matlab с помощью утилиты fdatool. Построение структурной схемы во вкладке Realize model. Общий вид линейного разностного уравнения.
курсовая работа [2,9 M], добавлен 19.03.2012Постановка задачи линейного программирования и формы ее записи. Понятие и методика нахождения оптимального решения. Порядок приведения задач к каноническому виду. Механизмы решения задач линейного программирования аналитическим и графическим способами.
методичка [366,8 K], добавлен 16.01.2010Математическое программирование. Линейное программирование. Задачи линейного программирования. Графический метод решения задачи линейного программирования. Экономическая постановка задачи линейного программирования. Построение математической модели.
курсовая работа [581,5 K], добавлен 13.10.2008Описание подхода по использованию методов оптимального управления для задачи следящих систем. Сопровождающая линейно-квадратичная задача оптимального управления. Свойства и алгоритм построения оптимальной стартовой обратной связи и дискретного управления.
дипломная работа [871,4 K], добавлен 20.08.2013Рассмотрение модели механизма перемотки ленты в форме структурной схемы и передаточной функции. Определение характеристического уравнения и коэффициентов обратных связей. Проверка результатов синтеза моделированием в программном пакете Classic 3.01.
курсовая работа [1,1 M], добавлен 18.08.2013Строение системы уравнений-ограничений и ее переменных, графический способ решения задач линейного программирования на плоскости. Выражение неизвестных через две независимые переменные, являющиеся координатными осями графика. Значение целевой функции.
лабораторная работа [61,4 K], добавлен 07.01.2011Постановка задачи линейного программирования. Решение системы уравнений симплекс-методом. Разработка программы для использования симплекс-метода. Блок-схемы основных алгоритмов. Создание интерфейса, инструкция пользователя по применению программы.
курсовая работа [1,7 M], добавлен 05.01.2015Программа вычисления интеграла методом прямоугольников. Решение задачи Коши для дифференциальных уравнений. Модифицированный метод Эйлера. Методы решения краевой задачи для обыкновенного дифференциального уравнения. Задачи линейного программирования.
методичка [85,2 K], добавлен 18.12.2014