单频信号Prony法谱估计与参数提取工具包

资 源 简 介

本项目提供了一套完整的基于Prony法的单频信号谱估计例程与MATLAB工具包，旨在通过参数化建模方法精确提取信号的频率、幅值、相位和衰减因子。Prony法通过将等间距采样的观测数据拟合为一组指数函数的线性组合，有效解决了传统傅里叶变换在短数据记录下频率分辨率受限的问题。功能涵盖了从基础的线性预测方程构建、扩展Prony法的最小二乘求解，到多项式求根及模型参数映射的全过程。工具包内包含多个精心设计的MATLAB脚本与函数，展示了在不同信噪比条件下单频信号的参数估计稳定性。本资源不仅详细演示了数学推导的算法

详 情 说 明

单频信号谱估计Prony法MATLAB工具包

项目介绍

本项目提供了一套完整的基于Prony法的信号谱估计与参数提取方案。Prony法是一种经典的高分辨率谱估计技术，其核心思想利用指数函数的线性组合来拟合等间距采样的观测数据。相比于传统的傅里叶变换（FFT），Prony法在处理短数据记录时具有显著的频率分辨率优势，能够同时精确估计信号的频率、幅值、相位以及衰减因子。本项目通过一个高度集成的MATLAB主程序，完整演示了从数学理论到算法实现的转化过程。

功能特性

1. 多参数联合估计：不仅限于探测信号频率，还能够提取信号的指数衰减因子、初始相位和精确幅值。
2. 稳健的算法实现：采用扩展Prony法，通过设置高于理论最小值的模型阶数，增强了算法在噪声环境下的稳定性。
3. 闭环验证流程：包含从合成信号、添加噪声、参数求解到信号重构的完整闭环流程，方便用户通过拟合残差评估估计精度。
4. 直观的图形化结果：通过四项对比图表（时域重构、拟合残差、极坐标根分布、功率谱对比）提供详实的分析维度。
5. 自主性强：内置自定义加噪函数，不依赖特定的MATLAB通信或信号处理工具箱即可独立运行。

系统要求

• 运行环境：MATLAB R2016a 或更高版本。
• 硬件要求：标准个人计算机，具备基础的计算性能即可。
• 依赖项：无第三方库依赖。

核心功能逻辑说明

本程序在主逻辑中依次实现了以下关键步骤：

• 模拟信号建模：构建一个频率为50Hz、含负衰减因子（指数减弱）、指定初始相位的单频正弦信号。为了模拟真实环境，程序根据预设信噪比引入了高斯白噪声。

• 线性预测方程构建：利用观测到的带噪信号构建Hankel矩阵。通过建立线性预测模型，将信号样本与其历史值建立联系，从而将非线性估计问题转化为线性最小二乘求解问题，推导出预测系数。

• 特征多项式求根：将预测系数作为多项式系数，在复数域内求解特征方程的根。这些根（极点）直接反映了信号的频率和衰减特性，其在Z平面上的位置决定了信号的稳定性。

• 物理参数映射：根据采样周期，通过复对数运算将Z域的根映射回S域，从中提取出各分量的实时频率（Hz）和衰减因子（Alpha）。

• 复振幅求解：在已知极点的前提下，通过构建范德蒙德（Vandermonde）类型的矩阵，再次运用最小二乘法解出每个指数分量的复振幅，进而换算出信号的物理幅值和初始相位。

• 信号重构与评价：利用估计出的所有参数重新合成时间序列，计算重构信号与理想信号之间的残差，并对比Prony法生成的参数化谱与传统FFT周期图的性能差异。

关键函数与算法细节分析

• 模型阶数选择：程序中设置阶数p等于4。对于单频实信号，理论上仅需2个复指数分量（一对共轭根）即可描述，增加到4阶有助于吸收部分噪声干扰，提升主成分提取的准确度。

• 最小二乘求解：在求解预测系数和复振幅时，程序使用了左除算子（），这在内部执行的是基于矩阵分解的最小二乘拟合，确保了在过定方程组下的最优求解。

• 参数筛选策略：考虑到模型阶数可能高于信号分量数，程序通过寻找具有最大幅值的复分量来确定主频率，有效过滤了由于过建模产生的伪根干扰。

• 功率谱密度计算：不同于FFT直接对时域数据做变换，程序基于估计的预测系数a构建了类似自回归（AR）模型的传递函数，通过频率响应函数计算出平滑的连续谱，展现了参数化模型在频率分辨率上的优势。

• 自定义辅助功能：程序末尾包含一个名为custom_awgn的子函数，它通过能量计算在指定信噪比下生成并叠加分量，保证了代码在基础MATLAB环境下的可移植性。

使用方法

1. 启动MATLAB软件。
2. 将项目文件夹设置为当前工作目录或将其添加到MATLAB搜索路径中。
3. 在命令行窗口直接输入程序名称并回车。
4. 程序将自动运行并在命令行输出估计出的频率、幅值、相位及衰减参数。
5. 同时会弹出图形窗口，展示时域拟合效果、残差分布、Z平面极点轨迹以及功率谱对比。