声光频谱分析仪对瞬变信号参数的提取方法技术

本发明专利技术公开了一种声光频谱分析仪对瞬变信号参数的提取方法，1）入射光经过由被测电信号驱动的声光偏转器后，在透镜的后焦平面上形成一维光点分布，线阵图像传感器按串行帧扫描方式对该光点分布进行读取，输出数据代表光点强度值；2）扫描形成的多帧数据便形成按时间顺序排列的多行数组矩阵A；3）矩阵A中任意j点对应的被测信号频率值为f=f0+j*δf，该信号出现时间点t=t0+iT；4）对矩阵A中的每个数据Aij进行i和j的相关运算即建立起被测信号的瞬态二维时频函数A（t，f）。本发明专利技术拓展了声光频谱分析仪识别信号瞬变信息的功能，数据运算简便且数据量小，多分量信号检测精准，不存在交叉项干扰。

【技术实现步骤摘要】
声光频谱分析仪对瞬变信号参数的提取方法
本专利技术涉及声光频谱分析仪，主要涉及声光频谱分析仪对接收信号的时频幅参数进行有效提取的一种方法。
技术介绍
对常见的不平稳信号，如语音信号、雷达信号以及核医学的图像信号等，它们的频域特性是随时间变化的，人们需要了解某些局部时段上所对应的主要频率特性是什么，也需要了解某些频率的信息出现在哪些时段上。对于这种时-频局部化要求，傅里叶变换是无能为力的，它只能获取一段信号总体上包含哪些频率成分，而对各成分出现的时刻并无所知。因此直接采用傅里叶变换处理，时域相差很大的两个信号，可能频谱图却一样。为得到信号随时间变化的局部化信息，最简单的方法就是对信号进行加窗处理，进行短时傅里叶变换。声光频谱分析仪利用声光调制技术和空间傅里叶变换原理，实现信号的频谱测量。传统的声光频谱分析仪只是对信号积累一段时间的功率谱进行输出，并未对信号的时间瞬变信息进行识别。随着高速光电检测阵列的发展，已经可以实现将图像信号进行准实时输出(特别是图像传感器可直接集成A/D转换进行ns量级的数字信号输出)，由此声光频谱分析仪也可进行短时傅里叶变换处理：通过将被测信号的空域傅里叶变换数据按瞬态时-空关联输出并进行一定运算处理就可以获取信号瞬态时频域全参数信息。
技术实现思路
针对传统声光频谱分析仪只能检测信号频谱，而无时域参数识别功能的不足，本专利技术的目的是提供一种声光频谱分析仪对瞬变信号参数的提取方法——瞬态时空级联算法，本方法结合短时傅里叶变换(STFT)理论对信号的时频幅参数进行实时有效提取。为了实现上述目的，本专利技术采用的技术方案是这样的：声光频谱分析仪对瞬变信号参数的提取方法，其特征在于：步骤如下，1)入射光经过由被测电信号驱动的声光偏转器后，在透镜L的后焦平面P上形成一维光点分布，设计线阵图像传感器按串行帧扫描方式对该光点分布进行读取，帧扫描时间T＝D/v，D为声光偏转器的有效孔径，v为声光介质中的声速；扫描输出的一帧数据代表线阵图像传感器光敏区上依次探测到的光点强度值；线阵图像传感器的光敏区长度应满足系统总工作带宽Δf要求，像素个数应满足系统频率分辨率δf要求，设像素个数为n，则n＝Δf/δf；2)线阵图像传感器输出的一帧数据建立一行数组，反复读出的多帧数据便形成按时间顺序排列的多行数组矩阵A，此为对信号进行时频分析的时-空原始数据矩阵；矩阵中的数据Aij表示第i行、第j列的信号幅值；其中i为行序号,对应时间点；j为列序号，对应频率点；3)任意j点对应的被测信号频率值为f＝f0+j*δf，该信号出现时间点t＝t0+iT；其中，f0为起点频率修正值；t0为起点时间修正值；4)对时-空原始数据矩阵A中的每个数据Aij按步骤3)进行i和j的相关运算即可建立起被测信号的瞬态二维时频函数A(t，f)。相比现有技术，本专利技术具有以下优点：1拓展了声光频谱分析仪识别信号瞬变信息的功能。2数据运算简便且数据量小，有效避免系统数据冗余。3多分量信号检测精准，不存在交叉项干扰。4数据处理模块只需采用常规器件，成本低、功耗低，便于推广应用。附图说明图1-声光频谱分析系统空域傅里叶变换示意图。图2-声光频谱分析系统建立时-空原始数据矩阵示意图。图3-声光频谱分析仪提取信号时频幅参数流程图。图4-声光频谱分析仪实测数据的时频窗口显示图。具体实施方式由声光频谱分析仪的工作原理可知，当输入电信号转换为超声波在介质中建立衍射光栅使入射光波被调制时具有固定孔径时间。超声波陆续通过声光偏转器的有效孔径发生声光互作用，输出衍射光再通过透镜进行空域傅里叶变换，即相当于在时域上对信号进行滑动窗口的短时空域傅立叶变换。所以，后续采用恰当的光电转换时序和数据处理，使光信号空域与电信号时域数据关联输出，便可以利用短时傅立叶变换理论对瞬变信号进行时频幅参数的分析提取。下面结合实施例及附图对本专利技术作进一步详细的描述，但本专利技术的实施方式不限于此。1、分析建立声光频谱分析仪短时空域傅里叶变换的等效数学模型。如图1所示，在声光频谱分析系统中，当入射光入射到由被测信号s(t)＝A(t)cos2πfCt驱动的声光器件时，产生声光互作用。若产生沿声光器件长度x方向的空间调制光分布为f(x，t)，则根据傅里叶光学可知，经过声光器件后的光信号在透镜L的后焦平面P上实现一维傅里叶变换，其傅里叶变换表达式为式中u为空间频率。介质中的超声波陆续通过孔径D(0≤x≤D)发生声光互作用，离开声光器件的一级布喇格衍射光振幅分布可表示为(只需考虑衍射复振幅的相对分布，入射光与光轴倾斜角的影响可不计入表达式)：f+(x,t)＝ma(x)s(t-x/v)(2)式中m是调制度(常数)，v是声光介质中的声速(常数)，a(x)为声光孔径内振幅权函数。将(2)式代入(1)式得众所周知，信号s(t)时间域的短时傅里叶变换定义式为：式中，h(t)为移动窗函数。对比公式(3)和(4)可见运算形式基本相同，差异仅在公式(3)的积分变量为空间窗参数x(被测电信号陆续转换成超声波进入声光介质的有效孔径完成分段空域傅里叶变换)，公式(4)的积分变量为时间窗参数τ(人为将连续信号按时间窗截断进行分段傅里叶运算)。由声光偏转器工作原理可知，空间窗参数x与时间窗参数τ具有正比对应关系：x＝vτ，将其代入公式(3)可得只考虑衍射复振幅的相对分布时，可忽略常数因子v和mv。由此公式(4)和(5)是完全等效的。又声光介质中空间频率u与其驱动电信号(被测信号)的瞬时频率f之间存在正比对应关系：u＝1/Λ＝f/v(Λ为声波波长)。由此，对于驱动电信号之中的任一瞬时频率，声光器件都能产生与之对应的空间频率。所以，声光频谱分析装置实质通过模拟光学技术对被测信号进行了短时傅里叶变换处理，不过此时完成傅里叶变换的信号是以透镜后焦平面P上的光幅分布存在，还需后续的光电转换和数据采集运算配合才能最终完成信号时频幅参数的识别提取。2、以短时空域傅里叶变换理论为基础对变换后的光信号按时序控制读出，并建立时-空原始数据矩阵；由声光频谱分析仪工作原理可知：透镜后焦平面P上的空间频率坐标ξ与空间频率u以及驱动电信号频率f存在如下关系：(λ为入射光波长，F为透镜焦距)。所以可根据检测信号频率范围设计图像传感器的光敏区长度并确定其在焦平面P上的摆放位置。声光器件的有效孔径设计为D时，其渡越时间为T＝D/v，也即系统的短时空域傅里叶变换时间窗为T。设计线阵图像传感器按串行帧扫描模式输出，像素动态自复位，设置帧扫描时间(即各像素积分时间)也为T，此时系统对信号参数的检测精度达到最佳状态，在T时间内线阵中的每个像素依次读出数据，周而复始。如图2所示，按线阵图像传感器输出的一帧数据(对应T段时长内信号空域傅里叶变换的光幅分布)建立一行数组，反复读出的多帧数据便形成按时间顺序排列的多行数组矩阵A，此为对信号进行时频分析的时-空原始数据矩阵(此数据矩阵可以根据应用需求对数据筛选后再建立，判断为背景噪声或不关心的数据可以不保留)。按矩阵位置记录好有用数据Aij，其表示第i行、第j列的信号幅值；其中i为行序号,对应时间点；j为列序号，对应频率点。3、对时-空原始数据矩阵进行信号时-频关联处理，建立信号的瞬态时频函数。系统中，线阵图像传感器的光敏区长度与本文档来自技高网...

