谱相关密度

✍ dations ◷ 2025-07-11 16:50:00 #谱相关密度

谱相关密度 (spectral correlation density , SCD),有时也称为循环谱密度(cyclic spectral density)或频谱相关函数(spectral correlation function),是描述时间序列的所有频移版本对的交叉频谱密度的函数。谱相关密度仅适用于周期平稳过程,或称为循环平稳过程,普通平稳过程不具备谱相关性。 谱相关被广泛用于信号检测和信号分类。 谱相关密度与每个双线性时频分布密切相关,但不被认为是 Cohen 类分布。

时间序列的循环自相关函数 x ( t ) {textstyle x(t)} 计算如下:

对数字信号而言,SCD 可按任意频率和时间分辨率进行估计。由于直接计算SCD具有较高的计算复杂性,为满足信号实时分析的需求,有几类较为有效的信号谱相关估计方法被提出。

目前常用的算法是 FFT 累加法 (FFT Accumulation Method, FAM) 和带状谱相关法 (Strip-Spectral Correlation Algorithm), 近日,又有一种新的快速频谱相关 (fast-spectral-correlation, FSC) 算法被提出 。

在本节中,我们将介绍实际在计算机上估计SCD的方法。如使用MATLAB或Python中的NumPy库,以下步骤的实现将相当简单。

FFT累加法 (FAM) 是一种计算 SCD 的数字方法。它的输入是一组 IQ 样本矩阵,输出是复值图像(或者说是一复值矩阵),即目标 SCD。FAM输入的信号、或说是 IQ 样本矩阵 x {displaystyle x} ,应为复值张量的形式,或者是尺寸为 ( N , ) {displaystyle (N,)} 的多维数组的形式 ,其中数组中的每个元素都是一个 IQ 样本点。

FAM的第一步,是将所输入信号 x {displaystyle x} 分为多个相互重叠且长度为 N {displaystyle N'} 的数据帧,并将其组合成矩阵形式,记为 X {displaystyle X}

X = x x x . . . ] , {displaystyle X={begin{bmatrix}x\x\x\x\.\.\.end{bmatrix}},}

其中, L {displaystyle L} 两数据帧间起始位置相距的长度。为实现重叠,应有 L < N {displaystyle L<N'} X {displaystyle X} 是形状为 ( P , N ) {displaystyle (P,N')} 的张量, P {displaystyle P} 取决于 x {displaystyle x} 能够容纳多少帧 。

随后,将一形状为 ( N , ) {displaystyle (N',)} 的窗函数 a ( N ) {displaystyle a(N')} ,应用于 X {displaystyle X} 的每一行 (如汉明窗等),得到 X {displaystyle X'}

X = X , {displaystyle X'={begin{bmatrix}a(N')\a(N')\a(N')\.\.\.\end{bmatrix}}otimes X,}

其中 {displaystyle otimes } 是逐元素乘法,也就是将矩阵中的每个元素分别与对应位置的窗函数相乘。接下来,要对中的每一行进行 FFT ,得到 W {displaystyle W}

W = ) F F T ( X ) F F T ( X ) . . . ] . {displaystyle W={begin{bmatrix}FFT(X)\FFT(X)\FFT(X)\.\.\.\end{bmatrix}}.}

W {displaystyle W} 就是通常称为瀑布图或频谱图的矩阵。 FAM 的下一步是校正FFT后数据帧的相位延迟。

