Octubre 19, (más) ejercicios de laboratorio

En el laboratorio virtual de hoy quiero que trabajen con el comando

 polyfit

disponible en Octave y Matlab. De acuerdo a la documentación, polyfit da los valores de los coeficientes de un polinomio P(X) de grado N que minimiza el error por mínimos cuadrados del ajuste. La función se emplea con la instrucción

  P = polyfit (X, Y, N)

donde (X,Y) son los datos que se desean ajustar y N el grado del polinomio. Usaremos esta función para ajustar líneas rectas con el siguiente ejercicio:

x = -5:5;
y = [0,0,0,1,1,1,0,0,0,0,0];
plot(x,y,'*k');
set(gca, 'ylim', [-1,5])

Continue reading →

Lab: Básicos de Octave

En este laboratorio haremos algunos ejercicios para aprender a usar Octave de manera que puedan implementar los algoritmos que discutimos de la mejor manera posible. Recuerden que, como las matemáticas y la física, aprender a programar requiere de tiempo y paciencia, pero sobre todo que hagan los ejercicios y “metan las manos”. No sirve de nada ver como hace las cosas el compañero, hay que hacerlo uno mismo. Los ejercicios que veremos hoy no pretenden ser una guía completa de Octave sino una breve introducción a comandos importantes y útiles para que puedan superar la curva de aprendizaje. Cualquier pregunta no duden en acercarse.

Lab: Introducción Método de Euler

En el laboratorio de hoy discutimos el método de Euler para resolver ecuaciones diferenciales lineales de primer grado y una implementación en Octave/MATLAB. El método de Euler con diferencias finitas es

$y_{i+1}=y_i + \frac{dy_i}{dt} \Delta t$

y su interpretación es muy sencilla: el nuevo valor de la variable $y$ al tiempo $t+dt$ se obtiene del valor de la variable al tiempo anterior $y_i$ incrementada por la tasa de cambio de la variable multiplicada por el incremento del tiempo $\frac{dy}{dt}\Delta t$ . La tasa de cambio de la variable generalmente se obtiene del problema matemático, como en el caso estudiado en clase, donde Continue reading →