一种改进频域插值的动态功率分析方法技术

本发明专利技术公开了一种改进频域插值的动态功率分析方法，它包括：a、对实际电压信号和电流信号进行离散采样得到待分析的正弦信号序列X(n)；b、对待分析的正弦信号序列X(n)做快速傅里叶变换，计算频谱插值系数δ

【技术实现步骤摘要】
一种改进频域插值的动态功率分析方法
本专利技术属于动态功率分析，尤其涉及一种改进频域插值的动态功率分析方法。
技术介绍
我国经济的不断发展导致国内电力需求旺盛，远距离、大规模输电网络的建设以及各种非线性负载的大量使用，电能质量问题日益严重。我国电力系统标准频率为50Hz，电压、电流均为50Hz的标准正弦波，但是由于受到各种干扰影响，在实际电力系统中，电压电流波形总在不同程度上产生畸变。有效的检测基波电能成为专家和学者越来越关注的问题。目前，电力系统信号频率分析应用最为普遍的算法是快速傅里叶变换(FFT),然而当被检测信号频率发生变换导致非同步采样时，FFT存在栅栏效应和泄露现象，使得计算的电力信号参数——频率、幅值和相位等不准确，尤其是相位误差很大，无法满足准确测量要求。为了消除电网频率变化对测量的影响，可以采用硬件锁相环电路实现同步采样；也可以对FFT算法进行改进，对于频谱泄露通常采用加窗函数的方法解决，对于栅栏效应通常采用插值算法来改善。常用的窗函数有矩形窗、三角窗、Hanning窗、Blackman窗、Blackman-Harris窗、矩形卷积窗和最大旁瓣衰减速率窗等。常用的插值算法有单谱线差值、应用最普遍的双谱线差值和近几年比较流行的三谱线差值。但传统的加窗插值FFT算法计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。
技术实现思路
本专利技术要解决的技术问题是：提供一种改进频域插值的动态功率分析方法，以解决现有技术采用加窗插值FFT算法计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。本专利技术的技术方案是：一种改进频域插值的动态功率分析方法，它包括：步骤a，对实际电压信号u(t)和电流信号i(t)进行离散采样，采样频率为fs，采样点数为N，得到离散电压U0(n)和电流序列I0(n)，n＝0,1,2,…,N-1，采用待分析的正弦信号序列X(n)表示所采样的离散电压或离散电流序列；步骤b，对待分析的正弦信号序列X(n)做快速傅里叶变换，变换后的离散频谱序列X(k)，记X(k)的峰值谱线的序号k为k0，记X(k0)＝a+jb，X(k0+1)＝c+jd，X(k0－1)＝e+jf，计算频谱插值系数δ+和δ－；步骤c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数步骤d，计算待分析的正弦信号序列X(n)的频率、幅值和相位：步骤e，将离散电压、离散电流序列U0(n)、I0(n)分别代入步骤b至步骤d的公式中，得到电压序列的幅值和相角分别记为Au和θu，电流序列的幅值和相角分别记为Ai和θi，动态功率为P＝AiAucos(θu－θi)。步骤b所述频谱插值系数δ+和δ－的计算公式为：步骤c所述根据峰值谱线邻近的两根谱线幅度来计算加权后的频谱插值系数的表达式为：步骤d所述计算待分析的正弦信号序列X(n)的频率幅值A和相位θ的表达式分别为：本专利技术有益效果：本专利技术采集正弦电压信号、电流信号，用一种改进的频域插值方法计算两个方向的频谱插值系数，再经过加权处理得到最终的频谱插值系数，将最终的频率插值系数带入频率估计式中计算正弦信号的频率，然后计算电压、电流的幅值和相位，最后计算动态功率，从而实现对动态功率实时、准确的测量。本专利技术功率计量方法所采用的改进FFT频域插值频率估计算法比传统插值算法数学模型简单，计算量小、抗干扰能力强，特别当峰值谱线邻近的两个峰值谱线收到严重干扰，产生错位时，本专利技术公开的方法由于采用加权插值方法，仍能够得到较为准确的频率、幅值和相角计算结果，从而确保动态功率计算的准确性。传统的加窗插值FFT算法来计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用目前广泛应用的双谱线、三谱线插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题；本专利技术公开的方法在步骤b中，巧妙地运用三根相邻的幅度最大谱线X(k0－1)、X(k0)和X(k0+1)的实部和虚部，构建两个方向的插值多项式δ+和δ－，通过引入负频率部分的参数，使频域插值结果的准确度得到明显提升。本专利技术提出的电功率计量方法数学模型简单，用简单的MATLAB语句即可完成计算，便于在工程应用中实现。也可用C语言编程本专利所提出的电功率计量方法数据处理程序，将程序下载到移动设备中，便于生产生活中动态电功率的计量。解决了现有技术采用加窗插值FFT算法计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。附图说明图1为本专利技术方法程序流程图；图2为具体实施例电压采样图形；图3为具体实施例电流采样图形；图4为具体实施例电压序列频谱图；图5为具体实施例电流序列频谱图。具体实施方式一种改进频域插值的动态功率分析方法，采集正弦电压信号、电流信号，用一种改进的频域插值方法计算两个方向的频谱插值系数，再经过加权处理得到最终的频谱插值系数，将最终的频率插值系数带入频率估计式中计算正弦信号的频率，然后计算电压、电流的幅值和相位，最后计算动态功率，从而实现对动态功率实时、准确的测量，具体包括以下步骤：步骤a，对实际电压信号u(t)和电流信号i(t)进行离散采样，采样频率为fs，采样点数为N，得到离散电压、电流序列U0(n)、I0(n)，其中n＝0,1,2,…,N-1，为便于叙述，后续步骤用待分析的正弦信号序列x(n)表示所采样的离散电压或离散电流序列；步骤b，对待分析的正弦信号序列x(n)做快速傅里叶变换，变换后的离散频谱序列X(k)，记X(k)的峰值谱线的序号k为k0，记X(k0)＝a+jb，X(k0+1)＝c+jd，X(k0－1)＝e+jf，计算频谱插值系数δ+和δ－；所述方法，其特征在于步骤b中，巧妙地运用三根相邻的幅度最大谱线：X(k0－1)、X(k0)和X(k0+1)的实部和虚部，构建两个方向的插值系数δ+和δ－，且频谱插值系数δ+和δ－的数学模型简单，计算量小只需要8次加/减法、8次乘法和2次除法。它们的表达式为：步骤c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数所述方法，其特征在于步骤c中，根据峰值谱线邻近的两根谱线幅度来加权后的频谱插值系数的表达式为：步骤d，计算待分析的正弦信号序列x(n)的频率、幅值和相位：所述方法，其特征在于步骤d中，采用加权后的频谱插值系数来计算相应的待分析的正弦信号序列x(n)的频率幅值A和相位θ，它们的表达式分别为：传统的加窗插值FFT算法来计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用目前广泛应用的双谱线、三谱线插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。步骤e，将离散电压、离散电流序列U0(n)、I0(n)分别代入步骤b至步骤d，得到电压序列的幅值和相角分别记为Au和θu，电流序列的幅值和相角分别记为Ai和θi，动态功率为P＝AiAucos(θu－θi)。具体实施例：本实施例利用电压、电流互感器引入受到5次谐波干扰的普通居民生活用电，这里电压、电流互感器就是将居民所用的电压和电流按照一定比例，转本文档来自技高网...

