agua

Esta es una vista previa del archivo. Inicie sesión para ver el archivo original

% Solucion a problema opcional 2 en Test 1, 521230, S1 2020
% Datos
x = 0:10:100;
densidad = [1000 1000 998 996 992 988 983 978 972 965 958];
visc = [1.788 1.307 1.003 0.799 0.657 0.548 0.467 0.405 0.355 0.316 0.283];%*1e-3;
kvisc = [1.788 1.307 1.005 0.802 0.662 0.555 0.475 0.414 0.365 0.327 0.295];%*1e-5
% aproximar densidad con spline c'ubica, pero tambi'en es v'alido usar
% polinomios de interpolaci'on.
pdensidad = polyfit(x,densidad,length(x)-1);
sdensidad = spline(x,densidad);
ppval(sdensidad,[34 68 86 91])
% para aproximar viscosidad y densidad cin'etica s'i deben utilizar spline
pvisc = polyfit(x,visc,length(x)-1);
svisc = spline(x,visc);
ppval(svisc,[34 68 86 91])
% una manera de determinar para qu'e valor de temperatura la densidad cin'etica
% es 1e-5 es con ayuda de la spline c'ubica que interpola la densidad cin'etica
% En flab4_ej2 se calcula y eval'ua f(x) = spline(x) - 1
% Dado que en 20 la densidad cin'etica es mayor que 1e-5 y en 30 es menor, el
% intervalo [20,30] sirve como intervalo inicial para bisecci'on
biseccion(@flab5_ej2,20,30,1e-8,100)
% Dado que en el video del laboratorio se utiliz'o fzero, los estudiantes tambi'en
% podr'ian utilizar este comando
skvisc = spline(x,kvisc); 
f = @(x) ppval(skvisc,x)-1;
fzero(f,20)
xx = linspace(x(1),x(end),10000);
plot(x,densidad,'ro')
hold on
plot(xx,ppval(sdensidad,xx))
plot(xx,polyval(pdensidad,xx),'r');
figure(2)
plot(x,visc,'ro')
hold on
plot(xx,ppval(svisc,xx))
plot(xx,polyval(pvisc,xx),'c')