UNIVERSIDAD TÉCNICA DE MANABÍ

INSTITUTO DE CIENCIAS BÁSICAS

DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICAS

NOMBRE: ANDREA CORRALES GARCÍA DOCENTE: LIC. FELIPE RUMBAUT

ASIGNATURA: MÉTODOS NUMÉRICOS CA.: OCTUBRE 2016-FEBRERO 2017

NIVEL: CUARTO FECHA: 24/01/2017

PARCIAL: SEGUNDO DEBER: 13

 Dado el conjunto de datos, 𝒙 = [−𝟑 − 𝟐 𝟎 𝟒 𝟓 𝟔]; 𝒚 = [𝟓 𝟒 𝟐 𝟕 𝟔 𝟑]

a) Obtenga f(x) para cada valor, f(-2.5); f(0.6); f(2); f(3); f(5.12)
b) Realice la gráfica de la interpolación lineal
c) Obtenga el polinomio de interpolación de los datos utilizando interpolación de potencia
ayudado por Matlab
d) Obtenga los resultados para el inciso a utilizando el polinomio de interpolación encontrado.
e) Grafique el resultado con el polinomio encontrado

>> % INTERPOLACIÓN LINEAL

7
>> x=[-3 -2 0 4 5 6]; y=[5 4 2 7 6 3];
>> x1=-3:0.2:6; y1=interp1(x,y,x1); 6

>> interp1(x,y,[-2.5,0.6,2,3,5.12])
5

ans =
4
4.5000 2.7500 4.5000 5.7500 5.6400
>> plot(x,y,x1,y1,'o');grid on 3

>> hold on
2

>> % INTERPOLACIÓN POR POTENCIAS

1
>> n=length(x); -3 -2 -1 0 1 2 3 4 5 6

>> a=polyfit(x,y,5)
a=
0.0027 -0.0324 -0.0200 0.6941 0.1651 2.0000
>> y2=polyval(a,x1);
>> plot(x1,y2,'or');grid on
>> polyval(a,[-2.5,0.6,2,3,5.12])
ans =
4.7042 2.3406 4.5159 6.2440 5.7426
 UNIVERSIDAD TÉCNICA DE MANABÍ
INSTITUTO DE CIENCIAS BÁSICAS
DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICAS

NOMBRE: ANDREA CORRALES GARCÍA DOCENTE: LIC. FELIPE RUMBAUT

ASIGNATURA: MÉTODOS NUMÉRICOS CA.: OCTUBRE 2016-FEBRERO 2017

NIVEL: CUARTO FECHA: 24/01/2017

PARCIAL: SEGUNDO DEBER: 14

 Dados los valores obtenidos en una medición determine un polinomio de interpolación por
el método de Lagrange. Evalúe para x=0.27. Represéntelo.

x -1.5 2.4 3.18 4.21

F(x) -12 -5.24 -2.37 1.12

>> p=collect(sym('(x-2.4)*(x-3.18)*(x-4.21)*-12/((-1.5-2.4)*(-1.5-3.18)*(-1.5-4.21))+(x+1.5)*(x-3.18)*(x-
4.21)*-5.24/((2.4+1.5)*(2.4-3.18)*(2.4-4.21))+(x+1.5)*(x-2.4)*(x-4.21)*-2.37/((3.18+1.5)*(3.18-
2.4)*(3.18-4.21))+(x+1.5)*(x-2.4)*(x-3.18)*1.12/((4.21+1.5)*(4.21-2.4)*(4.21-3.18))'),'x')

p=

- 0.1009972640199380662676673835642*x^3 + 0.82791367612310930905715853309445*x^2 +
1.4336089591504618274223038150664*x - 12.053258098618594177898528396392

>> a=[-
0.1009972640199380662676673835642,0.82791367612310930905715853309445,1.433608959150
4618274223038150664,-12.053258098618594177898528396392];

>> polyval(a,0.27)

ans =

-11.6078 -...- 12.053258098618594177898528396392

2
>> ezplot('-
0.1009972640199380662676673835642*x^3 + 0
0.82791367612310930905715853309445*x^2 +
1.4336089591504618274223038150664*x - -2

12.053258098618594177898528396392',[-
-4
1.5,4.21]);grid on
-6