【技术保护点】
1.一种改进频域插值的动态功率分析方法，它包括：步骤a，对实际电压信号u(t)和电流信号i(t)进行离散采样，采样频率为fs，采样点数为N，得到离散电压U0(n)和电流序列I0(n)，n＝0,1,2,…,N‑1，采用待分析的正弦信号序列X(n)表示所采样的离散电压或离散电流序列；步骤b，对待分析的正弦信号序列X(n)做快速傅里叶变换，变换后的离散频谱序列X(k)，记X(k)的峰值谱线的序号k为k0，记X(k0)＝a+jb，X(k0+1)＝c+jd，X(k0－1)＝e+jf，计算频谱插值系数δ

【技术特征摘要】
1.一种改进频域插值的动态功率分析方法，它包括：步骤a，对实际电压信号u(t)和电流信号i(t)进行离散采样，采样频率为fs，采样点数为N，得到离散电压U0(n)和电流序列I0(n)，n＝0,1,2,…,N-1，采用待分析的正弦信号序列X(n)表示所采样的离散电压或离散电流序列；步骤b，对待分析的正弦信号序列X(n)做快速傅里叶变换，变换后的离散频谱序列X(k)，记X(k)的峰值谱线的序号k为k0，记X(k0)＝a+jb，X(k0+1)＝c+jd，X(k0－1)＝e+jf，计算频谱插值系数δ+和δ－；步骤c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数步骤d，计算待分析的正弦信号序列X(n)的频率、幅值和相位：...

【专利技术属性】
技术研发人员：张俊玮，温和，张秋雁，欧家祥，单铉昇，张博，古庭赟，牛唯，吴冕之，安甦，丁超，孙航，王扬，王蓝苓，胡厚鹏，李航锋，王泽贤，熊德智，
申请(专利权)人：贵州电网有限责任公司，
类型：发明
国别省市：贵州,52