CN103454495B - 自适应高精度快速频谱分析方法 - Google Patents
自适应高精度快速频谱分析方法 Download PDFInfo
- Publication number
- CN103454495B CN103454495B CN201310418656.0A CN201310418656A CN103454495B CN 103454495 B CN103454495 B CN 103454495B CN 201310418656 A CN201310418656 A CN 201310418656A CN 103454495 B CN103454495 B CN 103454495B
- Authority
- CN
- China
- Prior art keywords
- frequency
- correction
- spectrum
- spectral line
- amplitude
- 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.)
- Active
Links
- 230000003595 spectral effect Effects 0.000 title claims abstract description 87
- 238000012937 correction Methods 0.000 claims abstract description 83
- 238000001228 spectrum Methods 0.000 claims abstract description 65
- 238000012545 processing Methods 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 44
- 238000001914 filtration Methods 0.000 claims abstract description 29
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 24
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 238000004458 analytical method Methods 0.000 claims description 25
- 238000005070 sampling Methods 0.000 claims description 14
- 238000007781 pre-processing Methods 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
Abstract
本发明涉及数字信号处理技术。本发明解决了现有频谱分析方法误差较大精度不佳的问题,提供了一种自适应高精度快速频谱分析方法,其技术方案可概括为:首先采集模拟信号,然后进行预处理滤波,再进行模数转换,得到采样数据,然后选取28个数据进行FFT处理,得到初步的信号频谱信息,再进行峰值搜索,然后再对初步的信号频谱信息进行谱线混叠识别,若未混叠则对初步的信号频谱信息进行频谱校正处理,得到信号的频率、幅值及初相信息,若混叠则对初步的信号频谱信息进行复解析带通滤波细化处理,再进行频谱校正处理,分别得到信号的频率、幅值及初相信息。本发明的有益效果是,在保障计算精度的前提下,大幅度降低了计算量,适用于频谱分析。
Description
技术领域
本发明涉及数字信号处理技术,特别涉及数字信号的频谱分析技术。
背景技术
信号的频谱分析一直是数字信号处理的热点之一,获取信号的频谱信息,是诸多领域的关键点所在。从1965年Cooely-Tukey在《计算数学》(MathematicsofComputation)杂志上首次发表了FFT(FastFourierTransform,快速傅里叶变换)算法,FFT和频谱分析很快发展成为了机械设备故障诊断、电力***、无线电通信和信息图像处理等多种学科重要的理论基础,是一种应用极为广泛的动态信号处理方法。在FFT算法的基础上,衍生出了很多处理方法,诸如ZoomFFT(Zoom-FastFourierTransform,变焦快速傅里叶变换)、MEM(MaximumEntropyMethod,最大熵分析)等频谱分析手段。
MOROZOVJA、DongHui、谭思炜、温和等人在理论上***地分析了FFT及改进FFT在离散频谱分析在频率,幅值,初相均存在较大误差,且对于频率间距小于五个频谱分辨率的频谱分析存在严重谱线混叠问题。夏均忠、LiangkaiLiu、KuiWang、WANGK实现并且改进了ZoomFFT算法,虽然能解决密集频谱分析存在的局部谱线混叠问题,但是处理结果在频率,幅值,初相同样存在较大误差,并且无法处理频率间距较大和频率间距很小并存的频谱分析问题。同样AbheyRamBansa、柴山、王俊等人提出分析的MEM,在短序列处理上有很大优势,但是在精度方面仍有欠缺。
发明内容
本发明的目的是要克服目前频谱分析方法误差较大精度不佳的缺点,提供一种自适应高精度快速频谱分析方法。
本发明解决其技术问题,采用的技术方案是,自适应高精度快速频谱分析方法,其特征在于,包括如下步骤:
步骤1、采集需处理的模拟信号;
步骤2、对采集的模拟信号进行预处理滤波;
步骤3、对预处理滤波后的模拟信号进行模数转换,得到采样数据;
步骤4、从采样数据中选取28个数据进行FFT处理,得到初步的信号频谱信息;
步骤5、对初步的信号频谱信息进行峰值搜索,得到峰值;
步骤6、对初步的信号频谱信息进行谱线混叠识别,判断谱线是否被混叠,若是则进入步骤7,若不是则进入步骤9;
步骤7、对初步的信号频谱信息进行复解析带通滤波细化处理;
步骤8、对进行了复解析带通滤波细化处理的信号进行频谱校正处理,分别得到信号的频率、幅值及初相信息;
步骤9、对初步的信号频谱信息进行频谱校正处理,得到信号的频率、幅值及初相信息。
具体的,步骤5中,所述峰值处理搜索的方法为:设Y(n)为n点FFT处理后的谱线幅值序列,设第i条谱线对应于峰值,则搜索公式为:Y(i)>Y(i-1)&Y(i)>Y(i+1)&Y(i)>A,其中,A为幅值处理精度要求,具体为大于0的常数。
进一步的,步骤6中,所述谱线混叠识别的方法为:
步骤61、选择靠近峰值的三条谱线;
步骤62、对该三条谱线分别进行比值法校正,判断校正结果是否相同,若相同则认为无谱线混叠,进入步骤9,若不相同则认为谱线混叠,识别出谱线混叠的谱线号,进入步骤7。
再进一步的,步骤62中,所述比值法校正是指:根据谐波信号加hanning窗进行离散频谱校正的重心定理,若频率不完全相同,则为发生谱线混叠,此时任意两条谱线不完全重叠,则从代表频率大小的轴的从左到右或者从右到左方向来进行校正,校正结果不同,此时即为谱线发生混叠,识别出其谱线号,若校正结果相同,则无谱线混叠。
具体的,步骤7中,所述进行复解析带通滤波细化处理的方法为:
步骤71、根据识别出来的存在谱线混叠的谱线号,记为k,得到需要细化的频现段,记为[fl,fh];
步骤72、针对D=128,得到细化方案为:[fl=(k-1)▽f,fh=(k+1)▽f],其中,fs指代信号采样频率,N为使用到的信号个数;
步骤73、根据细化方案,得到该频现段的频率中心为:fe=(fl+fh)/2,对其归一化为
步骤74、确定低通数字滤波器的截止频率为fc=fs/2D,对其归一化得:
步骤75、通过窗函数法设计线性相位FIR低通数字滤波器为:a=(M-1)/2,
其中,c为小于等于M的正整数,M为数字滤波器阶数,a为数字滤波器半阶数;
步骤76、设计复解析带通滤波器为:其中,we为细化区间归一化频率;
步骤77、用复解析带通滤波器对获得的采样信号进行选抽滤波,其具体公式为:
其中,b为计算变量,x为采集信号,t为小于等于N的正整数,N为使用到的信号个数;
步骤78、选抽滤波后进行复调制频移,采用如下公式:
其中,
步骤79、对得到的复调制后的序列进行加hanning窗N点FFT运算,得到复解析带通细化的谱线信息,针对谱线号为s的谱线,其未进行校正处理的真实频率为
再进一步的,步骤8中,所述频谱校正处理包括如下步骤:
步骤81、设yk为某一主瓣内的最大值,k为该最大值对应的谱线号,Δf1为最大离散采样值的频率和该主瓣对称中心处频率的差值,W1(Δf1)为归一化变量之后的窗函数幅值函数;
步骤82、采用频率校正公式、幅值校正公式及相位校正公式进行频谱校正处理,所述频率校正公式为:f0=(k-Δf1)▽f,
所述幅值校正公式为:
所述相位校正公式为:
具体的,步骤9中,所述频谱校正处理包括如下步骤:
步骤91、设yk为某一主瓣内的最大值,k为该最大值对应的谱线号,Δf1为最大离散采样值的频率和该主瓣对称中心处频率的差值,W1(Δf1)为归一化变量之后的窗函数幅值函数;
步骤92、采用频率校正公式、幅值校正公式及相位校正公式进行频谱校正处理,所述频率校正公式为:f0=(k-Δf1)▽f,
所述幅值校正公式为:
所述相位校正公式为:
本发明的有益效果是,在本发明方案中的自适应高精度快速频谱分析方法,能够在不提高频率分辨率的情况下,自动识别出是否存在频率混叠,自动对存在频率交叠区域进行频谱细化,自动校正,在保障计算精度的前提下,大幅度降低了计算量。
附图说明
图1为本发明自适应高精度快速频谱分析方法的流程图。
具体实施方式
下面结合实施例及附图,详细描述本发明的技术方案。
本发明的自适应高精度快速频谱分析方法为:首先采集需处理的模拟信号,然后对采集的模拟信号进行预处理滤波,再对预处理滤波后的模拟信号进行模数转换,得到采样数据,然后从采样数据中选取28个数据进行FFT处理,得到初步的信号频谱信息,再对初步的信号频谱信息进行峰值搜索,得到峰值,然后再对初步的信号频谱信息进行谱线混叠识别,判断谱线是否被混叠,若不是则对初步的信号频谱信息进行频谱校正处理,得到信号的频率、幅值及初相信息,若是则对初步的信号频谱信息进行复解析带通滤波细化处理,并对进行了复解析带通滤波细化处理的信号进行频谱校正处理,分别得到信号的频率、幅值及初相信息。
实施例
本发明实施例的自适应高精度快速频谱分析方法为:首先采集需处理的模拟信号,然后对采集的模拟信号进行预处理滤波,再对预处理滤波后的模拟信号进行模数转换,得到采样数据,然后从采样数据中选取28个数据进行FFT处理,得到初步的信号频谱信息,再对初步的信号频谱信息进行峰值搜索,得到峰值,然后再对初步的信号频谱信息进行谱线混叠识别,判断谱线是否被混叠,若不是则对初步的信号频谱信息进行频谱校正处理,得到信号的频率、幅值及初相信息,若是则对初步的信号频谱信息进行复解析带通滤波细化处理,并对进行了复解析带通滤波细化处理的信号进行频谱校正处理,分别得到信号的频率、幅值及初相信息。
其中,本方法的核心在于峰值搜索、谱线混叠识别、复解析带通滤波细化及频谱校正,其具体如下:
1、峰值搜索:
峰值搜索的目的是找到FFT处理后的频谱的峰值,为后面的谱线混叠识别及频谱校正做准备,设Y(n)为n点FFT处理后的谱线幅值序列,设第i条谱线对应于峰值,则搜索公式为:Y(i)>Y(i-1)&Y(i)>Y(i+1)。对于窗函数,主瓣区间有若干条谱线,实际信号处理中,广泛被使用的窗函数为hanning窗,本例仅对hanning窗加以说明,通过hanning窗性质可知,其主瓣宽度是四个频率分辨率▽f,其中,fs指代信号采样频率,N为使用到的信号个数,因而hanning窗与单频信号卷积后,在主瓣内最多有四条谱线。同时由于窗函数旁瓣的影响,卷积处理之后,会在主瓣边缘有实际并不存在的旁瓣峰值,同时为了应对在实际处理的信号中混有随机噪声,所以在峰值搜索时,增加一个常数A,A为大于0的常数,以排除旁瓣峰值和随机噪声导致的虚假峰值,增加不必要的搜索,所以,搜索公式为:Y(i)>Y(i-1)&Y(i)>Y(i+1)&Y(i)>A,其中,A为幅值处理精度要求,具体为大于0的常数。
2、谱线混叠识别:
(1)、选择靠近峰值的三条谱线。
(2)、对该三条谱线分别进行比值法校正,判断校正结果是否相同,若相同则认为无谱线混叠,若不相同则认为谱线混叠,识别出谱线混叠的谱线号。
本步骤中,比值法校正是指:所述比值法校正是指:根据谐波信号加hanning窗进行离散频谱校正的重心定理,若频率不完全相同,则为发生谱线混叠,此时任意两条谱线不完全重叠,则从代表频率大小的轴的从左到右或者从右到左方向来进行校正,校正结果不同,此时即为谱线发生混叠,识别出其谱线号,若校正结果相同,则无谱线混叠。
3、复解析带通滤波细化处理:
(1)、根据识别出来的存在谱线混叠的谱线号,记为k,得到需要细化的频现段,记为[fl,fh]。
(2)、针对D=128,得到细化方案为:[fl=(k-1)▽f,fh=(k+1)▽f]。
(3)、根据细化方案,得到该频现段的频率中心为:fe=(fl+fh)/2,对其归一化为
(4)、确定低通数字滤波器的截止频率为fc=fs/2D,对其归一化得:
(5)、设计线性相位FIR低通数字滤波器为:a=(M-1)/2,
其中,c为小于等于M的正整数,M为数字滤波器阶数,a为数字滤波器半阶数;
步骤76、设计复解析带通滤波器为:其中,we为细化区间归一化频率;
步骤77、用复解析带通滤波器对获得的采样信号进行选抽滤波,其具体公式为:
其中,b为计算变量,x为采集信号,t为小于等于N的正整数,N为使用到的信号个数;
步骤78、选抽滤波后进行复调制频移,采用如下公式:
其中,
步骤79、对得到的复调制后的序列进行加hanning窗N点FFT运算,得到复解析带通细化的谱线信息,针对谱线号为s的谱线,其未进行校正处理的真实频率为
3、频谱校正处理:
频谱校正处理分为无混叠频谱校正处理及细化频谱校正处理。
无混叠频谱校正处理:
在时域对信号进行截断,实际上是把原来的信号频谱函数和窗函数频谱函数进行了卷积,因而得到的结果是卷积后的数据,而不是原始的频谱数据。再者,计算机在进行数据处理时,不可能采用连续的变量值,而只可能是采用密集离散的自变量数据,进而导致得到频谱信息进一步变成了前面所叙的原始频谱函数与窗函数频谱函数的卷积后的离散采样值,(而采样的精细程度看频率分辨率),同时在读取信息的时候,是选择针对的主瓣内最大的离散采样值的频率是信号频率的,其大小对应于该频率的特征信息的,包括幅值和相位。这样就有可能导致误差,即如果选择的离散点并不是对窗函数的对称中心的离散采样值,那么就引入了误差。
通过上面的误差来源分析,可以知道,如果能精确得求解出作为某一主瓣内的最大离散采样值的频率和该主瓣对称中心处频率的差值Δf1,然后在结合窗函数的具体表达式,就能求解得到该主瓣对应的真实频率,进而精确求解出幅值,相位等信息。
具体步骤如下:
(1)、设yk为某一主瓣内的最大值,k为该最大值对应的谱线号,用频率归一化之后差值为1的主瓣峰顶附近的两条谱线(可通过分析窗函数的性质得到)的窗函数比值为:其中,W1(Δf1)为归一化变量之后的窗函数幅值函数;这样就可以精确的求解出采用窗函数截断,频域离散化之后得到的频率与真实理论上频率之间的差值。
现针对hanning窗进行数学分析,hanning窗的序列表达式为:
同时通过理论分析,可以得到其理论上推导的频域幅值函数为:
(2)、设yk+1为主瓣内离散次最大值,则:
(3)、计算得到最后频率误差公式:
yk+1≥yk-1指离散次最大值在最大值后面,yk+1≤yk-1指离散次最大值在最大值前面。
所以某一主瓣内的真实频率为:f0=(k-Δf1)▽f。设A为原来单频信号的幅值,则有:
因此得到原始谱线对应的真实频谱信息如下:
频率校正公式为:f0=(k-Δf1)▽f,
幅值校正公式为:
相位校正公式为:
细化频谱校正处理:
(1)、设yk为某一主瓣内的最大值,k为该最大值对应的谱线号,用频率归一化之后差值为1的主瓣峰顶附近的两条谱线(可通过分析窗函数的性质得到)的窗函数比值为:其中,W1(Δf1)为归一化变量之后的窗函数幅值函数;这样就可以精确的求解出采用窗函数截断,频域离散化之后得到的频率与真实理论上频率之间的差值。
现针对hanning窗进行数学分析,hanning窗的序列表达式为:
同时通过理论分析,可以得到其理论上推导的频域幅值函数为:
(2)、设yk+1为主瓣内离散次最大值,则:
(3)、计算得到最后频率误差公式:
yk+1≥yk-1指离散次最大值在最大值后面,yk+1≤yk-1指离散次最大值在最大值前面。。
(4)、采用窗函数设计方法设计低通滤波器,具体为:h′(c)=h(c)×w(c),其中,c为小于等于M的正整数,得到该窗函数设计的低通滤波器的频谱函数为:
进而可得:H′(ejω)=H′(ω)ejθ(ω)
(5)、设计复解析带通滤波器,其由低通滤波器乘以所得,即:
(6)、由傅里叶变换的频移性质可得:
(7)、由于数字滤波丢失了在前的M-1个数据,通过傅里叶变换的性质可得结果经过复解析带通滤波器滤波之后,信号频谱函数变为:
因此得:
进行复调制后得:
最后的频谱信息为:
由此可知,经过复解析带通滤波处理的序列,其幅值信息是不会变的,而相位信息发生了变化,所以还可以用上面的校正方法进行校正,对于频率和幅值,其校正公式不变,但是相位校正公式发生了变化,通过进一步的分析,可得如下具体公式:
频率校正公式为:f0=(k-Δf1)▽f,
幅值校正公式为:
相位校正公式为:
Claims (6)
1.自适应高精度快速频谱分析方法,其特征在于,包括如下步骤:
步骤1、采集需处理的模拟信号;
步骤2、对采集的模拟信号进行预处理滤波;
步骤3、对预处理滤波后的模拟信号进行模数转换,得到采样数据;
步骤4、从采样数据中选取28个数据进行FFT处理,得到初步的信号频谱信息;
步骤5、对初步的信号频谱信息进行峰值搜索,得到峰值,所述峰值搜索的方法为:设Y(n)为n点FFT处理后的谱线幅值序列,设第i条谱线对应于峰值,则搜索公式为:Y(i)>Y(i-1)&Y(i)>Y(i+1)&Y(i)>A,其中,A为幅值处理精度要求,具体为大于0的常数;
步骤6、对初步的信号频谱信息进行谱线混叠识别,判断谱线是否被混叠,若是则进入步骤7,若不是则进入步骤9;
步骤7、对初步的信号频谱信息进行复解析带通滤波细化处理;
步骤8、对进行了复解析带通滤波细化处理的信号进行频谱校正处理,分别得到信号的频率、幅值及初相信息;
步骤9、对初步的信号频谱信息进行频谱校正处理,得到信号的频率、幅值及初相信息。
2.如权利要求1所述的自适应高精度快速频谱分析方法,其特征在于,步骤6中,所述谱线混叠识别的方法为:
步骤61、选择靠近峰值的三条谱线;
步骤62、对该三条谱线分别进行比值法校正,判断校正结果是否相同,若相同则认为无谱线混叠,进入步骤9,若不相同则认为谱线混叠,识别出谱线混叠的谱线号,进入步骤7。
3.如权利要求2所述的自适应高精度快速频谱分析方法,其特征在于,步骤62中,所述比值法校正是指:根据谐波信号加hanning窗进行离散频谱校正的重心定理,若频率不完全相同,则为发生谱线混叠,此时任意两条谱线不完全重叠,则从代表频率大小的轴的从左到右或者从右到左方向来进行校正,校正结果不同,此时即为谱线发生混叠,识别出其谱线号,若校正结果相同,则无谱线混叠。
4.如权利要求1所述的自适应高精度快速频谱分析方法,其特征在于,步骤7中,所述进行复解析带通滤波细化处理的方法为:
步骤71、根据识别出来的存在谱线混叠的谱线号,记为k,得到需要细化的频现段,记为[fl,fh];
步骤72、针对D=128,得到细化方案为:[fl=(k-1)▽f,fh=(k+1)▽f],其中,fs指代信号采样频率,N为使用到的信号个数;
步骤73、根据细化方案,得到该频现段的频率中心为:fe=(fl+fh)/2,对其归一化为
步骤74、确定低通数字滤波器的截止频率为fc=fs/2D,对其归一化得:
步骤75、通过窗函数法设计线性相位FIR低通数字滤波器为:a=(M-1)/2,
其中,c为小于等于M的正整数,M为数字滤波器阶数,a为数字滤波器半阶数;
步骤76、设计复解析带通滤波器为:其中,we为细化区间归一化频率;
步骤77、用复解析带通滤波器对获得的采样信号进行选抽滤波,其具体公式为:
其中,b为计算变量,x为采集信号,t为小于等于N的正整数;
步骤78、选抽滤波后进行复调制频移,采用如下公式:
其中,
步骤79、对得到的复调制频移后的序列进行加hanning窗N点FFT运算,得到复解析带通细化的谱线信息,针对谱线号为s的谱线,其未进行校正处理的真实频率为
5.如权利要求4所述的自适应高精度快速频谱分析方法,其特征在于,步骤8中,所述频谱校正处理包括如下步骤:
步骤81、设yk为某一主瓣内的最大值,k为该最大值对应的谱线号,Δf1为最大离散采样值的频率和该主瓣对称中心处频率的差值,W1(Δf1)为归一化变量之后的窗函数幅值函数;
步骤82、采用频率校正公式、幅值校正公式及相位校正公式进行频谱校正处理,所述频率校正公式为:
所述幅值校正公式为:
所述相位校正公式为:
6.如权利要求5所述的自适应高精度快速频谱分析方法,其特征在于,步骤9中,所述频谱校正处理包括如下步骤:
步骤91、设yk为某一主瓣内的最大值,k为该最大值对应的谱线号,Δf1为最大离散采样值的频率和该主瓣对称中心处频率的差值,W1(Δf1)为归一化变量之后的窗函数幅值函数;
步骤92、采用频率校正公式、幅值校正公式及相位校正公式进行频谱校正处理,所述频率校正公式为:
所述幅值校正公式为:
所述相位校正公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310418656.0A CN103454495B (zh) | 2013-09-13 | 2013-09-13 | 自适应高精度快速频谱分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310418656.0A CN103454495B (zh) | 2013-09-13 | 2013-09-13 | 自适应高精度快速频谱分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103454495A CN103454495A (zh) | 2013-12-18 |
CN103454495B true CN103454495B (zh) | 2016-01-20 |
Family
ID=49737086
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310418656.0A Active CN103454495B (zh) | 2013-09-13 | 2013-09-13 | 自适应高精度快速频谱分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103454495B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110068730A (zh) * | 2019-05-07 | 2019-07-30 | 中国科学院电子学研究所 | 双频段频谱数据采集方法及装置 |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375111B (zh) * | 2014-11-16 | 2017-03-29 | 甘肃省机械科学研究院 | 对密集频谱进行快速高精度细化校正的方法 |
CN105004913A (zh) * | 2015-07-21 | 2015-10-28 | 黎文安 | 一种电力***交流采样频率实时跟踪方法 |
CN105510812A (zh) * | 2015-10-12 | 2016-04-20 | 国家电网公司 | 一种换流变分接开关状态监测***及非平稳信号频谱估计方法 |
CN105865616B (zh) * | 2016-03-31 | 2018-10-12 | 湖南科技大学 | 基于fft的调制谱快速细化方法 |
CN106018955B (zh) * | 2016-05-07 | 2018-11-13 | 天津大学 | 快速卷积可调滤波器组的低速率端频率估计方法 |
CN107645364B (zh) * | 2016-07-22 | 2022-12-09 | 山东壹号方舟机器人智能科技有限公司 | 互补编码方法及装置、互补译码方法及装置、OvXDM*** |
CN108241091B (zh) * | 2016-12-27 | 2021-09-03 | 北京普源精电科技有限公司 | 利用频谱仪进行2fsk信号峰值搜索的方法及频谱仪 |
CN106855787A (zh) * | 2016-12-30 | 2017-06-16 | 深圳市易联网络技术有限公司 | 一种物联网模拟信号的采样方法 |
CN106990288A (zh) * | 2017-03-20 | 2017-07-28 | 成都米风通信技术有限公司 | 基于stm32f4的高速信号频谱分析方法 |
CN108334682B (zh) * | 2018-01-24 | 2021-08-17 | 北京机械设备研究所 | 一种频谱自动细化的批处理方法 |
CN108761117B (zh) * | 2018-07-27 | 2023-07-28 | 江苏省特种设备安全监督检验研究院 | 便携式电流检测转速测试仪 |
CN108982964A (zh) * | 2018-07-28 | 2018-12-11 | 华中科技大学 | 一种基于细化傅里叶变换的信号分析方法及设备 |
CN109580221A (zh) * | 2018-11-30 | 2019-04-05 | 电子科技大学 | 一种机械传动振动频谱去混叠技术 |
CN109738697B (zh) * | 2019-01-24 | 2021-02-26 | 中国电子科技集团公司第二十九研究所 | 一种基于有限点离散频谱校正的频率测量方法 |
CN110022181B (zh) * | 2019-05-22 | 2021-06-04 | 成都天奥测控技术有限公司 | 一种提升频谱能量测量精度的方法和装置 |
CN110471018B (zh) * | 2019-09-19 | 2021-12-24 | 中南大学 | 一种频谱校正方法 |
CN111312291B (zh) * | 2020-02-24 | 2023-01-03 | 厦门快商通科技股份有限公司 | 信噪比检测方法、***、移动终端及存储介质 |
CN114258495B (zh) * | 2020-07-23 | 2023-08-15 | 刘保国 | 一种有限复数信号测量***与高精度分解方法 |
CN112213560A (zh) * | 2020-09-01 | 2021-01-12 | 东南大学 | 一种基于z-adaline的高精度电网宽频信号测量方法 |
CN115333603B (zh) * | 2022-07-14 | 2023-06-27 | 航天恒星科技有限公司 | 载波同步方法、装置、电子设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0996244A2 (de) * | 1998-10-23 | 2000-04-26 | Abb Research Ltd. | Vorrichtung und Verfahren zum Orten von Störstellen in elektrischen Verteilnetzen |
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN101718816A (zh) * | 2009-11-23 | 2010-06-02 | 浙江大学 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
CN101915874A (zh) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | 一种基于傅立叶变换的谐波检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5011340B2 (ja) * | 2009-03-31 | 2012-08-29 | アンリツ株式会社 | Apd測定装置 |
-
2013
- 2013-09-13 CN CN201310418656.0A patent/CN103454495B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0996244A2 (de) * | 1998-10-23 | 2000-04-26 | Abb Research Ltd. | Vorrichtung und Verfahren zum Orten von Störstellen in elektrischen Verteilnetzen |
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN101718816A (zh) * | 2009-11-23 | 2010-06-02 | 浙江大学 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
CN101915874A (zh) * | 2010-07-20 | 2010-12-15 | 北海市深蓝科技发展有限责任公司 | 一种基于傅立叶变换的谐波检测方法 |
Non-Patent Citations (2)
Title |
---|
基于LabWindows/CVI和DSP的高精度频谱分析***;焦新涛;《计量技术》;20081231(第5期);第17-20页 * |
离散频谱分析中的连续谱校正技术;焦新涛;《中国优秀博硕士学位论文全文数据库 (硕士) 工程科技Ⅱ辑》;20050315;第7页2.2.1节、第9页2.3节,第13页2.3.3节 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110068730A (zh) * | 2019-05-07 | 2019-07-30 | 中国科学院电子学研究所 | 双频段频谱数据采集方法及装置 |
CN110068730B (zh) * | 2019-05-07 | 2020-08-07 | 中国科学院电子学研究所 | 双频段频谱数据采集方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103454495A (zh) | 2013-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454495B (zh) | 自适应高精度快速频谱分析方法 | |
CN107505135B (zh) | 一种滚动轴承复合故障提取方法及*** | |
Gryllias et al. | Estimation of the instantaneous rotation speed using complex shifted Morlet wavelets | |
Chen et al. | Time-varying frequency-modulated component extraction based on parameterized demodulation and singular value decomposition | |
CN107329932B (zh) | 基于非线性调频分量分解的时频域模态参数辨识方法 | |
CN104375111B (zh) | 对密集频谱进行快速高精度细化校正的方法 | |
CN104006962A (zh) | 一种齿轮故障特征提取方法及*** | |
CN109765055B (zh) | 基于ewt、谱有效值和knn的滚动轴承故障检测方法及*** | |
CN103837740A (zh) | 一种高精度数字瞬时测频方法及装置 | |
CN103823177B (zh) | 基于窗函数设计的滤波器性能检测方法和*** | |
CN103743470A (zh) | 一种汽车噪声频谱分析方法 | |
Zhang et al. | Improved local cepstrum and its applications for gearbox and rolling bearing fault detection | |
CN104809702A (zh) | 基于频率域处理的脉诊曲线图像网格消除方法 | |
CN104849221B (zh) | 一种基于光学相干层析成像技术的相位定标方法 | |
CN104991119B (zh) | 一种消除伪峰、谱泄漏效应的互素谱分析方法及其装置 | |
CN107367640B (zh) | 基于相位信息的非均匀包络未知脉冲边沿估计方法及设备 | |
Yi et al. | Time-varying fault feature extraction of rolling bearing via time–frequency sparsity | |
CN107345834B (zh) | 基于共振加强奇异值分解的线性窄带滤波方法 | |
CN107688167B (zh) | 一种多时宽线性调频脉冲压缩信号幅度包络曲线生成方法 | |
CN108181486A (zh) | 加速度信号的处理方法及装置 | |
CN106027454A (zh) | 基于频率补偿的qam信号的频偏估计方法 | |
CN111611686A (zh) | 一种通信信号时频域的检测方法 | |
CN112345248A (zh) | 一种滚动轴承的故障诊断方法及装置 | |
CN109374969B (zh) | 一种基于相干积累的dft的下采样正弦信号相位差测量方法 | |
CN106646348A (zh) | 适用于短时多频率信号的干涉仪相位差测量电路及测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |