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

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

Метод последовательных приближений - student2.ru , (6)

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

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

Метод последовательных приближений - student2.ru , (7)

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

Метод последовательных приближений - student2.ru (8)

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

Метод последовательных приближений - student2.ru , (9)

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

Задание 1

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

Решение:

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

F(x) = Метод последовательных приближений - student2.ru + Метод последовательных приближений - student2.ru ;

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ruМетод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru .

Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru ≤ 1 => Метод последовательных приближений - student2.ruМетод последовательных приближений - student2.ru .

Т.е. при Метод последовательных приближений - student2.ruМетод последовательных приближений - student2.ru применим метод сжимающих отображений.

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

Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru ≤ 0,001 => Метод последовательных приближений - student2.ruМетод последовательных приближений - student2.ru .

Следовательно, Метод последовательных приближений - student2.ru будет решением с требуемой точностью.

С помощью программы посчитаем Метод последовательных приближений - student2.ru (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 = Метод последовательных приближений - student2.ru , тогда Метод последовательных приближений - student2.ru .

Подставив Метод последовательных приближений - student2.ru в исходное уравнение, получим:

Метод последовательных приближений - student2.ru (1.1)

Тогда Метод последовательных приближений - student2.ru = 4.6E-5 ≤ 0,001.

Рассмотрим пространство Метод последовательных приближений - student2.ru [0, 1]. Оценим ядро:

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru ≤ 1.

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

Метод последовательных приближений - student2.ru < 0,001 (1.2)

Откуда получаем: Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru < 0,001.

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

Задание 2

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

Метод последовательных приближений - student2.ru . (2.1)

Приведём это уравнение к виду Метод последовательных приближений - student2.ru . Сделаем это, выразив Метод последовательных приближений - student2.ru :

Метод последовательных приближений - student2.ru (2.2)

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

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

Метод последовательных приближений - student2.ru (2.3)

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru , Метод последовательных приближений - student2.ru . Тогда:

Метод последовательных приближений - student2.ru (2.4)

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

Метод последовательных приближений - student2.ru (2.5)

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

С помощью программы посчитаем Метод последовательных приближений - student2.ru :

Листинг 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];

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

Ответ: Метод последовательных приближений - student2.ru , Метод последовательных приближений - student2.ru .

Задание 3

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

Е = Метод последовательных приближений - student2.ru [-1, 1], Метод последовательных приближений - student2.ru .

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

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ruМетод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru Метод последовательных приближений - student2.ru .

То есть f является сжимающим отображением с коэффициентом Метод последовательных приближений - student2.ru .

Тогда:

Метод последовательных приближений - student2.ru = 0,

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru ,

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru ,

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru .

Оценим расстояние от Метод последовательных приближений - student2.ru до неподвижной точки:

Метод последовательных приближений - student2.ru = Метод последовательных приближений - student2.ru = 0.3363.

Задание 4

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

X = C[-2, 4], Y = C[-2, 4], F(x) = Метод последовательных приближений - student2.ru .

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

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

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

Метод последовательных приближений - student2.ru , что при Метод последовательных приближений - student2.ru => Метод последовательных приближений - student2.ru ≤ ε.

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

Наши рекомендации