1 Aplicaciones de interpolaciónProgramación Numérica

2 Caso 1: Modelo de poblaciónEl crecimiento de población de bacterias puede modelarse mediante Donde p es la población y k es la velocidad de crecimiento específico. La solución de esta ecuación diferencial es p(t) = p0 ekt Obviamente k no puede ser constante. Un modelo de k es suponerlo un modelo de crecimiento de saturación kmax = velocidad de crecimiento máxima K = constante de saturación media

3 Resolviendo en MatLab se obtiene: kmax = 1.230431 dia–1 Evaluar kmax y K para los siguientes datos usando ajuste de mínimos cuadrados de 1/k y 1/f. f (mg/L) k (dia–1) o Solución: Resolviendo en MatLab se obtiene: kmax = dia–1 K = mg/L

4 Gráfico del ajuste

5 Trazadores para transferencia de calorLos lagos de zona templada se dividen en estratos térmicos durante el verano: epilimnion y hipolimnion, separadas por un plano llamado termoclina. La termoclina se caracteriza por ser el punto donde la curva de temperatura tiene un punto de inflexión d2T/dt = 0 y la primera derivada tiene un máximo. Utilizar trazadores cúbicos para determinar la profundidad de la termoclina para los siguientes datos: T (ºC) z (m)

6 El siguiente guión calcula los trazadores y los grafica en pasos uniformes. Calcula también la primera y segundas derivadas y las grafica. function termoclina2 prof = [ ]; temp = [ ]; %calcula trazadores cúbicos pt = spline3(prof,temp); [r c] = size(pt); x = []; y = []; y1 = []; y2 = [];

7 %evalua los trazadores desde 0 a 28 m de profundidad en pasos uniformesfor p=0:28 for j=1:r if p>=prof(j) & p k = j; end x = [x p]; pol = pt(k,1)*p^3+pt(k,2)*p^2+pt(k,3)*p +pt(k,4); dpol = 3*pt(k,1)*p^2+2*pt(k,2)*p+pt(k,3); d2pol = 6*pt(k,1)*p+2*pt(k,2); y = [y pol]; y1 = [y1 dpol]; y2 = [y2 d2pol];

8 subplot(3,1,1) plot(x,y) subplot(3,1,2) plot(x,y1) subplot(3,1,3) plot(x,y2) De las gráficas puede verse que el valor de profundidad donde la segunda derivada es cero o la primera derivada es máxima es m con un gradiente de temperatura de -1.61ºC/m, aproximadamente.