基于改进的SmartDFT算法的非平衡系统频率估计方法技术方案

在非平衡系统中，估计非圆信号的频率、幅度和相位是非常重要的非线性问题。本专利中，我们扩展了原始的Smart DFT技术（SDFT），使此项技术不仅可以应用于实值正弦信号，也可以处理复值非圆信号。基于连续DFT基频分量间的线性预测（LP）性质，应用最小二乘框架可以减小模型的均方误差，得到改进的复值最小二乘算法（CLS）。同时，专利还提出了一种复值改进的Pisarenko谐波分解算法（CRPHD）,此方法可以移除噪声的干扰获得精确的频率估计，可以有效应用在含有噪声的非平衡三相电力系统中。

Frequency estimation method of unbalanced system based on improved SmartDFT algorithm

It is very important to estimate the frequency, amplitude and phase of non circular signals in unbalanced systems. This patent, we extended the original Smart DFT Technology (SDFT), so that this technology not only can be applied to real valued sinusoidal signal, can also handle complex noncircular signals. Based on the linear prediction (LP) property of the baseband components of continuous DFT, the mean square error of the model can be reduced by using the least square frame, and the improved complex least squares (CLS) algorithm is obtained. At the same time, the patent also presents an improved complex valued Pisarenko harmonic decomposition algorithm (CRPHD), this method can remove the interference noise to obtain accurate frequency estimation can be effectively applied in unbalanced three-phase power system with noise in.

【技术实现步骤摘要】
基于改进的SmartDFT算法的非平衡系统频率估计方法
本专利技术涉及非平衡系统频率估计
，特别是涉及基于改进的SmartDFT算法的非平衡系统频率估计方法。
技术介绍
在附加白高斯噪声下，估计正弦信号或复指数信号的频率、幅度、相位是非常重要的非线性问题，在电力系统分析、无线通信、雷达信号监测和语音分析中有着广泛的应用。在不同的应用中，众多文献也提出了相对应的估计方法，如基于DFT算法，最小二乘算法，自适应陷波器，卡尔曼滤波器和它的拓展，最大后验概率算法和基于子空间投影的算法。典型的，标准的N点DFT操作在频谱上可以得到2π/N的分辨率，但是在异步采样下是不准确的。传统方法上如果想要提高基于DFT的频率估计算法的精确度，额外的计算复杂度不可避免。在已有文献中，提出过了一种SmartDFT(SDFT)技术改善实时正弦信号频率估计的准确度。也有文献进一步应用了复值最小二乘框架改善了SDFT的估计效果。基于窗口DFT的插值算法可以获得频率估计的重要改善，因为信号序列的窗口可以减少频谱泄漏和频谱插值。DFT算法的相角误差可以估计相应模拟信号的真实频率,由非相干DFT采样引起的泄漏效果也可以用于频率估计。进一步，有文献提出采用两个阶段去改善估计效果，第一阶段实行N点DFT的粗估计，接着在频谱峰值附近进行第二阶段的精确估计。非圆信号广泛应用于描述非平衡系统的动态状态，例如非平衡三相电力系统、I-Q不平衡通信系统等等。基于恰当的非圆信号统计分析，利用非平衡电压的非圆统计特性，众多文献提出了一系列的估计模型，标准的频率估计算法也可以进一步扩展为通用的形式。无论是频域的SDFT技术还是时域的RPHD技术还无法处理非平衡系统中的复值非圆信号，本专利意义在于基于复值非圆信号的特性，实现与以上技术的结合，完成在有噪条件下，对复值非圆信号频率精确的估计和追踪。
技术实现思路
为了解决上述存在的问题，本专利技术提供基于改进的SmartDFT算法的非平衡系统频率估计方法，扩展了原始的SmartDFT技术(SDFT)，使此项技术不仅可以应用于实值正弦信号，也可以处理复值非圆信号。基于连续DFT基频分量间的线性预测(LP)性质，应用最小二乘框架可以减小模型的均方误差，得到改进的复值最小二乘算法(CLS)。同时，专利还提出了一种复值改进的Pisarenko谐波分解算法(CRPHD),此方法可以移除噪声的干扰获得精确的频率估计，可以有效应用在含有噪声的非平衡三相电力系统中，为达此目的，本专利技术提供基于改进的SmartDFT算法的非平衡系统频率估计方法，具体估计方法如下：步骤1计算非圆信号第k点频率分量：没有噪声和谐波干扰的非平衡系统中的非圆信号可以表示为：参数α和β分别定义为这里A,B和是未知的确定性常量分别表示当前信号的幅度和相位，定义信号频率这里k∈{1,2,…,N-1}，δ∈(0,1)分别是未知系统频率的整数和小数部分，目标是根据N点DFT变换后的第k频点序列{Xk(m)}找到ω0；信号x(n)第k点DFT频率分量如下：步骤2无噪条件下的SDFT算法：现在提出频率估计算法：用定义参数的估计值；定义公式(1)中的指数内核为；如此，根据公式(1)和(2)，找到如下的关系：ak(m+1)＝rak(m)和bk(m+1)＝r-1bk(m)，以及：对公式(3)进行代数变换后，可以找到连续DFT分量的线性预测关系：Xk(m+1)+Xk(m-1)＝μXk(m)，这里可以通过下式估计因此，SDFT算法可以通过计算来估计系统频率这里表示取复数的实部；步骤3有噪条件下CLS优化处理：现在考虑在噪声环境下的非平衡系统：噪声q(n)是均值为零，方差为的高斯白过程，继续使用复值最小二乘CLS框架改善了智能离散傅里叶变换SDFT算法的性能，同样，将此方法应用于非平衡系统信号上，基于CLS的算法利用了序列的线性预测性质：然后最小化误差e(n)的平方和：这里Qk(m)是噪声q(n)的DFT变换；组合L点序列定义向量相应的误差向量为此误差向量的共轭转置形式可表示为CLS框架作用在于找到一个最优的值以最小化均方误差：令则提出的CLS改进SDFT的算法可以表示如下：步骤4有噪条件下CRPHD实现进一步抗噪优化：另一方面，重新考虑均方误差值继续推导如下：这里是Q(m)的方差，由于公式(6)的第二项是关于的噪声项，一般情况下并不在取得最小值，这也是CLS算法估计含噪信号的劣势所在，在非平衡系统中使用复值RPHD方法CRPHD，如公式(6)所示，为了移除噪声对频率估计误差的影响，依据来最小化重新定义一个新的代价函数为：步骤5寻求最优的值最小化计算求解频率估计值：为了找到最优的值最小化求偏导令其等于零：这里：尽管公式(8)的最后推导式有两个根，但只有一个根可以表示频率估计量，如此推导得到了相对于CLS更加抗噪的，适用于复值非圆信号的CRPHD频率估计方法，如下所示：本专利技术的进一步改进，步骤二，通过下面的方法根据Xk(m+1),Xk(m)和Xk(m-1)的值可以估计k和δ；1)基于DFT最大点的位置进行频率粗估计：2)估计δ：定义公式(1)前一项为a(m)，第二项为b(m)，这样公式(1)中的Xk(m)表示为：Xk(m)＝ak(m)+bk(m)。本专利技术基于改进的SmartDFT算法的非平衡系统频率估计方法，本申请扩展了原始的SmartDFT技术(SDFT)，使此项技术不仅可以应用于实值正弦信号，也可以处理复值非圆信号。基于连续DFT基频分量间的线性预测(LP)性质，应用最小二乘框架可以减小模型的均方误差。同时，提出了一种复值改进的Pisarenko谐波分解算法(CRPHD)，此方法可以移除噪声的干扰获得精确的频率估计。本专利技术在含有噪声的非平衡三相电力系统中，经过仿真显示本申请所提出算法的优越的抗噪性能。附图说明图1为不同信噪比下频率估计算法的均方误差对比图；图2为两种算法在频率波动下的频率跟踪性能图(a)系统频率经历频率调制(b)系统频率经历了上升和下降；图3真实环境非平衡三相电力系统中的频率估计效果图(a)三相电压幅值波动(b)频率估计效果。具体实施方式下面结合附图与具体实施方式对本专利技术作进一步详细描述：本专利技术提供基于改进的SmartDFT算法的非平衡系统频率估计方法，扩展了原始的SmartDFT技术(SDFT)，使此项技术不仅可以应用于实值正弦信号，也可以处理复值非圆信号。基于连续DFT基频分量间的线性预测(LP)性质，应用最小二乘框架可以减小模型的均方误差，得到改进的复值最小二乘算法(CLS)。同时，专利还提出了一种复值改进的Pisarenko谐波分解算法(CRPHD),此方法可以移除噪声的干扰获得精确的频率估计，可以有效应用在含有噪声的非平衡三相电力系统中。具体实施例1如下：在这部分，我们将利用MATLAB平台展示专利算法的仿真结果。我们分别设置采样频率为fs＝1600，基点频率为f0＝50，这样，为了计算DFT的基础频率，系统电压的采样点需要设置为N＝fs/f0＝32。同时，我们将估计算法的窗长设置为L＝15。首先，我们评估在噪声环境下CRPHD算法相对于CLS优越的估计性能。非平衡系统电压真实频率f＝f0+Δf设置在51Hz，将不同信噪比条件下的噪声加入系统中，图本文档来自技高网...

