一种基于快速傅里叶变换的相位判读方法技术

本发明专利技术涉及一种基于快速傅里叶变换的相位判读方法，包含：S1、根据采样定理确定采样频率；S2、分别对输入信号和输出信号进行采集，并对采集得到的时间序列数据开展有效性检查；S3、确定用于进行快速傅里叶变换的输入信号和输出信号的时间序列数据长度；S4、分别对选定的输入信号和输出信号进行快速傅里叶变换并提取相位信息，将计算得到的输出信号和输入信号的相位信息相减，得到相位滞后信息。本发明专利技术采用快速傅里叶变换算法判读输出信号和输入信号在目标频率点的相位差，通过对采样频率、输入信号采集数据和输出信号采集数据进行合适的选取及设定，提高目标频率点的相位判读的精度。

【技术实现步骤摘要】
一种基于快速傅里叶变换的相位判读方法
本专利技术涉及一种相位判读方法，具体是指一种基于快速傅里叶变换的硬件系统相位特性判读方法，属于工程
中的试验与测试技术。
技术介绍
现有技术中，对于一些硬件系统，通过对其施加正弦形式的物理激励输入信号或数字激励输入信号，并且比较系统输出信号相对输入信号的相位滞后，从而达到考核系统的相位特性的目的。在传统的时域判读方法中，按下述公式计算输入信号与输出信号之间的相位差单位为度(°)，即：其中，f是输入信号的频率，单位为赫兹(Hz)，Δt是输入正弦信号与输出正弦信号达到极大值(或极小值)的最小时间差。具体如图1所示，图中示出的Δt是输入信号与输出信号达到极大值的最小时间差。由于系统输出信号可能为模拟信号，其中包含有随机噪声。上述传统的时域判读方法中，需判读系统输出信号的极大值(或极小值)的对应时刻，由于受噪声影响较大，导致判读误差较大。如果引入均值滤波对输出的模拟信号进行预处理，仅能提高对其峰峰值的判读精度，却又引入了相位滞后，导致相位误差较大，仍然不能准确反映被测硬件系统的相位滞后特性。进一步，快速傅里叶变换方法被广泛应用于工程
内的时间序列信号的频谱分析。但快速傅里叶变换方法存在截断效应、栅栏效应、信号混叠引起的误差。在不同的应用场合下，采用快速傅里叶变换方法时，如果盲目增大数据长度和采样频率，却没有合理选取采样频率和数据长度等参数，可能导致相位判读结果仍无法满足判读精度的要求。基于上述，目前亟需提出一种基于快速傅里叶变换的相位判读方法，采用快速傅里叶变换算法判读输入、输出信号某一频率点的相位差，通过对输入、输出数据的选取和设定，提高该频率点的相位判读精度。
技术实现思路
本专利技术的目的是提供一种基于快速傅里叶变换的相位判读方法，可计算得到输出信号和输入信号在目标频率点的相位差，通过对采样频率、输入信号采集数据和输出信号采集数据进行合适的选取及设定，提高目标频率点的相位判读的精度。为实现上述目的，本专利技术提供一种基于快速傅里叶变换的相位判读方法，包含以下步骤：S1、根据采样定理确定采样频率，且确定的采样频率必须避免信号混叠引起的误差；S2、根据确定的采样频率分别对输入信号和输出信号进行采集，并对采集得到的时间序列数据开展有效性检查，确保进行相位判读的原始数据有效，否则需重新进行数据采集；S3、确定用于进行快速傅里叶变换的输入信号和输出信号的时间序列数据长度，且确定的时间序列数据长度必须保证快速傅里叶变换的频率分辨率满足要求，避免栅栏效应引起的误差；S4、分别对选定的输入信号和输出信号进行快速傅里叶变换并提取相位信息，将计算得到的输出信号和输入信号的相位信息相减，得到相位滞后信息。所述的S1中，为避免信号混叠引起的误差，确定的采样频率Fs与输入信号频率f满足公式：Fs＞2f。所述的S2中，对时间序列数据的有效性检查包括漏帧检查、等间隔采样检查以及时间基准一致性检查。所述的S2中，具体包含以下步骤：S21、漏帧检查；输入信号和输出信号在时间段[Tj(k)，Tj(k+n-1)]内采集到的数据个数n和采样频率Fs满足公式：其中，j＝1时代表输入信号，j＝2时代表输出信号；ρ表示漏帧率；ρ0表示漏帧率要求；即应使得漏帧率满足要求，否则当漏帧率ρ超过ρ0时，需重新采集数据或选取未发生漏帧的时间段内的数据；S22、等间隔采样检查；每相邻两个采集数据分别对应的时间信号满足公式：其中，j＝1时代表输入信号，j＝2时代表输出信号；即应确保采样频率始终不变，否则需重新采集数据；S23、时间基准一致性检查；输入信号进行数据采集对应的时间段[T1(k)，T1(k+n-1)]与输出信号进行数据采集对应的时间段[T2(k)，T2(k+n-1)]满足公式：T1(k)＝T2(k)；T1(k+n-1)＝T2(k+n-1)；即应选用同一时间段对输入信号和输出信号进行相位判读，否则需重新采集数据。所述的S21中，ρ0是根据相位判读精度的要求预先选取的，其取值范围在5％以内。所述的S3中，为保证快速傅里叶变换的频率分辨率满足要求，避免栅栏效应引起的误差，确定的输入信号和输出信号的时间序列数据长度Nfft与采样频率Fs、输入信号频率f满足公式：其中，N为整数；此时频率分辨率为：目标频率f的频域信息为傅里叶变换第N+1个点的值。所述的S4中，具体包含以下步骤：S41、分别对选定的输入信号和输出信号进行快速傅里叶变换；Yj＝fft(xj，Nfft)；Yj(N+1)＝a+ib；其中，j＝1时代表输入信号，j＝2时代表输出信号；xj为快速傅里叶变换输入的时间序列数据；Yj为复数，是快速傅里叶变换的输出；S42、分别提取输入信号和输出信号在目标频率f处的相位其中，j＝1时代表输入信号，j＝2时代表输出信号；S43、将输出信号在目标频率f处的相位与输入信号在目标频率f处的相位相减，得到相位滞后信息：当时，当时，当时，其中，k为整数，且使得综上所述，本专利技术所提供的基于快速傅里叶变换的相位判读方法，对采样频率、输入信号采集数据和输出信号采集数据进行合适的选取及设定，保证输入信号和输出信号满足使用要求，保证采用快速傅里叶变换方法得到的频率分辨率满足要求，提高目标频率点的相位判读的精度。附图说明图1为系统输出信号相对于输入信号的相位滞后示意图；图2为本专利技术中的基于快速傅里叶变换的相位判读方法的流程图；图3为本专利技术中的数据有效性检查方法的流程图。具体实施方式以下结合图1～图3，详细说明本专利技术的一个优选实施例。如图2所示，为本专利技术所提供的基于快速傅里叶变换的相位判读方法，包含以下步骤：S1、根据采样定理确定采样频率，且确定的采样频率必须避免信号混叠引起的误差；S2、根据确定的采样频率分别对输入信号和输出信号进行采集，并对采集得到的时间序列数据开展有效性检查，确保进行相位判读的原始数据有效，否则需重新进行数据采集；S3、确定用于进行快速傅里叶变换的输入信号和输出信号的时间序列数据长度，且确定的时间序列数据长度必须保证快速傅里叶变换的频率分辨率满足要求，避免栅栏效应引起的误差；S4、分别对选定的输入信号和输出信号进行快速傅里叶变换并提取相位信息，将计算得到的输出信号和输入信号的相位信息相减，得到相位滞后信息。所述的S1中，为避免信号混叠引起的误差，确定的采样频率Fs与输入信号频率f满足公式：Fs＞2f。所述的S2中，对时间序列数据的有效性检查包括漏帧检查、等间隔采样检查以及时间基准一致性检查，以确保进行相位判读的原始数据有效，避免原始数据引入相位判读误差。如图3所示，所述的S2中，具体包含以下步骤：S21、漏帧检查；输入信号和输出信号在时间段[Tj(k)，Tj(k+n-1)]内采集到的数据个数n和采样频率Fs满足公式：其中，j＝1时代表输入信号，j＝2时代表输出信号；ρ表示漏帧率；ρ0表示漏帧率要求；即应使得漏帧率满足要求，否则当漏帧率ρ超过ρ0时，需重新采集数据或选取未发生漏帧的时间段内的数据，避免原始数据引入相位判读误差；S22、等间隔采样检查；每相邻两个采集数据分别对应的时间信号满足公式：其中，j＝1时代表输入信号，j＝2时代表输出信号；即应确保采样频率始终不变，否则需重新采集数据，避免原始数据引入相位判读误差；S2本文档来自技高网...

