💀
Второй курс РПО
Численные методы
Численные методы
  • Погрешности вычислений
  • Вычисление абсолютной и относительной погрешностей
    • Вариант 2
    • Вариант 3
    • Вариант 4
  • Локализация корня уравнения
  • Метод хорд
  • Метод Ньютона
  • Метод деления отрезка пополам
  • Практическая № X
  • Приближённые методы
    • Задача: Вариант 3 #1
    • Задача: Вариант 3 #2
    • Задача: Вариант 3 #3
  • Метод гауса: С чем его едят?
  • Интерполяционный многочлен Лагранжа
    • Вариант 1
  • Интерполяция: Формулы Ньютона
  • Полигональное моделирование
    • Интерполяция кубическими сплайнами
  • Самостоятельная В3
  • Распределение Пуассона
  • Суть интерполяции Ньютона
    • Объяснение
  • Суть интерполяции сплайнами
  • Численное интегрирование
  • Метод Гауса
  • Метод Эйлера
Powered by GitBook
On this page
  • Формула Лагранжа
  • Этап 1: Вычисление базисных многочленов Лагранжа
  • Этап 2: Построение интерполяционной полинома

Самостоятельная В3

04.12.2024

PreviousИнтерполяция кубическими сплайнамиNextРаспределение Пуассона

Last updated 5 months ago


Формула Лагранжа

Интерполяционная формула Лагранжа для функции f(x)f(x)f(x) на узлах x0,x1,…,xnx_0, x_1, \dots, x_n x0​,x1​,…,xn​выглядит следующим образом:

Pn(x)=∑i=0nf(xi)⋅Li(x)P_n(x) = \sum_{i=0}^{n} f(x_i) \cdot L_i(x)Pn​(x)=i=0∑n​f(xi​)⋅Li​(x)

где Li(x)L_i(x)Li​(x) — это базисные многочлены Лагранжа, которые вычисляются по формуле:

Li(x)=∏0≤j≤nj≠ix−xjxi−xjL_i(x) = \prod_{\substack{0 \leq j \leq n \\ j \neq i}} \frac{x - x_j}{x_i - x_j}Li​(x)=0≤j≤nj=i​∏​xi​−xj​x−xj​​

В данной задаче у нас n=2n = 2n=2, то есть интерполяция на 3 узлах 000, π4\frac{\pi}{4} 4π​, fracπ2frac{\pi}{2} fracπ2.

Этап 1: Вычисление базисных многочленов Лагранжа

Для каждого узла xix_ixi​ вычислим Li(x)L_i(x)Li​(x).

Для L0(x)L_0(x)L0​(x):

L0(x)=(x−x1)(x−x2)(x0−x1)(x0−x2)L_0(x) = \frac{(x - x_1)(x - x_2)}{(x_0 - x_1)(x_0 - x_2)}L0​(x)=(x0​−x1​)(x0​−x2​)(x−x1​)(x−x2​)​

Подставляем значения x0=0x_0 = 0x0​=0, x1=π4x_1 = \frac{\pi}{4} x1​=4π​, x2=π2x_2 = \frac{\pi}{2}x2​=2π​:

L0(x)=(x−π4)(x−π2)(0−π4)(0−π2)L_0(x) = \frac{(x - \frac{\pi}{4})(x - \frac{\pi}{2})}{(0 - \frac{\pi}{4})(0 - \frac{\pi}{2})}L0​(x)=(0−4π​)(0−2π​)(x−4π​)(x−2π​)​

Упростим:

L0(x)=(x−π4)(x−π2)(−π4)(−π2)=(x−π4)(x−π2)π28L_0(x) = \frac{(x - \frac{\pi}{4})(x - \frac{\pi}{2})}{\left(-\frac{\pi}{4}\right) \left(-\frac{\pi}{2}\right)} = \frac{(x - \frac{\pi}{4})(x - \frac{\pi}{2})}{\frac{\pi^2}{8}}L0​(x)=(−4π​)(−2π​)(x−4π​)(x−2π​)​=8π2​(x−4π​)(x−2π​)​

Для L1(x)L_1(x)L1​(x):

L1(x)=(x−x0)(x−x2)(x1−x0)(x1−x2)L_1(x) = \frac{(x - x_0)(x - x_2)}{(x_1 - x_0)(x_1 - x_2)}L1​(x)=(x1​−x0​)(x1​−x2​)(x−x0​)(x−x2​)​

Подставляем значения:

L1(x)=(x−0)(x−π2)(π4−0)(π4−π2)L_1(x) = \frac{(x - 0)(x - \frac{\pi}{2})}{(\frac{\pi}{4} - 0)(\frac{\pi}{4} - \frac{\pi}{2})}L1​(x)=(4π​−0)(4π​−2π​)(x−0)(x−2π​)​