【技术保护点】
声光频谱分析仪对瞬变信号参数的提取方法，其特征在于：步骤如下，1）入射光经过由被测电信号驱动的声光偏转器后，在透镜L的后焦平面P上形成一维光点分布，设计线阵图像传感器按串行帧扫描方式对该光点分布进行读取，扫描输出的一帧数据代表线阵图像传感器光敏区上依次探测到的光点强度值；线阵图像传感器的光敏区长度应满足系统总工作带宽Δf要求，像素个数应满足系统频率分辨率δf要求，设像素个数为n，则n=Δf/δf；2）线阵图像传感器输出的一帧数据建立一行数组，反复读出的多帧数据便形成按时间顺序排列的多行数组矩阵A，此为对信号进行时频分析的时‑空原始数据矩阵；矩阵中的数据Aij表示第i行、第j列的信号幅值；其中i为行序号,对应时间点；j为列序号，对应频率点；3）任意j点对应的被测信号频率值为f=f0+j*δf，该信号出现时间点t=t0+iT；其中，f0为起点频率修正值；t0为起点时间修正值；4）对时‑空原始数据矩阵A中的每个数据Aij按步骤3）进行i和j的相关运算即可建立起被测信号的瞬态二维时频函数A（t，f）。

【技术特征摘要】
1.声光频谱分析仪对瞬变信号参数的提取方法，其特征在于：步骤如下，1）入射光经过由被测电信号驱动的声光偏转器后，在透镜L的后焦平面P上形成一维光点分布，设计线阵图像传感器按串行帧扫描方式对该光点分布进行读取，帧扫描时间为T，扫描输出的一帧数据代表线阵图像传感器光敏区上依次探测到的光点强度值；线阵图像传感器的光敏区长度应满足系统总工作带宽Δf要求，像素个数应满足系统频率分辨率δf要求，设像素个数为n，则n=Δf/δf；2）线阵图像传感器输出的一帧数据建立一行数组，反复读出的多帧数据便形成按时间顺序排列的多行数组矩阵A，此为对信号进行时频分析的时-空原始数据矩阵；矩阵...

【专利技术属性】
技术研发人员：吴泽，吴冉，陈永峰，
申请(专利权)人：中国电子科技集团公司第二十六研究所，
类型：发明
国别省市：重庆;85