示例#1
0
计算 phi_c ,b_c ,Phi_c, flux_c
'''
ssn_in = range_info.loc[:, 'ssn_delay'].values

plt.figure(figsize=(8, 4))
plt.plot(range(len(delay_ssn)),
         delay_ssn,
         label='Data',
         c='b',
         lw=2,
         zorder=12)
plt.plot(range(len(delay_ssn)), ssn_in, label='Smooth', c='r', lw=2)
plt.legend()
plt.show()

phi_c = FUNC.ssn_phi(ssn_in)
b_c = FUNC.ssn_b(ssn_in)
Phi_c = FUNC.fit_obj_vary(0.274, phi_c, b_c)

flux_c = np.array([])
for i in range(len(Phi_c)):
    flux_sig = FUNC.ffm_fun(0.274, Phi_c[i])
    flux_c = np.append(flux_c, flux_sig)
# print(range_info)
range_info['phi_c'] = phi_c
range_info['b_c'] = b_c
range_info['Phi_c'] = Phi_c
range_info['flux_c'] = flux_c

range_info.to_csv(r'./output/info_test.csv')