NIST: Chwirut2 — Calibración ultrasónica [Baja]

[ ]

Xv = [0.5,1,1.75,3.75,5.75,0.875,2.25,3.25,5.25,0.75,1.75,2.75,4.75,0.625,1.25,2.25,4.25,0.5,3,0.75,3,1.5,6,3,6,1.5,3,0.5,2,4,0.75,2,5,0.75,2.25,3.75,5.75,3,0.75,2.5,4,0.75,2.5,4,0.75,2.5,4,0.5,6,3,0.5,2.75,0.5,1.75]
Yv = [92.9,57.1,31.05,11.5875,8.025,63.6,21.4,14.25,8.475,63.8,26.8,16.4625,7.125,67.3,41,21.15,8.175,81.5,13.12,59.9,14.62,32.9,5.44,12.56,5.44,32,13.95,75.8,20,10.42,59.5,21.67,8.55,62,20.2,7.76,3.75,11.81,54.7,23.7,11.55,61.3,17.7,8.74,59.2,16.3,8.62,81,4.87,14.62,81.7,17.17,81.3,28.9]
b0 = [0.1, 0.01, 0.02]

res = []
ajuste_no_lineal(Xv, Yv, 3, b0, res)
b1 = res[1]
b2 = res[2]
b3 = res[3]

b2 = res[2]
b3 = res[3]
mostrar "b1 (cert 0.1666):"
mostrar res[1]
mostrar "RSS (cert 513.0):"
mostrar res[4]
n_graf = 54.0
Y_fit = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
i = 1
para i de 1 hasta 54
    Y_fit[i] = exp(-b1 * Xv[i]) / (b2 + b3 * Xv[i])
finpara
graficar_dispersion(Xv, Yv, "Datos observados")
graficar_xy(Xv, Y_fit, "Ajuste NIST")

Salida Jota: b1 (cert 0.1666): Salida Jota: 0.167 Salida Jota: RSS (cert 513.0): Salida Jota: 513.048

Shift + Enter ejecuta la celda actual

Ayuda de Jota