用于计算幂律拟合

Question

与斯隆数字天空频谱工作类星体光谱的斜率，我创建了类星体的合成频谱。频谱是波长（x轴）之间的Angstrom和通量（y轴）测量在测量ergs/cm^2/s/Angstrom的曲线图。

对于幂定律的计算，我转换成波长这导致THz规模频率的频率。在Y轴上，我首先通过平移this到MATLAB和比jansky单元改变了单位从ergs/cm^2/s/Angstrom到ergs/cm^2/s/Hz。 This是结果图。

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

Answer 1

使用polyfit用1度，以适应线性模型。

P = polyfit(Xtrain, Ytrain, 1); 

P将是两个分量的向量：

• 的偏移（或0度分量）和
• 斜率（或第一度分量）。

然后可以适合通过施加拟合线性模型

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

新数据（测试数据）