---
title: "复算真实SMA图像与红外谱，保留多种解释"
description: "用固定HD48682材料独立运行Python条件拟合，核对束单位并比较半径与颗粒退化。"
canonical_url: "https://sandboxphysics.com/zh/stories/debris-disks/fit-real-sma-and-infrared-data/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/debris-disks/fit-real-sma-and-infrared-data/"
last_modified: "2026-10-28"
---

[从颗粒碰撞走到碎屑盘辐射证据](https://sandboxphysics.com/zh/stories/debris-disks/)

AC16 · Tutorial

# 复算真实SMA图像与红外谱，保留多种解释

Sandbox Physics · 2026-10-28

一幅暗弱恢复图与六个流量点，在指定尘埃模型下究竟要求什么？这个练习拟合辐射总体，不拟合隐藏行星。

## 开始前需要什么

需要上一章的热辐射映射与加权最小二乘。Python3.9及以上即可，无需私有仓库、账号或额外包。

## 先拿到数据，再读拟合

下载下方公开真实数据压缩包，解压到空目录。里面有未改动的74880字节HD48682图像、REASONS原始拟合注记、六个转录的观测流量点、五个Python文件及八个预期输出。出处文件记录精确版本与SHA256。恢复图来自Zenodo12583396v0.0，按CC BY4.0分发；单独发表的流量数值保留各自引用，不把归档许可扩展到论文。

文件名虽然带有calibratedvis，它仍是一幅图像，不是可见度表。连续谱参考频率为229.4830981104吉赫兹，RESTFRQ另外记录230.538吉赫兹。脚本保留区别，屏蔽4799个非有限边缘像素，保留赤经负增量，使用实际椭圆恢复束；这些细节会改变物理预测。

选用论文自身PACS约70、100、160微米的处理结果，再加入已发表250、350、850微米探测。同一PACS观测的另一种处理不作为额外数据重复统计。未经核实置信约定的上限，以及相关的Spitzer合成谱箱不进入拟合。2020年勘误修订观测与标定描述，这里没有重做标定流水线。

对应阅读：Matrà · REASONS archive 12583396, version v0.0, CC BY 4.0 · Hengst et al. 2020 · HD 48682, Tables 1–2 · Hengst et al. 2020 · published erratum

## 在独立目录里运行

先运行教程的七个独立检查，覆盖Planck有限积分、质量积分、平衡温度、真实束归一化与图像和SED的流量闭合。再运行分析，与expected目录比较。相同平目录包已在CPython3.9.6和3.13.11中断网复算。

分析比较448个条件族。尘埃位于试探半径的0.9、1.0、1.1倍处，质量比例固定；24个对数粒径箱保留给定数量律的质量积分，每个粒径与环带分别求辐射平衡。同一总体生成六波段光谱和SMA频率处的图像。恒星黑体半径由指定光度与温度推得，同时加入两种预测。

拟合中刻意没有行星参数。动力学Lab是另一颗合成恒星周围的独立实验；环带辐射拟合成功，不能被改名为用Lab行星拟合HD48682。材料也明确记录了未核实的交付主束响应、缺失标定协方差和粗几何网格。

对应阅读：Matrà · REASONS archive 12583396, version v0.0, CC BY 4.0 · Hengst et al. 2020 · HD 48682, Tables 1–2

## 推导那个被精确求解的参数

每个固定网格点先预测单位参考质量的尘埃信号，再加入固定恒星贡献。只有非负尘埃质量振幅被连续优化。对加权残差平方和求振幅导数，一条线性方程给出加权乘积之比；比值为负时，受约束最优值在零边界。

流量点按发表方差倒数加权。图像采用15角秒内每隔一个的有限像素，以外环稳健估计噪声，再乘采样面积与束面积之比。这避免了把过采样全部当成独立信息，却没有提供完整CLEAN协方差。光谱和图像得分分别保存，总和只用于条件例子的排序，不是校准后的概率或置信区间。

$$
S(A)=\sum_k w_k[y_k-s_k-A d_k]^2,\qquad \frac{dS}{dA}=-2\sum_k w_kd_k(y_k-s_k-A d_k)
$$

观测、固定恒星信号与单位质量尘埃预测必须采用相同单位。

$$
A_{\rm best}=\max\!\left(0,\frac{\sum_k w_kd_k(y_k-s_k)}{\sum_k w_kd_k^2}\right)
$$

只要试探尘埃模型不是处处零信号，分母就是正的。

对应阅读：

## 看替代解，而不只看第一行

保存网格选中半径40天文单位、最小粒径5微米、光学指数二的条件族，总诊断得分为14.3188。另一组半径25天文单位、最小粒径15微米，得分14.3248；半径60天文单位、最小粒径0.2微米、指数一，得分16.1689。这些物理假设差异很大，在受限计算中的排序却相近。

把conditional-profile.csv、geometry-profile.csv与model-sed.csv放在一起读，检查哪些解落在网格边界，再看原始注记怎样描述方向与宽度的弱约束。最小得分不能直接变成唯一半径或行星质量；扩大网格、改变颗粒组成，都可能移动最小值。

刻意受限的黑体替代模型还在更宽的连续半径区间内优化，避免仅因网格漏掉最佳半径就将它拒绝。其光谱得分约28.7303，半径约6.715天文单位。在给定误差约定下，它比灵活颗粒谱更差；这个比较仍依赖简化发射律与尚未建模的标定相关性。

对应阅读：

## 只改变一个条件，并保存过程

先固定几何，改变最小粒径。读新质量振幅之前，先解释温度与长波效率为什么变化。再只改变图像倾角：光学薄总体谱不应因此直接改变。如果投影后总流量改变，先检查有限视场损失与归一化。

选中的辐射模型还用48个粒径箱、四倍积分加密复算。保存波段变化低于0.234%，仅积分加密低于0.00731%。这约束的是这个特定辐射计算，不是所有光学律，更不是随机Lab动力学的误差。

迁移练习可以去掉图像项，重新比较半径剖面。说明图像实际减少了哪一种歧义，还有哪些留着。把变更假设、分开的得分与残差列成表，比只保留一张最佳拟合图更有信息。

对应阅读：

## [区分合成动力学与真实观测](https://sandboxphysics.com/zh/labs/debris-disk-sculpting/)

1. 先查看Lab合成记录及恒星条件，再运行独立真实数据包。
2. 比较两类预测的假设，不把模型行星搬到HD48682。

## 练习与解释性答案

1. 两组得分为14.3188和14.3248，能从差值直接报告概率比吗？

   展开解释

   不能。图像项是近似束面积诊断，没有校准似然与完整协方差。接近的分数说明本计算没有清楚区分两组条件族，不是后验赔率。

2. 要把盘中心偏移解释成偏心率，还缺哪些信息？

   展开解释

   需要观测时刻可靠恒星位置、天体测量标定、投影几何与真实偏心盘模型。本教程只将恒星点源固定于FITS参考像素，探索有条件的偏移。

## 公开科学材料

- [下载真实SMA与SED输入、脚本及预期输出](https://sandboxphysics.com/stories/debris-disks/materials/real-data-materials.zip)
- [阅读独立目录运行说明与边界](https://sandboxphysics.com/stories/debris-disks/materials/real-data-README.md)
- [检查条件替代解](https://sandboxphysics.com/stories/debris-disks/materials/real-data-summary.json)
- [输入哈希、单位、许可与精确来源版本](https://sandboxphysics.com/stories/debris-disks/materials/real-data-provenance.json)

```
python3 test_tutorial.py .
python3 analyze_debris.py --input . --output rerun
python3 verify_outputs.py --expected expected --actual rerun
```

[上一篇：用同一份颗粒记录生成热辐射图与光谱](https://sandboxphysics.com/zh/stories/debris-disks/turn-one-grain-record-into-an-image-and-sed/) [下一篇：沿ARKS的尘埃弧追问行星证据](https://sandboxphysics.com/zh/stories/debris-disks/read-arks-without-turning-every-arc-into-a-planet/)

## 来源与计算范围

1. [Matrà · REASONS archive 12583396, version v0.0, CC BY 4.0](https://doi.org/10.5281/zenodo.12583396)

   未改动的HD48682 SMA恢复FITS与拟合注记，保留图像单位和实测束；虽然文件名含 calibratedvis，交付的是恢复图像，不能按可见度表读取。论文背景固定arXiv2501.09058v1，2025年1月15日，A&A693A151。

2. [Hengst et al. 2020 · HD 48682, Tables 1–2](https://doi.org/10.1093/mnras/staa1972)

   已发表恒星参数和六个冷尘流量探测数值；PACS只保留一种处理，不纳入置信约定未核实的上限与相关合成谱箱。

3. [Hengst et al. 2020 · published erratum](https://doi.org/10.1093/mnras/staa2874)

   2020年10月9日勘误修订并扩展SMA观测与标定描述；交付REASONS产品另行固定版本，没有重跑标定流水线。