Instituto Politécnico Nacional

Unidad Profesional Interdisciplinaria de Biotecnología

MÉTODOS NUMÉRICOS (TALLER)

PROBLEMARIO
Grupo 4LM5

Integrantes
Arenas Villalba Daniela
Benítez Martínez Yaremi Arlet
Del Angel Hernandez Dominga Isabel
Morales Hernández Iris Jeanette
Rodríguez Alaníz Cristian Alejandro
Ruiz Martínez Aylin Brenda

Profesores
Dra. Alvarado Palacios Quetzaliztli
Rios Albarran Sergio

Fecha de entrega: 23/10/2018

INSTRUCCIONES: Resuelva detalladamente cada uno de los siguientes ejercicios contestando con
claridad las preguntas que se le hagan. Escriba el sistema de ecuaciones que se deba resolver y
aplique el método de solución que Ud. guste para encontrar los coeficientes de la curva de ajuste.
Calcule el coeficiente de determinación en cada caso.¡En todos los ejercicios grafique el (los) modelo(s)
de ajuste y los puntos etiquetándolos!
Ejercicio 1

La temperatura de ebullición del agua, Te, registrada a diferentes alturas se muestra en la

h(pies) 0 2000 5000 7500 10 000 20 000 26 000

Ajuste un modelo de regresión lineal a los datos. Use la ecuación para estimar la temperatura de
ebullición a 16,000 pies.

Código:
clc;clear all;close all
h=[0 2000 5000 7500 10000 20000 26000]
T=[212 210 203 198 194 178 168]
A=[length(h) sum(h)
sum(h) sum(h.^2)]
b=[sum(T)
sum(T.*h)]
x=inv(A)*b
syms l
p=x(2)*l+x(1)
plot(h,T,'mo','Linewidth',3,'Markersize',2);grid on;hold on
ezplot(p,[min(h),max(h)])
polyval(x,16000)

Comando:
La temperatura de ebullición a 16,000 pies es de 3.3900e+06

Grafico:

Ejercicio 2.

Los siguientes datos experimentales representan el calor especifico del agua como función de la
temperatura
T (K) 273 283 293 303 313 323 333 343 353 363
Cp(cal/gK) 1.00783 1.00129 0.99880 0.99800 0.99800 0.99850 0.99940 1.00060 1.00229 1.00437

Realice los ajustes polinomiales de 2o. y 4o. grado. ¿Qué diferencia existe entre uno y otro ajuste?,
¿es significativa (observe los coeficientes) ?, ¿cuál ajuste recomendaría?

Estime, con los dos modelos, el calor especifico del agua a 20°C.

Segundo orden.
P2(x)= ao+ a1x +a2x2

273= ao+1.00783 a1+1.015721309 a2

283= ao+1.00129 a1+1.002581664 a2

293= ao+0.99880 a1+0.99760144 a2

303= ao+0.99800 a1+0.996004 a2

313= ao+0.99800 a1+0.996004 a2

323= ao+0.99850 a1+0.99700225 a2

333= ao+0.99940a1+0.99880036 a2

343= ao+1.00060 a1+1.00120036 a2

353= ao+1.00229 a1+1.004585244 a2

363= ao+1.00437 a1+1.008759097 a2

Codigo:
clc;clear all;close all
T=[273 283 293 303 313 323 333 343 353 363]
Cp=[1.00783 1.00129 0.99880 0.99800 0.99800 0.99850 0.99940 1.00060
1.00229 1.00437]
A=[length(T) sum(T) sum(T.^2)
sum(T) sum(T.^2) sum(T.^3)
sum(T.^2) sum(T.^3) sum(T.^4)]
b=[sum(Cp)
sum(Cp.*T)
sum((Cp.^2).*T)]
x=inv(A)*b
syms l
p=x(3)*l^2+x(2)*l+x(1)
plot(T,Cp,'bo','Linewidth',3,'Markersize',2);grid on;hold on
ezplot(p,[min(T),max(T)])

código Ajuste.:
clc;clear all;close all
T=[273 283 293 303 313 323 333 343 353 363]
Cp=[1.00783 1.00129 0.99880 0.99800 0.99800 0.99850 0.99940 1.00060
1.00229 1.00437]
x=T;y=log(Cp)
A=[length(T) sum(T) sum(T.^2)
sum(T) sum(T.^2) sum(T.^3)
sum(T.^2) sum(T.^3) sum(T.^4)]
b=[sum(Cp)
sum(Cp.*T)
sum((Cp.^2).*T)]
X=inv(A)*b
alpha=exp(x(1));
beta=x(2);
gama=x(3);
syms V
D=vpa(alpha*V^beta,3)
plot(x, y, 'or', 'LineWidth', 4, 'Markersize', 4); grid on;
hold on
xlabel('T'); ylabel('Cp')
St=sum((T-mean(T)).^2);
Sr=sum((T-(X(1)+X(2)*Cp)).^2)
r=sqrt((St-Sr)/St)
polyval(X,20)

Comando:

El Cp del agua a 20°C es de -7.7262e+06 para el ajuste de polinomio de grado 2

Grafico:
Cuarto orden.

P2(x)= ao+ a1x +a2x2

273= ao+1.00783 a1+1.015721309 a2

