NIST: Hahn1 — Viscosidad bronce [Alta]

[ ]

mostrar "Hahn1 — start cerca de cert para convergencia"
Xv = [24.41,34.82,44.09,45.07,54.98,65.51,70.53,75.70,89.57,91.14,96.40,97.19,114.26,120.25,127.08,133.55,151.84,168.10,204.16,212.13,257.10,261.34,328.65,342.62,362.17,382.44,400.37,429.04,452.50,472.44,513.13,524.01,539.70,566.10,599.25,626.75]
Yv = [130.21,131.48,132.92,133.59,133.14,131.84,130.96,130.01,127.11,126.82,125.51,125.26,121.68,120.12,118.42,116.83,112.48,108.44,99.39,97.28,83.27,79.69,66.46,64.15,61.38,58.40,55.75,51.30,46.28,44.42,40.37,38.89,37.15,34.38,31.09,28.43]
b0 = [1.077,0.03977,-0.0001,0.0000001,-0.000438,0.0000296,-0.000000571]

res = []
ajuste_no_lineal(Xv, Yv, 11, b0, res)
b1 = res[1]
b2 = res[2]
b3 = res[3]
b4 = res[4]
b5 = res[5]
b6 = res[6]
b7 = res[7]

b2 = res[2]
b3 = res[3]
b4 = res[4]
b5 = res[5]
b6 = res[6]
b7 = res[7]
mostrar "b1 (cert 1.0776):"
mostrar res[1]
mostrar "RSS (cert 1.5324):"
mostrar res[8]
n_graf = 36.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]
i = 1
para i de 1 hasta 36
    Y_fit[i] = (b1 + b2*Xv[i] + b3*Xv[i]*Xv[i] + b4*Xv[i]*Xv[i]*Xv[i]) / (1.0 + b5*Xv[i] + b6*Xv[i]*Xv[i] + b7*Xv[i]*Xv[i]*Xv[i])
finpara
graficar_dispersion(Xv, Yv, "Datos observados")
graficar_xy(Xv, Y_fit, "Ajuste NIST")

Salida Jota: Hahn1 — start cerca de cert para convergencia Salida Jota: b1 (cert 1.0776): Salida Jota: 103.605 Salida Jota: RSS (cert 1.5324): Salida Jota: 49.954

Shift + Enter ejecuta la celda actual

Ayuda de Jota