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

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);