(x); = length (y); (nx ~ = ny), ('(x, y) do not have the same # values') k = 1: m = 0; i = 1: nx (i) = 1; j = 1: nx (j ~ = i), (i) = delt (i) * (x (k) - xt (j))/(xt (i ) - xt (j)); = sum + ft (i) * delt (i); (k) = sum; (x, y, 'o', r)
Керівництво програміста:
Призначення програми: Побудова полінома Лагранжа.
Використовувані функції:
plot (x, y) - команда plot (x, y) відповідає побудові звичайної функції, коли одновимірний масив x відповідає значенням аргументу, а одновимірний масив y - значенням функції.
length (x) - вбудована функція, обчислює кількість елементів в одновимірному масиві.
interp1 (x, y, xi, text ) - наближення функції однієї змінної сплайнами, де x, y - це табличні дані вихідної функції, xi - значення аргументу сітки, а text < span align = "justify"> = nearest - інтерполяція по сусідніх елементам, linear - лінійна ітерполяція, Spline - інтерполяція кубічними сплайнів.
Керівництво користувача :
Призначення програми: Побудова полінома Лагранжа c допомогою інтерполяційної формули Лагранжа.
Виконання програми: У віконці Command Window введемо табличні дані вихідної функції x span> і y , а також X i> абсцис точок, в яких обчислюються значення інтерполяційного полінома. Після цього введемо Y = LagrangeP (x, y, X) для побудови полінома Лагранжа. Виводиться графіки вихідної функції і інтерполянти.
Роздруківка серії тестів
Для прикладу візьмемо вихідну функцію:, побудуємо графік цієї функції:
В
Знаходимо точки дискретизації, де крок дорівнює 0.5:
В В
x22, 533,544,555,566,57 y-0 ,68-0 ,57-0 ,03-0 ,23-0 ,74-0, 400,520,700,140,090,65
Інтервал обчислення інтерполянти: від 2 до 7 з кроком 0.1
Використовуємо вбудовану функцію: - інтерполяція по сусідніх елементам.
Будуємо точки дискретизації і графік інтерполянти:
В
- лінійна ітерполяція.
Будуємо точки дискретизації і графік інтерполянти:
В
- інтерполяція кубічними сплайнів.
Будуємо точки дискретизації і графік інтерполянти:
В
Скористаємося написаної функцією:
- інтерполяці...