Численное решение уравнения f(x)=0 методом хорд и касательных. Метод монотонной прогонки

Учебное издание

Татьяна Владимировна

ЛАБОРАТОРНЫЕ РАБОТЫ ПО ЧИСЛЕННЫМ МЕТОДАМ

Методические указания и задания для cтудентов математических cпециальноcтей

В авторской редакции

Технический редактор

Компьютерный набор и верстка Л. А. Молчановой

Подписано в печать 28.11.2007

Формат 60×84 1/16. Усл. печ. л. 1,3. Уч.-изд. л. 1,5 Тираж 50 экз.

Издательство Дальневосточного университета 690950, Владивосток, ул. Октябрьская, 27.

Отпечатано в лаборатории кафедры компьютерных наук ИМКН ДВГУ

690950, Владивосток, ул. Октябрьская, 27, к. 132.

Министерство образования и науки Российской Федерации

Федеральное агентство по образованию

Дальневосточный государственный университет

ЛАБОРАТОРНЫЕ РАБОТЫ ПО ЧИСЛЕННЫМ МЕТОДАМ

Методические указания и задания для cтудентов математических cпециальноcтей

Владивосток

Издательство Дальневосточного университета 2007

.

Варианты заданий

Таблица 3. Построение сплайна через наклоны

Таблица 4. Построение сплайна через моменты

В полях таблиц указан номер варианта.

Тема 7. Численное решение уравнения f(x)=0 методом хорд и касательных

Найти действительные корни уравнения

                                                                 f(x) = 0                                                            (1)

20

R₂(x) = f^′′′(ξ)ω₃(x)/6, ξ ∈ [x_i−₁,x_i+1], ω₃(x) = (x−x_i−₁)(x−x_i)(x−x_i+1), на [x_i−₁,x_i+1] оценить минимальное и максимальное значения f^′′′(x), а затем минимальное и максимальное значения остаточного члена R₂(x).

•  Проверить, выполняется ли неравенство minR₂ < R₂(x^∗) < maxR₂, R₂(x^∗) = L₂(x^∗) − f(x^∗). Ответить на вопрос п.3.

•  Построить таблицу разделенных разностей по узлам x_i−₁, x_i,x_i+1. Вычислить интерполяционные многочлены Ньютона:

L₁(x^∗) = f(x_i) + f(x_i,x_i+1)(x^∗ − x_i) и

L₂(x^∗) = f(x_i−1)+f(x_i−1,x_i)(x^∗−x_i−1)+f(x_i−1,x_i,x_i+1)(x^∗−x_i−1)·(x^∗−x_i). Сравнить с соответствующими результатами, полученными по формулам Лагранжа.

Таблица 1.

5

Умножая первое уравнение на λ_N−₁ и складывая со вторым, получают: µ       λ             γ             m              γ             λ             m

.

С учетом того, что

                          h2N                      hN                                   h2N

µN−1+λN−1h2_N−1 = h_N + h_N−₁ + (h_N + h_N−₁)h_N−₁ = µN−1(1+γN) = γN,

2 − λ_N−1 + λ_N−1γ_N² = 1 + µ_N−1 + µ_N−1γ_N = 1 + µ_N−1(1 + γ_N) = 1 + γ_N,

имеют

.

Итак, в случае краевых условий 4-ого типа имеют:

;

2;

.

4. Построение кубического сплайна через моменты

Моментами называются;

 на [x_i, x_i+1],

               ′                                                                (x − x_i)²

S (x) = c_i + a_i(x − x_i) + b_i,

2 S^′′(x) = a_i + b_i(x − x_i),

S′′′(x) = bi, S′′(xi) = ai = Mi; S′′(xi+1) = ai + bihi+1 ⇒ bi = Mi+1 − Mi. hi+1

Отсюда:, тогда

.

Из условия непрерывности первой производной для S(x) в узлах x_i, имеют:

,

16

16. n = 4,      t_1,4 = ±0,861136,           c₁ = c₄ = 0,347855, t_2,3 = ±0,339981