专利名称：：电气信号频率的数字化实时检测方法
技术领域：
：本发明属于电工
技术领域：
，涉及一种电气信号频率的数字化实时检测方法，该方法特别适应于电网电压工作频率和电气设备运行频率的检测。
背景技术：
：频率是电力系统的主要技术指标之一，也是电机等电气设备运行的主要监视参数。电频率测量是电力系统及电气设备运行与多种控制、调节的基础。用电气信号交流采样值的数值运算来估计信号频率的软件测频技术，是目前电气频率测量的主流技术。已有的测频方法有周期法、解析法、误差最小化原理算法、DFT类算法、正交去调制法等。这些方法在测量速度、精度、计算量、实现难易程度、抑制谐波及噪声能力等方面各有优劣，适合于周期信号及频率变化缓慢的准周期信号的频率测量；对于频率变化较快的动态信号，这些方法并不适用。电力系统正常运行时，电网频率随电力负荷变化而不断波动，变化速率较慢；电网发生事故出现有功功率不平衡时，系统频率变化速率加快，这时电网频率及其变化速率的快速跟踪测量对事故处理和电网安全十分重要。电机等旋转电气设备在起动、运行及事故发生时的频率变化也可能相当快，其快速跟踪测量是其并列控制、转速调节和事故处理的重要依据。因此，对变化较快的电频率进行快速跟踪测量十分重要。可满足电网和电气设备在不同运行方式下的频率测量要求的频率测量方法是迫切需要的，极具工程实用价值。
发明内容本发明要解决的技术问题是提出一种电气信号频率的数字化实时检测方法，以满足电力系统和电气设备在各种运行方式下对频率测量精度和响应速度的要求。本发明的技术解决方案如下一种电气信号频率的数字化实时检测方法，其特征在于，包括以下步给定信号频率为50Hz时电压信号每周波的采样点数M值和p值，p取2或3，窗函数长度W-;M，以固定采样周期7^=0.02/^秒对电网或电气设备的电压进行不间断采样，得到电压采样值序列{};步骤2):从"(")序列中第A点开始取W点得，从第A:+1点起取TV点得步骤3):对W|(n)和"2(")分别加p阶矩形自巻积窗以求取初相位角计算的加窗离散傅立叶变换所得A(")变换后实部为<formula>formulaseeoriginaldocumentpage5</formula>所得变换后虚部为.2;r"":一5^(")'w(").sin(—.")，得^(")初相位角0)'=tan—"二01/计算^()的加窗离散傅立叶变换所得^(w)变换后实部为22;r仏《=:》2(")'w(")'cos(;'");TV=07V所得"2(")变换后虚部为仏,=丄y"V(w),.sin(^.w);2/2、7、/、；v'得W2(")初相位角①^tan22/其中2阶矩形自巻积窗—w)的表达式为;2"阪J、TV22--(—<"S7V—1)TV23阶矩形自巻积窗的时域表达式为:121—w+—w22-"2—6-)111TV(os"——1)3(—SwS--1);27V步骤4):计算信号的频率,=1+①2—^AT;2;t7;再更新k,即令l^k+l，返回步骤2)不断实时求取信号的频率,作为改进，根据当前频率估计值/自适应调整窗函数长度N，所述步骤l)中的采样点数M取32128间的整数。该检测方法的基本思想(原理)是，把电气信号划分成一些小的时间间隔，加适当的窗函数后，用傅立叶变换分析每一个时间间隔内采样数据，以便确定在那个时间间隔内的信号频率；同时，根据信号频率的变化，自适应地调整采样数据长度。方法的具体工作原理如下设电气信号x(O在较短时间间隔内近似为简单正弦信号，其数学表达式为=^cos(叫+a)(1)式中4、a、a分别为信号幅值、频率和初相位。40的傅立叶变换为=￡x(y)e—声^=;r单—^)e州+5—+^](2)设窗函数vKr)为偶对称的实函数，窗宽r,在区间[-772,7V2)外的值为0，其傅立叶变换为一)=r(3)Loo根据傅立叶变换的性质，『(w)亦为偶对称的实函数。对x(/)加滑动窗vv(卜r)(r为滑动窗的中心位置)，再作傅立叶变换，即得到x(/)的短时傅立叶变换(ShortTimeFourierTransform,STFT):^y,r)=5TF7;(6;,r)=「十〃、一(卜r)e-,(4)由式(2)、(3)及巻积定理0=—^(—附)+;r4『("+q)e"(肌+竿,)(5)由于一般窗函数相当于低通滤波器，只考虑^>0时，若忽略负频率分量的影响，则X(w,r)=;r单ff(w—^)(—肌+卬,)(6)A^,r)可认为是在[r-772,r+772;)范围内的信号的局部频谱。实际计算中，需将信号和窗函数都进行离散化。设在窗宽r内对x(/)均匀采样的点数为W，则采样周期7;^77/V，STFT的角频率间隔A"二2;r/r。4=mA6;、r=Ars、Z=w7;，由式(4)得STFT的离散化形式为A72-l根据STFT性质，X(w,A:)是式(6)在时、频域分别以Aw、7;抽样离散2r,,、-/——w("十A)化的结果，即设时间窗宽度r与信号周期7;(7^=2^/^)的关系为rat.7:7:=/+CT其中/^row^(7V7;/r》，为最接近M;/7;的整数；a=，&—roww^(i^)为取整后的余数。于是/9+O"当,/时，根据式(6)、(8)，有其中"AW为％(AW的相位角8(10)(11)于是<formula>formulaseeoriginaldocumentpage9</formula>所以有<formula>formulaseeoriginaldocumentpage9</formula>(13)信号频率为<formula>formulaseeoriginaldocumentpage9</formula>(14)基于上述检测原理，电气信号频率的数字检测技术方案如下1)以固定采样频率对电网(或电气设备)电压进行等间隔采样，得到电压采样值序列；2)设M为当信号频率为50Hz时电压信号每周波的采样点数，/7为矩形自巻积窗阶数(p取值为2或3，这样既可达到很高的测量精度，又保持了算法的简洁性)。从电压采样值序列中从第^点开始取长度为pM的采样数据，加p阶矩形自巻积窗，并对加窗数据进行离散傅立叶变换(DFT)，得到信号初相位角<formula>formulaseeoriginaldocumentpage9</formula>3)在电压采样值序列中从第A+1点开始同样取长度为PM的采样数据，加p阶矩形自巻积窗，并对加窗数据进行离散傅立叶变换(DFT)，得到信号初相位角<formula>formulaseeoriginaldocumentpage9</formula>4)根据式(14)求取信号频率，并根据信号频率自适应调整采样数据窗长度，实现电气信号频率的自适应测量。本发明的优点与效果-1)适用范围广。即可用于电网频率的实时测量，也可用于电机等电气设备的实时测量；即适应于变化较慢的信号频率测量，也适应于快速变化的频率测量。因此，在电力系统和电气设备的监测、控制和保护中，凡需要实时测量电压、电流信号频率的场合，均可采用本专利。2)采样周期固定，采样频率不需随信号频率的变化而调整，实现简单。3)算法不需要迭代运算，可一次得到信号频率的估计值，响应时间快。综上所述，该方法实现简单、计算量小，测频范围大，在信号频率缓慢变化和快速变化时，均具有较好的测量精度和跟踪速度，因此，本发明在电网和电气设备的频率实时检测方面，具有显著的理论和技术优势，具有极高的应用价值。具体实施例方式下面结合附图和具体实施例对本发明作进一步的说明。实施例1:给定M值(M—般取32128间的整数)和p值(p取2或3)。以固定采样周期^=0.02/^秒对电网(或电气设备)电压进行不间断采样，得到电压采样值序列{"(")};步骤2):令A^pM，从M(")序列中第A点开始取iV点得A(w)，从第A:+1点起取W点得^(");步骤3):<formula>formulaseeoriginaldocumentpage11</formula>由下式计算得到信号的频率<formula>formulaseeoriginaldocumentpage12</formula>步骤6):当信号频率变化较大时，根据当前频率估计值/，自适应调整窗函数长度<formula>formulaseeoriginaldocumentpage12</formula>其中ra"m/[函数的含义为，取最接近的整数。步骤7):令k-k+l，更新序列m,(")和^(")，重复步骤2)步骤6),不断实时求取电气信号的频率。权利要求1.一种电气信号频率的数字化实时检测方法，其特征在于，包括以下步骤步骤1)给定信号频率为50Hz时电压信号每周波的采样点数M值和p值，p取2或3，窗函数长度N＝pM，以固定采样周期TS＝0.02/N秒对电网或电气设备的电压进行不间断采样，得到电压采样值序列{u(n)}；步骤2)从u(n)序列中第k点开始取N点得u1(n)，从第k+1点起取N点得u2(n)；步骤3)对u1(n)和u2(n)分别加p阶矩形自卷积窗w(n)以求取初相位角计算u1(n)的加窗离散傅立叶变换所得u1(n)变换后实部为<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>U</mi><mrow><mn>1</mn><mi>R</mi></mrow></msub><mo>=</mo><mfrac><mn>2</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>u</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mo>·</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math>id="icf0001"file="A2009100438330002C1.tif"wi="73"he="11"top="150"left="23"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>所得u1(n)变换后虚部为<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>U</mi><mrow><mn>1</mn><mi>I</mi></mrow></msub><mo>=</mo><mfrac><mn>2</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>u</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>sin</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mo>·</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math>id="icf0002"file="A2009100438330002C2.tif"wi="69"he="11"top="166"left="23"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>得u1(n)初相位角<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>Φ</mi><mn>1</mn></msub><mo>=</mo><msup><mi>tan</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mfrac><msub><mi>U</mi><mrow><mn>1</mn><mi>R</mi></mrow></msub><msub><mi>U</mi><mrow><mn>1</mn><mi>I</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math>id="icf0003"file="A2009100438330002C3.tif"wi="35"he="13"top="165"left="132"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>计算u2(n)的加窗离散傅立叶变换所得u2(n)变换后实部为<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>U</mi><mrow><mn>2</mn><mi>R</mi></mrow></msub><mo>=</mo><mfrac><mn>2</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>u</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mo>·</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math></maths>所得u2(n)变换后虚部为<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msub><mi>U</mi><mrow><mn>2</mn><mi>I</mi></mrow></msub><mo>=</mo><mfrac><mn>2</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>u</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>·</mo><mi>sin</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>N</mi></mfrac><mo>·</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math>id="icf0005"file="A2009100438330002C5.tif"wi="71"he="12"top="210"left="84"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>得u2(n)初相位角<mathsid="math0006"num="0006"><math><![CDATA[<mrow><msub><mi>Φ</mi><mn>2</mn></msub><mo>=</mo><msup><mi>tan</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mfrac><msub><mi>U</mi><mrow><mn>2</mn><mi>R</mi></mrow></msub><msub><mi>U</mi><mrow><mn>2</mn><mi>I</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math>id="icf0006"file="A2009100438330002C6.tif"wi="36"he="13"top="225"left="69"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中2阶矩形自卷积窗w(n)的表达式为<mathsid="math0007"num="0007"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><mfrac><mrow><mn>2</mn><mi>n</mi></mrow><mi>N</mi></mfrac></mtd><mtd><mrow><mo>(</mo><mn>0</mn><mo>≤</mo><mi>n</mi><mo>≤</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>2</mn><mo>-</mo><mfrac><mrow><mn>2</mn><mi>n</mi></mrow><mi>N</mi></mfrac></mtd><mtd><mrow><mo>(</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo><</mo><mi>n</mi><mo>≤</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math></maths>3阶矩形自卷积窗的时域表达式为<mathsid="math0008"num="0008"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>n</mi><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>n</mi></mtd><mtd><mrow><mo>(</mo><mn>0</mn><mo>≤</mo><mi>n</mi><mo>≤</mo><mfrac><mi>N</mi><mn>3</mn></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><msup><mi>n</mi><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>n</mi><mo>-</mo><mrow><mo>(</mo><mfrac><msup><mi>N</mi><mn>2</mn></msup><mn>6</mn></mfrac><mo>-</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><mfrac><mi>N</mi><mn>3</mn></mfrac><mo>≤</mo><mi>n</mi><mo>≤</mo><mfrac><mrow><mn>2</mn><mi>N</mi></mrow><mn>3</mn></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>n</mi><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>-</mo><mi>N</mi><mo>)</mo></mrow><mi>n</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msup><mi>N</mi><mn>2</mn></msup><mo>-</mo><mi>N</mi><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>N</mi></mrow><mn>3</mn></mfrac><mo>≤</mo><mi>n</mi><mo>≤</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math></maths>步骤4)计算信号的频率<mathsid="math0009"num="0009"><math><![CDATA[<mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><mfrac><mi>p</mi><msub><mi>NT</mi><mi>s</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>Φ</mi><mn>2</mn></msub><mo>-</mo><msub><mi>Φ</mi><mn>1</mn></msub></mrow><mrow><mn>2</mn><mi>π</mi><msub><mi>T</mi><mi>s</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math>id="icf0009"file="A2009100438330003C3.tif"wi="40"he="11"top="127"left="67"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>再更新k，即令k＝k+1，返回步骤2)不断实时求取信号的频率。2.根据权利要求1所述的电气信号频率的数字化实时检测方法，其特征在于，根据当前频率估计值/自适应调整窗函数长度N，3.根据权利要求1所述的电气信号频率的数字化实时检测方法，其特征在于，所述步骤l)中的采样点数M取32128间的整数。全文摘要针对电网和电气设备各种运行工况对频率测量性能的要求，本发明提供一种电气信号频率的数字化实时检测方法，该方法基于短时傅里叶变换，其主要步骤为1)对电气信号进行等间隔采样，并取采样序列中的2小段加矩形自卷积窗，然后进行离散傅里叶变换，求得每段信号的初相角；2)根据2小段信号的相角差，估计信号的频率；3)保持采样频率不变，根据信号估计频率，自适应地确定分析信号段的长度；4)按相同步骤再次自适应检测电气信号实时频率。该方法实现简单、计算量小，测频范围大，在信号频率缓慢变化和快速变化时，均有较好的测量精度和跟踪速度，极具工程实用价值。文档编号G01R23/00GK101603984SQ20091004383公开日2009年12月16日申请日期2009年7月1日优先权日2009年7月1日发明者何怡刚,江亚群,纯黄申请人:湖南大学