【发布时间】:2016-02-18 16:19:37
【问题描述】:
简介
我正在根据给定的数据计算和绘制围绕脉冲星运行的行星的光谱能量。
我之前已将列表变体中的所有数据排序,维度为 [172, 2](172 行和 2 列)。
首先,我必须相应地计算预设模型的参数和光谱能量(根据这些参数)。
为此,我定义了一个函数,在该函数下定义了预设模型和获取模型和变化数据的 find_fit 函数。
代码
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
然后我调用函数并打印值:
c, v, energy = spectrum(20) #enter arbitray angular frequency here
print "Values for the given angular frequency : \n \n a = %f, b = %f, spectral_energy = %f " % (c, v, energy)
现在我必须绘制 only sp_en 输出。
“半解”
如果 spectrum 函数返回 仅 sp_en,这很容易。写就够了:
var('t')
plot(spectrum(t), (t, 1, 100))
问题是:如果我想打印所有三个输出,如何绘制这个函数?
【问题讨论】:
标签: python matplotlib plot physics multipleoutputs