Добавил:

Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.

Вуз:

Белорусский национальный технический университет

Предмет:

[НЕСОРТИРОВАННОЕ]

Файл:

курсовой_ДубатовкаСВ_ч2.doc

Скачиваний:

Добавлен:

17.12.2018

Размер:

465.41 Кб

Скачать

☆

<<< < Предыдущая 12 / 32 3 > Следующая >>>

4. Приближённое решение краевых задач и задач о собственных значениях метода конечных разностей.

Основная мысль этого метода состоит в следующем: пусть дана краевая задача

Для того чтобы получить приближенное решение этой задачи, которую мы предполагаем разрешимой, разделим отрезок [а,b] на равные части длины и для каждой точки деления напишем уравнение, которое получается из дифференциального уравнения (5), если в нем производные заменить их приближенными выражениями с помощью соответствующим образом выбранных разностных отношений. Такую же замену выполним и для производных, входящих в краевые условия. Тогда для нахождения приближенной величины значения искомого решения мы получим систему линейных уравнений, и тем самым исходную задачу сведем к задаче конечной.

Если дана задача о собственных значениях

то коэффициенты получающейся однородной системы линейных уравнений зависят еще от . Значение следует выбрать так, чтобы система имела нетривиальное решение. Таким образом, для нахождения приближенной величины первого собственного значения мы получаем алгебраическое уравнение.

Применение выражений более высоких порядков дает при той же длине шага увеличение точности. Здесь Rр означает некоторое число, меньшее, чем максимум Ι Ι на некотором интервале, содержащем все абсциссы, встречающиеся в приведенных выше выражениях.

Если пользоваться конечными разностями более высоких порядков так же и для краевых условий, то при этом могут получиться члены, содержащие значения Y, соответствующее точкам, лежащим вне интервала (а,b). Этого можно избежать, если на концах данного интервала пользоваться приближениями более низких порядков.

С теоретической точки зрения интересно то, что этот метод даже и при более грубом способе расчетов дает приближенные решения, которые при стремится к точному решению. Именно, согласно Планшерелю, имеет место следующее: пусть в краевой задаче

(5)

и дважды непрерывно дифференцируема, а и непрерывны на отрезке . Если производные заменить через

то соответствующая система алгебраических уравнений будет иметь вид

кроме того, краевые условия дает Если не есть собственное значение соответствующей однородной краевой задачи, т.е. (8) имеет единственное решение , тогда где

Если есть -е собственное значение задачи и если - упорядоченные по возрастанию нули детерминанта написанной выше линейной системы, то при .

5. Решение систем нелинейных уравнений

В отличие от систем линейных уравнений для систем нелинейных уравнений не известны прямые методы решения. Лишь в отдельных случаях систему можно решить непосредственно. Например, для системы из двух уравнений иногда удается выразить одно неизвестное через другое и таким образом свести задачу к решению одного нелинейного уравнения относительно одного неизвестного. Поэтому итерационные методы для нелинейных систем приобретают особую актуальность.

Метод Ньютона.

Рассмотрим нелинейную систему уравнений

(1)

или в векторной форме

f (x) = 0,

(1 ')

где

f x

Для решения системы (23 ) будем пользоваться методом последовательных приближений.

Предположим, известно k-е приближение

x⁽^k⁾=

одного из изолированных корней x= векторного уравнения (23 '). Тогда точный корень уравнения (1') можно представить в виде

х = x⁽^k⁾ +  x⁽^k⁾,

(2)

где  x⁽^k⁾ = - поправка (погрешность корня).

Подставляя выражение (24) в (23'), будем иметь

f (x⁽^k⁾ +  x⁽^k⁾) = 0.

(3)

Предполагая, что функция f (x) непрерывно дифференцируема в некоторой выпуклой области, содержащей x и x⁽^k⁾, разложим левую часть уравнения (3) по степеням малого вектора  x⁽^k⁾ , ограничиваясь линейными членами,

f (x⁽^k⁾ +  x⁽^k⁾) = f (x⁽^k⁾) + f  (x⁽^k⁾)  x⁽^k⁾= 0

(4)

или, в развернутом виде,

(4')

Из формул (4) и (4') вытекает, что под производной f '(x) следует понимать матрицу Якоби системы функций f₁, f₂, ..., f_nотносительно переменных x₁, x₂, ..., x_n, т. е.

f ' (x) = W(x) =,

или в краткой записи

f ' (x) = W(x) = (i, j = 1, 2, …, n).

Поэтому формула (4) может быть записана в следующем виде:

f (x⁽^k⁾ ) + W (x⁽^k⁾ )  x⁽^k⁾ = 0

Если det W (х) =, то  x⁽^k⁾ = - W ^-1(x⁽^k⁾) f (x⁽^k⁾).

Отсюда видно, что метод Ньютона решения системы (1) состоит в построении итерационной последовательности:

x⁽^{k
+ 1}⁾ = x⁽^k⁾ - W ^-1(x⁽^k⁾) f (x⁽^k⁾) (k = 0, 1, 2, …).

(5)

Если все поправки становятся достаточно малыми, счет прекращается. Иначе новые значения x_i используются как приближенные значения корней, и процесс повторяется до тех пор, пока не будет найдено решение или не станет ясно, что получить его не удастся.

Пример 1. Методом Ньютона приближенно найти положительное решение системы уравнений