计算类星体谱的斜率以进行幂律拟合

Question

使用 Sloan Digital Sky Spectrum，我创建了类星体的复合光谱。光谱是在Angstrom 中测量的波长（x 轴）和在ergs/cm^2/s/Angstrom 中测量的通量（y 轴）之间的图。

为了计算幂律，我将波长转换为频率，得到THz 比例频率。在 y 轴上，我首先将单位从 ergs/cm^2/s/Angstrom 更改为 ergs/cm^2/s/Hz，方法是将 this 转换为 MATLAB，然后再转换为 jansky 单位。 This 是结果图。

现在我想计算这个图的斜率，我应该使用基本的拟合工具并从那里获取m 的值吗？幂律f = f_o nu^(-slope)计算类星体谱斜率的其他方法有哪些。

Answer 1

使用 1 阶的 polyfit 来拟合线性模型。

P = polyfit(Xtrain, Ytrain, 1);

P 将是两个分量的向量：

• 偏移量（或 0 度分量）和
• 坡度（或一级分量）。

然后您可以通过应用拟合的线性模型来拟合新数据（测试数据）

Ytest = P(1) * Xtest + P(2);