def I_ds_mod_skr(k):
return d / np.pi * np.cos(k * g / 2) * np.sin(
k * d / 2) / (k * d / 2) * np.cos(ds_y * k)
def I_es_mod_skr(k):
return d / np.pi * np.sin(k * d / 2) / (k * d / 2) * np.cos(es_y * k)
def ds_breite_nst(b):
alpha = np.arctan2(max1[0], schirmabstand)
x_skr = np.pi * b * np.sin(alpha) / laser
y = 2. * np.pi * ds_abstand * np.sin(alpha) / laser
return (np.sin(x_skr) / x_skr)**2 * (np.sin(y) / y)**2 - max1[1]
def I_ds_theo(x):
alpha = np.arctan2(x, schirmabstand)
x_skr = np.pi * ds_breite * np.sin(alpha) / laser
y = 2. * np.pi * ds_abstand * np.sin(alpha) / laser
return (np.sin(x_skr) / x_skr)**2 * (np.sin(y) / y)**2
def I_es_theo(x):
alpha = np.arctan2(x, schirmabstand)
x_skr = np.pi * es_breite * np.sin(alpha) / laser
return (np.sin(x_skr) / x_skr)**2