发明名称 汉宁窗函数连续频谱内插电力谐波参数获取方法
摘要 汉宁(Hanning)窗函数连续频谱内插电力谐波参数获取方法适用于对电力网电压、电流的谐波分析与监测。首先使用线性调频Z变换(Chirp-ZTransform)或称CZT从含有谐波的电力信号中高精度提取基波信号参数(幅值、频率和相位)。然后从被分析电力信号中减去基波信号后加Hanning窗函数截断电力信号,并用FFT计算出剩余信号的频谱。再依据基波频率精确计算出各谐波的频率值。最后依据各谐波频率对Hanning窗函数在频域内插值,精确计算出各电力谐波的参数。本发明与加Hanning窗FFT双谱线插值拟合分析电力谐波方法有基本同等的估值精度,而计算量约为1/2。
申请公布号 CN102331526A 申请公布日期 2012.01.25
申请号 CN201110154993.4 申请日期 2011.06.10
申请人 中国矿业大学 发明人 唐轶;刘昊;秦媛倩;李辉;周斌;陈雷;杨洛;尹远;张跃;渐伟
分类号 G01R23/16(2006.01)I 主分类号 G01R23/16(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 1. 一种汉宁窗函数连续频谱内插电力谐波参数获取方法,其特征是该方法采用如下步骤:步骤a.采样被分析电力信号电压或电流,并按快速CZT算法流程计算其线性调频Z变换CZT值<img file="2011101549934100001DEST_PATH_IMAGE002.GIF" wi="108" he="22" />,取自然正整数,再由式1、式2和式3分别计算出电力信号的基波参数,幅值、频率和相位;估计基波幅值:<img file="2011101549934100001DEST_PATH_IMAGE004.GIF" wi="256" he="29" />式1估计基波频率值:<i>f</i><sub><i>1</i></sub><i>=(θ+k’Ф)/2</i><img file="2011101549934100001DEST_PATH_IMAGE006.GIF" wi="16" he="16" />式2估计基波相位值:<img file="2011101549934100001DEST_PATH_IMAGE008.GIF" wi="277" he="28" />式3这里,<i>M</i>为线性调频Z变换时在频域内的抽样点数;<i>k</i><sup><i>’</i></sup>为<i>M</i>个<i>X(k)</i>中取得最大值的<i>k</i>值;<img file="2011101549934100001DEST_PATH_IMAGE010.GIF" wi="14" he="20" />为起始采样点的角频率;<img file="2011101549934100001DEST_PATH_IMAGE012.GIF" wi="14" he="22" />为相邻两采样点之间的角频率差;<img file="2011101549934100001DEST_PATH_IMAGE014.GIF" wi="72" he="25" />为<img file="2011101549934100001DEST_PATH_IMAGE016.GIF" wi="42" he="25" />的虚部;<img file="2011101549934100001DEST_PATH_IMAGE018.GIF" wi="72" he="25" />为<img file="372547DEST_PATH_IMAGE016.GIF" wi="42" he="25" />的实部;步骤b. 对Hanning窗函数连续频谱在频域内抽样求得各次电力谐波的校正系数<img file="2011101549934100001DEST_PATH_IMAGE020.GIF" wi="22" he="26" />;<img file="2011101549934100001DEST_PATH_IMAGE022.GIF" wi="194" he="46" />式4式4中:<img file="2011101549934100001DEST_PATH_IMAGE024.GIF" wi="17" he="16" />为数字角频率;<img file="2011101549934100001DEST_PATH_IMAGE026.GIF" wi="92" he="28" />为一预定值,<img file="2011101549934100001DEST_PATH_IMAGE028.GIF" wi="112" he="28" />,<img file="2011101549934100001DEST_PATH_IMAGE030.GIF" wi="90" he="28" />,<img file="2011101549934100001DEST_PATH_IMAGE032.GIF" wi="72" he="26" />为电网基波额定频率;<i>Ts</i>为采样周期,采样频率<img file="2011101549934100001DEST_PATH_IMAGE034.GIF" wi="62" he="25" />等于电网基波额定频率的2<sup>i</sup>倍,i取自然正整数,i=1、2、…;<img file="2011101549934100001DEST_PATH_IMAGE036.GIF" wi="92" he="26" />为以电网p次谐波信号的实际数字角频率<img file="2011101549934100001DEST_PATH_IMAGE038.GIF" wi="97" he="29" />在Hanning窗函数连续频谱上在频域内抽样值,p为自然正整数;步骤c.从被分析电力信号采样值中减去基波信号采样值,并加Hanning窗<img file="2011101549934100001DEST_PATH_IMAGE040.GIF" wi="54" he="25" />截断,得序列<img file="2011101549934100001DEST_PATH_IMAGE042.GIF" wi="41" he="25" />,再对<img file="855088DEST_PATH_IMAGE042.GIF" wi="41" he="25" />进行快速傅立叶变换(FFT),得<img file="2011101549934100001DEST_PATH_IMAGE044.GIF" wi="68" he="25" />,最终由式5和式6分别计算出各次电力谐波的幅值和相位;P次谐波的幅值:<img file="2011101549934100001DEST_PATH_IMAGE046.GIF" wi="378" he="30" />式5P次谐波的相位为:<img file="2011101549934100001DEST_PATH_IMAGE048.GIF" wi="185" he="42" />式6式5和式6中:<img file="2011101549934100001DEST_PATH_IMAGE050.GIF" wi="26" he="26" />是p次谐波FFT离散主谱角频率<img file="2011101549934100001DEST_PATH_IMAGE052.GIF" wi="42" he="26" />与<img file="2011101549934100001DEST_PATH_IMAGE054.GIF" wi="22" he="28" />的数字角频率差;<img file="2011101549934100001DEST_PATH_IMAGE056.GIF" wi="21" he="28" />是电力p次谐波信号的FFT离散频谱的主谱<img file="DEST_PATH_IMAGE058.GIF" wi="74" he="26" />值的相位;<img file="DEST_PATH_IMAGE060.GIF" wi="86" he="20" />为频域采样间隔;k<sub>1</sub>为基波主谱谱线;k<sub>P</sub>=pk<sub>1</sub>为p次谐波主谱谱线。
地址 221116 江苏省徐州市大学路南1号中国矿业大学科技处