wandering1
0f7b5a87fb
登记簿、文档与代码注释里的存量字母加数字编号（N-x / S-x / T-x / E-x / A-x / 卡号） 一律改写为语义短语，净删 2,028 处引用；删除基线与待确认、并行分道方案、排期详述、 采集保护收尾 4 篇（内容已并入 docs/任务规划/README 与登记簿）。 门禁机器与三张账同批重打：check_ids.py 新增变更记录行豁免与第二根护栏 undecidable_total_ceiling；check_test_traceability.py 只加 stdout/stderr UTF-8 reconfigure（TC_RE 未动）；check_cubemx_survival.py 的 --write 固定 newline=LF； id_registry / id_aliases（+19 条 retired-to-phrase）/ id_known_backlog 重打。 test/ 56 个文件同步 TC 标签语义化，其中 18 个 TU 改了 CHECK / printf 里的标签文字， 是本批唯一改变构建产物的文件集合。 为什么合成一笔、不拆四笔（实测撞闸，非偏好）：pre-commit 跑钩子前会 stash 未 staged 部分，拆开后中间提交必然是「新规则 + 旧语料」——实跑可判定 85.15% < 基线 86.00%、 不可判定 1,453 > 上限 1,166，编号门禁直接拦下。编号识别面与它描述的语料只能同码移动。 实跑基线（工作区全量）：754 文件 / 8,230 处编号字面量 / 可判定 7,079 = 86.01% （基线 86.0%，容差 0.1pp）/ 不可判定绝对数 1,151（上限 1,166）/ 校验 1-6 错误级 0 / 告警 = 悬空 2（App/algo/freq.c 的 S2，属 §4.4 同号不同义）+ 孤儿 50。 基线从 84.0 抬到 86.0 的理由见 docs/门禁说明/id-check.md §7.6：本批真删了 2,028 处 引用，可判定面随真实内容上升，不是把分母做小。同批新增排除 docs/硬件平台/ 语料，使 识别面少 264 处——芯片手册位号与网表料号不属编号空间（registry 有 excluded_paths_note）。 ⚠ 余量只剩 0.11pp（约 10 处）与绝对数 15 处：把「可判定的引用」改成散文会立刻撞第二 根护栏。另遗留一项：id_aliases 里 retired 的 N / S / T / A 别名没有 scope（U1~U5 有）， 今后有人新写一个 A-4，校验 3 不会响，只能靠人，建议下批补 scope。 固件不变证据：115 个 .c/.h 剥注释并归一字符串后与 HEAD 逐字节全等，残留代码差异 0； App / Src / Drivers 的字符串字面量零改动。@build 0 错误 0 警告；40/40 主机套件、2,889 断言全通过；覆盖率 line 89.2% / func 92.6%；pre-commit 48 钩子 --all-files 全绿； check_doc_links --all 通过（176 文档，章节引用无法判定 0）；逐字节卫生核过：14 个带 BOM 文件与 HEAD 一致、0 个 CR、0 个 U+FFFD，字符串字面量纯 ASCII。 .cubemx_survival_baseline.json 的 watchdog-module 哈希 56bf27 变 2bb051 属合法： Src/freertos.c 的 USER CODE 区内只有注释文字随编号改写，代码未变。 同批带上体积同码凭证与 RAM 台账行（@build 之后盖戳）：RW_IRAM1 79,008 B / 60.3%， RW_CCM 57,376 B / 87.9%，阈值 92%。
120 lines
4.7 KiB
C
120 lines
4.7 KiB
C
/*
|
||
* 模块名称：相量计算层（L2 相量计算）
|
||
* 模块功能：全周傅氏（全周波傅里叶）算法——对一周波 20 点采样做单频 DFT，
|
||
* 产出基波/谐波相量。系数表**预乘 √2/N**（口径见
|
||
* `docs/算法/傅氏算法详解.md` §2.3），故 `fAmSQ = √(fRe²+fIm²)`
|
||
* **直接就是有效值**，无需再除 √2。
|
||
* 覆盖规划 §3.2：L2-1 系数表 / L2-2 实例 / L2-3 有效值输出 /
|
||
* L2-6 每 1 ms 滑动窗 / L2-7 双区指针原子交换。
|
||
* ⚠ 本层**纯算法**：不碰寄存器、不碰任务调度、不自己取数（输入由
|
||
* `sample_ring` 提供）→ 主机侧可单测（L0.5）。
|
||
* 适用平台：STM32F407ZGTx（Cortex-M4F，硬件 FPU）
|
||
* 作者：王建锋
|
||
* 创建日期：2026-10-10
|
||
* 修改记录：
|
||
* v1.0 2026-10-10 王建锋 创建：L2-1/2/3/6/7（基波 + 谐波相量）
|
||
*/
|
||
|
||
#ifndef PHASOR_H
|
||
#define PHASOR_H
|
||
|
||
#include <stdint.h>
|
||
|
||
#include "adc_input.h" /* adc_semantic_t：实例按语义通道绑定 */
|
||
|
||
#ifdef __cplusplus
|
||
extern "C" {
|
||
#endif
|
||
|
||
/* 每周波采样点数（1000 Hz @ 50 Hz，见 docs/子系统设计/ADC采样节拍设计说明.md） */
|
||
#define PHASOR_N 20U
|
||
|
||
/* 相量结果（**有效值口径**）。字段名沿用参考系统 `{fRe, fIm, fAm, fAmSQ}`
|
||
* （规划 §10.2 术语表：实部/虚部/幅值/有效值）：
|
||
* fRe - 实部 = a1/√2（有效值相量实部）
|
||
* fIm - 虚部 = b1/√2（有效值相量虚部）
|
||
* fAm - 幅值 = √2·√(fRe²+fIm²)（正弦峰值 A1）
|
||
* fAmSQ - 有效值 = √(fRe²+fIm²)（系数已折算 √2/N，故直接即有效值）
|
||
* 相位（L2-4）按需由 `atan2(fIm, fRe)` 导出，不常驻本结构。 */
|
||
typedef struct {
|
||
float fRe;
|
||
float fIm;
|
||
float fAm;
|
||
float fAmSQ;
|
||
} phasor_t;
|
||
|
||
/* 相量实例（L2-2）：**数据驱动**——新增实例只需在实例表尾加一行 */
|
||
typedef struct {
|
||
const char *name; /* 实例名（日志 / UI 用） */
|
||
adc_semantic_t sem; /* 语义通道 */
|
||
uint8_t harm; /* 谐波次数（1=基波，2=二次，5=五次…） */
|
||
} phasor_inst_t;
|
||
|
||
/*
|
||
* 函数功能：相量层初始化
|
||
* 入口参数：无
|
||
* 返 回 值：0 - 成功
|
||
* 限定条件：内核启动前调用（无 RTOS 依赖）
|
||
* 函数说明：清双区结果与有效标志，并复位采样环缓冲；可重复调用（幂等）。
|
||
*/
|
||
int phasor_init(void);
|
||
|
||
/*
|
||
* 函数功能：取相量实例数（L2-2）
|
||
* 入口参数：无
|
||
* 返 回 值：实例个数
|
||
* 限定条件：无
|
||
* 函数说明：供遍历 / 日志打印。
|
||
*/
|
||
uint8_t phasor_inst_count(void);
|
||
|
||
/*
|
||
* 函数功能：取相量实例描述
|
||
* 入口参数：idx - 实例下标
|
||
* 返 回 值：指向静态实例表的指针；idx 越界返回 NULL
|
||
* 限定条件：无
|
||
* 函数说明：返回静态存储地址，调用方不得修改。
|
||
*/
|
||
const phasor_inst_t *phasor_inst(uint8_t idx);
|
||
|
||
/*
|
||
* 函数功能：重算全部实例并从计算区原子切换到使用区（L2-6 滑动窗 + L2-7 双区交换）
|
||
* 入口参数：无
|
||
* 返 回 值：无
|
||
* 限定条件：仅任务上下文调用；须先经 `sample_ring_push_frame()` 压入本拍数据
|
||
* 函数说明：以"当前时刻为终止的前一周波"为窗，**20 点全部重算**（暴力法，
|
||
* 无累积误差）。**任一实例数据不足一周波则本拍整体不产出**（不切换
|
||
* 使用区），避免部分新部分旧的混合结果。每拍更新，故相量刷新周期 = 1 ms。
|
||
*/
|
||
void phasor_run(void);
|
||
|
||
/*
|
||
* 函数功能：读实例结果（使用区，无锁）
|
||
* 入口参数：idx - 实例下标
|
||
* out - 输出相量
|
||
* 返 回 值：0 - 成功；-1 - 参数错；-2 - 尚无有效结果（未攒满一周波）
|
||
* 限定条件：无
|
||
* 函数说明：读的是"使用区"，与正在写的"计算区"物理隔离 → 无撕裂、无锁。
|
||
*/
|
||
int phasor_get(uint8_t idx, phasor_t *out);
|
||
|
||
/*
|
||
* 函数功能：对 n 点窗口算 harm 次相量（L2-1 系数表 + L2-3 有效值）
|
||
* 入口参数：x - 采样窗（按时间序，旧 → 新；长度须 = PHASOR_N）
|
||
* n - 窗长（必须等于 PHASOR_N）
|
||
* harm - 谐波次数（>=1）
|
||
* out - 输出相量
|
||
* 返 回 值：无（参数非法时 out 保持不动）
|
||
* 限定条件：无
|
||
* 函数说明：**纯函数**（无状态），可脱离采样环独立调用，便于单测与复用。
|
||
* 索引步进 = 谐波次数即可取任意次谐波（sin 以 2π 为周期，
|
||
* 故 `sin(2π·((n·h) mod N)/N) ≡ sin(2π·n·h/N)`）。
|
||
*/
|
||
void phasor_calc(const float *x, uint8_t n, uint8_t harm, phasor_t *out);
|
||
|
||
#ifdef __cplusplus
|
||
}
|
||
#endif
|
||
|
||
#endif /* PHASOR_H */
|