声光频谱分析仪对瞬变信号参数的提取方法
【技术领域】
[0001] 本发明涉及声光频谱分析仪，主要涉及声光频谱分析仪对接收信号的时频幅参数 进行有效提取的一种方法。
【背景技术】
[0002] 对常见的不平稳信号，如语音信号、雷达信号以及核医学的图像信号等，它们的频 域特性是随时间变化的，人们需要了解某些局部时段上所对应的主要频率特性是什么，也 需要了解某些频率的信息出现在哪些时段上。对于这种时-频局部化要求，傅里叶变换是 无能为力的，它只能获取一段信号总体上包含哪些频率成分，而对各成分出现的时刻并无 所知。因此直接采用傅里叶变换处理，时域相差很大的两个信号，可能频谱图却一样。为得 到信号随时间变化的局部化信息，最简单的方法就是对信号进行加窗处理，进行短时傅里 叶变换。
[0003] 声光频谱分析仪利用声光调制技术和空间傅里叶变换原理，实现信号的频谱测 量。传统的声光频谱分析仪只是对信号积累一段时间的功率谱进行输出，并未对信号的时 间瞬变信息进行识别。随着高速光电检测阵列的发展，已经可以实现将图像信号进行准实 时输出（特别是图像传感器可直接集成A/D转换进行ns量级的数字信号输出），由此声光 频谱分析仪也可进行短时傅里叶变换处理：通过将被测信号的空域傅里叶变换数据按瞬态 时-空关联输出并进行一定运算处理就可以获取信号瞬态时频域全参数信息。
【发明内容】
[0004] 针对传统声光频谱分析仪只能检测信号频谱，而无时域参数识别功能的不足，本 发明的目的是提供一种声光频谱分析仪对瞬变信号参数的提取方法一一瞬态时空级联算 法，本方法结合短时傅里叶变换（STFT)理论对信号的时频幅参数进行实时有效提取。
[0005] 为了实现上述目的，本发明采用的技术方案是这样的：
[0006] 声光频谱分析仪对瞬变信号参数的提取方法，其特征在于：步骤如下，
[0007] 1)入射光经过由被测电信号驱动的声光偏转器后，在透镜L的后焦平面P上形成 一维光点分布，设计线阵图像传感器按串行帧扫描方式对该光点分布进行读取，帧扫描时 间T=D/v，D为声光偏转器的有效孔径，v为声光介质中的声速；扫描输出的一帧数据代表 线阵图像传感器光敏区上依次探测到的光点强度值；线阵图像传感器的光敏区长度应满足 系统总工作带宽Af?要求，像素个数应满足系统频率分辨率Sf?要求，设像素个数为n，则n =Af/ 8f；
[0008] 2)线阵图像传感器输出的一帧数据建立一行数组，反复读出的多帧数据便形成按 时间顺序排列的多行数组矩阵A，此为对信号进行时频分析的时-空原始数据矩阵；矩阵中 的数据'表示第i行、第j列的信号幅值；其中i为行序号，对应时间点；j为列序号，对 应频率点；
[0009] 3)任意j点对应的被测信号频率值为f=心+j*Sf，该信号出现时间点t=tfiT; 其中，fc为起点频率修正值；t。为起点时间修正值；
[0010] 4)对时-空原始数据矩阵A中的每个数据Ax]按步骤3)进行i和j的相关运算 即可建立起被测信号的瞬态二维时频函数A(t，f)。
[0011] 相比现有技术，本发明具有以下优点：
[0012] 1拓展了声光频谱分析仪识别信号瞬变信息的功能。
[0013] 2数据运算简便且数据量小，有效避免系统数据冗余。
[0014] 3多分量信号检测精准，不存在交叉项干扰。
[0015] 4数据处理模块只需采用常规器件，成本低、功耗低，便于推广应用。
【附图说明】
[0016]图1-声光频谱分析系统空域傅里叶变换示意图。
[0017]图2-声光频谱分析系统建立时-空原始数据矩阵示意图。
[0018] 图3-声光频谱分析仪提取信号时频幅参数流程图。
[0019] 图4-声光频谱分析仪实测数据的时频窗口显示图。
【具体实施方式】
[0020] 由声光频谱分析仪的工作原理可知，当输入电信号转换为超声波在介质中建立衍 射光栅使入射光波被调制时具有固定孔径时间。超声波陆续通过声光偏转器的有效孔径发 生声光互作用，输出衍射光再通过透镜进行空域傅里叶变换，即相当于在时域上对信号进 行滑动窗口的短时空域傅立叶变换。所以，后续采用恰当的光电转换时序和数据处理，使光 信号空域与电信号时域数据关联输出，便可以利用短时傅立叶变换理论对瞬变信号进行时 频幅参数的分析提取。
[0021] 下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限 于此。
[0022] 1、分析建立声光频谱分析仪短时空域傅里叶变换的等效数学模型。
[0023] 如图1所示，在声光频谱分析系统中，当入射光入射到由被测信号s(t) =A(t) C〇s2 31fet驱动的声光器件时，产生声光互作用。若产生沿声光器件长度x方向的空间调制 光分布为f(x，t)，则根据傅里叶光学可知，经过声光器件后的光信号在透镜L的后焦平面P 上实现一维傅里叶变换，其傅里叶变换表达式为
[0024]
(13
[0025] 式中u为空间频率。介质中的超声波陆续通过孔径D(0<x<D)发生声光互作 用，离开声光器件的一级布喇格衍射光振幅分布可表示为（只需考虑衍射复振幅的相对分 布，入射光与光轴倾斜角的影响可不计入表达式）：
[0026]f+(x,t)=ma(x)s(t-x/v) (2)
[0027] 式中m是调制度（常数），v是声光介质中的声速（常数），a(x)为声光孔径内振 幅权函数。将（2)式代入（1)式得
[0028]
[0029] 众所周知，信号s(t)时间域的短时傅里叶变换定义式为：
[0030]
[0031] 式中，h(t)为移动窗函数。
[0032] 对比公式（3)和（4)可见运算形式基本相同，差异仅在公式（3)的积分变量为空 间窗参数x(被测电信号陆续转换成超声波进入声光介质的有效孔径完成分段空域傅里叶 变换），公式（4)的积分变量为时间窗参数t(人为将连续信号按时间窗截断进行分段傅 里叶运算）。由声光偏转器工作原理可知，空间窗参数x与时间窗参数t具有正比对应关 系：x=vt，将其代入公式⑶可得
[0033]
[0034] 只考虑衍射复振幅的相对分布时，