基于全相位频谱分析的正弦信号高精度参数估计方法.docx
《基于全相位频谱分析的正弦信号高精度参数估计方法.docx》由会员分享，可在线阅读，更多相关《基于全相位频谱分析的正弦信号高精度参数估计方法.docx（12页珍藏版）》请在冰豆网上搜索。
基于全相位频谱分析的正弦信号高精度参数估计方法
电力系统中,对电网电压和电流基波参数的测量分析,通常采用传统的离散傅里叶变换（DFT进行频谱分析来实现,当信号的频率不是DFT频率分辨率的整数倍的时候,即对信号非整周期采样,会产生
频谱泄漏,使测得的幅值、
频率和相位偏离实际值,尤其相位测量误差更大[1-4]。
现提出一种基于全相位频谱分析的正弦信号高精度参数估计算法,在精度与实时性与算法复杂度上均优于现有算法,软件编写简单,在无噪情况下参数估计近似为无偏估计,尤其对相位的估计,误差可达到0.001%。
1全相位频谱分析
文献[5]提出一种新型的频谱估计算法,对传统DFT频谱分析时数据的截断方式进行了改进,可以很大程度地减小频谱泄漏。
若将N阶中心对称窗和N阶矩形窗卷积产生的一个2N-1阶窗作为窗函数,则是全相位单窗频谱分析,若将N阶中心对称窗和自身卷积产生的一个2N-1阶窗作为窗函数,则是全相位双窗频谱分析,其框图如图1所示。
首先,推导对具有单一频率f0的单频信号进行全相位频谱分析得到的幅度。
设单频信号:
x=ej2πf0
fs
n
其中,f0为信号频率,fs为采样频率。
对于时间序列中的一点x（N,存在也只存在N个包含该点的N维向量:
X0=[x（Nx（N+1…x（2N-1]TX1=[x（N-1x（N…x（2N-2]T
…
XN-1=[x（1x（2…x（N]T
将每个向量进行循环移位,把样本点x（N移到首位,则得到另外的N个N维向量:
X0′=[x（Nx（N+1…x（2N-1]TX1′=[x（Nx（N+1…x（N-1]T
…
XN-1′=[x（Nx（1…x（N-1]T
对准x（N相加得到全相位数据向量:
XAP=1N
[Nx（N（N-1x（N+1+x（1…
x（2N-1+（N-1x（N-1]T
基于全相位频谱分析的正弦信号
高精度参数估计方法
黄晓红1,王兆华2
（1.河北理工大学信息学院,河北唐山063000;
2.天津大学电子信息工程学院,天津300072
摘要:
对传统频谱分析的输入数据截断方式进行了改进,提出了全相位频谱分析,理论证明全相位频谱分析幅度谱是传统FFT频谱分析幅度谱的平方,可以很大程度上减小频谱泄露,同时,给出了其实现框图。
另外,全相位频谱分析方法具有良好的相位分析性能,不受频偏影响,在信号是非整数倍频率采样情况下,由该方法分析的相位与信号真实相位误差极小,无需校正。
提出基于全相位相位差法估计正弦信号幅值、频率和相位的新算法,对传统相位差法其序列的取法进行了改进,对单频余弦信号进行非整周期采样,分别用全相位相位差法和离散频谱综合相位差校正法对信号的频率、幅值和相位进行了校正。
仿真结果证明该算法参数估计精度高于现有算法,具有较好的实用价值,在无噪情况下参数估计近似为无偏估计,尤其对相位的估计,误差达到0.001%。
关键词:
全相位频谱分析;参数估计;相位差;校正中图分类号:
TN911
文献标识码:
A
文章编号:
1006-6047（200807-0054-03
收稿日期:
2007-04-06;修回日期:
2007-12-05
电力自动化设备
ElectricPowerAutomationEquipment
Vol.28No.7Jul.2008
第28卷第7期2008年7月
图1传统频谱分析与全相位频谱分析框图
Fig.1BlockdiagramsoftraditionalandapFFTspectrumanalysis
x（n
Z-1
Z-1Z-1
…
N阶对称窗WN阶FFT…
…（aN阶传统FFT频谱分析框图
（bN阶全相位频谱分析框图
Z
-1Z-1
…
N阶对称窗W1与N阶对称窗W2卷积
N阶FFT
…
…
x（n
Z-1+
Z-1
+
+
根据DFT移位性质,X′i（i=0,1,…N-1的离散傅里叶变换X′i（k和Xi（i=0,1,…,N-1的离散傅里叶变换Xi（k之间有很明确的关系:
X′i（k=Xi（kej
2πki
N
（1全相位频谱是由X′i（k之和组成,所以有
XAP（k=1N!
i=0N-1
X′i
（k=1N!
i=0
N
-1Xi（kej2πkiN=1N!
i=0N
-1!
n=0
N-1
ej2π
f0fs
（N-i+ne-j2πNknej2πkiN=1N
ej2π
f0fs
N!
i=0N-1e-j2πf0fs
-kN"#i!
n=0N-1ej2πf0fs
-k
N"$n
=
1Nej2πf0fs
Nsin2πN（f0/fs-k/Nsin2π（f0/fs-k/N
（2由式（2可见全相位频谱分析幅度谱为
1N
sinπN（f0/fs-k/Nsinπ（f0/fs-k/N
2
其为传统DFT频谱分析幅度谱的平方,这对减
小频谱泄漏很有益。
全相位频谱分析另一个重要的特点是其相位不用校正,恒定不变,不受频偏影响,即当信号未做整周期截断时,全相位频谱分析仍能求出信号的真实相位,得到的相位与信号的真实相位误差极小,近似相等,这是传统DFT频谱分析所不具备的,也是用全相位频谱分析进行参数校正的基础。
现以余弦信号cos（1.2×2πt/6+100π/180为例进行N=6阶全相位频谱分析,分析全相位方法求相位误差小的原因。
11个取样信号为:
-0.1736,-0.9903,-0.4384,0.7193,0.8829,-0.1736,-0.9903,-0.4384,0.7193,0.8829,-0.1736。
全相位输入信号由6组N=6的取样信号组成,第1组是由11个取样信号的最后6个组成,第2组由向左移1位的6个数组成,但-0.1736要循环移到首位,其他类同。
信号排列如表1所示。
这6组N=6的取样信号的相位φ如表2所示。
因为频率是1.2Hz,从每组第2个相位可见,6个相位中3个偏离100°增大,3个偏离100°减小。
全相位输入数据是上面6组信号之和的平均,相位互相抵消,使得与原始信号的相位差为零,所以全相位频谱分析后求出的相位即为信号真实相位。
实验结果表明:
全相位无窗频谱分析时,求整数倍频率相位准确,如果求偏离整数倍频率相位时,全相位加kaiser（N,9.5双窗误差最小。
此例使用全相位kaiser（N,9.5双窗求出的信号相位如下:
180.0000°
100.0069°,100.0004°180.0000°
259.9996°,259.9931°与信号真实的相位误差仅为0.69%,可认为相等。
从具体的数据分析,可清楚地看出全相位频谱分析具有良好的相位分析性能,在信号是非整数倍频率情况下相位不变。
而传统谱分析方法相位会发生很大误差。
用计算机生成如下信号,采样频率128Hz,作谱点数为128的分析,频率分辨率为1Hz,有
s1（t=cos（9×2πt/N+50π/180+
cos（19.1×2πt/N+50π/180+cos（29.2×2πt/N+50π/180+cos（39.3×2πt/N+50π/180+
cos（49.4×
2πt/N+50π/180全相位频谱分析采用hanning（N双窗,传统DFT频谱分析采用hanning（N窗,分析结果如图2（a（b所示。
可见全相位双窗频谱分析的各频率成分信号的相位均为信号真实相位50°,而传统加窗分析的信号相位除了第一个没有偏离整数倍频率信号的相位为50°,其余频率成分信号的相位分析均与真实相位有很大误差,且频率偏离整数倍越远,相位分析结果误差越大。
若将信号相位发生变化,即产生信号:
s2（t=cos（9×2πt/N+10π/180+
cos（19.1×
2πt/N+30π/180+cos（29.2×2πt/N+50×π/180+cos（39.3×2πt/N+70π/180+
cos（49.4×
2πt/N+90π/180分析结果如图2（c（d所示。
由图2可明显看出,全相位频谱分析得到的相位和信号的真实相位误差极小,而传统DFT方法分析得到相位误差则较大。
2全相位相位差法估计正弦信号参数
文献[2]中给出了关于相位差法校正的原理,全相位相位差法对其序列的取法进行了改进。
具体算法步骤如下:
a.对信号做一定长度的非整周期采样,从采样序列的第i点开始起,取2N-1点,得到一个序列u1（t;再从i+N点开始起,也取2N-1点,得到另一个序列u2（t;
b.对序列u1（t进行全相位双窗（汉宁与汉恩窗的卷积窗N点频谱分析,得到频率表示为f=kfs/N,频率分辨率F0=fs/N（k为主谱线的序号,频率校正
第1组
第2组
第3组
第4组
第5组
第6组
-0.1736-0.1736-0.1736-0.1736-0.1736-0.1736-0.9903-0.9903-0.9903-0.9903-0.9903-0.1736-0.4384-0.4384-0.4384-0.4384-0.9903-0.99030.71930.71930.7193-0.4384-
0.4384-
0.43840.88290.88290.71930.71930.71930.7193-0.17360.88290.88290.8829
0.8829
0.8829
表1取样信号
Tab.1Sampledsignals
第1组
第2组
第3组
第4组
第5组
第6组
180.0000
00180.0000180.0000180.0000132.5107112.9053112.497888.459883.627572.5107-25.769951.729853.4788138.4375161.4769214.23010180.0000180.00000
0180.000025.7699308.2702306.5212-138.4375-161.4769145.7699-132.5107247.0947247.5022-88.4598-83.6275287.4893
表2取样信号相位（N=6
Tab.2Phasesofsampledsignals（N=6
（°黄晓红,等:
基于全相位频谱分析的正弦信号高精度参数估计方法
第7
期
图2传统DFT与全相位频谱分析多信号的相位谱
Fig.2Phasespectrumsanalyzedbytraditional
DFTandapFFTspectrumanalysis
100
50
20
40
6080
（c双窗频谱分析相位谱（10°30°50°70°90°
200100
20
406080（d传统加窗DFT相位谱（10°30°50°70°90°
φ/（°
φ/（°
f/Hz
f/Hz
10
305070
90
10
47.885.7123.6161.5
200100
20
406080
（b传统加窗DFT相位谱（50°50°50°50°50°
φ/（
°
100500
20
406080
（a双窗频谱分析相位谱（50°50°50°50°50°
φ/（°
f/Hz
f/Hz
50
67.8
85.7103.5121.4
量Δf=dF0（d为泄露误差系数;
c.信号校正泄漏误差后的频率为
f!
=f+Δf=f+dF0=（k+dF0
同时得到相位φ1;
d.对序列u2（t进行全相位双窗（汉宁与汉恩窗的卷积窗N点频谱分析,得到主谱线处相位φ2,根据相位差Δφ=φ2-φ1=2πd,计算出泄漏误差系数d;
e.校正其频率、
幅值、相位,将d代入步骤c后得到校正后的频率。
令主谱线处的粗估幅值为Ak,则幅值校正公式为
A=2πΔ
φ（1-Δφ2
Ak!
2sin（πΔφ
"#
2
信号的初始相位公式为
φ0=2φ1-φ2
3仿真与结果
单频余弦信号为
s（t=cos（4.39×2πt/fs+40π/180
对其进行非整周期采样,采样频率fs=1024Hz,频谱分析点数为N=1024点,分别用全相位相位差法和文献[6]相位差法进行参数估计,结果如表3所示（表中幅值为归一化数值。
4结语
在传统频谱分析的基础上,对输入数据截断方
式进行了改进,得到了全相位频谱分析,相对于传统
频谱分析,具有频谱泄露小、参数估计精度高的优点。
该算法无需对信号进行严格的整周期采样,可有效抑制DFT算法中的频谱泄露和栅栏效应,且无需相位校正,算法简单,可用于电力自动化中的向量测量以及正弦交流电谐波分析等工程,有良好的实用价值。
参考文献:
[1]JAINVK,COLLINSWL,DAVIDDC.High-accuracyanalog
measurementsviainterpolatedFFT[J].IEEETransonInstrum
Meas,1979,28（2:
113-122.
[2]曹国剑,黄纯,谢雁鹰,等.基于时移相位差校正的准周期算法的研究
[J].继电器,2004,32（13:
1-5.
CAOGuojian,HUANGChun,XIEYanying,etal.Researchonquasi-synchronizationalgorithmbasedontime-shiftingandphasedifferencecorrecting[J].Relay,2004,32（13:
1-5.
[3]RENSERSH,SCHOUKENSJ,VILAING.High-accuracyspectrum
analysisofsampleddiscretefrequencysignalsbyanalyticalleakagecompensation[J].IEEETransactionsonInstrumentationandMeasurement,1984,33（4:
287-292.
[4]丁仁杰.电力系统同步相量动态测量技术的研究与实现[D].北
京:
清华大学,1995.
[5]王兆华,侯正信,苏飞.全相位FFT频谱分析[J].通信学报,
2003,24（11:
16-19.
WANGZhaohua,HOUZhengxin,SUfei.AllphaseFFTspectrumanalysis[J].JounalofCommunications,2003,24（11:
16-19.[6]丁康,朱小勇,谢明,等.离散频谱综合相位差校正法[J].振动
工程学报,2002,15（1:
114-117.
DINGKang,ZHUXiaoyong,XIEMing,etal.Synthesizedcorrectingmethodofphasedifferenceondiscretespectrum[J].JournalofVibrationEngineering,2002,15（1:
114-117.
[7]朱小勇,丁康.离散频谱校正方法的综合比较[J].信号处理,2001,
17（1:
91-97.
ZHUXiaoyong,DINGKang.Thesyntheticalcomparisionofcor-rectingmethodsondiscretespectrum[J].SignalProcessing,2001,17（1:
91-97.
[8]吴国乔,王兆华,黄晓红.离散频谱的全相位校正法[J].数据采
集与处理,2005,20（3:
286-290.
WUGuoqiao,WANGZhaohua,HUANGXiaohong.Allphasecorrectionmethodfordiscretespectrum[J].JournalofDataAcqui-sition&Processing,2005,20（3:
286-290.
[9]RIFEDC,VINCENTGA.UseofthediscreteFouriertransform
inthemeasurementoffrequenciesandlevelsoftones[J].BellSystTechJ,1970,49（2:
197-228.
[10]PALMERLC.Coarsefrequencyestimationusingthediscrete
Fouriertransform[J].IEEETransInformTheory,1974,20（1:
104-109.
[11]齐国清,贾欣乐.基于DFT相位的正弦波频率和初相的高精度
估计算法[J].电子学报,2001,29（9:
1164-1167.
QIGuoqing,JIAXinle.High-accuracyfrequencyandphaseestimationofsignal-tonebasedonphaseofDFT[J].ActaElectronicSinica,2001,29（9:
1164-1167.
（责任编辑:
康鲁豫
作者简介:
黄晓红（1973-,女,黑龙江海伦人,副教授,博士,研究方向为数字信号处理及CPLD（E-mail:
tshxh@;
王兆华（1937-,男,上海人,教授,博士研究生导师,研究方向为数字信号处理及数字图像处理。
方法幅值
频率/Hz
相位/（°
全相位相位差法
0.9999999
4.3900000
39.9999999文献[6]中相位差法0.9988161
4.386835040.5486832理想值1.00000004.390000040.0000000
表3无噪情况下全相位相位差法与相位
差法参数校正结果比较
Tab.3ComparisonofparametercorrectionsbetweenapFFTphasedifferencemethodandphasedifference
methodunderno-noiseconditions
第28卷
电力自动化
设备
Highaccuracyparameterestimationofsinusoidalsignal
basedonall-phasespectrumanalysis
HUANGXiaohong1,WANGZhaohua2
（1.HebeiPolytechnicUniversity,Tangshan063000,China;
2.TianjinUniversity,Tianjin300072,China
Abstract:
TheapFFT（allphaseFFTspectrumanalysisanditsimplementationisputforwardtoim-provetheinputdatatruncationoftraditionalspectrumanalysis.ItistheoreticallyprovedthattheamplitudespectrumofapFFTspectrumanalysisisthesquareofthatoftraditionalspectrumanalysis,whichreducesleakagegreatly.TheapFFThasperfectphaseanalysisperformancewithlessimpactoffrequencyshift.ThephaseanalyzedbyapFFThasextremelysmallerrorevenifthesamplingfrequencyisnotthemultiplicationofsignalfrequency.AnapFFTphasedifferencealgorithmisproposedtoestimatetheamplitude,frequencyandphaseofsinusoidalsignal.Simulationresultsprovethatthealgorithmismoreaccurateandpracticablethantheexistingones.Theestimationisnearlyerror-freeunderno-noiseconditions,andthephaseerrorislessthan0.001%.
Keywords:
all-phasespectrumanalysis;parameterestimation;phasedifference;correction
电力期货市场的价格发现功能
刘思东1,2,杨洪明1,童小娇1
（1.长沙理工大学电气与信息工程学院,湖南长沙410076;
2.五邑大学数理系,广东江门529020
摘要:
价格发现功能是期货市场的基本功能之一,它在期货市场的发挥程度直接反映了期货市场的有效性。
为检测电力期货市场的价格发现能力,借助协整理论和误差修正模型分析了北欧电力期货市场的“简单效率”,并利用Granger因果检验、冲击反应分析、方差分解等方法对北欧电力期货市场与现货市场的领先-滞后关系进行了研究。
研究结果显示,期货电价和现货电价存在协整关系,且期货电价是现货电价的无偏估计量,北欧电力期货市场满足“简单效率”;期货市场的价格变化领先于现货市场的价格变化,且期货市场在价格发现功能中处于主导作用,北欧电力期货市场具有良好的价格发现功能。
关键词:
电力期货市场;价格发现;协整关系中图分类号:
TM73;F224.0
文献标识码:
A
文章编号:
1006-6047（200807-0057-050引言
电力期货市场概念自20世纪80年代被提出后,引起各国重视。
为了提高电力系统的运行效率,降低用电价格,各国分别建立起电力期货市场。
世界上第1份电力期货于1995年出现在北欧电力交
易所（NordicPowerExchange,美国、
澳大利亚和德国等国家的交易所也相继推出了电力期货[1]。
实践表明,电力期货市场的功能为欧美电力工业提供了较好的风险管理工具,降低了因现货市场电价高度波动带来的经营风险。
2004年初国家电监会在下发的《
关于建立东北区域电力市场的意见》中勾勒出东北区域电力市场
的远期目标:
实现所有市场主体参与的全面竞争,建立电力期货、期权等电力金融市场。
今年国家电监会指
升级会员 