---
{
  "path": "/research/code",
  "title": "计算与复现代码",
  "category": "研究资料",
  "summary": "汇集实际 LVE 演示、方向钟几何试运行与实验复现模板，说明输入、运行、输出和未完成部分。",
  "related": [
    "/theory/lve",
    "/research/data",
    "/collaborate/how-to-contribute"
  ],
  "id": "ASD-DOC-RESEARCH-CODE",
  "template": "document",
  "version": "1.1",
  "updated": "2026-10-03",
  "canonical": false,
  "source_tier": "website-guide",
  "depends_on": [],
  "open_problems": [],
  "equations": [],
  "sources": [
    {
      "label": "ASD 建站蓝图与页面模板（用户提供）",
      "kind": "supplied-blueprint"
    }
  ],
  "theory_status": "FOUNDATION",
  "evidence_status": "DOCUMENT",
  "title_en": "Code & Computation",
  "locale": "zh",
  "translation_path": "/en/research/code",
  "section_anchors": [
    "s0",
    "clock-pilot",
    "comparison-kit",
    "s1",
    "s2"
  ]
}
---

## 可下载的演示计算

本目录汇集已经存在的脚本、复现模板和研究试运行。先确定要核对的物理关系，再选择相应材料；这些项目使用的输入和完成的计算各不相同。

[LVE 演示脚本 ASD_LVE_demo.py](/downloads/ASD_LVE_demo.py)只使用 Python 标准库。它计算 $\beta=v/c$、$\gamma=(1-\beta^2)^{-1/2}$、横向尺度比 $L'_\perp/L_\perp=\gamma$，以及同速平行、横向分离构型中的力比 $F/F_0=1-\beta^2$ 和磁作用差额 $F_m/F_0=\beta^2$。长度关系来自 [LVE](/theory/lve)；力比还使用核外真值作用距离和逆平方本构输入，见[库仑作用](/theory/electromagnetism/coulomb-interaction)与[磁力推导](/theory/electromagnetism/magnetic-force)。脚本中的 `force_ratio` 对应这里的真值物理力 $F_T/F_0$，不能直接当作参考坐标广义力或探针读数。速度域为 $0\leq\beta<1$；这个速度范围不把给定构型推广为任意源与探针。

下载脚本后运行：

```sh
python ASD_LVE_demo.py --beta 0.6 --csv lve-model.csv
python ASD_LVE_demo.py --beta 0
```

第一条命令的终端结果为 $\gamma=1.25$、$F/F_0=0.64$、$F_m/F_0=0.36$；第二条给出静止极限 $1,1,0$。`--beta` 选择终端显示的单点，CSV 则固定扫描 0.00 至 0.95、步长 0.01，共 96 行，字段为 `beta,gamma,force_ratio,magnetic_fraction`。这些数值由给定公式生成，不含测量记录、误差估计或一般动力学求解。

## 方向钟试运行：几何与可识别性

[方向钟试运行](/research/clock-response-pilot)对应 `ASD-RUN-CLOCK-001`。复算代码已经包含在 [ASD_AI_Research_Pack.zip](/downloads/ASD_AI_Research_Pack.zip)中：`tasks/clock-pilot.py`、装置记录 `tasks/clock-pilot.json` 和冻结结果 `tasks/clock-pilot-expected.json`。这里的 JSON 是输入与未知量的记录；脚本在代码内固定下面的几何例子，并不读取 JSON 来设置装置参数。

设纬度 $40^\circ$，用户选定的背景轴赤纬 $30^\circ$，A、B 轴分别水平指北、指东，一个名义恒星日内采样 96 次。脚本使用未由 ASD 原子响应推导的几何模板 $Q=(\mathbf n\cdot\mathbf b)^2-1/3$，计算候选读数中三个系数对应的设计矩阵 $X=(Q_A,-Q_B,1)$。解压资料包后，在包的根目录运行：

```sh
python tasks/clock-pilot.py --verify tasks/clock-pilot-expected.json --output clock-pilot-output.json
```

输出 JSON 给出默认构型秩 3、赤纬 $0^\circ$ 时秩 2、两轴相同时秩 2；同时检查轴反向、一个周期的重复和同轴等系数的抵消。比较的背景轴和采样均为假设输入，秩只说明给定模板中参数能否分辨。代码目录中的方向钟结果只给出几何可识别性，尚未计算 ASD 频移振幅，也未使用观测数据。所选 $^{87}\mathrm{Sr}$ 跃迁、参考钟及传输链的响应系数仍是[钟响应任务](/open-problems/clock-response)的下一份输入；不能用满秩代替这些物理计算或检测功效。

## 实验比较基线与复现模板

[ASD_Experimental_Reproducibility_Kit.zip](/downloads/ASD_Experimental_Reproducibility_Kit.zip)包含 `README.md`、`equations.md`、`models.py`、`run.py`、`verify.py`、`notebook.ipynb`、`input/`、`schemas/`、`expected/` 和 `checksums.json`。它用 Python 3.10 及以上的标准库运行比较模型，输入文件与单位见[数据目录](/research/data)及包内说明。解压后运行：

```sh
python run.py --output results
python verify.py
```

`results/summary.json` 汇总结果，`runtime.json`记录实际 Python 版本；`clock_geometry.csv`、`muon_bins.csv` 和 `charge_response.csv` 分别保存几何因子、合成 cohort 的参考存活计数、标准点电荷场与假定高斯响应。这一包的钟几何有 97 个半小时采样点；它与上面的 96 点双钟秩试运行是不同计算，不能混用结果文件。`expected/` 提供冻结数值，`verify.py` 核对分发文件的 SHA-256、数值和静止/周期/对称等极限。

μ 子部分使用恒能量、平直大气的标准存活公式和三个明确标为 `SYNTHETIC` 的 bin；高 γ 部分使用匀速点电荷的标准归一化横向场和单位面积高斯卷积，不含实际束团、晶体与电光传递函数。四份 `schemas/*.csv` 只有表头，是接收未来输入的格式模板；S301 没有轨道积分或拟合。详细前提、输出含义和缺项见[复现资源](/research/reproducibility)与[数据目录](/research/data)。包内 `asd_output`、`predicted_difference` 和 `fit` 均保留 `null`。

包内已有诊断图。需要重绘时，另行安装说明中指定的 `matplotlib==3.10.8`，运行 `python figure.py --output results`；数值计算不依赖它。Notebook 调用相同函数，需要读者自己的 Jupyter 环境。本包不包含真实观测比较或 Ota 实验复现。

## 验证哪些结果

上述命令分别核对公式演示、几何秩与比较基线。2026-10-03 使用 Python 3.12.14 复算了实际分发材料：LVE 的 96 行 CSV、双钟秩和三个比较模型的 JSON/CSV 与所附冻结数值相符。运行结果应结合输入、代码、冻结输出和校验清单一起保存；数值一致只说明这份计算可核对，不会补出尚未求得的 ASD 装置响应。

## 完整复现仓库的要求

本目录沿用现有下载包和登记，当前没有经确认的公共计算仓库链接。进一步计算应明确问题、方程、边界条件、真实输入来源、环境、输出和核对方法，并把数值误差与测量不确定度分开。方向钟先求响应系数，μ 子先闭合路径到衰变概率的映射，高 γ 场先闭合源与探针的联合响应，S301 先核实对象、历元、参考系与协方差；相应入口保留在[预测登记](/tests/predictions)与[开放问题](/open-problems)。这些明确缺项完成后，再登记真实数据比较。