W = e j ω L W e j ω 2 L W e j ω 3 L W . . . ] , {displaystyle W'={begin{bmatrix}W\e^{jomega L}otimes W\e^{jomega 2L}otimes W\e^{jomega 3L}otimes W\.\.\.\end{bmatrix}},}

其中 ω {displaystyle omega } 对应于 FFT 结果中的每个数字频率,是形状为 ( N , ) {displaystyle (N',)} 张量。

ω = . {displaystyle omega ={bigg }.}

随后,通过求经 FFT 后结果的自相关,得到形状为 ( P , N , N ) {displaystyle (P,N',N')} 张量 S {displaystyle S}

S = W W , {displaystyle S=W'W'^{*},}

其中 {displaystyle ^{*}} 表示复共轭。换言之,若记 W i = W {displaystyle W_{i}=W'} ( 1 , N ) {displaystyle (1,N')} 的矩阵, S {displaystyle S} 可改写为

S = W i H W i , {displaystyle S=W_{i}^{H}W_{i},}

其中 H 表示矩阵的Hermitian (共轭转置)矩阵。接下来的一步,是将 S {displaystyle S} 沿着第一维分别进行 FFT。

S = F F T ( S ) . {displaystyle S'=FFT(S).}

S {displaystyle S'} 是一个包含完整 SCD 信息的三维张量,但我们的目标是构建形状为 ( P N , N ) {displaystyle (PN',N')} 的二维张量,即矩阵或着图像的形式,张量的两个维度分别对应特定频率 f {displaystyle f} 和循环频率 α {displaystyle alpha } S {displaystyle S'} 中所有 α {displaystyle alpha } 的值可以通过张量 A {displaystyle A} 的到,而所有的频率值 f {displaystyle f} 则记录在张量 F {displaystyle F} 中。这里的 F = j + k 2 N .5 . {displaystyle F={frac {j+k}{2N'}}-.5.}

A = j k N + ( i P 2 ) Δ α . {displaystyle A={frac {j-k}{N'}}+{bigg (}i-{frac {P}{2}}{bigg )}Delta alpha .}

上式中, Δ α = 1 / N {displaystyle Delta alpha =1/N} 。至此,SCD 可以退化为一个二位的图像或矩阵 s {displaystyle s} S {displaystyle S'} 中的 ( f , α ) {displaystyle (f,alpha )} 对都可以赋为0,有效值可以通过 A {displaystyle A} F {displaystyle F} 获取 。

完整计算一次 SCD 具有相当大的复杂度,复杂度的主要来源是第二轮 FFT。幸运的是,从 S {displaystyle S'} 估计 s {displaystyle s'} SCD 的计算公式为

s = 1 P i = 0 P 1 S . {displaystyle s'={frac {1}{P}}sum _{i=0}^{P-1}S'.}

为了更小的计算复杂度,我们可以通过下式,直接从 S {displaystyle S} 计算 s {displaystyle s'} ,因为在 FFT 前或后计算FFT中所有数值的均值是等效的。

s = 1 P i = 0 P 1 S , {displaystyle s'={frac {1}{P}}sum _{i=0}^{P-1}S,}

需要注意的是, s {displaystyle s'} 将看起来像真正 SCD 的 s {displaystyle s} 旋转45 度的版本 。

相关

  • 乌拉尔山乌拉尔山脉(俄语:Ура́льские го́ры,简称Ура́л)是俄罗斯境内大约南北走向的一座山脉,它位于俄罗斯的中西部。山脉位于乌拉尔地理区域内,与乌拉尔联邦管区和乌
  • 搜寻来自近地外智慧生命群落的无线电波计划搜寻来自近地外智慧生命群落的无线电波计划(Search for Extraterrestrial Radio Emissions from Nearby Developed Intelligent Populations,缩写SERENDIP) 是一个搜索地外生命
  • E-6冲印处理E-6冲印处理(E-6 Process,简称E-6)是用作冲印柯达Ektachrome、富士Fujichrome及其他彩色正片(反转片、幻灯片)的一种冲印程序。与一些其他彩色正片冲印程序(如柯达Kodachrome的K-1
  • 黄祖洽黄祖洽(1924年10月2日-2014年9月7日),湖南省长沙市人,中国理论物理和核物理学家。1948年毕业于清华大学,1950年于该校理论物理研究生毕业。1980年当选为中国科学院学部委员。 北京
  • 司马光司马光(1019年11月17日-1086年10月11日),字君实,号迂叟,陕州夏县涑水乡(今山西省夏县)人,北宋文学家、史学家。西晋安平献王司马孚的后裔。世称涑水先生,身后称司马温公。历仕仁宗、英
  • 徐冠巨徐冠巨(1961年-),浙江萧山人,汉族,浙江传化集团有限公司董事长,中华人民共和国企业家、第十一届全国政协委员。父亲徐传化。担任全国工商联常委、浙江省工商联主席、浙江省政协副主
  • 迈克尔·克罗恩-德利迈克尔·克罗恩-德利(丹麦语:Michael Krohn-Dehli,1983年6月6日-)是一名丹麦边锋球员,现效力西乙俱乐部拉科鲁尼亚。克罗恩-德利在2006年10月11日作客瓦都兹对列支敦士登的2008年
  • 大同大学 (日本)大同大学(日语:だいどうだいがく,英语:Daido University),是一所位在日本爱知县名古屋市的私立大学。1939年创校。1964年设立大学部,简称大同大。
  • 唐摭言《唐摭言》,唐末五代南昌人王定保(870年—940年)撰,凡十五卷。《唐摭言》的书名含有谦意,摭字本身是“摘取”之意,是专记唐代科举之笔记小说。全书15卷,分103门,记述大量唐代科举取
  • 穆罗斯 (拉科鲁尼亚省)穆罗斯(加利西亚语:Muros) 是西班牙加利西亚自治区拉科鲁尼亚省的一个沿海港口市镇,下辖8个堂区。市政厅