我在计算和绘制行星围绕脉冲星运行的光谱能量。在
我以前对列表中的所有数据进行了排序,其维度为[172,2](172行,2列)。在
首先,我必须计算一个预设模型的参数和相应的光谱能量(从这些参数中)。在
定义了一个函数,在这个函数下,我们找到了一个函数。在
var('a, b, t')
def spectrum(omega):
model = a*sin(omega*t) + b*cos(omega*t)
fit = find_fit(variations, model, parameters= [a, b], variables = [t], solution_dict = True)
sp_en = ((fit[a])**2 + (fit[b])**2)/2
return fit[a], fit[b], sp_en
然后调用函数并打印值:
^{pr2}$现在我只需要绘制输出。在
如果一个谱函数返回只返回sp
var('t')
plot(spectrum(t), (t, 1, 100))
问题是:如果我想打印所有三个输出,如何绘制这个函数?在
只需使用
energy
变量调用plot函数plot(energy, (omega, 1, 100)
)只需对spectrum的返回值使用索引:
相关问题 更多 >
编程相关推荐