* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
510
Приближение функций и прогноз
Рис. 6.35. График полинома степени 40, иллюстрирующий возникновение числовой неустойчивости в правой части – при больших x
g1:=Plot[f[x],{x,0,10}] p[x_]=SetPrecision[Collect[InterpolatingPolynomial[data, x ],x],20] 1.0000000000000000000 + 0.50000000026808300024x - 2.1946079883927716025?10-9x2 0.062499992009412544403x3 - 1.7338956400433230831?10-8x4 + 0.0026041919664282057491x5 - 2.6520381936664922480?10-8x6 0.000054232658330779002056x7 - 1.2572612631091120302?10-8x8 + 6.8412544115568677809?10-7x9 + 2.2426515004772608460?10-9x10 4.97580739219249427092?10-9x11 - 1.6405048247087269161?10-10x12 + 6.5645938318937200218?10-11x13 - 5.0007652405250705692?10-12x14 + 4.6925914016242405225?10-13x15 - 5.9718776421275391973?10-14x16 + 4.8572376755578655131?10-15x17 - 2.2237845150851505265?10-16x18 + 5.4431086706688516340?10-18x19 - 5.6368826561019998150?10-20x20 g3:=Plot[p[x],{x,0,10}] Show[g3,g2]
В выводе этого фрагмента представлен полученный аппроксимирующий поли ном степени n=20. График полинома и точек функции f(x) представлен на рис. 6.36. От графика рис. 6.35 он отличается лишь большей вдвое густотой узловых точек.
Рис. 6.36. График полинома двадцатой степени