Формула Тейлора и приложения

Идея полиномиальной аппроксимации

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

«Согласование» означает: многочлен и функция совпадают в точке a по значению, по первой производной, по второй, ..., по n-й. Этот многочлен единственен и называется многочленом Тейлора.

Формула Тейлора

Pₙ(x) = Σₖ₌₀ⁿ f^(k)(a)/k! · (x-a)^k = f(a) + f'(a)(x-a) + f''(a)(x-a)²/2! + ... + f^(n)(a)(x-a)ⁿ/n!

Остаточный член в форме Лагранжа: Rₙ(x) = f^(n+1)(c)/(n+1)! · (x-a)^(n+1) для некоторого c между a и x.

Разложения стандартных функций (a = 0, ряд Маклорена)

eˣ = 1 + x + x²/2! + x³/3! + ... + xⁿ/n! + O(xⁿ⁺¹)

sin x = x - x³/3! + x⁵/5! - ... + O(x²ⁿ⁺¹)

cos x = 1 - x²/2! + x⁴/4! - ... + O(x²ⁿ)

ln(1+x) = x - x²/2 + x³/3 - ... + O(xⁿ⁺¹), |x| ≤ 1

(1+x)^α = 1 + αx + α(α-1)/2! x² + ... (биномиальный ряд)

Применения формулы Тейлора

Вычисление пределов: lim(x→0) (eˣ - 1 - x)/x² = lim(x→0) (x²/2 + O(x³))/x² = 1/2.

Приближённые вычисления с оценкой погрешности: sin(0.1) ≈ 0.1 - (0.1)³/6 = 0.1 - 0.000167 = 0.09983. Погрешность ≤ (0.1)⁵/120 ≈ 10⁻⁸.

Анализ экстремумов при f'(a) = f''(a) = 0: смотрим на первое ненулевое слагаемое разложения.

Применение в физике и технике

Формула Тейлора — инструмент линеаризации. Когда физик записывает «малые колебания», он разлагает потенциальную энергию по Тейлору и оставляет квадратичный член. Когда инженер анализирует устойчивость системы, он линеаризует уравнение в окрестности равновесия.

Формула eⁱˣ = cos x + i sin x (формула Эйлера) вытекает из разложений eˣ, sin x, cos x по Тейлору — это одна из красивейших формул математики, связывающая показательную функцию с тригонометрическими.

Погрешность аппроксимации и практика

Оценка остаточного члена через формулу Лагранжа |Rₙ(x)| ≤ M_{n+1}/(n+1)! · |x−a|^(n+1) критически важна в вычислительной практике. Рассмотрим задачу: вычислить e с точностью 10⁻⁶. Нужно найти n так, чтобы |Rₙ(1)| ≤ 1/n! · |e| ≤ 3/n! < 10⁻⁶. Из таблицы: 3/9! = 3/362880 ≈ 8.3·10⁻⁶ (мало), 3/10! ≈ 8.3·10⁻⁷ < 10⁻⁶. Значит, достаточно 10 членов ряда.

Метод Тейлора в численном анализе

Формула Тейлора — теоретическая основа численного интегрирования и дифференцирования. Метод Рунге–Кутта для ОДУ: Разложение решения по Тейлору, усечённое до степени 4, даёт ошибку O(h⁵) на каждом шаге. Численное дифференцирование: центральная разность (f(x+h) − f(x−h))/(2h) = f'(x) + O(h²) — ошибка второго порядка, потому что члены нечётных степеней взаимно уничтожаются при вычитании разложений Тейлора.

Трудность: при малых h вступает в силу погрешность округления в компьютерных вычислениях. Оптимальный шаг h ≈ √ε (ε — машинная точность), что составляет примерно 10⁻⁸ для 64-битной арифметики.

Многомерный ряд Тейлора и матрица Гессе

Для функции f(x₁, ..., xₙ) ряд Тейлора вокруг точки a до второго порядка:

f(a + h) ≈ f(a) + ∇f(a)ᵀh + (1/2)hᵀH(a)h,

где H — матрица Гессе (матрица вторых производных). Это квадратичная форма определяет, является ли точка a минимумом (H > 0), максимумом (H < 0) или седловой точкой (неопределённый знак). Именно этот факт лежит в основе методов оптимизации второго порядка (метод Ньютона в многомерном случае).

Оценка остатка и практическая точность

Остаток Лагранжа Rₙ(x) = f^(n+1)(c)/(n+1)! · (x−a)^(n+1) позволяет оценить точность приближения. Для sin x: |Rₙ(x)| ≤ |x|^(n+1)/(n+1)! — при x = 0.1 и n = 3 погрешность не превышает 0.1⁴/24 ≈ 4·10⁻⁶. Это делает полиномы Тейлора практическим инструментом вычислительной математики: инженерные библиотеки реализуют тригонометрические и показательные функции именно через усечённые ряды Тейлора с контролем остатка.

Вопрос для размышления: Почему разложение eˣ в ряд Тейлора сходится для всех вещественных x, а разложение 1/(1+x²) — только при |x| < 1? Как это связано с полюсами в комплексной плоскости?

Вычисление пределов методом разложения Тейлора

Разложение Тейлора — мощный инструмент вычисления пределов. Пример 1: lim(x→0) (e^x − 1 − x)/x² = lim (1 + x + x²/2 + ... − 1 − x)/x² = lim (x²/2 + O(x³))/x² = 1/2. Пример 2: lim(x→0) (sin x − x)/x³ = lim (x − x³/6 + ... − x)/x³ = −1/6. Пример 3: lim(x→0) (1 − cos x)/x² = lim (x²/2 − x⁴/24 + ...)/x² = 1/2. Во всех случаях метод одинаков: разложить числитель до нужного порядка и сократить с знаменателем. Этот подход значительно эффективнее многократного применения правила Лопиталя.