GUÍA : APROXIMACIÓN DE FUNCIONES POR

POLINOMIOS Y FUNCIONES RACIONALES.

Por : Alvaro H. Salas

INTRODUCCION.
Examinemos la siguiente gráfica :

7 x3
x-
60
In[13]:= PlotEvaluateSin[x], , {x, - 3 π / 4, 3 π / 4}, PlotStyle → {Red, {Dashed, Thick}}
x2
1+
20

1.0

0.5

Out[13]=
-2 -1 1 2

-0.5

-1.0

7 x3
x-
60
Según esto, sen (x) ≈ para x entre - 3 π / 4 y 3 π / 4
x2
1+
20

Esta aproximación viene dada como un cociente de dos polinomios,

es decir, constituye una aproximación racionals cerca de x = 0.
Esta aproximación es similar a las aproximaciones polinomiales. Por ejemplo, sabemos que
x3 x5 x7
sen (x) = x - + - + o (x)9 , de modo que
6 120 5040
x3 x5 x7
sen (x) ≈ x - + - .
6 120 5040

La gráfica de seno y su polinomio de McLaurin  polinomio de Taylor en x = 0  para - 3 π / 4 ≤

x ≤ 3 π / 4 se parecia a continuación :
2 APROXIMACIÓN DE FUNCIONES POR POLINOMIOS Y FUNCIONES RACIONALES.nb

x3 x5 x7
In[17]:= PlotEvaluateSin[x], x - + -,
6 120 5040
{x, - 3 π / 4, 3 π / 4}, PlotStyle → {Red, {Dashed, Thick}}

1.0

0.5

Out[17]=
-2 -1 1 2

-0.5

-1.0

Observamos que la aproximación polinomial es un tanto mejor que la racional. Sin embargo,
en otras situaciones ocurre que una aproximación racional es mucho mejor que una polinomial.
Para obtener aproximaciones polinomiales y racionales en Mathematica,
la librería << FunctionApproximations` eds la herramienta adecuada. Primeo se craga así :

In[49]:= << FunctionApproximations`(*Shift+Enter o simplemente Intro*)

El comando para aproximar es EconomizedRationalApproximation, el cual se usa así :

EconomizedRationalApproximation[f[x], {x, {a, b}, m, n}]

Esta inmstrucción nosd proporciona una aproximación racional económica para la función y =
a+b
f (x) en el intervalo [a, b]. El punto medio es el centro de la aproximación. Mathematica
2
P (x)
nos proporciona una aproximación  cuando existe en la forma ,
Q (x)
en donde P (x) es de grado m o menor , mientras que Q (x) es de grado n o
menor. Mientras más cerca esté x del centro del intervalo,
es decir, de a + b  2, mejor será la aproximación.
Otros formatos :

EconomizedRationalApproximation[f[x], {x, {a, b}, m, 0}] : Aproximación por un

polinomio de grado m o menor EN EL INTERVALO [A,B]

EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, m, n}] : Aproximación

racional económica alrededror de x = a.
EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, m, n}] //. ϵ -> 0 :
Aproximación racional económica muy cerca de x = a.
EJEMPLOS :
 APROXIMACIÓN DE FUNCIONES POR POLINOMIOS Y FUNCIONES RACIONALES.nb 3

In[23]:= EconomizedRationalApproximation[Exp[x] Sin[x], {x, {- ϵ, ϵ}, 3, 3}] //. ϵ → 0

x+ 3 x2 + 19 x3
Out[23]= 4 120
1- x + 3 x 2 + x3
4 40 120

3 x2 19 x3
x+ +
4 120
In[30]:= PlotEvaluateExp[x] Sin[x], ,
x 3 x2 x3
1- + +
4 40 120

{x, - 3, 3}, PlotStyle → {Red, {Dashed, Thick}}

12

10

Out[30]= 6

-3 -2 -1 1 2 3

In[86]:= era = EconomizedRationalApproximation[Exp[x] Sin[x], {x, {- ϵ, ϵ}, 3, 3}] //. ϵ → 2;

APROX ⩵ era
PlotEvaluate[{Exp[x] Sin[x], era}], x, - π  2, π  2, PlotStyle → {Red, {Dashed, Thick}}
6101 x + 1341 x2 + 67 x3
Out[87]= APROX ⩵ 8400 2800 840
2041 - x + 3 x2 + x3
2800 4 40 120

Out[88]=
2

-1.5 -1.0 -0.5 0.5 1.0 1.5

4 APROXIMACIÓN DE FUNCIONES POR POLINOMIOS Y FUNCIONES RACIONALES.nb

In[89]:= era = EconomizedRationalApproximation[Sin[x], {x, {- ϵ, ϵ}, 5, 5}] //. ϵ → 2;

Approx == era
Plot[Evaluate[{Sin[x], era}], {x, - 3, 3}, PlotStyle → {Red, {Dashed, Thick}}]

Out[90]= Approx ⩵ 99 912 227 x - 7 289 783 x3 + 3467 x5  13 321 631 + 59 327 x2 + 5 x4
114 307 200 62 868 960 1 247 400 15 240 960 1 995 840 11 088

1.0

0.5

Out[91]=
-3 -2 -1 1 2 3

-0.5

-1.0

In[92]:= era = EconomizedRationalApproximation[Sin[x], {x, {- ϵ, ϵ}, 7, 1}] //. ϵ → 0;

