CN105425037A - 一种频谱分析方法及装置 - Google Patents
一种频谱分析方法及装置 Download PDFInfo
- Publication number
- CN105425037A CN105425037A CN201510905737.2A CN201510905737A CN105425037A CN 105425037 A CN105425037 A CN 105425037A CN 201510905737 A CN201510905737 A CN 201510905737A CN 105425037 A CN105425037 A CN 105425037A
- Authority
- CN
- China
- Prior art keywords
- component
- centerdot
- increment
- sinusoidal
- rsqb
- 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
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
本发明涉及一种频谱分析方法及装置,依次包含如下步骤:步骤1:设定参数T,N,f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值;步骤2:设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量k 的初值;步骤3:依次读取一个采样数据,采用迭代法,对采样数据h(s)进行分解处理;步骤4:重复步骤3,直到完成对所有采样数据的分解处理,进行频谱分析。本发明结构简单、运算量小,减小了对内存的消耗。
Description
技术领域
本发明涉及信号分析领域,尤其涉及一种关于偏置正弦叠加信号的频谱分析方法。
背景技术
电能是现代社会中使用最为广泛的一种能源,随着社会的发展,各种复杂的、精密的、对电能质量敏感的用电设备越来越多,对电能质量的要求越来越高。然而,同时造成电能质量问题的因素也在不断增加,如以电力电子装置为代表的非线性负荷的使用,各种大型用电设备的启停等,使得电力信号发生严重畸形,给电力***带来了很大的“电网污染”。这些“电网污染”主要表现为电力***中的谐波与间谐波,因此为保证电能质量,对电力***中的谐波成分进行精确适时的分析和检测显得越来越紧迫和重要。
目前在谐波分析中常用分析方法是适于计算机分析的基于傅里叶和小波变换的分析方法,将原始谐波信号通过小波多分辨率分析分解为低频和高频两部分,然后利用傅里叶变换分析得到低频成分中稳态谐波的幅值、相位等参数,并利用小波包进行分析得到高频部分非稳态谐波的时域信息。然而这种方法需要利用傅里叶和小波变换两种方法分开进行计算,计算过程复杂,而且在计算过程中需要保存一个完整周期内的全部采样数据,会占用更多的内存空间,而且无法实现实时分析。
发明内容
本发明的目的是提供一种频谱分析方法,用以解决现有谐波分析方法中计算过程复杂的问题。
为实现上述目的,本发明的方案包括:
一种频谱分析方法,依次包含如下步骤:
步骤1:设定参数T,N,f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值;其中T是采样数据的采样周期,N是设定正弦分量的个数,f1f2…fN是采样数据各个正弦分量的频率数值,ε1ε2…εN依次是对应于指定频点f1f2…fN处的通频带的带宽,ε0是低通滤波的通频带;各正弦分量f1f2…fN和ε0ε1ε2…εN均为不大于2π/T的正数;
步骤2:设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值;
步骤3:依次读取一个采样数据,采用迭代法,对采样数据h(s)进行分解处理;
步骤4:重复步骤3,直到完成对所有采样数据的分解处理,进行频谱分析。
优选的,所述步骤1中,设定频率f1、f2、…、fN的初值分别等于被测信号正弦分量的实际频率,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值均为0;所述步骤2中,设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值均为0。
进一步的,所述步骤3中迭代法依次包含如下步骤:
S01:利用式(1)获得偏置分量增量k2(0),令n分别取值为1,2,…,N,循环执行式(2),获得正弦分量增量k2(n)和余弦分量增量
S02:依据式(3)对偏置分量s0进行迭代处理,令n分别取值为1,2,…,N,循环执行式(4),先对正弦分量sn、余弦分量cn进行迭代处理;
s0=s0+0.5T[k1(0)+k2(0)](3),
S03:依据迭代处理后的偏置分量s0和正弦分量s1、s2、…、sN的值,利用式(5)更新偏置分量增量k1(0),令n分别取值为1,2,…,N,循环执行式(6),更新正弦分量增量k1(n)和余弦分量增量
S04:利用式(7),将偏置分量s0输出至s0(m)、令n分别取值为1,2,…,N,循环执行式(8),正弦分量sn、余弦分量cn依次输出至sn(m)、cn(m)。其中m=1,2,3,…;
s0(m)=s0(7),
进一步的,根据所述(S02)中计算得到的正弦分量sn、余弦分量cn,计算得到估计幅值dn的值,按照下式(9)得到估计幅值dn,并最终输出到dn(m);
进一步的,步骤3中,依次读取离线采样数据,或者通过中断方式读取实时数据。
一种频谱分析装置,依次包含如下模块:
模块1:用于设定参数T,N,f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值的装置;其中T是采样数据的采样周期,N是设定正弦分量的个数,f1f2…fN是采样数据各个正弦分量的频率数值,ε1ε2…εN依次是对应于指定频点f1f2…fN处的通频带的带宽,ε0是低通滤波的通频带;各正弦分量f1f2…fN和ε0ε1ε2…εN均为不大于2π/T的正数;
模块2:用于设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值的装置;
模块3:用于依次读取一个采样数据,采用迭代法,对采样数据h(s)进行分解处理的装置;
模块4:用于重复执行模块3,完成对所有采样数据的分解处理,进行频谱分析的装置。
优选的,所述模块1中,设定频率f1、f2、…、fN的初值分别等于被测信号正弦分量的实际频率,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值均为0;所述模块2中,设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值均为0。
进一步的,所述模块3中迭代法依次包含如下子模块:
子模块1:用于利用式(1)获得偏置分量增量k2(0),令n分别取值为1,2,…,N,循环执行式(2),获得正弦分量增量k2(n)和余弦分量增量的装置;
子模块2:用于依据式(3)对偏置分量s0进行迭代处理的装置,令n分别取值为1,2,…,N,循环执行式(4),先对正弦分量sn、余弦分量cn进行迭代处理;
s0=s0+0.5T[k1(0)+k2(0)](3),
子模块3:用于更新正弦分量增量k1(n)和余弦分量增量的装置,依据迭代处理后的偏置分量s0和正弦分量s1、s2、…、sN的值,利用式(5)更新偏置分量增量k1(0),令n分别取值为1,2,…,N,循环执行式(6);
子模块4:用于将正弦分量sn、余弦分量cn依次输出至sn(m)、cn(m)的装置,利用式(7),将偏置分量s0输出至s0(m)、令n分别取值为1,2,…,N,循环执行式(8),其中m=1,2,3,…;
s0(m)=s0(7),
进一步的,根据所述子模块2中计算得到的正弦分量sn、余弦分量cn,计算得到估计幅值dn的值,按照下式(9)得到估计幅值dn,并最终输出到dn(m);
进一步的,模块3中,依次读取离线采样数据,或者通过中断方式读取实时数据。
与现有技术相比,本发明具有以下优点:
1、本发明针对经定时采样后得到的离散时间数据,采用迭代方法进行信号分析,便于利用计算机实现,迭代公式选定好即可根据设定的初始值进行迭代,得到频谱分析结果,并且不需要计算正弦函数和余弦函数,结构简单、运算量小,程序简单不复杂。
2、本发明在计算过程中不需要保存一个周期内的多个采样数据,减小了内存的消耗。
3、本发明可以对采集到的信号实时进行分析,而且具有两种模式选择:离线模式和在线模式。
4、本发明不要求所指定的正弦分量的频率在数值上保持特定关系,便于分析被测信号的谐波和间谐波成分,尤其适用于电力***领域。
附图说明
图1是离线模式下频谱分析的流程示意图;
图2是在线模式下频谱分析的流程示意图;
图3是被测信号中偏置分量和3个正弦分量分析前后的对比图;
图4是信号分量跟踪总误差的变化曲线;
图5是幅值跟踪总误差变化曲线。
具体实施方式
下面结合附图对本发明做进一步详细的说明。
实施例1
本实施例以机械振动的振动分量指定频率信号的频谱分析为例,将采集到的采样数据进行分解以实现频谱分析。具体为离线模式下的频谱分析,这种分析方法适用于分析保存在存储器件中的离散时间序列,这些离散时间序列是被测信号经周期采样得到的。
如图1所示,是本实施例进行频谱分析的流程图,首先对程序的运行进行初始化设置。具体的,首先,设置离散时间序列中数据的个数K,设定当前要处理的数据在离散时间序列中的位置号k为1,设定当前迭代处理后的数据在输出序列中的位置号m为1;
然后设定参数T,N,f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值;其中T是采样周期,N是设定正弦分量的个数,f1f2…fN是各个正弦分量的频率数值,ε1ε2…εN依次对应于指定频点f1f2…fN处的通频带的带宽,ε0是低通滤波的通频带。接着设定参数T、N、f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0、正弦分量s1s2…sN、余弦分量c1c2…cN的初值。
其中,采样周期T、正弦分量个数N、指定频率f1f2…fN、参数ε0ε1ε2…εN的初值,均可依据被测信号的先验知识与信号分析要求设定。
采样周期T首先要满足香农采样定理要求,在满足香农采样定理要求下,采样周期越小,分析精度越高。对于在线分析方式,受实时性的限制,正弦分量个数N和采样周期T还应该满足(N+1)Δt<T,其中Δt表示执行一遍定时中断步骤所需要的最大时间。
参数ε0的物理意义相当于低通滤波的通频带,ε0ε1ε2…εN的物理意义依次相当于指定频点f1f2…fN处的通频带的带宽,依据香农采样定理,限定其数值均不大于2π/T,也可以设为相同的数值ε以简化算法结构。ε数值大小对估计幅值的收敛速度具有主要影响,ε的值越大,估计幅值越快地收敛到实际值,但同时加大了干扰对幅值估计精度的不良影响。
指定频率f1、f2、…、fN均分别等于被测信号正弦分量的实际频率,迭代变量偏置分量s0、正弦分量s1、s2、…、sN余弦分量c1、c2、…、cN总能分别收敛到各自的实际值。因此对于这些迭代变量的初值没有特别限制。优选地,均设定为0。
当这些基本值都设定完毕后,则开始设定偏置分量增量k1(0),正弦分量增量k1(1)k1(2)…k1(N),余弦分量增量的初值;对于偏置分量增量k1(0),正弦分量增量k1(1)k1(2)…k1(N)和余弦分量增量 的初值没有特别限制。优选地,均设定为0。
上述步骤是对信号进行处理时的初始化步骤,给出了一定的初始化方式,属于现有技术,同样也可以采用其他方式来实现初始化;而步骤3则是本发明的创新所在,下面做出具体说明。
参数设定完毕后,则开始对采样数据通过迭代方法进行处理,从而实现对全部采样数据的分解,以便于后续的频谱分析。
首先读取离散时间序列中的第k个数据作为当前处理的采样数据h(s),依次执行步骤3中数据处理的过程;为提高分析方法的有效性,把迭代处理后的数据存储到输出序列中,并把输出序列中的位置号m增加1;然后进入循环控制步骤。
选用迭代法中的一种具体公式对采样信号进行处理,具体处理过程如下:
S01:利用式(1)获得偏置分量增量k2(0),令n分别取值为1,2,…,N,循环执行式(2),获得正弦分量增量k2(n)和余弦分量增量
S02:依据式(3)对偏置分量s0进行迭代处理,令n分别取值为1,2,…,N,循环执行式(4),先对正弦分量sn、余弦分量cn进行迭代处理,然后计算估计幅值dn的值;
s0=s0+0.5T[k1(0)+k2(0)](3)
S03:依据迭代处理后的偏置分量s0和正弦分量s1、s2、…、sN的值,利用式(5)更新偏置分量增量k1(0),令n分别取值为1,2,…,N,循环执行式(6),更新正弦分量增量k1(n)和余弦分量增量
S04:利用式(7),将偏置分量s0输出至s0(m)、令n分别取值为1,2,…,N,循环执行式(8),正弦分量sn、余弦分量cn和估计幅值dn依次输出至sn(m)、cn(m)、dn(m)。其中m=1,2,3,…。
s0(m)=s0(7)
在循环控制步骤中,先把位置号k增加1,再依据位置号k和离散时间序列中数据的个数K的值判断是否返回执行增量计算及数据迭代步骤。若k≤K,返回执行增量计算及数据迭代步骤;若k>K,表示离散时间序列中所有采样数据都已处理完毕,则终止运行,结束离线分析过程。
实施例2
上述实施例给出了机械振动的振动分量指定频率信号提取分析中本发明的使用,具体为离线模式下的频谱分析;作为其他实施方式,频谱分析还可以运行在非离线模式下,即在线模式。在本实施例中,具体为对电力***中谐波和间谐波的在线分析,在线模式下频谱分析的特点在于边采样边进行频谱分析,适用于需要对被测信号的每个采样数据都进行实时分析处理的情况。
如图2所示,是本实施例进行频谱分析的流程图,首先对程序的运行进行初始化设置。具体的,对应于步骤1,设定当前要处理的数据在离散时间序列中的位置号k为1,设定当前迭代处理后的数据在输出序列中的位置号m为1;接着设定参数T、N、f1、f2、…、fN、ε0、ε1、ε2、…、εN、的值,设定偏置分量s0、正弦分量s1、s2、…、sN、余弦分量c1、c2、…、cN的初值;
当这些基本值都设定完毕后,对应于步骤2的,再设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值;
由于本实施例与上述实施例1的不同仅在于步骤3的处理,步骤1和2完全相同,而步骤1和2中对这些参数的设置已经完全说明清楚,因此在此不再赘述。值得注意的是,中断定时器的定时时间为T,以定时开放***的定时中断。
对应于步骤3的,采用迭代法对采样数据进行处理。首先保存中断现场的各个寄存器的当前值;接着对被测信号进行采样获得采样数据h(s),执行处理方案中计算增量数值并进行迭代处理;处理完成后再恢复中断现场的各个寄存器的值,然后中断返回到主程序执行。其中,对采样数据h(s)的处理方法,与上述实施例1中的处理方法完全相同,在此不再赘述。
对应于步骤4的,循环执行步骤3,完成对所有采样数据的处理。
在线分析实施方式通过定时器循环产生定时中断事件,引起中断服务程序循环执行。
上述实施例给出了两种具体的实施方式,为对本发明的技术特征、目的、效果有更加清楚的理解,现以一个实际工程为例,做出更具体详细的说明。
设被测信号为s=S0+S1sin(w1t+θ1)+S2sin(w2t+θ2)+S3sin(w3t+θ3),其中三个交流频率w1、w2、w3数值分别等于22π、46π、74π,单位为弧度/秒,其余各个参数随时间t的变化如表1所示,其中π为圆周率。
为通过变化曲线说明本发明的偏置正弦叠加信号的频谱分析方法的有效性,定义信号分量总跟踪误差errs为:
errs=|S0-s0|+|S1sin(w1t+θ1)-s1|+|S2sin(w2t+θ2)-s2|+|S3sin(w3t+θ3)-s3|
定义幅值总跟踪误差errd为:
errd=|S1-d1|+|S2-d2|+|S3-d3|
表1被测信号各个参数随时间t变化
设采样周期T=0.1毫秒,对被测信号进行定时采样,获得采样数据形成离散时间序列,再按照图1所示的离线分析实施方式,编写程序在计算机中仿真运行。设定N=3,指定频率f1、f2、f3的值依次等于w1、w2、w3,偏置分量s0、正弦分量s1、s2、s3、余弦分量c1、c2、c3的初值均为0。设定偏置分量增量k1(0)、正弦分量增量k1(1)、k1(2)、k1(3)、余弦分量增量的初值均为0。设定ε0、ε1、ε2、ε3的值都等于ε,所得结果显示在图3~图5中。
图3给出了被测信号的偏置分量及3个正弦分量的精确值和运用本发明得出的偏置分量及3个正弦分量的跟踪曲线。用于说明本发明的偏置正弦叠加信号的频谱分析方法的正弦分量的跟踪情况。
图4用于说明本发明的偏置正弦叠加信号的频谱分析方法的信号分解与跟随性能以及参数ε对信号分量收敛速度的影响,显示了本发明的偏置正弦叠加信号的频谱分析方法的信号分析性能,当信号分量跟踪总误差errs等于零时,说明偏置分量与各个正弦分量都分别收敛到其实际值,实现对被测信号的精确分析和准确跟随。
图5用于说明本发明的偏置正弦叠加信号的频谱分析方法的幅值跟随性能以及参数ε对幅值收敛速度的影响,显示了本发明的偏置正弦叠加信号的频谱分析方法的幅值跟踪性能,如果幅值跟踪总误差errd等于零,说明所有分量的估计幅值都收敛到其实际值。
图4和图5都显示出,本发明的偏置正弦叠加信号的频谱分析方法的收敛速度均决定于参数ε的数值,参数ε的数值越大,收敛速度越快。
本发明提供的一种频谱分析方法,适用于电力***中的谐波与间谐波分析、机械振动的振动分量指定频率信号提取,但不限于此两种领域的应用,还可以适用于其他技术领域。
以上给出了本发明具体的实施方式,但本发明不局限于所描述的实施方式。在本发明给出的思路下,采用对本领域技术人员而言容易想到的方式对上述实施例中的技术手段进行变换、替换、修改,并且起到的作用与本发明中的相应技术手段基本相同、实现的发明目的也基本相同,这样形成的技术方案是对上述实施例进行微调形成的,这种技术方案仍落入本发明的保护范围内。
Claims (10)
1.一种频谱分析方法,其特征在于,依次包含如下步骤:
步骤1:设定参数T,N,f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值;其中T是采样数据的采样周期,N是设定正弦分量的个数,f1f2…fN是采样数据各个正弦分量的频率数值,ε1ε2…εN依次是对应于指定频点f1f2…fN处的通频带的带宽,ε0是低通滤波的通频带;各正弦分量f1f2…fN和ε0ε1ε2…εN均为不大于2π/T的正数;
步骤2:设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值;
步骤3:依次读取一个采样数据,采用迭代法,对采样数据h(s)进行分解处理;
步骤4:重复步骤3,直到完成对所有采样数据的分解处理,进行频谱分析。
2.根据权利要求1所述的一种频谱分析方法,其特征在于,所述步骤1中,设定频率f1、f2、…、fN的初值分别等于被测信号正弦分量的实际频率,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值均为0;所述步骤2中,设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值均为0。
3.根据权利要求1所述的一种频谱分析方法,其特征在于,所述步骤3中迭代法依次包含如下步骤:
S01:利用式(1)获得偏置分量增量k2(0),令n分别取值为1,2,…,N,循环执行式(2),获得正弦分量增量k2(n)和余弦分量增量
S02:依据式(3)对偏置分量s0进行迭代处理,令n分别取值为1,2,…,N,循环执行式(4),先对正弦分量sn、余弦分量cn进行迭代处理;
s0=s0+0.5T[k1(0)+k2(0)](3),
S03:依据迭代处理后的偏置分量s0和正弦分量s1、s2、…、sN的值,利用式(5)更新偏置分量增量k1(0),令n分别取值为1,2,…,N,循环执行式(6),更新正弦分量增量k1(n)和余弦分量增量
S04:利用式(7),将偏置分量s0输出至s0(m)、令n分别取值为1,2,…,N,循环执行式(8),正弦分量sn、余弦分量cn依次输出至sn(m)、cn(m)。其中m=1,2,3,…;
s0(m)=s0(7),
4.根据权利要求3所述的一种频谱分析方法,其特征在于,根据所述(S02)中计算得到的正弦分量sn、余弦分量cn,计算得到估计幅值dn的值,按照下式(9)得到估计幅值dn,并最终输出到dn(m);
5.根据权利要求1所述的一种频谱分析方法,其特征在于,步骤3中,依次读取离线采样数据,或者通过中断方式读取实时数据。
6.一种频谱分析装置,其特征在于,依次包含如下模块:
模块1:用于设定参数T,N,f1f2…fN,ε0ε1ε2…εN的值,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值;其中T是采样数据的采样周期,N是设定正弦分量的个数,f1f2…fN是采样数据各个正弦分量的频率数值,ε1ε2…εN依次是对应于指定频点f1f2…fN处的通频带的带宽,ε0是低通滤波的通频带;各正弦分量f1f2…fN和ε0ε1ε2…εN均为不大于2π/T的正数;
模块2:用于设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值;
模块3:用于依次读取一个采样数据,采用迭代法,对采样数据h(s)进行分解处理;
模块4:用于重复执行模块3,完成对所有采样数据的分解处理,进行频谱分析。
7.根据权利要求6所述的一种频谱分析装置,其特征在于,所述模块1中,设定频率f1、f2、…、fN的初值分别等于被测信号正弦分量的实际频率,设定偏置分量s0,正弦分量s1、s2、…、sN和余弦分量c1、c2、…、cN的初值均为0;所述模块2中,设定偏置分量增量k1(0),正弦分量增量k1(1)、k1(2)、…、k1(N),余弦分量增量的初值均为0。
8.根据权利要求6所述的一种频谱分析装置,其特征在于,所述模块3中迭代法依次包含如下子模块:
子模块1:用于利用式(1)获得偏置分量增量k2(0),令n分别取值为1,2,…,N,循环执行式(2),获得正弦分量增量k2(n)和余弦分量增量
子模块2:用于依据式(3)对偏置分量s0进行迭代处理,令n分别取值为1,2,…,N,循环执行式(4),先对正弦分量sn、余弦分量cn进行迭代处理;
s0=s0+0.5T[k1(0)+k2(0)](3),
子模块3:依据迭代处理后的偏置分量s0和正弦分量s1、s2、…、sN的值,利用式(5)更新偏置分量增量k1(0),令n分别取值为1,2,…,N,循环执行式(6)更新正弦分量增量k1(n)和余弦分量增量
子模块4:用于将正弦分量sn、余弦分量cn依次输出至sn(m)、cn(m),利用式(7),将偏置分量s0输出至s0(m)、令n分别取值为1,2,…,N,循环执行式(8),其中m=1,2,3,…;
s0(m)=s0(7),
9.根据权利要求8所述的一种频谱分析装置,其特征在于,根据所述子模块2中计算得到的正弦分量sn、余弦分量cn,计算得到估计幅值dn的值,按照下式(9)得到估计幅值dn,并最终输出到dn(m);
10.根据权利要求6所述的一种频谱分析装置,其特征在于,模块3中,依次读取离线采样数据,或者通过中断方式读取实时数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510905737.2A CN105425037A (zh) | 2015-12-09 | 2015-12-09 | 一种频谱分析方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510905737.2A CN105425037A (zh) | 2015-12-09 | 2015-12-09 | 一种频谱分析方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105425037A true CN105425037A (zh) | 2016-03-23 |
Family
ID=55503369
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510905737.2A Pending CN105425037A (zh) | 2015-12-09 | 2015-12-09 | 一种频谱分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105425037A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111077369A (zh) * | 2019-12-28 | 2020-04-28 | 苏州摩联通信技术有限公司 | 频谱分析方法、装置、电子设备及存储介质 |
CN111487463A (zh) * | 2020-04-21 | 2020-08-04 | 全球能源互联网研究院有限公司 | 一种间谐波检测***及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070242489A1 (en) * | 2006-04-13 | 2007-10-18 | Tatung Company | Method of designing an RPWM inverter with unwanted harmonic elimination |
RU2435168C1 (ru) * | 2010-04-09 | 2011-11-27 | ООО предприятие "КОНТАКТ-1" | Способ гармонического анализа периодического многочастотного сигнала |
CN102967760A (zh) * | 2012-11-30 | 2013-03-13 | 合肥工业大学 | 指定频率的信号分析方法 |
CN102967761A (zh) * | 2012-11-30 | 2013-03-13 | 合肥工业大学 | 频率自适应信号分析方法 |
CN102981045A (zh) * | 2012-11-30 | 2013-03-20 | 合肥工业大学 | 归一化自适应电功率计量方法 |
-
2015
- 2015-12-09 CN CN201510905737.2A patent/CN105425037A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070242489A1 (en) * | 2006-04-13 | 2007-10-18 | Tatung Company | Method of designing an RPWM inverter with unwanted harmonic elimination |
RU2435168C1 (ru) * | 2010-04-09 | 2011-11-27 | ООО предприятие "КОНТАКТ-1" | Способ гармонического анализа периодического многочастотного сигнала |
CN102967760A (zh) * | 2012-11-30 | 2013-03-13 | 合肥工业大学 | 指定频率的信号分析方法 |
CN102967761A (zh) * | 2012-11-30 | 2013-03-13 | 合肥工业大学 | 频率自适应信号分析方法 |
CN102981045A (zh) * | 2012-11-30 | 2013-03-20 | 合肥工业大学 | 归一化自适应电功率计量方法 |
Non-Patent Citations (1)
Title |
---|
储昭碧等: "基于离散频率自适应梳状滤波器的时变谐波分析", 《电测与仪表》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111077369A (zh) * | 2019-12-28 | 2020-04-28 | 苏州摩联通信技术有限公司 | 频谱分析方法、装置、电子设备及存储介质 |
CN111077369B (zh) * | 2019-12-28 | 2022-04-15 | 苏州摩联通信技术有限公司 | 频谱分析方法、装置、电子设备及存储介质 |
CN111487463A (zh) * | 2020-04-21 | 2020-08-04 | 全球能源互联网研究院有限公司 | 一种间谐波检测***及方法 |
CN111487463B (zh) * | 2020-04-21 | 2022-03-01 | 全球能源互联网研究院有限公司 | 一种间谐波检测***及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102967760B (zh) | 指定频率的信号分析方法 | |
CN102967761B (zh) | 频率自适应信号分析方法 | |
CN102981045B (zh) | 归一化自适应电功率计量方法 | |
CN103235518A (zh) | 一种对光电跟踪稳定***的试验方法及平台 | |
CN103018546B (zh) | 指定频率的电功率计量方法 | |
CN104849545A (zh) | 一种电力***频率测量方法及测量装置 | |
CN104376368A (zh) | 一种基于频域分解的风力发电短期负荷预测方法及装置 | |
CN103516279A (zh) | 一种基于fpga的永磁同步电机控制芯片 | |
CN116306798A (zh) | 一种超短时风速预测方法及*** | |
CN105114189A (zh) | 基于fpga实现的电子节气门模型预测控制*** | |
CN105425037A (zh) | 一种频谱分析方法及装置 | |
CN104682387A (zh) | 一种基于多区域交互迭代的概率潮流计算方法 | |
Bai et al. | Application of Time‐Frequency Analysis in Rotating Machinery Fault Diagnosis | |
Lv et al. | Longitudinal synchroextracting transform: A useful tool for characterizing signals with strong frequency modulation and application to machine fault diagnosis | |
Zhang et al. | Absolute position acquisition for linear synchronous motor with passive mover | |
CN101806834B (zh) | 基于卡尔曼滤波器的信号实时时频谱仪 | |
CN105510698A (zh) | 一种基于指定交流成分的用电特性计量方法 | |
Lu et al. | Real-time chatter detection based on fast recursive variational mode decomposition | |
CN105550516A (zh) | 基于Adams算法的频率自适应谐波间谐波分析方法 | |
Liu et al. | Sparse coefficient fast solution algorithm based on the circulant structure of a shift-invariant dictionary and its applications for machine fault diagnosis | |
CN109617051B (zh) | 一种新能源电力***低频振荡参数辨识方法 | |
CN103018558A (zh) | 主从式多处理器实时信号分析方法 | |
Zou | Motion automation control system of ship assembly robot based on plc | |
CN105510708A (zh) | 一种指定频率的谐波或间谐波的计量方法 | |
CN105425034A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160323 |