---
title: "从出射强度走到能拟合的光谱"
description: "把连续谱、仪器卷积和误差放入同一个观测模型，用完整算例区分随机误差、波长零点与大气模型误差。"
canonical_url: "https://sandboxphysics.com/zh/stories/stellar-spectroscopy/from-intensity-to-a-fitted-spectrum/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/stellar-spectroscopy/from-intensity-to-a-fitted-spectrum/"
last_modified: "2026-10-28"
---

[从恒星暗线走到丰度模型](https://sandboxphysics.com/zh/stories/stellar-spectroscopy/)

AC03 · L3

# 从出射强度走到能拟合的光谱

Sandbox Physics · 2026-10-28

先预测一条谱线，再问仪器会把它记录成什么样，最后才能拿预测和保存的样本比较。连续谱的位置、仪器响应和波长标定都要写进这个过程；它们也决定了拟合结果该相信到什么程度。

## 开始前需要什么

先读前两章。观测值减去预测值得到残差；把每个残差除以随机标准差，平方后相加，得到加权平方和。卷积在这里指仪器把一个波长的光分散到周围波长。

## 先把仪器的模糊写清楚

用 $q(\lambda)$ 表示上一章给出的归一化强度。理想仪器会直接记录它；实际仪器把一个波长的光分散到邻近像素。响应函数 $G$ 告诉我们各处收到的比例。把所有输入波长的贡献加起来，就是卷积。

仪器响应的积分必须等于一。这样一段平坦连续谱经过仪器后仍平坦，完整谱线的缺损面积也保持不变。数值计算先对 $q-1$ 作卷积，采样域外按零亏损延伸，相当于域外恢复单位连续谱，避免在边缘凭空出现暗带。

$$
q_{\rm inst}(\lambda)=\int G(\lambda-\lambda^\prime)q(\lambda^\prime)d\lambda^\prime,\qquad\int G(u)du=1
$$

响应函数的量纲是波长的倒数。若积分窗口太窄，展宽会把面积散到窗口外，此时有限窗口内的等值宽度仍会改变。

对应阅读：Condon & Ransom, Essential Radio Astronomy, Appendix E · Krumholz, ASTR 4003/8003 (2019), Class 1, Radiative transfer I

## 连续谱也要和谱线一起估计

真实标定后的光谱仍可能带着斜坡。短波一端的背景高一点，长波一端低一点，若硬把连续谱当成水平线，拟合常会挪动线心或改变线深来补偿。局部连续谱先用直线 $C(\lambda)=a+b x$ 表示，$x$ 是以固定中心为原点的无量纲波长坐标。

观测预测取 $M_i=C(\lambda_i)q_{\rm inst}(\lambda_i)$。给定速度、宽度和光学深度后，连续谱的两个系数仍是线性的，可直接解两个联立方程。外层只搜索非线性的三个谱线参数，计算量小，也方便逐项检查。

令 $u_i=q_i$、$v_i=x_iq_i$，其中 $q_i=q_{\rm inst}(\lambda_i)$，再令权重 $w_i=1/\sigma_i^2$。预测便是 $a u_i+b v_i$。平方和对连续谱截距和斜率分别求导，最小值处两个导数都为零；把含系数的项移到左边，就得到下面的两个方程。

矩阵只是把两个联立方程排在一起的写法。先算五个加权和，再解截距和斜率。只有两个基函数提供不同信息时，行列式才不为零；这也是为什么连续谱不能只靠同一个波长点同时决定高度和斜率。

$$
\chi^2=\sum_i\left[\frac{y_i-(a+b x_i)q_{\rm inst}(\lambda_i)}{\sigma_i}\right]^2
$$

观测值为 y，随机标准差为 sigma。只有误差独立且估计合理时，这个平方和才具有熟悉的卡方解释。

$$
\begin{pmatrix}\sum_i w_i u_i^2&\sum_i w_i u_i v_i\\\sum_i w_i u_i v_i&\sum_i w_i v_i^2\end{pmatrix}\begin{pmatrix}a\\b\end{pmatrix}=\begin{pmatrix}\sum_i w_i u_i y_i\\\sum_i w_i v_i y_i\end{pmatrix}
$$

这是给定谱线模型后的连续谱正规方程。解出 a、b 后，外层才比较下一组速度、宽度和光学深度。

## 完整算例：很好看的残差仍可带错速度

教程先固定一组可控的合成样本。注入速度为 $18\,\mathrm{km\,s^{-1}}$，仪器分辨本领为 $60000$。同时拟合倾斜连续谱、速度、Doppler 参数和参考光学深度，冻结脚本得到速度约 $18.0\,\mathrm{km\,s^{-1}}$，约化平方和约 $1.04$。它验证的是保存样本的参数回收。

把全部波长故意加上 $0.002\,\mathrm{nm}$，强度和噪声都不动，再拟合一次。速度变成约 $19.2\,\mathrm{km\,s^{-1}}$，残差几乎同样好。全局波长零点错误，恰好能伪装成整层气体的速度。

如果零点标准差给定为 $0.0002\,\mathrm{nm}$，在五百纳米处对应约 $0.120\,\mathrm{km\,s^{-1}}$ 的共同速度不确定度。再多测几个同一标定下的像素，也消不掉这一项。

$$
\delta v\simeq c\frac{\delta\lambda_0}{\lambda_0}
$$

这是共同零点误差的传播。减小这项共同误差，需要独立的波长标定信息。

## 真实光谱会把简化模型的缺点暴露出来

配套实测教程固定 CALSPEC 的织女星光谱，只取表头注明来自 STIS 的光学段。氢线宽而有翼，简单高斯函数在很小的统计误差面前留下巨大平方和。换积分窗口、改连续谱、换线型，表观线心与宽度都会改变。

这正是从教学谱走向天文资料时应该遇到的检查。一个能输出数字的优化程序，离可靠的恒星参数还差一个经得住残差和模型比较的物理描述。先报告线心、窗口和拟合形状的敏感性；要解释元素丰度，再进入三维非 LTE 大气。

对应阅读：STScI CALSPEC, alpha\_lyr\_stis\_011.fits (25 October 2022) · MAST reference atlases and CALSPEC releases · Amarsi, Liljegren & Nissen (2022), 3D non-LTE iron abundances in FG-type dwarfs, arXiv:2209.13449v3

教学薄层默认预测的固定结果

均匀薄层、固定电离阶段与人为跃迁，没有随机噪声。固定计算采用七百二十一个波长样本和归一化高斯仪器响应。

| 薄层温度 | $4500\,\mathrm K$ |
| --- | --- |
| 柱密度与参考值之比 | $N/N_0=1$ |
| 远离速度 | $12\,\mathrm{km\,s^{-1}}$ |
| 光谱分辨本领 | $R=60000$ |
| 强线中心 | $500.0200\,\mathrm{nm}$ |
| 总等值宽度 | $11.952\,\mathrm{pm}$ |
| 薄层与背景的源亮度比 | $0.201$ |

## [检查面积和深度](https://sandboxphysics.com/zh/labs/stellar-absorption-lines/)

1. 保持薄层参数不变，改变仪器分辨本领。
2. 比较中心深度与总等值宽度，并解释差异。
3. 导出两份模型谱，确认它们的输入和样本一起保存，再读教程中的真实谱对比。

## 练习与解释性答案

1. 为什么同一波长零点错误不会被大量像素平均掉？

   展开解释

   它把所有像素朝同一方向挪动，像同一把尺子的起点偏了。各像素不是独立地朝正负两边错，误差不会互相抵消。

2. 一个观测谱的约化平方和很大，先把误差棒放大到一就能解决吗？

   展开解释

   应先检查残差形状、连续谱、仪器响应和被遗漏的物理。任意放大误差棒会抹掉模型失配的信息，得不到更可靠的物理解读。

[上一篇：把吸收和发射放进同一个方程](https://sandboxphysics.com/zh/stories/stellar-spectroscopy/absorption-emission-and-growth/) [下一篇：亲手拟合一份固定的恒星光谱](https://sandboxphysics.com/zh/stories/stellar-spectroscopy/fit-a-frozen-calibrated-spectrum/)

## 来源与计算范围

1. [Krumholz, ASTR 4003/8003 (2019), Class 1, Radiative transfer I](https://www.mso.anu.edu.au/~krumholz/teaching/astr8003_s1_2019/notes1.pdf)

   第 24—28 式：积分因子、形式解与常源函数薄层；这里另作推导，并用数值输运检验。

2. [Condon & Ransom, Essential Radio Astronomy, Appendix E](https://www.cv.nrao.edu/~sransom/web/A5.html)

   辐射转移第 2.27 式、LTE Kirchhoff 关系第 2.30 式和 Planck 定律第 2.86 式；光学波段保留完整 Planck 函数。

3. [STScI CALSPEC, alpha\_lyr\_stis\_011.fits (25 October 2022)](https://ssb.stsci.edu/trds/calspec/alpha_lyr_stis_011.fits)

   固定版本 FITS 表头与四千七百五十至四千九百七十五真空埃之间的八十二行 STIS 观测值；教程附校验和及原表行号。

4. [MAST reference atlases and CALSPEC releases](https://archive.stsci.edu/hlsp/reference-atlases)

   用于说明数据来源与发布背景；练习固定第十一版，不追随会变化的 current 目录。

5. [Amarsi, Liljegren & Nissen (2022), 3D non-LTE iron abundances in FG-type dwarfs, arXiv:2209.13449v3](https://arxiv.org/pdf/2209.13449v3)

   第三版，2022 年 12 月 9 日：第 2.1—2.4 节、第 3.1 节及图 2、3、8、10。阅读指南检查约定和敏感性，没有重跑三维求解器。