-8

-10

-12

-14
-1 0 1 2 3 4
x
 UNIVERSIDAD TÉCNICA DE MANABÍ
INSTITUTO DE CIENCIAS BÁSICAS
DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICAS

NOMBRE: ANDREA CORRALES GARCÍA DOCENTE: LIC. FELIPE RUMBAUT

ASIGNATURA: MÉTODOS NUMÉRICOS CA.: OCTUBRE 2016-FEBRERO 2017

NIVEL: CUARTO FECHA: 31/01/2017

PARCIAL: SEGUNDO DEBER: 15

Dados los datos de la tabla. Encuentre un polinomio interpolador por el método de las diferencias
divididas de Newton. Obtenga el valor de f(2.2) y realice el gráfico del polinomio.
x 1 1.5 2 2.5
f(x) -14 -2 6 15

x f(x) 1 dif div 2 dif div 3 dif div

1 -14 24 -8 6,66666667
1,5 -2 16 2
2 6 18
2,5 15

MATLAB
>> p=collect(sym('-14+24*(x-1)-8*(x-1)*(x-1.5)+6.66666667*(x-1)*(x-1.5)*(x-2)'),'x')
p=
6.66666667*x^3 - 38.000000015*x^2 + 87.333333355*x - 70.00000001
>> a=[6.66666667,-38.000000015,87.333333355,-70.00000001];
>> polyval(a,2.2)
ans =
9.2000
>> ezplot('6.66666667*x^3 - 38.000000015*x^2 + 87.333333355*x - 70.00000001',[1,2.5]);grid on
 6.66666667 x 3 - 38.000000015 x 2 +...- 70.00000001

15

X: 2.199
Y: 9.177
10

-5

-10

-15

1 1.5 2 2.5
x
 UNIVERSIDAD TÉCNICA DE MANABÍ
INSTITUTO DE CIENCIAS BÁSICAS
DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICAS

NOMBRE: ANDREA CORRALES GARCÍA DOCENTE: LIC. FELIPE RUMBAUT

ASIGNATURA: MÉTODOS NUMÉRICOS CA.: OCTUBRE 2016-FEBRERO 2017

NIVEL: CUARTO FECHA: 31/01/2017

PARCIAL: SEGUNDO DEBER: 16

 Utilice la interpolación inversa para obtener el valor de x para el cual f(x)= 7.94. Utilizando
el método que elija para obtener un polinomio interpolador de tercer grado.
X 0 1 2.5 3 4.5 5 6
Y 2 5.4375 7.3516 7.5625 8.4453 9.1875 12

MÉTODO DE LAGRANGE

>> p=collect(sym('(x-3)*(x-4.5)*(x-5)*7.3516/((2.5-3)*(2.5-4.5)*(2.5-5))+(x-2.5)*(x-4.5)*(x-5)*7.5625/((3-
2.5)*(3-4.5)*(3-5))+(x-2.5)*(x-3)*(x-5)*8.4453/((4.5-2.5)*(4.5-3)*(4.5-5))+(x-2.5)*(x-3)*(x-4.5)*9.1875/((5-
2.5)*(5-3)*(5-4.5))'),'x')

p=

0.14582666666666666666666666666667*x^3 - 1.3749*x^2 + 4.6661933333333333333333333333333*x

+ 2.0007

