CN113295923A - 一种基于改进s变换的VFTO信号频谱分析方法 - Google Patents
一种基于改进s变换的VFTO信号频谱分析方法 Download PDFInfo
- Publication number
- CN113295923A CN113295923A CN202110409569.3A CN202110409569A CN113295923A CN 113295923 A CN113295923 A CN 113295923A CN 202110409569 A CN202110409569 A CN 202110409569A CN 113295923 A CN113295923 A CN 113295923A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- signal
- frequency
- window function
- vfto
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000010183 spectrum analysis Methods 0.000 title claims abstract description 19
- 238000001228 spectrum Methods 0.000 claims abstract description 64
- 238000005070 sampling Methods 0.000 claims abstract description 19
- 238000002789 length control Methods 0.000 claims abstract description 12
- 230000009466 transformation Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims abstract description 3
- 230000000737 periodic effect Effects 0.000 claims abstract description 3
- 238000004458 analytical method Methods 0.000 claims description 12
- 229940035637 spectrum-4 Drugs 0.000 claims 1
- 230000015556 catabolic process Effects 0.000 description 11
- 239000002131 composite material Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于改进s变换的VFTO信号频谱分析方法,对输入信号进行时域采样,获得离散序列;对离散序列进行FFT变换,获得信号频谱;根据信号频谱和分辨率要求,确定窗长控制函数以及窗函数;对窗函数进行FFT得到窗函数频谱;对信号频谱进行周期延拓,并将扩维后的信号频谱跟窗函数频谱相乘;对相乘的结果进行傅里叶逆变换,得到单个频率点的时间分布结果;完成所有频率点的计算,最终获得二维矩阵时频谱。本发明,VFTO波形具有频带宽、频率分量幅值跨度大等特性,一种基于改进s变换的VFTO信号频谱分析方法能够较好地适应VFTO频谱特性,并反映频率分量随时间变化的局部特征,更适合于VFTO的频谱分析过程。
Description
技术领域
本发明属于暂态电磁信号处理技术领域,具体涉及一种基于改进s变换的VFTO信号频谱分析方法。
背景技术
在各种快速暂态现象中,由DS操作引起的VFTO的危害最为严重,也最令人关注。VFTO频谱分析依赖于性能良好的信号处理方法,将VFTO信号从时域转换到频域,进而对其频谱特征进行分析。从频域的角度分析,VFTO是一种非平稳信号,其包含的各频率成分的幅值随时间变化。傅里叶变换是一种整体变换,对于信号的表征可以是时域,也可以是频域,但不能表征信号中频率随时间的变化情况,所以只能用于VFTO频谱的粗略分析。
发明内容
本发明的目的在于提供一种基于改进s变换的VFTO信号频谱分析方法,以解决上述问题。
为实现上述目的,本发明采用以下技术方案:
一种基于改进s变换的VFTO信号频谱分析方法,包括以下步骤:
步骤1:对输入信号进行时域采样,获得离散序列;
步骤2:对离散序列进行FFT变换,获得信号频谱;
步骤3:根据信号频谱和分辨率要求,确定窗长控制函数以及窗函数;
步骤4:对窗函数进行FFT得到窗函数频谱;
步骤5:对信号频谱进行周期延拓,并将扩维后的信号频谱跟窗函数频谱相乘;
步骤6:对步骤5相乘的结果进行傅里叶逆变换,得到单个频率点的时间分布结果;
步骤7:对窗函数引入调节因子λ,利用λ替代f,从而通过调节因子λ的作用来控制f的变化速度;
步骤8:重复步骤4-7,直至完成所有频率点的计算,最终获得二维矩阵时频谱。
进一步的,步骤3中,确定窗长控制函数以及窗函数的具体方法为:
步骤3-1、根据信号频谱和分辨率要求,确定参数a、b、c;
步骤3-2、在取值范围内确定a和c的值,将各个参数值代入窗长控制函数中:
步骤3-3、将窗长控制函数带入到高斯窗函数中,得到改进窗函数表达式:
其中,α(f)为窗函数尺度因子。
其中,n从0开始取值,N为频率点个数,T为采样周期。
进一步的,步骤5中,将扩维后的信号频谱跟窗函数频谱相乘,具体为:
进一步的,步骤7中,通过引入调节因子λ对窗函数进行改进,具体表达式为:
改进S变换定义为:
其中λ为调节因子,且λ>0。
进一步的,步骤8中,判断所有频率点是否都已计算具体判断方式是:判断n≥N-1是否成立,若不成立,n加1后重复步骤4、步骤5、步骤6;若成立,输出时频谱结果。
与现有技术相比,本发明有以下技术效果:
本发明根据信号频率自适应调整窗函数,能够实现高分辨率时频分析,具有实时性;
本发明由于该方法不同于WVD变换属于二次型时频变换,不存在交叉项干扰问题;
时频分辨率更佳,能够从时域上区分两段较近的频率分量。综上所述,VFTO波形具有频带宽、频率分量幅值跨度大等特性,一种基于改进s变换的VFTO信号频谱分析方法能够较好地适应VFTO频谱特性,并反映频率分量随时间变化的局部特征,更适合于 VFTO的频谱分析过程。
附图说明
图1为本发明基于改进s变换的VFTO信号频谱分析方法的流程图;
图2为VFTO合成信号波形图;
图3为VFTO合成信号频谱图;
图4为VFTO合成信号S变换灰度图;
图5为VFTO合成信号改进S变换灰度图;
图6为快速隔离开关合闸时负载侧VFTO波形;
图7为快速隔离开关合闸时负载侧首次击穿VFTO波形;
图8为快速隔离开关合闸时负载侧首次击穿VFTO频谱图;
图9为首次预击VFTO波形的改进S变换时频灰度图;
图10为VFTO频率分量的幅值变化曲线。
具体实施方式
以下结合附图对本发明进一步说明:
请参阅图1至图10,一种基于改进s变换的VFTO信号频谱分析方法,包括以下步骤:
步骤3、根据频谱特征,确定参数a、b、c的值,确定窗长控制函数以及窗函数的具体方法为:
步骤3-1、根据信号频谱和分辨率要求,确定参数a、b、c;
步骤3-2、在取值范围内确定a和c的值,将各个参数值代入窗长控制函数中:
步骤3-3、将窗长控制函数带入到高斯窗函数中,得到改进窗函数表达式:
步骤5、将扩维后的信号频谱跟窗函数频谱相乘,具体为:
步骤7、通过引入调节因子λ对窗函数进行改进,具体表达式为:
步骤8、判断所有频率点是否都已计算具体判断方式是:判断n≥N-1是否成立,若不成立,n加1后重复步骤4、步骤5、步骤6;若成立,输出时频谱结果。
下面以三种情况下的实例来说明本发明的时频分析方法。
实施例一:VFTO合成测试信号仿真分析
仿照VFTO波形特点,利用MATLAB软件合成具有单次击穿VFTO波形特征的简单信号。VFTO合成信号各频率分量及持续时间如表1所示,相应波形如图2所示,其中为电压幅值,可看出VFTO合成信号整体上是一个多频率衰减振荡波。VFTO合成信号 Fourier频谱如图3所示,包含有0.00005、1、10、40、80和100MHz共6个频率分量,但Fourier频谱无法给出该信号的局部特征,即无法反映信号的频谱随时间变化的特征。
表1VFTO合成信号频率分量
图4是合成信号的S变换时频分析,可看出其时频分辨率随频率而变化,低频有较好的Δf,高频有较好的Δt,但高频Δf的较差,不能从频域上区分80MHz和100MHz两个分量;图5是合成信号的改进S变换频谱分析,取λ=0.3,其时频分辨率性能主要区别在于使得时频分辨率随频率变化的速度减慢,改进S变换通过λ弥补了S变换在高频上的不足,能从频域上区分80MHz和100MHz两个频率分量。
实施例二:VFTO实测波形分析
图6为GIS变电站快速隔离开关合闸时负载侧的VFTO实测波形。其纵轴为测量电压标幺值pu,基准为峰–峰值电压。
由图6可知,隔离开关在合闸过程中,触头间隙发生多次预击穿,每次击穿形成一个台阶。一方面,因采样间隔极短(1.6ns),图6波形的数据量非常大,若对全部数据进行频谱分析,工作量很大且缺乏针对性;另一方面,单次击穿波形持续时间约为几到几十μs,而相邻两次击穿的时间间隔为级,相邻单次击穿波形不重叠。因此研究VFTO波形频谱特性时,可分解为对各次击穿波形进行频谱分析。
图7为图6中首次击穿前2μs的波形图(20ms处放大波形),可看出单次击穿VFTO 波形是一个多频率衰减振荡波。
由图8可知,波形包含有丰富的频率成分,具体包括0.31、0.94、1.25、4.7、6.9、11.3、 14.1、34.4、42、44及60MHz频率分量,此外还包括直流分量。其中6.9MHz频率分量与其他频率分量相比幅值相对较高,为单次击穿VFTO主要频率分量,对VFTO波形形状影响较大,为主导频率。傅里叶频谱只能从整体上描述VFTO波形所包含的频率成分,无法反映波形频率分量随时间变化的特性。
图9为隔离开关合闸时负载侧首次击穿VFTO波形时频分析结果的时频灰度图。由图9可知,VFTO实测波形具有非平稳特性,其频率成份丰富,且各频率分量持续时间不尽相同。
其中,低频成分自始至终都存在,因其由基本电气振荡而成:高频分量则发速衰减,这是因其产生原因行波的折反射所决定的,分析结果前边所述VFTO的理论分析比较一致。
为深入分析各频率分量的非平稳特性,图10进一步给出了各频率分量的频率切片图,即各频率分量的幅值–时间曲线。
由图10可知,在VFTO波形的2内,直流、310kHz、940kHz及1250kHz等4个低频分量的幅值基本保持不变,而4.7MHz、69MHz、113MH、14.1MHz及344MHz等5个频率分量幅值则会随时间发生衰减,但衰减速度较慢,其中6.9MHz为VFTO主导频率分量,其幅值最大,衰减也相对较慢;而42MHz、44MHz及60MHz3个频率分量的幅值的衰减速度很快,很短时间内衰减至零,可见VFTO高频分量的持续时间非常有限。因此,VFTO频率分量的频率越大、幅值衰减速度越快、持续时间越短。
Claims (9)
1.一种基于改进s变换的VFTO信号频谱分析方法,其特征在于,包括以下步骤:
步骤1:对输入信号进行时域采样,获得离散序列;
步骤2:对离散序列进行FFT变换,获得信号频谱;
步骤3:根据信号频谱和分辨率要求,确定窗长控制函数以及窗函数;
步骤4:对窗函数进行FFT得到窗函数频谱;
步骤5:对信号频谱进行周期延拓,并将扩维后的信号频谱跟窗函数频谱相乘;
步骤6:对步骤5相乘的结果进行傅里叶逆变换,得到单个频率点的时间分布结果;
步骤7:对窗函数引入调节因子λ,利用λ替代f,从而通过调节因子λ的作用来控制f的变化速度;
步骤8:重复步骤4-7,直至完成所有频率点的计算,最终获得二维矩阵时频谱。
9.根据权利要求1所述的一种基于改进s变换的VFTO信号频谱分析方法,其特征在于,步骤8中,判断所有频率点是否都已计算具体判断方式是:判断n≥N-1是否成立,若不成立,n加1后重复步骤4、步骤5、步骤6;若成立,输出时频谱结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110409569.3A CN113295923A (zh) | 2021-04-16 | 2021-04-16 | 一种基于改进s变换的VFTO信号频谱分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110409569.3A CN113295923A (zh) | 2021-04-16 | 2021-04-16 | 一种基于改进s变换的VFTO信号频谱分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113295923A true CN113295923A (zh) | 2021-08-24 |
Family
ID=77319821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110409569.3A Pending CN113295923A (zh) | 2021-04-16 | 2021-04-16 | 一种基于改进s变换的VFTO信号频谱分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113295923A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115420949A (zh) * | 2022-11-04 | 2022-12-02 | 中国电力科学研究院有限公司 | 一种vfto时频分析方法、装置、介质及设备 |
CN116449077A (zh) * | 2023-04-23 | 2023-07-18 | 国网江苏省电力有限公司 | 基于Wigner-Ville分布算法对PT二次侧骚扰电压进行时频分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111198357A (zh) * | 2019-12-19 | 2020-05-26 | 南京理工大学 | 一种基于可调窗函数的s变换时频分析方法 |
CN111257619A (zh) * | 2020-02-17 | 2020-06-09 | 南京工程学院 | 一种基于多属性决策和改进s变换方法的电压暂降检测方法 |
-
2021
- 2021-04-16 CN CN202110409569.3A patent/CN113295923A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111198357A (zh) * | 2019-12-19 | 2020-05-26 | 南京理工大学 | 一种基于可调窗函数的s变换时频分析方法 |
CN111257619A (zh) * | 2020-02-17 | 2020-06-09 | 南京工程学院 | 一种基于多属性决策和改进s变换方法的电压暂降检测方法 |
Non-Patent Citations (1)
Title |
---|
康兵;关向雨;舒乃秋;侯铁信;卜正良;: "几种时频分析方法在VFTO频谱分析中的应用", 高电压技术, no. 02, 28 February 2015 (2015-02-28) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115420949A (zh) * | 2022-11-04 | 2022-12-02 | 中国电力科学研究院有限公司 | 一种vfto时频分析方法、装置、介质及设备 |
CN115420949B (zh) * | 2022-11-04 | 2022-12-30 | 中国电力科学研究院有限公司 | 一种vfto时频分析方法、装置、介质及设备 |
CN116449077A (zh) * | 2023-04-23 | 2023-07-18 | 国网江苏省电力有限公司 | 基于Wigner-Ville分布算法对PT二次侧骚扰电压进行时频分析方法 |
CN116449077B (zh) * | 2023-04-23 | 2024-03-26 | 国网江苏省电力有限公司 | 基于Wigner-Ville分布算法对PT二次侧骚扰电压进行时频分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ohtomo et al. | New method of time series analysis and its application to Wolf's sunspot number data | |
CN113295923A (zh) | 一种基于改进s变换的VFTO信号频谱分析方法 | |
CN109343020B (zh) | 一种基于改进窗函数的s变换时频分析方法 | |
CN107608935B (zh) | 基于时间重排压缩变换的冲击类信号时频分析与重构方法 | |
CN108918964B (zh) | 一种稀疏性增强的谐波分析方法 | |
CN110224394B (zh) | 适用于非平稳功率振荡信号特征提取的傅里叶分解算法 | |
Sedláček et al. | Interpolations in frequency and time domains used in FFT spectrum analysis | |
CN106404386A (zh) | 一种用于采集、提取及诊断齿轮箱早期故障特征信号的方法 | |
CN111198357A (zh) | 一种基于可调窗函数的s变换时频分析方法 | |
CN109682492B (zh) | 基于频域高斯拟合的频率估计方法 | |
CN114881072A (zh) | 一种基于峰值包络谱的傅里叶分解信号降噪方法 | |
CN110687595A (zh) | 一种基于时间重采样和同步挤压变换的地震数据处理方法 | |
Chen et al. | Surface-wave dispersion measurements using Hilbert-Huang transform | |
CN110112757B (zh) | 基于sure小波消噪和改进hht的低频振荡分析方法 | |
CN112328956A (zh) | 一种强频变信号时频分析方法 | |
Abdullah et al. | Bilinear time-frequency analysis techniques for power quality signals | |
CN112507280A (zh) | 信号瞬时频率估计方法 | |
CN113050043A (zh) | 基于非理想lpf的探地雷达超宽带高斯脉冲fri采样方法 | |
Zhang et al. | Noise Reduction of Speech Signal Based on Wavelet Transform with Improved Threshold Function | |
Yu et al. | Performance improvement of wavelet noise reduction based on new threshold function | |
Xu et al. | A noise suppression method for partial discharge based on two-dimensional SVD | |
CN117435864A (zh) | 一种基于分数阶傅里叶变换的多分量线性调频信号高分辨率时频分析方法 | |
Han | The analysis of signal based on the S-transform | |
Zhi-bin | Time frequency analysis of multi-component non-stationary signal with filter bank decomposition | |
CN105631133A (zh) | 基于频域平滑算法的chirp扫频信号处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |