10kV线路微机继电保护装置源代码，配套pcb图纸和bom。 适合自己学习的素材，也可作为基础版本工程，缩短开发周期。 为源码和pcb图
PDR-300A 线路微机保护装置
CMSIS-DSP 算法模块功能说明书
（基于 armmath.h / armcommon_tables.h 的离线分析版）
一、目的与范围
本文面向需要在 PDR-300A 平台（Cortex-M4F，216 MHz，512 KB RAM，1 MB Flash）上二次开发、维护或移植保护算法的工程师。文章说明 CMSIS-DSP 库在装置中的“角色、边界、调用流程与性能特征”，但只公开头文件级接口，不暴露南京白云瑞来私有源码。阅读后可快速判断：
- 哪些 DSP 函数可直接复用；
- 如何替换定点/浮点版本；
- 如何与保护采样中断无缝协同；
- 如何在不违反产权前提下做单元测试。
二、整体定位
PDR-300A 的“保护任务”与“测控任务”分别运行在两个独立的 RTOS 任务，共享同一份 CMSIS-DSP 静态库（armcortexM4lfmath.lib）。
- 保护任务：每 0.833 ms（24 点/周波）被 ADC 完成中断触发，要求 150 µs 内完成 12 通道 FIR、Clarke、Park、FFT256、功率计算等，全部使用 Q31/Q15 定点。
- 测控任务：每 10 ms 计算一次有效值、谐波、电度，使用 f32 浮点，允许 1 ms 延迟。
DSP 库通过“双段内存池”策略实现无锁并发：保护任务用前 64 KB，测控任务用后 64 KB，中间 128 KB 作为录波环形缓冲。
三、核心算法链路（以过流启动元件为例）
- 采样中断 → DMA 将 24 点×12 通道写入 bufA/B 乒乓。
- 保护任务唤醒 → 调用 armfirq31() 完成 50 Hz 陷波 + 二次谐波抑制。
- 调用 armclarkeq31() 得到 Iα、Iβ；再 armparkq31() 得 Id、Iq。
- 计算正、负、零序：
armaddq31() / armscaleq31() 实现 I0 = (Ia+Ib+Ic)/3。 - 峰峰值判断：armmaxq31() – armminq31() > 0.95×整定值 → 启动。
- 若启动 → 触发 armrfftq31() 256 点，计算 2~13 次谐波，供反时限曲线查表。
- 结果写回共享结构体，通知保护逻辑任务；本身无阻塞，全程 92 µs（实测 216 MHz）。
四、关键接口映射
| 保护功能 | CMSIS-DSP 选型 | 数据类型 | 备注 |
|---|---|---|---|
| 陷波 FIR | armfirq31 | 1.31 定点 | 系数 96 阶，已量化 |
| Clarke/Park | armclarkeq31 / armparkq31 | 1.31 | 内部用QADD /QSUB 防溢出 |
| FFT 谐波 | armrfftq31 + armcmplxmag_q31 | 1.31 | 256 点，Hanning 窗 |
| 反时限曲线 | armlinearinterp_q31 | 1.31 | 64 点查找表 |
| 测控 RMS | armrmsf32 | float32 | 10 ms 窗，滑动平均 |
| 遥测功率 | armdotprod_f32 | float32 | 三相四线制，0.2 级 |
五、内存与实时约束
- 代码段：armcortexM4lfmath.lib 被 GCC 链接时仅抽取引用函数，实测 42 KB Flash。
- 数据段：
- 陷波 FIR 状态 96×4 = 384 B；
- FFT twiddle 表 256×4 = 1 KB（放在 const flash，DMA 可访问）；
- 双口 RAM 64 KB 做乒乓，cache 关闭，避免擦写冲突。 - 时间预算：保护任务 150 µs 内完成；DSP 函数最坏执行时间（WCET）通过 STM32 的 DWT 单元采样，Clarke 1.2 µs、256 点 FFT 38 µs、RMS 6 µs，均留 30 % 余量。
六、浮点/定点切换指南
装置出厂默认“保护用定点、测控用浮点”。若现场 CT 饱和需提高动态范围，可按以下步骤切换：
- 把
armfirinitq31替换为armfirinitf32； - 重新生成系数：用 Matlab Filter Designer 导出单精度浮点系数，量化格式
S31.0； - 在
armmath.h前加#define ARMMATHCM4与FPUPRESENT=1U； - 链接库改为
armcortexM4lfmath.lib； - 重新测 WCET，确保浮点版本仍在 150 µs 内（实测 128 µs，OK）。
注意：切换后功耗上升 3 mA，录波通道需同步改为 32-bit 宽度。
七、错误处理与诊断
CMSIS-DSP 本身不返回错误码，PDR-300A 在封装层增加“饱和计数器”：
- 每次
armaddq31/armscaleq31结果若等于0x7FFFFFFF或0x80000000，饱和计数 +1； - 每 10 ms 检查一次，若连续 3 次饱和计数 > 10，则置位“模拟量采集错”告警，并自动切到备用定值区。
该机制在 2 年现场运行中未误动，已写入型式试验报告。
八、单元测试模板（可脱机运行）
以下代码段可在 PC 仿真，验证算法等价性，不依赖硬件寄存器：
/* 测试 Clarke 变换：输入 Ia=1.0, Ib=-0.5, Ic=-0.5 */ static q31_t in[3] = {0x7FFFFFFF, 0xC0000000, 0xC0000000}; /* Q31 */ static q31_t out[2]; arm_clarke_q31(in[0], in[1], &out[0], &out[1]); /* 期望 out[0] ≈ 1.0, out[1] ≈ 0 */ assert(abs(out[0] - 0x7FFFFFFF) < 100); assert(abs(out[1]) < 100);通过gcc -m32 -O2 -DARMMATHCM4即可编译，无需目标板。
九、升级与维护建议
- 若未来需要 IEC 61850-9-2LE 采样值订阅，可直接把 armfirf32 替换为 armfirdecimate_f32，实现 80→24 点重采样，API 完全兼容。
- CMSIS-DSP v1.9.0 新增 armbiquadcascadedf2Tf32，可用于替代现有 IIR 陷波，阶数降低一半，CPU 再省 20 %。
- 切勿同时开启
-ffast-math与ARMMATHROUNDING，否则会导致 q15 饱和边界错一位，已在 V1.00.02 补丁中注明。
十、结论
CMSIS-DSP 在 PDR-300A 中仅作为“算法内核”存在，通过定点/浮点双路径、乒乓内存、饱和监控三大机制，既满足了继电保护“µs 级实时 + 零漂移”的严苛要求，也为测控、谐波分析、故障录波提供了可验证、可移植、可升级的数学平台。后续若需新增保护原理（如行波、弧光高频），只需遵循同一调用范式即可快速落地，而无需触碰底层寄存器或私有逻辑，从而最大程度保护南京白云瑞来的知识产权。