283= ao+1.00129 a1+1.002581664 a2

293= ao+0.99880 a1+0.99760144 a2

303= ao+0.99800 a1+0.996004 a2

313= ao+0.99800 a1+0.996004 a2

323= ao+0.99850 a1+0.99700225 a2

333= ao+0.99940a1+0.99880036 a2

343= ao+1.00060 a1+1.00120036 a2

353= ao+1.00229 a1+1.004585244 a2

363= ao+1.00437 a1+1.008759097 a2

P2(x)= ao+ a1x +a2x2

273= ao+1.00783 a1+1.015721309 a2

283= ao+1.00129 a1+1.002581664 a2

293= ao+0.99880 a1+0.99760144 a2

303= ao+0.99800 a1+0.996004 a2

313= ao+0.99800 a1+0.996004 a2

333= ao+0.99940a1+0.99880036 a2

343= ao+1.00060 a1+1.00120036 a2

353= ao+1.00229 a1+1.004585244 a2

363= ao+1.00437 a1+1.008759097 a2

Codigo:
clc;clear all;close all
T=[273 283 293 303 313 323 333 343 353 363]
Cp=[1.00783 1.00129 0.99880 0.99800 0.99800 0.99850 0.99940 1.00060
1.00229 1.00437]
A=[length(T) sum(T) sum(T.^2) sum(T.^3) sum(T.^4)
sum(T) sum(T.^2) sum(T.^3) sum(T.^4) sum(T.^5)
sum(T.^2) sum(T.^3) sum(T.^4) sum(T.^5) sum(T.^6)
sum(T.^3) sum(T.^4) sum(T.^5) sum(T.^6) sum(T.^7)
sum(T.^4) sum(T.^5) sum(T.^6) sum(T.^7) sum(T.^8)]
b=[sum(Cp)
sum(Cp.*T)
sum((Cp.^2).*T)
sum((Cp.^3).*T)
sum((Cp.^4).*T)]
x=inv(A)*b
syms l
p=x(5)*l^4+x(4)*l^3+x(3)*l^2+x(2)*l+x(1)
plot(T,Cp,'bo','Linewidth',3,'Markersize',2);grid on;hold on
ezplot(p,[min(T),max(T)])
polyval(x,20)
xlabel('T'); ylabel('Cp')

Código ajuste:
clc;clear all;close all
T=[273 283 293 303 313 323 333 343 353 363]
Cp=[1.00783 1.00129 0.99880 0.99800 0.99800 0.99850 0.99940 1.00060
1.00229 1.00437]
x=T;y=log(Cp)
A=[length(T) sum(T) sum(T.^2) sum(T.^3) sum(T.^4)
sum(T) sum(T.^2) sum(T.^3) sum(T.^4) sum(T.^5)
sum(T.^2) sum(T.^3) sum(T.^4) sum(T.^5) sum(T.^6)
sum(T.^3) sum(T.^4) sum(T.^5) sum(T.^6) sum(T.^7)
sum(T.^4) sum(T.^5) sum(T.^6) sum(T.^7) sum(T.^8)]
b=[sum(Cp)
sum(Cp.*T)
sum((Cp.^2).*T)
sum((Cp.^3).*T)
sum((Cp.^4).*T)]
X=inv(A)*b
alpha=exp(x(1));
beta=x(2);
gama=x(3);
syms V
D=vpa(alpha*V^beta,3)
plot(x,y,'bo','Linewidth',4,'Markersize',4);grid on;hold on
st=sum((T-mean(T)).^2)
sr=sum((T-(X(1)+X(2)*Cp)).^2)
r=sqrt((st-sr)/st)
polyval(X,20)
xlabel('T'); ylabel('Cp')

Comando:

El Cp del agua a 20°C es de -2.9577e+14

Grafico:
Se recomienda ajuste de cuarto grado.

Ejercicio 3.

En un artículo de J. Agricultural Eng. Research, 1975 (pp. 353-361) se reportan los siguientes
datos con el numero días después de la floración , d, y el rendimiento de la cosecha
en kg/ha, R:
d(dı́as) 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46
R(kg/ha) 2508 2518 3000 3423 3507 3400 3500 3883 3823 3646 3708 3333 3517 3241 3103 2776

Ajuste el mejor modelo de regresión y estime el rendimiento de la cosecha 50 días

Codigo:
clc;clear all;close all
d=[16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46]
R=[2508 2518 3000 3423 3507 3400 3500 3883 3823 3646 3708 3333 3517 3241
3103 2776]
A=[length(d) sum(d) sum(d.^2)
sum(d) sum(d.^2) sum(d.^3)
sum(d.^2) sum(d.^3) sum(d.^4)]
b=[sum(R)
sum(R.*d)
sum((R.^2).*d)]
x=inv(A)*b
syms l
p=x(3)*l^2+x(2)*l+x(1)
plot(d,R,'or','LineWidth',3,'Markersize',2);grid on;
hold on
ezplot(p,[min(d),max(d)])
st=sum((d-mean(d)).^2)
sr=sum((d-(x(1)+x(2)*R)).^2)
r=sqrt((st-sr)/st)
polyval(x,50)

Comando:

el rendimiento de la cosecha 50 días después de la floración es de 1.3253e+11