>> x=solve('0.14582666666666666666666666666667*x^3 - 1.3749*x^2 +

4.6661933333333333333333333333333*x + 2.0007=7.94')

x=

3.9051575510541471755717623425169

1.6742426662094151173988647261885*i + 2.7615795859980247024216711739915

2.7615795859980247024216711739915 - 1.6742426662094151173988647261885*i

RESPUESTA:

La solución es x= 3.9051575510541471755717623425169
 UNIVERSIDAD TÉCNICA DE MANABÍ
INSTITUTO DE CIENCIAS BÁSICAS
DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICAS

NOMBRE: ANDREA CORRALES GARCÍA DOCENTE: LIC. FELIPE RUMBAUT

ASIGNATURA: MÉTODOS NUMÉRICOS CA.: OCTUBRE 2016-FEBRERO 2017

NIVEL: CUARTO FECHA: 07/02/2017

PARCIAL: SEGUNDO DEBER: 17

 Realizar un ajuste potencial al modelo siguiente.

 Realizar un ajuste parabólico (polinomio de grado 2).
 Realizar ajuste lineal.
 Decida cual modelo utilizar de acuerdo al método de los mínimos cuadrados y al
coeficiente de determinación.

x 0.75 2 3 4 6 8 8.5
y 1.2 1.95 2 2.4 2.5 2.7 2.6

x y modelo lineal modelo polinómico modelo potencial

0,75 1,2 1,58565 0,14872592 1,2909875 0,00827873 1,29926958 0,00985445
2 1,95 1,7824 0,02808976 1,7748 0,03069504 1,77170163 0,03179031
3 2 1,9398 0,00362404 2,0867 0,00751689 2,01405182 0,00019745
4 2,4 2,0972 0,09168784 2,3318 0,00465124 2,20585149 0,03769365
6 2,5 2,412 0,007744 2,6216 0,01478656 2,50758882 5,759E-05
8 2,7 2,7268 0,00071824 2,6442 0,00311364 2,74638839 0,00215188
8,5 2,6 2,8055 0,04223025 2,6081 6,561E-05 2,79954314 0,03981747
0,32282005 0,06910771 0,1215628

r 0,896325834
0,896325833611862
El mejor modelo es el método polinómico debido a que
y = 0.1574x + 1.4676 y = -0.0334x2 + 0.4789x + 0.9506 y = 1.423x0.3162
R² = 0.8034 R² = 0.9579 R² = 0.9404
3

2.5

2 y

1.5

0.5

0
0 1 2 3 4 5 6 7 8 9

y Linear (y) Poly. (y) Power (y) Linear (y)

 UNIVERSIDAD TÉCNICA DE MANABÍ
INSTITUTO DE CIENCIAS BÁSICAS
DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICAS

NOMBRE: ANDREA CORRALES GARCÍA DOCENTE: LIC. FELIPE RUMBAUT

ASIGNATURA: MÉTODOS NUMÉRICOS CA.: OCTUBRE 2016-FEBRERO 2017

NIVEL: CUARTO FECHA: 14/02/2017

PARCIAL: SEGUNDO DEBER: 18

La corriente eléctrica que circula por un conductor varía con el tiempo según la expresión
𝑰(𝒕) = 𝟐𝟎𝒆𝟓/(𝒕+𝟏) . Determine el valor aproximado de la carga que circula por el conductor en
el intervalo de 0 a 20 segundos. Divida el tiempo en intervalos de 1 segundos, utilice una hoja
de cálculo. En intervalos de 0.01 s, utilice Matlab. Encuentre el valor real de la integral.
MÉTODO DE LOS TRAPECIOS
t I A trapecios
0 2968,26318 1605,95653
1 243,649879 174,76984
2 105,889801 87,8483301
3 69,8068591 62,0862479
4 54,3656366 50,1925772
5 46,0195178 43,4370296
6 40,8545414 39,1097303
7 37,3649191 36,1115496
8 34,85818 33,9163027
9 32,9744254 32,2417837
10 31,5091421 30,923539
11 30,3379359 29,8594599
12 29,3809839 28,9828923
13 28,5848006 28,2485246
14 27,9122485 27,6245037
15 27,3367588 27,0877971
16 26,8388354 26,6213456
17 26,4038558 26,2122492
18 26,0206426 25,8505755
19 25,6805083 25,5285544
20 25,3766006

2442,60936
Matlab
>> t=0:0.01:20;
>> i=20*exp(5./(t+1));
>> I=trapz(t,i)
I=
1.6438e+003
>> syms t
>> i=20*exp(5/(t+1));
>> I=double(int(i,0,20))
I=
1.6437e+003