Упростим:

L1(x)=x(x−π2)π4⋅(−π4)=x(x−π2)−π216L_1(x) = \frac{x(x - \frac{\pi}{2})}{\frac{\pi}{4} \cdot \left(-\frac{\pi}{4}\right)} = \frac{x(x - \frac{\pi}{2})}{-\frac{\pi^2}{16}}L1​(x)=4π​⋅(−4π​)x(x−2π​)​=−16π2​x(x−2π​)​

Для L2(x)L_2(x)L2​(x):

L2(x)=(x−x0)(x−x1)(x2−x0)(x2−x1)L_2(x) = \frac{(x - x_0)(x - x_1)}{(x_2 - x_0)(x_2 - x_1)}L2​(x)=(x2​−x0​)(x2​−x1​)(x−x0​)(x−x1​)​

Подставляем значения:

L2(x)=(x−0)(x−π4)(π2−0)(π2−π4)L_2(x) = \frac{(x - 0)(x - \frac{\pi}{4})}{(\frac{\pi}{2} - 0)(\frac{\pi}{2} - \frac{\pi}{4})}L2​(x)=(2π​−0)(2π​−4π​)(x−0)(x−4π​)​

Упростим:

L2(x)=x(x−π4)π2⋅π4=x(x−π4)π28L_2(x) = \frac{x(x - \frac{\pi}{4})}{\frac{\pi}{2} \cdot \frac{\pi}{4}} = \frac{x(x - \frac{\pi}{4})}{\frac{\pi^2}{8}}L2​(x)=2π​⋅4π​x(x−4π​)​=8π2​x(x−4π​)​

Этап 2: Построение интерполяционной полинома

Теперь, используя значения f(x0)=sin⁡(0)=0f(x_0) = \sin(0) = 0 f(x0​)=sin(0)=0, f(x1)=sin⁡(π4)=22f(x_1) = \sin(\frac{\pi}{4}) = \frac{\sqrt{2}}{2}f(x1​)=sin(4π​)=22​​, f(x2)=sin⁡(π2)=1f(x_2) = \sin(\frac{\pi}{2}) = 1f(x2​)=sin(2π​)=1, мы можем записать интерполяционную формулу Лагранжа:

P2(x)=f(x0)⋅L0(x)+f(x1)⋅L1(x)+f(x2)⋅L2(x)P_2(x) = f(x_0) \cdot L_0(x) + f(x_1) \cdot L_1(x) + f(x_2) \cdot L_2(x)P2​(x)=f(x0​)⋅L0​(x)+f(x1​)⋅L1​(x)+f(x2​)⋅L2​(x)

Подставляем все значения:

P2(x)=0⋅L0(x)+22⋅L1(x)+1⋅L2(x)P_2(x) = 0 \cdot L_0(x) + \frac{\sqrt{2}}{2} \cdot L_1(x) + 1 \cdot L_2(x)P2​(x)=0⋅L0​(x)+22​​⋅L1​(x)+1⋅L2​(x)

Таким образом, формула для интерполяционного полинома:

P2(x)=22⋅x(x−π2)−π216+x(x−π4)π28P_2(x) = \frac{\sqrt{2}}{2} \cdot \frac{x(x - \frac{\pi}{2})}{-\frac{\pi^2}{16}} + \frac{x(x - \frac{\pi}{4})}{\frac{\pi^2}{8}}P2​(x)=22​​⋅−16π2​x(x−2π​)​+8π2​x(x−4π​)​

Этап 3: Упрощение

Теперь можно упростить полином:

P2(x)=22⋅16x(x−π2)−π2+8x(x−π4)π2P_2(x) = \frac{\sqrt{2}}{2} \cdot \frac{16x(x - \frac{\pi}{2})}{-\pi^2} + \frac{8x(x - \frac{\pi}{4})}{\pi^2}P2​(x)=22​​⋅−π216x(x−2π​)​+π28x(x−4π​)​
P2(x)=−82x(x−π2)π2+8x(x−π4)π2P_2(x) = \frac{-8\sqrt{2}x(x - \frac{\pi}{2})}{\pi^2} + \frac{8x(x - \frac{\pi}{4})}{\pi^2}P2​(x)=π2−82​x(x−2π​)​+π28x(x−4π​)​

Это и есть искомая интерполяционная формула Лагранжа для функции f(x)=sin⁡(x)f(x) = \sin(x)f(x)=sin(x) на узлах x0=0x_0 = 0x0​=0, x1=π4x_1 = \frac{\pi}{4}x1​=4π​, x2=π2x_2 = \frac{\pi}{2}x2​=2π​.