Приближенные решения задач математической физики

Автор работы: Пользователь скрыл имя, 07 Октября 2009 в 14:32, Не определен

Описание работы

Приближенные решения задач математической физики

Файлы: 1 файл

Вариационные методы - текст2008.doc

— 1.04 Мб (Скачать файл)

Если  оператор положительно определенный, то он положительный, обратное не всегда верно. 

Пример  4:  , .

 – положительный оператор. Докажем его положительную определенность.

      (с учетом  )

Неравенство Буняковского: . Примем , :

 
.

  Энергетическое пространство положительно определенного оператора

Оператор  - положительно определенный оператор, – линеал (область определения линейного оператора). Введем энергетическое пространство оператора (полное гильбертово пространство, совпадающее с ) с обозначениями:

.    (1) - энергетическое произведение

,   (4) - энергетическая норма .

.    (5) 

Теорема 1.

Все элементы пространства принадлежат также к пространству .

(Точнее: каждому элементу из  можно привести в соответствие один и только один элемент из , причем разным элементам из соответствуют разные элементы из .)

Сходимость  в энергетическом пространстве –  сходимость по энергии ( ).

Теорема 2.

Если  – положительно определенный оператор и по энергии, то одновременно в метрике исходного пространства :

 

  Энергетическое пространство только положительного оператора

Оператор  положительный, но не положительно определенный, называется только положительным. Не все элементы энергетического пространства только положительного оператора  принадлежат исходному пространству.

Элемент принадлежит исходному пространству тогда и только тогда, когда , что и 0. При этом последовательность стремится к тому же элементу в пространстве :

.

Теорема.

Для того, чтобы энергетическое пространство положительного оператора было сепарабельным, необходимо и достаточно, чтобы было сепарабельным исходное гильбертово пространство. 

(Сепарабельное пространство - плотное, счетное, – сепарабельное пространство).

Итак:

Положительный оператор

Положительно определенный оператор , что

  Главные и естественные краевые условия

, , ,    (1)

,    (2)

 – оператор задачи (1), (2) в  пространстве 

 – множество функций из  и удовлетворяет (2) ( –порядок уравнения (1)).

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

  1. , , , – главные (геометрические, кинематические)
  2. , , , – естественные (динамические)
 

Пример  5

  (3)

,    (4)

Докажем, что для уравнения (3) краевое условие (4) – естественное.

 – симметрично относительно , .

 – имеет смысл для любых , необязательно удовлетворяющих (3).

Построим  функционал .

Покажем, что точное решение задачи (3), (4) – реализует . Используем принцип виртуальных перемещений с параметром .

 – имеет  при и фиксированном .

.

При      (5)

По формуле  Грина

, для

,
  
,

Метод Ритца.

  Пусть – положительно определенный оператор на линеале в сепарабельном пространстве , и .  Пусть – гильбертово пространство. Рассмотрим в базис

  

.  (1)

  Обобщенное  решение уравнения  – это элемент , который минимизирует в функционал

  

,  (2)

т.е. элемент , для которого

       .   (3)

  Выберем целое положительное число , и будем искать аппроксимацию элемента в виде

  

,  (4)

где элементы базиса (1) , а – неизвестные пока вещественные постоянные. Эти постоянные определяются из условия

  Fun = min,  (5)

которое означает, что среди всех аппроксимаций  вида

  

 (6)

где bk – произвольные вещественные постоянные (т.е. в n-мерном подпространстве, порождаемом элементами j1, …, jk,), функционал F принимает минимальное значение в точности на аппроксимации (4). По предположению, (1) образует базис в HA , так что обобщенное решение u0 можно с произвольной точностью аппроксимировать соответствующей линейной комбинацией его элементов. Кроме того, условие (5) аналогично (3). Поэтому приближение (4) с постоянными, определенными в соответствии с (5), будет достаточно мало отличаться от искомого решения u0 в HA, если n будет достаточно велико.

  Постоянные  bk в (4) определяются подстановкой (6) вместо u в (2):

  F vn = (b1 j1 +…+ bnjn, b1j1 +…+ bnjn)A2(f, b1j1 +…+ bnjn)= 

  =(j1, j1)A b12 + (j1, j2)A b1 b2 +…+ (j1, jn)A b1 bn +  

  +(j2, j1)A b2 b1 + (j2, j2)A b22 +…+ (j2, jn)A b2 bn +…  

  +(jn, j1)A bn b1 + (jn, j2)A bn b2 +…+ (jn, jn)A bn2–  

  –2(f, j1) b1 2(f, j2) b2 …– 2(f, jn) bn ,  (7)

  а с учетом симметрии скалярного произведения:

  F vn =(j1, j1)A b12 + 2(j1, j2)A b1 b2 +…+ 2(j1, jn)A b1 bn +   

  + (j2, j2)A b22 +…+ 2(j2, jn)A b2 bn +…+ (jn, jn)A bn2–  

  – 2(f, j1)b1 2(f, j2)b2 …– 2(f, jn)bn .   (8)

  Скалярные произведения (ji,,jk)A – это фиксированные числа, определяемые элементами заданного базиса. Следовательно, в результате подстановки функционал F становится квадратичной функцией переменных b1,…,bn. Необходимые условия существования минимума этой функции в точке (a1,…,an) есть:

  

  (9)

т.е. должны выполняться равенства

 (10)

После простых преобразований эти уравнения можно записать в виде

 (11)

  Система (11) – это система n уравнений для n искомых постоянных a1,…,an. Так как по предположению j1,…,jn линейно независимы и определитель системы (11), являющийся определителем Грама, отличен от нуля, то существует единственное решение системы(11).

  Замечание. Базис в HA можно выбрать из элементов линеала DA, так как это множество плотно в HA. Тогда мы имеем (ji, jk)A = (Aji, jk) для всех i, k=1,…,n, и систему (11) можно записать в виде

  

 (12)

Применение  собственных элементов  сходного оператора.

Понятие о невязке

Если дано уравнение 

  (1)

и - какое-либо его приближенное решение, то разность называется невязкой этого приближенного решения. Пусть - положительно определенный оператор, действующий в гильбертовом пространстве , и - приближенное решение уравнения (1) по Ритцу. Будем считать, что координатные элементы принадлежат области определения оператора , тогда и выражение - невязка – имеет смысл. В случае ограниченного оператора невязка стремится к нулю:

Информация о работе Приближенные решения задач математической физики