2.
La tabla presenta datos acerca de la resistencia (y) del papel Kraft y el
porcentaje de madera dura (x) en el lote de pulpa con el que se fabrica.
Concentración de madera Resistencia a la tensión
dura x (%) y (psi)
1 6.3
1.5 11.1
2 20
3 24
4 26.1
4.5 30
5 33.8
5.5 34
6 38.1
6.5 39.9
7 42
8 46.1
9 53.1
10 52
11 52.5
12 48
13 42.8
14 27.8
15 21.9
a) Trazar un diagrama de dispersión con los datos de la tabla. ¿Parece que
será un polinomio cuadrático?
b) Con base en el diagrama de dispersión y el conocimiento del proceso de
producción parecen indicar que un modelo cuadrático puede describir en
forma adecuada la relación entre la resistencia a la tensión y la
concentración de fibra corta (es decir de madera dura). Determine el
modelo cuadrático
y=a x2 +bx +c
por el método de mínimos cuadrados, para los datos de la tabla.
c) Graficar en el mismo papel el modelo ajustado de segundo grado, junto con
el diagrama de dispersión.
d) Emplee el polinomio cuadrático ajustado para estimar la resistencia a la
tensión del papel Kraft si la concentración de fibra estuviera en 8.5%
e) Determine el coeficiente de correlación de Karl Pearson.
Codificación:
%Ajuste cuadrático por mínimos cuadrados
clear all
clc
x=[1;1.5;2;3;4;4.5;5;5.5;6;6.5;7;8;9;10;11;12;13;14;15]
sum(x)
sum(x.^2)
sum(x.^3)
sum(x.^4)
y=[6.3;11.1;20;24;26.1;30;33.8;34;38.1;39.9;42;46.1;53.1;52;52.5;48;42.8;
27.8;21.9]
sum(y)
sum(x.*y)
sum(x.^2.*y)
plot(x,y,'db');xlabel('Concentración de madera dura'),ylabel('Resistencia
a la tensión'), title('Grafico de dispersion')
axis([0,16,0 60])
n=length(x)
m=[n sum(x) sum(x.^2) sum(y); sum(x) sum(x.^2) sum(x.^3) sum(x.*y);
sum(x.^2) sum(x.^3) sum(x.^4) sum(x.^2.*y)]
m(1,:)=(1/m(1,1))*m(1,:)
m(2,:)=-m(2,1)*m(1,:)+m(2,:)
m(3,:)=-m(3,1)*m(1,:)+m(3,:)
m(2,:)=(1/m(2,2))*m(2,:)
m(3,:)=-m(3,2)*m(2,:)+m(3,:)
m(3,:)=(1/m(3,3))*m(3,:)
m(2,:)=-m(2,3)*m(3,:)+m(2,:)
m(1,:)=-m(1,2)*m(2,:)+m(1,:)
m(1,:)=-m(1,3)*m(3,:)+m(1,:)
a=m(1,4)
b=m(2,4)
c=m(3,4)
syms f
syms g
g=a+b*x+c*(x.^2);
f=0:0.01:8.5;
yg=subs(g,f);
%Estimar la resistencia de tensión si la concentración de la fibra
%estuviera en 8.5%
z=a+b*8.5+c*(8.5^2)
figure
hold on
plot(f,z)
plot(x,y,'db');xlabel('Concentración de madera dura'),ylabel('Resistencia
a la tensión'), title('Grafico de dispersion')
axis([0 16, 0 60])
plot (x,g,'g')
%Coeficiente de correlación Karl Pearson
r=sum(x.*y)/sqrt(sum(x.^2)*sum(y.^2))
Ejecución:
Gráfica de dispersión:
Curva de ajuste:
Resultados:
0.000 2 0.0014 0.0134 0.0065
m=[1 .0e+05]∗0. 0014 0. 0134 0.1494 0.0531
0. 0134 0.1494 1.8143 0.5167
0.000 0 0.00 01 0.0 007 0.0003
m=[1.0e+05 ]∗ 0.0014 0.0134 0.1494 0.0531
0.0134 0.1494 1.8143 0.5167
0.0000 0.0001 0.0007 0.0003
m=[1.0e+05 ]∗ 0 0. 0 0 3 3 0.0524 0.0 05 9
0.0134 0.1494 1.8143 0.5167
0.000 1 0.000 7 0.007 0 0.003 4
m=[1.0e+0 4]∗ 0 0.0 3 33 0.52 39 0.05 89
0 0. 5239 8.7626 0. 6031
0.0001 0.0007 0.0070 0.0034
m=[1.0e+04 ]∗ 0 0.0 001 0.0016 0.0 002
0 0.5239 8.7626 0.6031
0.0001 0.007 3 0.070 3 0.034 2
m=[1.0e+0 3 ]∗ 0 0.001 0 0.01 57 0.00 18
0 0 8.1181 −3.2477
1 .000 0 7 .26 3 2 70.2632 34.1842
m= 0 1.0000 15.7482 1.7710
0 0 1.0000 −0.6345
1.0000 7.2632 70.2632 34.1842
m= 0 1.0000 0.0000 1 1.7 640
0 0 1.0000 −0.6345
1.0000 7.2632 0.0000 78.7696
m= 0 1.0000 0.0000 11.7640
0 0 1.0000 −0.6345
1.0000 0 0.0000 −6.6742
m= 0 1.0000 0.0000 11.7640
0 0 1.0000 −0.6345
a=−6.674 2
b=11.764 0
c=−0.634 5
y=−6.6742 x 2 +11.7640 x−0.654 5
Resistencia de tensión si la fibra estuviera en 8.5%
Z=47.4736
Coeficiente de correlación Karl Pearson
R=0.9074