一种频率自适应的滑窗DFT谐波检测方法技术

本发明专利技术公开了一种频率自适应的滑窗DFT谐波检测方法，首先采集电网电压信号，再利用数字锁相环技术同步电网电压相位及频率；接着采样负载电流信号并通过频率自适应的滑窗迭代方法对采集的负载电流信号x(t)进行DFT正变换，分别得到各次谐波的x轴和y轴分量，进而计算出第h次谐波的指令信号，最后累加各次谐波指令得到总谐波指令。本发明专利技术解决了电网电压频率波动时，有源电力滤波器实时谐波指令提取的难题，同时还能准确、快速的进行单次谐波指令提取，以保证有源电力滤波器的实时补偿效果。

【技术实现步骤摘要】
一种频率自适应的滑窗DFT谐波检测方法
本专利技术属于数字信号处理
，更为具体地讲，涉及一种频率自适应的滑窗DFT谐波检测方法。
技术介绍
随着电力电子技术的飞速发展，越来越多的大功率开关器件应用到电力电子设备中，以及非线性负载在电网上的广泛应用，其产生的大量无功电流和谐波电流对电网造成了严重污染。为了保证电网的安全可靠运行，保证较高的电能质量，对电网进行高效的谐波抑制显得尤为重要。目前，关于电力系统中谐波检测方法大体可分为基于频域的检测方法和基于时域的检测方法两大类。基于频域的检测方法以傅里叶变换为基础，包括快速傅里叶变换法(FFT)、离散傅里叶变换法(DFT)和迭代傅里叶变换方法(RDFT)。基于傅里叶变换的谐波检测方法可以实现谐波分次检测，使用方便，但存在一个周期的固有延时，实时性较差；迭代傅里叶变换法，动态响应速度方面则优于傅里叶变换法，但对于每次谐波分量的检测需要为其额外增加存储空间。基于时域谐波检测方法，包括基于时域瞬时无功功率理论(p-q)检测方法、基于基波同步旋转变换(d-q)检测方法、基于谐波同步旋转变换检测方法。基于时域的谐波检测方法属于实时检测方法，一般动态响应速度快，但存在分次谐波检测和计算量大的矛盾。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种频率自适应的滑窗DFT谐波检测方法，解决了电网频率波动时有源电力滤波器(APF)对实时谐波指令提取的难题，同时还能准确、快速的检测各次谐波分量，以保证APF的实时补偿效果。为实现上述专利技术目的，本专利技术一种频率自适应的滑窗DFT谐波检测方法，其特征在于，包括以下步骤：(1)、在固定采样周期Ts下，采集电网电压信号，通过数字锁相环技术，同步电网电压相位θk及频率fk；(2)、根据步骤(1)所述的电网电压频率fk，确定第k次采样的迭代窗口的大小Nk；其中，符号[]表示取整，f0为电网电压理论频率50Hz，N0为工频时迭代窗口大小；(3)、采集负载电流信号x(t)，再通过频率自适应的滑窗迭代方法对负载电流信号x(t)进行DFT正变换，分别得到负载电流的第h次谐波的x轴分量Ahx(k)与y轴分量Ahy(k)；Ahx(k)＝Ahx(k-1)+Ahxin-Ahxout(2)Ahy(k)＝Ahy(k-1)+Ahyin-Ahyout(3)其中，Ahxin、Ahxout分别为当前进窗量、出窗量的周期平均值的x轴分量，Ahyin、Ahyout分别为当前进窗量、出窗量的周期平均值的y轴分量；(4)、根据步骤(1)所述的电网电压相位θk与步骤(3)所述的x轴分量Ahx(k)与y轴分量Ahy(k)，计算第k拍时的第h次谐波指令xh(k)；xh(k)＝Ahx(k)·cos(h·θk)+Ahy(k)·sin(h·θk)(4)(5)、将步骤(4)中得到的h次谐波指令信号xh(k)进行累加得到第k拍的总谐波指令x'(k)；其中，所述的步骤(3)中，进窗量和出窗量的周期平均值的计算方法为：(3.1)、根据最小电网电压频率fmin确定迭代窗口大小Nmax；fmin＝min{fk}(7)其中，min{}表示取最小运算；(3.2)、在负载电流的每相电流中设置三个大小为Nmax的数组，即x[Nmax]、Table[Nmax]、table[Nmax]，分别储存迭代窗口内负载电流信号采样值的幅值x(k)、幅值相对应的相位θk、及对应迭代窗口大小Nk；(3.3)、设置迭代计数变量n，其初始值为0，当其累加到(Nmax-1)后置0；设置变量p，变量p表示迭代窗口内末尾数据在数组x[Nmax]、Table[Nmax]、table[Nmax]中的位置；(3.4)、根据步骤(3.3)所述的迭代计数变量n，以及第k-1拍的迭代窗口大小Nk-1，计算出第k-1拍的变量p的值；(3.5)、根据第k拍的迭代窗口大小Nk和第k-1拍的迭代窗口大小Nk-1，分别计算出第k拍滑入迭代窗口的数据个数ni以及滑出迭代窗口的数据个数no；(3.6)、根据步骤(3.5)计算的第k拍滑入和滑出迭代窗口的数据个数ni和no，以及步骤(3.4)中的不同大小的变量p，按照以下公式，分别计算出Ahxin、Ahxout、Ahyin、Ahyout；其中，p1为滑入窗口的数据在迭代数组中的位置；其中，p2为滑出窗口的数据在迭代数组中的位置。本专利技术的专利技术目的是这样实现的：一种频率自适应的滑窗DFT谐波检测方法，首先采集电网电压信号，再利用数字锁相环技术同步电网电压相位及频率；接着采样负载电流信号并通过频率自适应的滑窗迭代方法对采集的负载电流信号x(t)进行DFT正变换，分别得到各次谐波的x轴和y轴分量，进而计算出第h次谐波的指令信号，最后累加各次谐波指令得到总谐波指令。本专利技术解决了电网电压频率波动时，有源电力滤波器实时谐波指令提取的难题，同时还能准确、快速的进行单次谐波指令提取，以保证有源电力滤波器的实时补偿效果。附图说明图1是本专利技术一种频率自适应的滑窗DFT谐波检测方法的框图图2是电压频率波动时采用非频率自适应谐波检测方法提取的7次谐波的时域仿真波形图及其频谱图；图3是电压频率波动时采用本算法提取的7次谐波的时域仿真图及其频谱图；图4是电压频率由50Hz变化为55Hz时采用本算法与未采用本算法提取各次谐波指令的误差图。具体实施方式下面就具体实施方式进行详细描述，以便本领域的技术人员更好地理解本专利技术。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本专利技术的主要内容时，这些描述在这里将被忽略。实施例为便于叙述，本实施例以提取单相负载电流的谐波检测为例，叙述当频率波动时提取单相负载电流x(t)的非3的倍数次的奇数次谐波波指令(即提取5、7、11、13、…、49次谐波指令)的具体方法：本实施例中，假设采样频率为fs＝10kHz，电网工频f0＝50Hz，周期T0＝0.02s，电网电压频率fk波动范围为(1±0.1)f0,即45Hz～55Hz，工频时一个迭代窗口内的大小下面进行详细说明：图1是本专利技术一种频率自适应的滑窗DFT谐波检测方法的框图。在本实施例中，如图1所示，一种频率自适应的滑窗DFT谐波检测方法，包括以下步骤：S1、在固定采样周期Ts下，采集电网电压信号，通过数字锁相环技术，同步电网电压相位θk及频率fk；在本实施例中，固定采样周期为Ts＝0.1ms；S2、根据步骤S1所述的电网电压频率fk，，确定第k次采样的迭代窗口的大小Nk；其中，符号[]表示取整，f0为电网电压理论频率50Hz，N0为工频时迭代窗口大小；S3、采集负载电流信号x(t)，再通过频率自适应的滑窗迭代方法对负载电流信号x(t)进行DFT正变换，分别得到负载电流的第h次谐波的x轴分量Ahx(k)与y轴分量Ahy(k)；Ahx(k)＝Ahx(k-1)+Ahxin-Ahxout(2)Ahy(k)＝Ahy(k-1)+Ahyin-Ahyout(3)其中，Ahxin、Ahxout分别为当前进窗量、出窗量的周期平均值的x轴分量，Ahyin、Ahyout分别为当前进窗量、出窗量的周期平均值的y轴分量；其中，对负载电流信号x(t)进行DFT正变换，从而得到进窗量、出窗量的周期平均值的过程的具体步骤包括：S3.1、设电网电压波动范围确本文档来自技高网...

