Дан ряд Тейлора функции ln(1+x) при x в (-1,1]. Объясните, как оценить погрешность при приближении функции третьим порядком разложения и приведите практические рекомендации
- Ряд и многочлен: ln(1+x)=∑n=1∞(−1) n−1xnn\displaystyle \ln(1+x)=\sum_{n=1}^\infty (-1)^{\,n-1}\frac{x^n}{n}ln(1+x)=n=1∑∞(−1)n−1nxn при −1<x≤1-1<x\le1−1<x≤1. Третьего порядка (Maclaurin) P3(x)=x−x22+x33\displaystyle P_3(x)=x-\frac{x^2}{2}+\frac{x^3}{3}P3(x)=x−2x2+3x3. Погрешность R3(x)=ln(1+x)−P3(x)R_3(x)=\ln(1+x)-P_3(x)R3(x)=ln(1+x)−P3(x). - Форма Лагранжа: для некоторого ξ\xiξ между 000 и xxxR3(x)=f(4)(ξ)4!x4,f(4)(t)=−6(1+t)4,
R_3(x)=\frac{f^{(4)}(\xi)}{4!}x^4,\qquad f^{(4)}(t)=-\frac{6}{(1+t)^4}, R3(x)=4!f(4)(ξ)x4,f(4)(t)=−(1+t)46,
откуда R3(x)=−x44(1+ξ)4,∣R3(x)∣=∣x∣44(1+ξ)4.
R_3(x)=-\frac{x^4}{4(1+\xi)^4},\qquad |R_3(x)|=\frac{|x|^4}{4(1+\xi)^4}. R3(x)=−4(1+ξ)4x4,∣R3(x)∣=4(1+ξ)4∣x∣4. - Практичные оценки: - Для 0≤x≤10\le x\le10≤x≤1: 1+ξ≥1⇒∣R3(x)∣≤x441+\xi\ge1\Rightarrow |R_3(x)|\le\frac{x^4}{4}1+ξ≥1⇒∣R3(x)∣≤4x4. (То же даёт признак знакоперемённого ряда: погрешность не превосходит модуля следующего члена.) Пример: для x=12x=\tfrac12x=21 получаем ∣R3∣≤(1/2)44=164=0.015625|R_3|\le\frac{(1/2)^4}{4}=\tfrac{1}{64}=0.015625∣R3∣≤4(1/2)4=641=0.015625. - Для −1<x<0-1<x<0−1<x<0: так как 1+ξ≥1+x1+\xi\ge1+x1+ξ≥1+x (и 1+x∈(0,1)1+x\in(0,1)1+x∈(0,1)), имеем более жёсткую оценку ∣R3(x)∣≤∣x∣44(1+x)4.
|R_3(x)|\le\frac{|x|^4}{4(1+x)^4}. ∣R3(x)∣≤4(1+x)4∣x∣4.
Пример: для x=−12x=-\tfrac12x=−21 получается ∣R3∣≤0.25|R_3|\le0.25∣R3∣≤0.25. - Заметьте: при x→−1+x\to-1^+x→−1+ оценка растёт неограниченно — приближение рядом в этом виде плохо близко к −1-1−1. - Практические рекомендации: 1. Для положительных xxx (особенно малого) пользуйтесь оценкой ∣R3∣≤x4/4|R_3|\le x^4/4∣R3∣≤x4/4; если нужна точность ε\varepsilonε, требуйте x4/4<εx^4/4<\varepsilonx4/4<ε или повышайте порядок. 2. Для xxx близких к −1-1−1 не используйте разложение в нуле — вместо этого: - делайте разложение вблизи другой точки aaa (например aaa такое, что 1+a1+a1+a далеко от нуля) и вычисляйте полином в точке xxx (т.н. развёртка в aaa); или - применяйте Padé-аппроксимации или стандартные численные функции логарифма. 3. Если нужен строгий численный верхний предел ошибки для данного xxx, используйте формулу Лагранжа и подставьте нижнюю границу для (1+ξ)4(1+\xi)^4(1+ξ)4 (для x≥0x\ge0x≥0 — 111, для x<0x<0x<0 — (1+x)4(1+x)^4(1+x)4).
ln(1+x)=∑n=1∞(−1) n−1xnn\displaystyle \ln(1+x)=\sum_{n=1}^\infty (-1)^{\,n-1}\frac{x^n}{n}ln(1+x)=n=1∑∞ (−1)n−1nxn при −1<x≤1-1<x\le1−1<x≤1. Третьего порядка (Maclaurin)
P3(x)=x−x22+x33\displaystyle P_3(x)=x-\frac{x^2}{2}+\frac{x^3}{3}P3 (x)=x−2x2 +3x3 . Погрешность R3(x)=ln(1+x)−P3(x)R_3(x)=\ln(1+x)-P_3(x)R3 (x)=ln(1+x)−P3 (x).
- Форма Лагранжа: для некоторого ξ\xiξ между 000 и xxx R3(x)=f(4)(ξ)4!x4,f(4)(t)=−6(1+t)4, R_3(x)=\frac{f^{(4)}(\xi)}{4!}x^4,\qquad f^{(4)}(t)=-\frac{6}{(1+t)^4},
R3 (x)=4!f(4)(ξ) x4,f(4)(t)=−(1+t)46 , откуда
R3(x)=−x44(1+ξ)4,∣R3(x)∣=∣x∣44(1+ξ)4. R_3(x)=-\frac{x^4}{4(1+\xi)^4},\qquad |R_3(x)|=\frac{|x|^4}{4(1+\xi)^4}.
R3 (x)=−4(1+ξ)4x4 ,∣R3 (x)∣=4(1+ξ)4∣x∣4 .
- Практичные оценки:
- Для 0≤x≤10\le x\le10≤x≤1: 1+ξ≥1⇒∣R3(x)∣≤x441+\xi\ge1\Rightarrow |R_3(x)|\le\frac{x^4}{4}1+ξ≥1⇒∣R3 (x)∣≤4x4 . (То же даёт признак знакоперемённого ряда: погрешность не превосходит модуля следующего члена.)
Пример: для x=12x=\tfrac12x=21 получаем ∣R3∣≤(1/2)44=164=0.015625|R_3|\le\frac{(1/2)^4}{4}=\tfrac{1}{64}=0.015625∣R3 ∣≤4(1/2)4 =641 =0.015625.
- Для −1<x<0-1<x<0−1<x<0: так как 1+ξ≥1+x1+\xi\ge1+x1+ξ≥1+x (и 1+x∈(0,1)1+x\in(0,1)1+x∈(0,1)), имеем более жёсткую оценку
∣R3(x)∣≤∣x∣44(1+x)4. |R_3(x)|\le\frac{|x|^4}{4(1+x)^4}.
∣R3 (x)∣≤4(1+x)4∣x∣4 . Пример: для x=−12x=-\tfrac12x=−21 получается ∣R3∣≤0.25|R_3|\le0.25∣R3 ∣≤0.25.
- Заметьте: при x→−1+x\to-1^+x→−1+ оценка растёт неограниченно — приближение рядом в этом виде плохо близко к −1-1−1.
- Практические рекомендации:
1. Для положительных xxx (особенно малого) пользуйтесь оценкой ∣R3∣≤x4/4|R_3|\le x^4/4∣R3 ∣≤x4/4; если нужна точность ε\varepsilonε, требуйте x4/4<εx^4/4<\varepsilonx4/4<ε или повышайте порядок.
2. Для xxx близких к −1-1−1 не используйте разложение в нуле — вместо этого:
- делайте разложение вблизи другой точки aaa (например aaa такое, что 1+a1+a1+a далеко от нуля) и вычисляйте полином в точке xxx (т.н. развёртка в aaa); или
- применяйте Padé-аппроксимации или стандартные численные функции логарифма.
3. Если нужен строгий численный верхний предел ошибки для данного xxx, используйте формулу Лагранжа и подставьте нижнюю границу для (1+ξ)4(1+\xi)^4(1+ξ)4 (для x≥0x\ge0x≥0 — 111, для x<0x<0x<0 — (1+x)4(1+x)^4(1+x)4).