基于FFT累积量的循环平稳相关熵谱密度计算方法

摘要

本发明公开了一种基于FFT累积量的循环平稳相关熵谱密度计算方法，包括以下步骤：给定信号x、核长σ、信号采样频率fs、循环频率分辨率Δα、谱频率分辨率Δf；将信号x进行分段，得到各分段信号xxi；计算各分段信号xxi的中心化相关熵Exx；创建中心化相关熵矩阵X；对X进行加窗处理，得到矩阵XW；按列计算矩阵XW的FFT，得到矩阵XF1及其共轭矩阵CXF1；将计算矩阵XF1的每一行向量的元素与CXF1的每一行向量的元素相乘，得到矩阵XM；按行计算矩阵XM的FFT并转置，得到循环平稳相关熵谱密度S。本发明提供了一种计算循环平稳相关熵谱密度的高效算法，提高了循环平稳相关熵谱密度的计算速度和频谱分辨率，实验证明本发明算法性能良好。