闪变值计算方法.doc

闪变值的计算方法:
假定为电压信号的采样值,并已经存放在一个时间长度为1min的数据文件中(,即1秒钟采6400个点,所以1min的数据文件共有60=384000个采样值)
求,为一个周波的电压均方根值序列。
(N=1,2, ……,256)
因为电压信号每个周波为20ms(),所以一个周波的采样点数为128,即m=128,
n1、n2为采样点的下标,因为是以一个周波为单位求,所以采样点的下标也是以128点为间隔变化的,例如:第一个周波的电压均方根值,第二个周波的电压均方根值,……,依次类推。
对数列进行FFT(256点FFT),求出其离散的频谱序列,即
(K=1,2,3, ……,256)
取出与下标K=5,13,18,23,28,36,41,45,54,59,77,82,87,113,118,123,128对应的,利用公式求出p:
,其中K就是上述列出的值,与K的对应关系如下表所示:(固定关系)

K
duK
5

13

18

23

28

K
duK
36

41

45

54

59

K
duK
77

82

87

113

118

K
duK
123

128

算1min所对应的所有P值,因为1min共有=3000个周波,因此一分钟需要计算
个p,取不到整数的原因是1min的周波数不是256的整数倍(),所以求最后一个p的时候,只能取得184个,即N=1,2,……,184,但在对这些进行FFT运算时,仍然用256点FFT,不够的用零代替,p的求法不变。
用界面将1min的p值显示出来。
注意:
只要是用FFT,都要首先对原始数据进行整序,然后再调用FFT算法进行处理。