雷达 · 通信 · 电子战 / 基础精学 156
普通幅平方相干比较同一频率的两路随机信号,频谱平均是使其具有统计意义的关键,单段谱构造会恒等于一。
实际记录只包含有限个采样值。由记录计算得到的数值与总体参数并不相同:前者受记录长度、噪声、边界约定及运算精度影响,后者属于给定模型的属性。建立可复核结论需要先固定输入含义、统计假设和归一化方式,再比较可计算量与理论量。
两路联合宽平稳;不同数据块近似独立;只考察普通同频相干,不涉及非零循环频率的谱相干。
Sxy:普通互功率谱(功率每赫兹)
Cxy:幅平方相干(无量纲)
K:独立平均次数(块)
两路互谱模平方除以各自功率谱乘积。
先平均复互谱和两路功率,再取模平方和比值。
仅用一块数据时,分子与分母完全相同,不能识别独立噪声。
输出由经过线性系统的输入和独立噪声叠加。
代入定义得到可由输入解释的输出功率比例。
独立圆对称复高斯频谱下,样本相干存在正偏差。
单个频点只用一段数据估计时,互谱模平方与自谱乘积完全相同:
在两项均非零处,无论两路真实统计关系如何均得到一。
对独立零相干的复高斯频谱样本,若使用 K>1 个独立段的普通等权平均,样本幅平方相干具有有限平均偏差,其均值为
该结果依赖独立段和相应高斯频点模型;
重叠窗口使有效独立平均次数小于名义段数。因而低但非零的测得相干不必意味着真实关联,比较结果时应同时报告平均方式、段数和重叠程度。
相干并不是单次FFT中两个数相乘后的相似性,而是跨实现保持稳定的复相位关系。独立噪声的互谱会在平均时相互抵消,但功率不会。普通相干的两个输入属于同一频率;循环相干比较相隔循环频率的谱分量,二者不能混用。
生成独立复高斯频谱样本,对1、2、8、32块平均计算零相干偏差,并对信号加噪模型验证总体相干为0.8的预测及单块恒等于一。数值脚本使用固定随机种子;随机实验的判据采用有限容差,代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差,关键数值写入独立结果文件。
| 量 | 数值 |
|---|---|
| repetitions | 6000 |
| single block identity error | 1.11022302463e-15 |
| null mean K2 | 0.505819441278 |
| null mean K32 | 0.0311979991262 |
| null theory K32 | 0.03125 |
| linked mean K32 | 0.802298550671 |
| linked population | 0.8 |
图 1 普通幅平方相干与平均次数的固定参数数值验证
表中数值和图中曲线由 MATLAB 实际运行得到,随机种子为 20261075;运行版本与完整结果另存为独立文件。
SystemVerilog实现复互谱实虚累加和两路谱能量累加;相干比值计算在软件中完成。端口为16位有符号码,内部64位,同步复位,valid空拍保持状态;
按外部首尾标志进行块累加并在块尾输出。自检覆盖多组数据、空拍、边界和复位,仅验证所述整数子核,不代表完整算法的硬件验证。
独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。
适用于线性系统辨识的输入输出一致性检查及通道噪声分离。平均次数决定估计波动,报告相干曲线时应同时标注分段长度、重叠与有效自由度。
重叠段并非独立,实际有效平均次数小于原始分段数。
非线性关系可能很强但线性相干仍较低。
有限精度实现必须区分输入量化、内部舍入和最终结果格式,不能仅依据浮点曲线估计硬件误差。
相干必须建立在复谱平均之上;单段结果恒等于一,有限平均还会产生非零偏差。
内容依据:每日精学第 156 课《普通幅平方相干与平均次数》。图表沿用原课固定参数数值结果,属于数值验证,不代表设备实测。