【技术保护点】
一种基于快速傅里叶变换的相位判读方法，其特征在于，包含以下步骤：S1、根据采样定理确定采样频率，且确定的采样频率必须避免信号混叠引起的误差；S2、根据确定的采样频率分别对输入信号和输出信号进行采集，并对采集得到的时间序列数据开展有效性检查，确保进行相位判读的原始数据有效，否则需重新进行数据采集；S3、确定用于进行快速傅里叶变换的输入信号和输出信号的时间序列数据长度，且确定的时间序列数据长度必须保证快速傅里叶变换的频率分辨率满足要求，避免栅栏效应引起的误差；S4、分别对选定的输入信号和输出信号进行快速傅里叶变换并提取相位信息，将计算得到的输出信号和输入信号的相位信息相减，得到相位滞后信息。

【技术特征摘要】
1.一种基于快速傅里叶变换的相位判读方法，其特征在于，包含以下步骤：S1、根据采样定理确定采样频率，且确定的采样频率必须避免信号混叠引起的误差；S2、根据确定的采样频率分别对输入信号和输出信号进行采集，并对采集得到的时间序列数据开展有效性检查，确保进行相位判读的原始数据有效，否则需重新进行数据采集；S3、确定用于进行快速傅里叶变换的输入信号和输出信号的时间序列数据长度，且确定的时间序列数据长度必须保证快速傅里叶变换的频率分辨率满足要求，避免栅栏效应引起的误差；S4、分别对选定的输入信号和输出信号进行快速傅里叶变换并提取相位信息，将计算得到的输出信号和输入信号的相位信息相减，得到相位滞后信息。2.如权利要求1所述的基于快速傅里叶变换的相位判读方法，其特征在于，所述的S1中，为避免信号混叠引起的误差，确定的采样频率Fs与输入信号频率f满足公式：Fs＞2f。3.如权利要求2所述的基于快速傅里叶变换的相位判读方法，其特征在于，所述的S2中，对时间序列数据的有效性检查包括漏帧检查、等间隔采样检查以及时间基准一致性检查。4.如权利要求3所述的基于快速傅里叶变换的相位判读方法，其特征在于，所述的S2中，具体包含以下步骤：S21、漏帧检查；输入信号和输出信号在时间段[Tj(k)，Tj(k+n-1)]内采集到的数据个数n和采样频率Fs满足公式：其中，j＝1时代表输入信号，j＝2时代表输出信号；ρ表示漏帧率；ρ0表示漏帧率要求；即应使得漏帧率满足要求，否则当漏帧率ρ超过ρ0时，需重新采集数据或选取未发生漏帧的时间段内的数据；S22、等间隔采样检查；每相邻两个采集数据分别对应的时间信号满足公式：

【专利技术属性】
技术研发人员：叶赛仙，奚勇，冯昊，陈光山，廖幻年，
申请(专利权)人：上海航天控制技术研究所，
类型：发明
国别省市：上海,31