【技术保护点】
基于改进的SmartDFT算法的非平衡系统频率估计方法，具体估计方法如下，其特征在于：步骤1计算非圆信号第k点频率分量：没有噪声和谐波干扰的非平衡系统中的非圆信号可以表示为：

【技术特征摘要】
1.基于改进的SmartDFT算法的非平衡系统频率估计方法，具体估计方法如下，其特征在于：步骤1计算非圆信号第k点频率分量：没有噪声和谐波干扰的非平衡系统中的非圆信号可以表示为：参数α和β分别定义为这里A,B和是未知的确定性常量分别表示当前信号的幅度和相位，定义信号频率这里k∈{1,2,…,N-1}，δ∈(0,1)分别是未知系统频率的整数和小数部分，目标是根据N点DFT变换后的第k频点序列{Xk(m)}找到ω0；信号x(n)第k点DFT频率分量如下：步骤2无噪条件下的SDFT算法：现在提出频率估计算法：用定义参数的估计值；定义公式(1)中的指数内核为；如此，根据公式(1)和(2)，找到如下的关系：ak(m+1)＝rak(m)和bk(m+1)＝r-1bk(m)，以及：对公式(3)进行代数变换后，可以找到连续DFT分量的线性预测关系：Xk(m+1)+Xk(m-1)＝μXk(m)，这里可以通过下式估计因此，SDFT算法可以通过计算来估计系统频率这里表示取复数的实部；步骤3有噪条件下CLS优化处理：现在考虑在噪声环境下的非平衡系统：噪声q(n)是均值为零，方差为的高斯白过程，继续使用复...

【专利技术属性】
技术研发人员：王开，柳旭，夏亦犁，裴文江，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32