Дан ряд Тейлора функции f(x) = arctan x в окрестности нуля. Сформулируйте задачу оценки остаточного члена при конечном числе членов ряда и обсудите разные способы получения верхних оценок остатка, укажите, какие методы даются более точные результаты на заданном отрезке
Ряд Тейлора в окрестности нуля: arctanx=∑n=0∞(−1)nx2n+12n+1,∣x∣≤1 (при x=±1 сходится условно).
\arctan x=\sum_{n=0}^\infty (-1)^n\frac{x^{2n+1}}{2n+1},\qquad |x|\le 1\ (\text{при }x=\pm1\ \text{сходится условно}). arctanx=n=0∑∞(−1)n2n+1x2n+1,∣x∣≤1(приx=±1сходитсяусловно). Остаточный член при конечном числе членов (приблизим суммой до NNN-го члена, т.е. степени 2N+12N+12N+1): обозначим SN(x)=∑n=0N(−1)nx2n+12n+1,RN(x)=arctanx−SN(x).
S_N(x)=\sum_{n=0}^N (-1)^n\frac{x^{2n+1}}{2n+1},\qquad R_N(x)=\arctan x-S_N(x). SN(x)=n=0∑N(−1)n2n+1x2n+1,RN(x)=arctanx−SN(x). Способы оценки остатка и результаты: 1) Оценка через знакочередующийся ряд (Лейбниц). Для ∣x∣≤1|x|\le1∣x∣≤1 члены по абсолютной величине убывают, поэтому ∣RN(x)∣≤∣x∣2N+32N+3,
|R_N(x)|\le \frac{|x|^{2N+3}}{2N+3}, ∣RN(x)∣≤2N+3∣x∣2N+3,
и при этом RN(x)R_N(x)RN(x) имеет знак (−1)N+1(-1)^{N+1}(−1)N+1 (т.е. следующий по порядку член задаёт направление ошибки). Этот метод прост и даёт точную по порядку (и часто почти точную по коэффициенту) оценку при ∣x∣≤1|x|\le1∣x∣≤1. 2) Интегральное представление остатка. Так как arctanx=∫0xdt1+t2\arctan x=\int_0^x\frac{dt}{1+t^2}arctanx=∫0x1+t2dt и для ∣t∣<1|t|<1∣t∣<111+t2=∑n=0∞(−1)nt2n\frac{1}{1+t^2}=\sum_{n=0}^\infty(-1)^n t^{2n}1+t21=∑n=0∞(−1)nt2n, получается точная формула RN(x)=∫0x(−1) N+1t2N+21+t2 dt.
R_N(x)=\int_0^x \frac{(-1)^{\,N+1} t^{2N+2}}{1+t^2}\,dt. RN(x)=∫0x1+t2(−1)N+1t2N+2dt.
Отсюда немедленно следует оценка ∣RN(x)∣≤∫0∣x∣t2N+21+t2 dt≤∣x∣2N+32N+3,
|R_N(x)|\le\int_0^{|x|}\frac{t^{2N+2}}{1+t^2}\,dt\le\frac{|x|^{2N+3}}{2N+3}, ∣RN(x)∣≤∫0∣x∣1+t2t2N+2dt≤2N+3∣x∣2N+3,
и при аккуратном учёте знаменателя можно иногда получить чуть лучшую (но по порядку тот же) оценку для конкретного xxx. 3) Форма Лагранжа (обобщённая оценка). Для многочлена степени mmmRm(x)=f(m+1)(ξ)(m+1)!xm+1для некоторого ξ∈(0,x).
R_m(x)=\frac{f^{(m+1)}(\xi)}{(m+1)!}x^{m+1}\quad\text{для некоторого }\xi\in(0,x). Rm(x)=(m+1)!f(m+1)(ξ)xm+1длянекоторогоξ∈(0,x).
Это даёт оценку, если можно явно ограничить ∣f(m+1)(t)∣|f^{(m+1)}(t)|∣f(m+1)(t)∣ на отрезке между 000 и xxx. Для arctan\arctanarctan производные имеют вид рациональных функций с знаменателем (1+t2)k(1+t^2)^k(1+t2)k, поэтому одна получает оценку вида ∣Rm(x)∣≤max∣t∣≤∣x∣∣f(m+1)(t)∣(m+1)!∣x∣m+1.
|R_m(x)|\le \frac{\max_{|t|\le|x|}|f^{(m+1)}(t)|}{(m+1)!}|x|^{m+1}. ∣Rm(x)∣≤(m+1)!max∣t∣≤∣x∣∣f(m+1)(t)∣∣x∣m+1.
Преимущество: применяется для любых xxx. Недостаток: при прямом грубом оценивании max∣f(m+1)∣\max|f^{(m+1)}|max∣f(m+1)∣ получается большая (иногда чрезмерно грубая) оценка по сравнению с методом Лейбинца на ∣x∣≤1|x|\le1∣x∣≤1. 4) Случай больших аргументов (∣x∣>1|x|>1∣x∣>1). Лучше сначала использовать тождество arctanx=π2sgn(x)−arctan1x,
\arctan x=\frac{\pi}{2}\operatorname{sgn}(x)-\arctan\frac{1}{x}, arctanx=2πsgn(x)−arctanx1,
и разложить arctan(1/x)\arctan(1/x)arctan(1/x) как ряд по малой переменной 1/x1/x1/x. Тогда оценка остатка даётся тем же правилом Лейбинца: ∣RN(x)∣≤∣1/x∣2N+32N+3,
|R_N(x)|\le\frac{|1/x|^{2N+3}}{2N+3}, ∣RN(x)∣≤2N+3∣1/x∣2N+3,
что даёт очень маленькую ошибку при больших ∣x∣|x|∣x∣. 5) Аппроксимации типа Паде, Чебышёва и минимакс. Для одновременной хорошей равномерной оценки на отрезке (особенно если нужен малый максимум ошибки на всём интервале) рациональные приближения (Паде) или полиномы с равномерной минимакс-ошибкой (Чебышёв) дают существенно лучшие результаты по сравнению с голым усечённым степенным рядом; для фиксированной степени они обычно дают намного меньшую максимальную ошибку. Эти методы сложнее в применении, но предпочтительны для практических вычислений с ограниченным числом членов на заданном отрезке. Краткие рекомендации: - Для точечной оценки при ∣x∣≤1|x|\le1∣x∣≤1 и простоты: используйте оценку Лейбинца (или интегральную форму) — она простая и почти оптимальная. - Для ∣x∣>1|x|>1∣x∣>1: сначала замените x↦1/xx\mapsto1/xx↦1/x через тождество, затем примените ту же оценку. - Для равномерной аппроксимации на отрезке: используйте Чебышёв- или Паде-аппроксимации — они дают более точные верхние оценки при том же числе параметров.
arctanx=∑n=0∞(−1)nx2n+12n+1,∣x∣≤1 (при x=±1 сходится условно). \arctan x=\sum_{n=0}^\infty (-1)^n\frac{x^{2n+1}}{2n+1},\qquad |x|\le 1\ (\text{при }x=\pm1\ \text{сходится условно}).
arctanx=n=0∑∞ (−1)n2n+1x2n+1 ,∣x∣≤1 (при x=±1 сходится условно).
Остаточный член при конечном числе членов (приблизим суммой до NNN-го члена, т.е. степени 2N+12N+12N+1):
обозначим
SN(x)=∑n=0N(−1)nx2n+12n+1,RN(x)=arctanx−SN(x). S_N(x)=\sum_{n=0}^N (-1)^n\frac{x^{2n+1}}{2n+1},\qquad R_N(x)=\arctan x-S_N(x).
SN (x)=n=0∑N (−1)n2n+1x2n+1 ,RN (x)=arctanx−SN (x).
Способы оценки остатка и результаты:
1) Оценка через знакочередующийся ряд (Лейбниц). Для ∣x∣≤1|x|\le1∣x∣≤1 члены по абсолютной величине убывают, поэтому
∣RN(x)∣≤∣x∣2N+32N+3, |R_N(x)|\le \frac{|x|^{2N+3}}{2N+3},
∣RN (x)∣≤2N+3∣x∣2N+3 , и при этом RN(x)R_N(x)RN (x) имеет знак (−1)N+1(-1)^{N+1}(−1)N+1 (т.е. следующий по порядку член задаёт направление ошибки). Этот метод прост и даёт точную по порядку (и часто почти точную по коэффициенту) оценку при ∣x∣≤1|x|\le1∣x∣≤1.
2) Интегральное представление остатка. Так как arctanx=∫0xdt1+t2\arctan x=\int_0^x\frac{dt}{1+t^2}arctanx=∫0x 1+t2dt и для ∣t∣<1|t|<1∣t∣<1 11+t2=∑n=0∞(−1)nt2n\frac{1}{1+t^2}=\sum_{n=0}^\infty(-1)^n t^{2n}1+t21 =∑n=0∞ (−1)nt2n, получается точная формула
RN(x)=∫0x(−1) N+1t2N+21+t2 dt. R_N(x)=\int_0^x \frac{(-1)^{\,N+1} t^{2N+2}}{1+t^2}\,dt.
RN (x)=∫0x 1+t2(−1)N+1t2N+2 dt. Отсюда немедленно следует оценка
∣RN(x)∣≤∫0∣x∣t2N+21+t2 dt≤∣x∣2N+32N+3, |R_N(x)|\le\int_0^{|x|}\frac{t^{2N+2}}{1+t^2}\,dt\le\frac{|x|^{2N+3}}{2N+3},
∣RN (x)∣≤∫0∣x∣ 1+t2t2N+2 dt≤2N+3∣x∣2N+3 , и при аккуратном учёте знаменателя можно иногда получить чуть лучшую (но по порядку тот же) оценку для конкретного xxx.
3) Форма Лагранжа (обобщённая оценка). Для многочлена степени mmm Rm(x)=f(m+1)(ξ)(m+1)!xm+1для некоторого ξ∈(0,x). R_m(x)=\frac{f^{(m+1)}(\xi)}{(m+1)!}x^{m+1}\quad\text{для некоторого }\xi\in(0,x).
Rm (x)=(m+1)!f(m+1)(ξ) xm+1для некоторого ξ∈(0,x). Это даёт оценку, если можно явно ограничить ∣f(m+1)(t)∣|f^{(m+1)}(t)|∣f(m+1)(t)∣ на отрезке между 000 и xxx. Для arctan\arctanarctan производные имеют вид рациональных функций с знаменателем (1+t2)k(1+t^2)^k(1+t2)k, поэтому одна получает оценку вида
∣Rm(x)∣≤max∣t∣≤∣x∣∣f(m+1)(t)∣(m+1)!∣x∣m+1. |R_m(x)|\le \frac{\max_{|t|\le|x|}|f^{(m+1)}(t)|}{(m+1)!}|x|^{m+1}.
∣Rm (x)∣≤(m+1)!max∣t∣≤∣x∣ ∣f(m+1)(t)∣ ∣x∣m+1. Преимущество: применяется для любых xxx. Недостаток: при прямом грубом оценивании max∣f(m+1)∣\max|f^{(m+1)}|max∣f(m+1)∣ получается большая (иногда чрезмерно грубая) оценка по сравнению с методом Лейбинца на ∣x∣≤1|x|\le1∣x∣≤1.
4) Случай больших аргументов (∣x∣>1|x|>1∣x∣>1). Лучше сначала использовать тождество
arctanx=π2sgn(x)−arctan1x, \arctan x=\frac{\pi}{2}\operatorname{sgn}(x)-\arctan\frac{1}{x},
arctanx=2π sgn(x)−arctanx1 , и разложить arctan(1/x)\arctan(1/x)arctan(1/x) как ряд по малой переменной 1/x1/x1/x. Тогда оценка остатка даётся тем же правилом Лейбинца:
∣RN(x)∣≤∣1/x∣2N+32N+3, |R_N(x)|\le\frac{|1/x|^{2N+3}}{2N+3},
∣RN (x)∣≤2N+3∣1/x∣2N+3 , что даёт очень маленькую ошибку при больших ∣x∣|x|∣x∣.
5) Аппроксимации типа Паде, Чебышёва и минимакс. Для одновременной хорошей равномерной оценки на отрезке (особенно если нужен малый максимум ошибки на всём интервале) рациональные приближения (Паде) или полиномы с равномерной минимакс-ошибкой (Чебышёв) дают существенно лучшие результаты по сравнению с голым усечённым степенным рядом; для фиксированной степени они обычно дают намного меньшую максимальную ошибку. Эти методы сложнее в применении, но предпочтительны для практических вычислений с ограниченным числом членов на заданном отрезке.
Краткие рекомендации:
- Для точечной оценки при ∣x∣≤1|x|\le1∣x∣≤1 и простоты: используйте оценку Лейбинца (или интегральную форму) — она простая и почти оптимальная.
- Для ∣x∣>1|x|>1∣x∣>1: сначала замените x↦1/xx\mapsto1/xx↦1/x через тождество, затем примените ту же оценку.
- Для равномерной аппроксимации на отрезке: используйте Чебышёв- или Паде-аппроксимации — они дают более точные верхние оценки при том же числе параметров.