【技术保护点】
一种频率自适应的滑窗DFT谐波检测方法，其特征在于，包括以下步骤：(1)、在固定采样周期Ts下，采集电网电压信号，通过数字锁相环技术，同步电网电压相位θk及频率fk；(2)、根据步骤(1)所述的电网电压频率fk，确定第k次采样的迭代窗口的大小Nk；其中，符号[]表示取整，f0为电网电压理论频率50Hz，N0为工频时迭代窗口大小；(3)、采集负载电电流信号x(t)，再通过频率自适应的滑窗迭代方法对负载电流信号x(t)进行DFT正变换，分别得到负载电流的第h次谐波的x轴分量Ahx(k)与y轴分量Ahy(k)；Ahx(k)＝Ahx(k‑1)+Ahxin‑Ahxout       (2) Ahy(k)＝Ahy(k‑1)+Ahyin‑Ahyout       (3) 其中，Ahxin、Ahxout分别为当前进窗量、出窗量的周期平均值的x轴分量，Ahyin、Ahyout分别为当前进窗量、出窗量的周期平均值的y轴分量；(4)、根据步骤(1)所述的电网电压相位θk与步骤(3)所述的x轴分量Ahx(k)与y轴分量Ahy(k)，计算第k拍时的第h次谐波指令xh(k)；xh(k)＝Ahx(k)·cos(h·θk)+Ahy(k)·sin(h·θk)     (4)(5)、将步骤(4)中得到的h次谐波指令信号xh(k)进行累加得到第k拍的总谐波指令x'(k)；...

