Лекции.Орг


Поиск:




Категории:

Астрономия
Биология
География
Другие языки
Интернет
Информатика
История
Культура
Литература
Логика
Математика
Медицина
Механика
Охрана труда
Педагогика
Политика
Право
Психология
Религия
Риторика
Социология
Спорт
Строительство
Технология
Транспорт
Физика
Философия
Финансы
Химия
Экология
Экономика
Электроника

 

 

 

 


Метод последовательных приближений




Задано линейное интегральное уравнение Фредгольма второго рода

, (6)

где y(t), K(t,s) – заданные функции, x(t) – искомая функция, λ – числовой параметр.

Запишем уравнение (6) в виде

, (7)

где отображение , X – банахово пространство. Предположим, что параметр λ, ядро K(t,s) таковы, что отображение F является сжимающим с коэффициентом сжатия λ. Тогда уравнение (6) имеет единственное решение, которое можно найти методом последовательных приближений с помощью следующих рекуррентных соотношений:

(8)

Приближения (8) сходятся к решению уравнения (6) с некоторой скоростью, причем величина погрешности n -го приближения определяется неравенством

, (9)

где - нулевое приближение, которое выбирается произвольным образом.

 

Задание 1

Задание: Определить, при каких для следующих интегральных уравнений Фредгольма 2-го рода в пространстве C[0, 1], [0, 1] можно применить метод сжимающих отображений. При найти приближенное решение методом последовательных приближений с точностью , сравнить его с точным решением: +

Решение:

Рассмотрим пространство C[0, 1]. Определим , при которых применим метод сжимающих отображений для данного уравнения.

F(x) = + ;

= = = .

≤ 1 => .

Т.е. при применим метод сжимающих отображений.

Положим , тогда . Оценим нужное количество итераций с помощью неравенства (положив = 0, тогда = , и = 1):

≤ 0,001 => .

Следовательно, будет решением с требуемой точностью.

С помощью программы посчитаем (t) = 0.214239 t1/3+t2

Листинг 1: Программа для вычисления приближённой функции x(t) (Wolfram Mathematica 7.0):

Clear[x,t];

x[t_]=0;

For[i=1,i<=8,i++,

x[t_] = 0.5*Integrate[(t*s)^(1/3)*x[s], {s, 0, 1}]+t^2;

];

Print[x[t]];

 

Сравним это решение с точным. Положим C = , тогда .

Подставив в исходное уравнение, получим:

(1.1)

Тогда = 4.6E-5 ≤ 0,001.

Рассмотрим пространство [0, 1]. Оценим ядро:

= = ≤ 1.

Т.е. F(x) отображает [0, 1] на себя и является сжимающим при . Поэтому при к данному уравнению можно применить метод сжимающих отображений. В данном случае понадобится число итераций, определяемое соотношением:

< 0,001 (1.2)

Откуда получаем: = < 0,001.

Т.е. при n = 6 достигается требуемая точность.

 

 

Задание 2

Задание: Вычислить приближенное решение уравнения с точностью 0,01:

. (2.1)

Приведём это уравнение к виду . Сделаем это, выразив :

(2.2)

Найдём и радиус такие, что шар B[ , ] инвариантен относительно отображения F и в этом шаре F – сжимающее. Т.к. F является дифференцируемой, то в качестве константы Липшица можно взять .

В данном случае = - . Выберем центр шара = -1, выберем из следующих условий:

(2.3)

 

= , . Тогда:

(2.4)

Выберем одно из решений системы, например, r = 1. Тогда отрезок [-2, 0] инвариантен относительно отображения F, на нём отображение F – сжимающее, и = . Оценим расстояние:

(2.5)

То есть при n = 6 будет достигнута требуемая точность.

С помощью программы посчитаем :

Листинг 2: Программа для вычисления приближённого значения одного из корней (Wolfram Mathematica 7.0):

Clear [F, x];

F[x_]=(1./8)*(-2*x^2-5);

x = -1;

For [i=1, i£6, i++,

x = F[x];

];

Print [x];

 

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

Ответ: , .

 

 

Задание 3

Задание: Определить, является ли отображение f нормированного пространства Е на себя сжимающим. Вычислить , где = , = 0, и оценить расстояние от до неподвижной точки.

Е = [-1, 1], .

Решение: Проверим, является ли f сжимающим.

= .

То есть f является сжимающим отображением с коэффициентом .

Тогда:

= 0,

 

= ,

= ,

= .

Оценим расстояние от до неподвижной точки:

= = 0.3363.

 

Задание 4

Задание: Выяснить, является ли отображение F: X -> Y непрерывным, равномерно непрерывным, удовлетворяющим условию Липшица.

X = C[-2, 4], Y = C[-2, 4], F(x) = .

Решение: Проверим, удовлетворяет ли отображение условию Липшица.

Т.е. удовлетворяет, и константа Липшица Т.к. отображение удовлетворяет условию Липшица, то оно равномерно непрерывно, т.е. для любого ε существует

, что при => ≤ ε.

Т.к. отображение равномерно непрерывное, то оно является и непрерывным.

 





Поделиться с друзьями:


Дата добавления: 2015-10-20; Мы поможем в написании ваших работ!; просмотров: 2516 | Нарушение авторских прав


Поиск на сайте:

Лучшие изречения:

Люди избавились бы от половины своих неприятностей, если бы договорились о значении слов. © Рене Декарт
==> читать все изречения...

2475 - | 2271 -


© 2015-2024 lektsii.org - Контакты - Последнее добавление

Ген: 0.008 с.