信号相干函数分析与频域评估系统

资 源 简 介

本项目主要用于在MATLAB开发环境中实现两个离散时间信号之间相干函数的计算与分析。相干函数（Magnitude Squared Coherence）是频域分析中的关键指标，用于评估两个随机信号在不同频率分量上的线性相关程度。该程序通过对输入的时域信号进行分段处理，采用经典的韦尔奇（Welch）改进周期图法，对信号施加特定的窗函数（如汉宁窗或汉明窗）以抑制频谱泄漏，并通过设置重叠区间来降低谱估计的方差。核心计算逻辑涉及对信号进行快速傅里叶变换（FFT）以获得频域表示，随后计算各信号的自功率谱密度（PSD）

详 情 说 明

MATLAB相干函数分析与实现系统

项目介绍

本项目是一个专门用于分析两个离散时间信号之间频域线性相关性的MATLAB系统。通过计算相干函数（Magnitude Squared Coherence），系统能够量化两个随机信号在各个频率分量上的互相关程度。该系统不仅提供了自动化的信号模拟环境，还完整实现了基于韦尔奇（Welch）法的手动功率谱估计逻辑，是研究信号传递关系、系统辨识及噪声干扰分析的有效工具。

功能特性

• 高精度谱估计：采用韦尔奇改进周期图法，通过分段加窗和重叠处理，有效降低了谱估计的方差，提升了频率分辨率。
• 线性相关性量化：计算结果严格限制在 [0, 1] 范围内，清晰直观地展示信号在不同频点的线性依赖程度。
• 全流程参数化：支持自定义采样频率、分段长度、重叠点数及FFT点数等关键谱分析参数。
• 多维度结果可视化：系统可自动生成时域对比图、自功率谱密度（PSD）对比图、相干函数谱、互功率谱（CSD）幅值谱及相位谱，提供全方位的频域特征描述。
• 自动化统计分析：自动识别并输出相干性最高的频率点及其对应系数值。

使用方法

1. 启动MATLAB开发环境。
2. 将核心程序脚本放置于当前路径下。
3. 直接运行该主程序函数。
4. 在命令行窗口查看最大相干点的统计数据。
5. 在弹出的两个图形窗口中分析时域信号、功率谱、相干函数以及互谱的物理含义。

系统要求

• MATLAB R2016b 或更高版本。
• 无需额外工具箱（程序内置了汉宁窗生成子函数，兼容性强）。

实现逻辑与算法细节

系统的核心执行流程严格遵循信号处理的标准规范：

1. 模拟信号仿真逻辑 系统首先构造两个具有物理关联的测试信号。信号1包含特定的谐波分量（50Hz和150Hz）及高斯白噪声。信号2则模拟了一个线性传递过程，通过对信号1施加样本延迟、增益缩减并叠加独立的背景噪声产生。这种设置确保了在特定频点存在高相干性，而在噪声占主导的频点呈现低相干性。

2. 预处理与分段参数 设定采样频率为1000Hz。通过定义256点的分段长度和128点的重叠区间，系统将长信号切分为多个重叠的子序列。这种重叠比例（50%）是减少估计偏差的经典配置。

3. 核心计算逻辑（Welch法手动实现） 相干函数的计算并非简单地调用库函数，而是通过以下底层步骤实现：

• 加窗处理：对每个分段应用自定义的汉宁窗（Hanning Window），旨在平滑信号边缘，抑制由于非周期截断引起的频谱泄漏。
• FFT变换：对每个加窗后的分段进行快速傅里叶变换，将时域能量映射到频域。
• 功率谱累加：在迭代循环中，分别累计算各个分段的自功率谱密度（信号各自的模平方）和互功率谱密度（信号1的频谱与信号2频谱共轭相乘）。
• 均值平滑：对累加后的结果进行段数平均，并进行能量归一化处理。
• 单边谱修正：为了遵循物理真实性，对单边谱（DC和奈奎斯特频率除外）进行了幅值倍增修正。

4. 相干函数定义应用 利用公式 $Cxy(f) = |Pxy(f)|^2 / (Pxx(f) * Pyy(f))$ 计算相干系数。该步骤通过互谱的模平方除以两个自谱的乘积，消除了信号绝对幅度的影响，仅保留两者之间的相位一致性和线性相关信息。

关键函数与细节说明

• 汉宁窗子函数：程序内部实现了 hanning(n) 函数，通过余弦公式生成窗函数序列，确保即便在缺少信号处理工具箱的环境下也能正常运行。
• 能量归一化：通过计算窗函数的平方和（window_norm）对FFT结果进行修正，保证了功率谱密度数值的物理标准性。
• 相位谱分析：系统不仅计算了相干幅值，还计算了互功率谱的相位（angle(Pxy)），这对于识别两个信号之间的时延和相位差至关重要。
• 交互式可视化：绘图部分使用了 subplot 布局，将时域、频域和统计域的信息整合对比，通过红色圆点标注出相干性峰值，通过分贝（dB/Hz）刻度展示功率谱动态范围。