【技术特征摘要】
1.一种频率自适应的滑窗DFT谐波检测方法，其特征在于，包括以下步骤：(1)、在固定采样周期Ts下，采集电网电压信号，通过数字锁相环技术，同步电网电压相位θk及频率fk；(2)、根据步骤(1)所述的电网电压频率fk，确定第k次采样的迭代窗口的大小Nk；其中，符号[]表示取整，f0为电网电压理论频率50Hz，N0为工频时迭代窗口大小；(3)、采集负载电流信号x(t)，再通过频率自适应的滑窗迭代方法对负载电流信号x(t)进行DFT正变换，分别得到负载电流的第h次谐波的x轴分量Ahx(k)与y轴分量Ahy(k)；Ahx(k)＝Ahx(k-1)+Ahxin-Ahxout(2)Ahy(k)＝Ahy(k-1)+Ahyin-Ahyout(3)其中，Ahxin、Ahxout分别为当前进窗量、出窗量的周期平均值的x轴分量，Ahyin、Ahyout分别为当前进窗量、出窗量的周期平均值的y轴分量；(4)、根据步骤(1)所述的电网电压相位θk与步骤(3)所述的x轴分量Ahx(k)与y轴分量Ahy(k)，计算第k拍时的第h次谐波指令xh(k)；xh(k)＝Ahx(k)·cos(h·θk)+Ahy(k)·sin(h·θk)(4)(5)、将步骤(4)中得到的h次谐波指令信号xh(k)进行累加得到第k拍的总谐波指令x'(k)；其中，所述步骤(3)中，进窗量和...

【专利技术属性】
技术研发人员：谢川，李凯，郑宏，邹见效，徐红兵，杨磊，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51