Ejemplo: uso de minerr para el ajuste no lineal de mínimos cuadrados
La función
minerr es similar a la función
find, excepto en que devuelve una solución aproximada en algunos casos en los que find indica que no existe ninguna solución.
1. Defina dos vectores.
2. Defina una función de ajuste (densidad Weibull con parámetros desconocidos).
3. Defina los valores de prueba inicial de los dos parámetros.
4. Utilice una ecuación para minimizar dentro de un bloque de resolución.
5. Añada un bloque de resolución y utilice minerr para resolver el problema. La función minerr utiliza el método de Levenberg-Marquardt para minimizar este problema. El método de Levenberg-Marquardt realiza sus propias operaciones de suma y cuadratura de los residuales.
Los parámetros del mejor ajuste son los valores calculados:
La función find no encuentra una solución para el problema anterior.
6. Calcule la suma de cuadrados que este método ha minimizado implícitamente.
7. Trace el mejor ajuste de Weibull con respecto a los datos x-y.
8. Evalúe el error cuadrático medio. Si la media es cero, existe una solución real:
Se puede minimizar directamente utilizando la ecuación SSE y la función minimize.