Approx == era
Plot[Evaluate[{Sin[x], era}], {x, - 2, 2}, PlotStyle → {Red, {Dashed, Thick}}]
3 5 7
Out[93]= Approx ⩵ x - x + x - x
6 120 5040

1.0

0.5

Out[94]=
-2 -1 1 2

-0.5

-1.0
 APROXIMACIÓN DE FUNCIONES POR POLINOMIOS Y FUNCIONES RACIONALES.nb 5

In[104]:= era = EconomizedRationalApproximation[Log[1 + x], {x, {- ϵ, ϵ}, 3, 3}] //. ϵ → 0;

Approx == era
Plot[Evaluate[{Log[1 + x], era}], {x, - 0.9, 2}, PlotStyle → {Red, {Dashed, Thick}}]
x + x2 + 11 x3
Out[105]= Approx ⩵ 60
1+ 3x + 3 x 2 + x3
2 5 20

1.0

0.5

-0.5 0.5 1.0 1.5 2.0

-0.5
Out[106]=

-1.0

-1.5

-2.0

FORMA SIMBÓLICA :

In[113]:= Clear[f];
poliap = era = EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 3, 0}] //. ϵ → 0 //
TraditionalForm
Out[114]//TraditionalForm=
1  3 1
f (a) (x - a)3 + (x - a)2 f ′′ (a) + (x - a) f ′ (a) + f (a)
6 2

Así, vemos que una aproximación polinomial de grado 3 MUY cerca de x =

a  o sea con ϵ = 0 es simplemente el polinomio de Taylor de grado 3 de y = f (x) cerca de x = a.

In[115]:= Clear[f];
poliap =
era = EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 3, 0}] // TraditionalForm
Out[116]//TraditionalForm=

1 4  4 1 1 ϵ2 f (a) f 4 (a) ϵ2 f 4 (a) (x - a) f ′ (a)

- ϵ f (a) + f 3 (a) (x - a)3 + (x - a)2 f ′′ (a) + (x - a) f ′ (a) - - + f (a) 
192 6 2 12 f ′′ (a) 12 f ′′ (a)
ϵ2 f 4 (a)
1-
12 f ′′ (a)

Obsérvese la diferencia que hay entre aproximación polinomial MUY cerca de x =

a ϵ = 0 y aproximación piolonomial cerca de x = a ϵ ≠ 0.
Lo mismo es válido apara aproximaciones racionales en general.
6 APROXIMACIÓN DE FUNCIONES POR POLINOMIOS Y FUNCIONES RACIONALES.nb

In[121]:= Clear[f];
poliap = era =
EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 1, 1}] //. ϵ → 0 // FullSimplify //
TraditionalForm
Out[122]//TraditionalForm=
2 (x - a) f ′ (a)2
+ f (a)
(a - x) f ′′ (a) + 2 f ′ (a)

In[123]:= Clear[f];
poliap =
era = EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 1, 1}] // FullSimplify //
TraditionalForm
Out[124]//TraditionalForm=

16 (x - a) f ′ (a)3   3 ϵ2 f ′′ (a)2 + 16 f ′ (a)2 + f ′ (a) 8 (a - x) f ′′ (a) - 2 ϵ2 f 3 (a) + f (a)

In[125]:= Clear[f];
poliap = era =
EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 2, 2}] //. ϵ → 0 // FullSimplify //
TraditionalForm
Out[126]//TraditionalForm=
2
 f (a) -4 f 3 (a) (a - x)2 - 36 f ′′ (a)2 + 3 f 4 (a) (a - x) ((a - x) f ′′ (a) + 2 f ′ (a)) + 12 f 3 (a) ((x - a) f ′′ (a) + 2 f ′ (a)) +
6 (a - x) 3 (x - a) f ′′ (a)3 +  f 4 (a) (x - a) - 4 f 3 (a) f ′ (a)2 + 2 f ′ (a) f ′′ (a) 2 f 3 (a) (a - x) + 3 f ′′ (a) 
2
-4 f 3 (a) (a - x)2 - 36 f ′′ (a)2 + 3 (a - x)  f 4 (a) (a - x) - 4 f 3 (a) f ′′ (a) + 6  f 4 (a) (a - x) + 4 f 3 (a) f ′ (a)

In[127]:= Clear[f];
poliap = era =
EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 1, 2}] //. ϵ → 0 // FullSimplify //
TraditionalForm
Out[128]//TraditionalForm=

2 f 3 (a) f (a)2 (x - a) + 12 f ′ (a)2 ((x - a) f ′ (a) + f (a)) - 6 f (a) f ′′ (a) (2 (x - a) f ′ (a) + f (a)) 
3 (a - x)2 f ′′ (a)2 + 12 f ′ (a)2 - 2 f 3 (a) (a - x) ((a - x) f ′ (a) + f (a)) - 6 f ′′ (a) ((x - a) f ′ (a) + f (a))

In[129]:= Clear[f];
poliap = era =
EconomizedRationalApproximation[f[x], {x, {a - ϵ, a + ϵ}, 2, 1}] //. ϵ → 0 // FullSimplify //
TraditionalForm
Out[130]//TraditionalForm=
3 (a - x)2 f ′′ (a)2
x f ′ (a) - a f ′ (a) + + f (a)
2 f 3 (a) (a - x) + 6 f ′′ (a)

NEXT : APROXIMACIÓN RACIONAL Y POLINÓMICA DE LA

SOLUCIOÓN DE UNA EDO CERCA DEL PUNTO INICIAL. : TO BE CONTINUED ... ..