NIST: Eckerle4 — Interferencia circular [Alta]

[ ]

Xv = [400,405,410,415,420,425,430,435,436.5,438,439.5,441,442.5,444,445.5,447,448.5,450,451.5,453,454.5,456,457.5,459,460.5,462,463.5,465,470,475,480,485,490,495,500]
Yv = [0.0001575,0.0001699,0.000235,0.0003102,0.0004917,0.000871,0.0017418,0.00464,0.0065895,0.0097302,0.0149002,0.023731,0.0401683,0.0712559,0.1264458,0.2073413,0.2902366,0.3445623,0.3698049,0.3668534,0.3106727,0.2078154,0.1164354,0.0616764,0.03372,0.0194023,0.0117831,0.0074357,0.0022732,0.00088,0.0004579,0.0002345,0.0001586,0.0001143,0.000071]
b0 = [1.5,5,450]

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

b2 = res[2]
b3 = res[3]
mostrar "b1 (cert 1.5544):"
mostrar res[1]
mostrar "b3 (cert 451.54):"
mostrar res[3]
mostrar "RSS (cert 1.464e-3):"
mostrar res[4]
n_graf = 35.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]
i = 1
para i de 1 hasta 35
    Y_fit[i] = (b1/res[2]) * exp(-0.5 * ((Xv[i]-res[3])/res[2]) * ((Xv[i]-res[3])/res[2]))
finpara
graficar_dispersion(Xv, Yv, "Datos observados")
graficar_xy(Xv, Y_fit, "Ajuste NIST")

Salida Jota: b1 (cert 1.5544): Salida Jota: 1.554 Salida Jota: b3 (cert 451.54): Salida Jota: 451.541 Salida Jota: RSS (cert 1.464e-3): Salida Jota: 0.001

Shift + Enter ejecuta la celda actual

Ayuda de Jota