一种改进频域插值的动态功率分析方法
阅读：1发布：2020-06-11
专利汇可以提供一种改进频域插值的动态功率分析方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种改进频域插值的动态功率分析方法，它包括：a、对实际 电压 信号 和 电流 信号进行离散 采样 得到待分析的正弦信号序列X(n)；b、对待分析的正弦信号序列X(n)做快速傅里叶变换，计算 频谱 插值系数δ+和δ－；c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数d，计算待分析的正弦信号序列X(n)的 频率 、幅值和 相位 ：步骤e，得到电压序列的幅值和相 角 分别记为Au和θu，电流序列的幅值、相角和动态功率；解决了 现有技术 采用 加窗 插值FFT 算法 计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用插值算法，也存在忽略负频率部分影响，导致插值结果 精度 不高的问题。，下面是一种改进频域插值的动态功率分析方法专利的具体信息内容。
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，计算频谱插值系数δ+和δ－；k0-1和k0+1分别代表峰值谱线前后紧邻的谱线；
+ －
步骤b所述频谱插值系数δ和δ的计算公式为：
步骤c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数 步骤c所述根据峰值谱线邻近的两根谱线幅度来计算加权后的频谱插值系数 的表达式为：
步骤d，计算待分析的正弦信号序列X(n)的频率、幅值和相位：
步骤d所述计算待分析的正弦信号序列X(n)的频率 幅值A和相位θ的表达式分别为：
步骤e，将离散电压、离散电流序列U0(n)、I0(n)分别代入步骤b至步骤d的公式中，得到电压序列的幅值和相角分别记为Au和θu，电流序列的幅值和相角分别记为Ai和θi，动态功率为P＝AiAucos(θu－θi)。
一种改进频域插值的动态功率分析方法
技术领域
[0001] 本发明属于动态功率分析，尤其涉及一种改进频域插值的动态功率分析方法。
背景技术
[0002] 我国经济的不断发展导致国内电力需求旺盛，远距离、大规模输电网络的建设以及各种非线性负载的大量使用，电能质量问题日益严重。我国电力系统标准频率为50Hz，电压、电流均为50Hz的标准正弦波，但是由于受到各种干扰影响，在实际电力系统中，电压电流波形总在不同程度上产生畸变。有效的检测基波电能成为专家和学者越来越关注的问题。
[0003] 目前，电力系统信号频率分析应用最为普遍的算法是快速傅里叶变换(FFT),然而当被检测信号频率发生变换导致非同步采样时，FFT存在栅栏效应和泄露现象，使得计算的电力信号参数——频率、幅值和相位等不准确，尤其是相位误差很大，无法满足准确测量要求。为了消除电网频率变化对测量的影响，可以采用硬件锁相环电路实现同步采样；也可以对FFT算法进行改进，对于频谱泄露通常采用加窗函数的方法解决，对于栅栏效应通常采用插值算法来改善。常用的窗函数有矩形窗、三角窗、Hanning窗、Blackman窗、Blackman-Harris窗、矩形卷积窗和最大旁瓣衰减速率窗等。常用的插值算法有单谱线差值、应用最普遍的双谱线差值和近几年比较流行的三谱线差值。
[0004] 但传统的加窗插值FFT算法计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。
发明内容
[0005] 本发明要解决的技术问题是：提供一种改进频域插值的动态功率分析方法，以解决现有技术采用加窗插值FFT算法计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。
[0006] 本发明的技术方案是：
[0007] 一种改进频域插值的动态功率分析方法，它包括：
[0008] 步骤a，对实际电压信号u(t)和电流信号i(t)进行离散采样，采样频率为fs，采样点数为N，得到离散电压U0(n)和电流序列I0(n)，n＝0,1,2,…,N-1，采用待分析的正弦信号序列X(n)表示所采样的离散电压或离散电流序列；
[0009] 步骤b，对待分析的正弦信号序列X(n)做快速傅里叶变换，变换后的离散频谱序列X(k)，记X(k)的峰值谱线的序号k为k0，记X(k0)＝a+jb，X(k0+1)＝c+jd，X(k0－1)＝e+jf，计算频谱插值系数δ+和δ－；
[0010] 步骤c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数[0011] 步骤d，计算待分析的正弦信号序列X(n)的频率、幅值和相位：
[0012] 步骤e，将离散电压、离散电流序列U0(n)、I0(n)分别代入步骤b至步骤d的公式中，得到电压序列的幅值和相角分别记为Au和θu，电流序列的幅值和相角分别记为Ai和θi，动态功率为P＝AiAucos(θu－θi)。
[0013] 步骤b所述频谱插值系数δ+和δ－的计算公式为：
[0014]
[0015]
[0016] 步骤c所述根据峰值谱线邻近的两根谱线幅度来计算加权后的频谱插值系数 的表达式为：
[0017]
[0018] 步骤d所述计算待分析的正弦信号序列X(n)的频率 幅值A和相位θ的表达式分别为：
[0019]
[0020]
[0021]
[0022] 本发明有益效果：
[0023] 本发明采集正弦电压信号、电流信号，用一种改进的频域插值方法计算两个方向的频谱插值系数，再经过加权处理得到最终的频谱插值系数，将最终的频率插值系数带入频率估计式中计算正弦信号的频率，然后计算电压、电流的幅值和相位，最后计算动态功率，从而实现对动态功率实时、准确的测量。
[0024] 本发明功率计量方法所采用的改进FFT频域插值频率估计算法比传统插值算法数学模型简单，计算量小、抗干扰能力强，特别当峰值谱线邻近的两个峰值谱线收到严重干扰，产生错位时，本发明公开的方法由于采用加权插值方法，仍能够得到较为准确的频率、幅值和相角计算结果，从而确保动态功率计算的准确性。
[0025] 传统的加窗插值FFT算法来计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用目前广泛应用的双谱线、三谱线插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题；本发明公开的方法在步骤b中，巧妙地运用三根相邻的幅度最大谱线X(k0－1)、X(k0)和X(k0+1)的实部和虚部，构建两个方向的插值多项式δ+和δ－，通过引入负频率部分的参数，使频域插值结果的准确度得到明显提升。
[0026] 本发明提出的电功率计量方法数学模型简单，用简单的MATLAB语句即可完成计算，便于在工程应用中实现。也可用C语言编程本专利所提出的电功率计量方法数据处理程序，将程序下载到移动设备中，便于生产生活中动态电功率的计量。
[0028] 图1为本发明方法程序流程图；
[0029] 图2为具体实施例电压采样图形；
[0030] 图3为具体实施例电流采样图形；
[0031] 图4为具体实施例电压序列频谱图；
[0032] 图5为具体实施例电流序列频谱图。
具体实施方式
[0033] 一种改进频域插值的动态功率分析方法，采集正弦电压信号、电流信号，用一种改进的频域插值方法计算两个方向的频谱插值系数，再经过加权处理得到最终的频谱插值系数，将最终的频率插值系数带入频率估计式中计算正弦信号的频率，然后计算电压、电流的幅值和相位，最后计算动态功率，从而实现对动态功率实时、准确的测量，具体包括以下步骤：
[0034] 步骤a，对实际电压信号u(t)和电流信号i(t)进行离散采样，采样频率为fs，采样点数为N，得到离散电压、电流序列U0(n)、I0(n)，其中n＝0,1,2,…,N-1，为便于叙述，后续步骤用待分析的正弦信号序列x(n)表示所采样的离散电压或离散电流序列；
[0035] 步骤b，对待分析的正弦信号序列x(n)做快速傅里叶变换，变换后的离散频谱序列X(k)，记X(k)的峰值谱线的序号k为k0，记X(k0)＝a+jb，X(k0+1)＝c+jd，X(k0－1)＝e+jf，计算频谱插值系数δ+和δ－；
[0036] 所述方法，其特征在于步骤b中，巧妙地运用三根相邻的幅度最大谱线：X(k0－1)、X(k0)和X(k0+1)的实部和虚部，构建两个方向的插值系数δ+和δ－，且频谱插值系数δ+和δ－的数学模型简单，计算量小只需要8次加/减法、8次乘法和2次除法。它们的表达式为：
[0037]
[0038]
[0039] 步骤c，根据峰值谱线邻近的两根谱线幅度，计算加权后的频谱插值系数[0040] 所述方法，其特征在于步骤c中，根据峰值谱线邻近的两根谱线幅度来加权后的频谱插值系数 的表达式为：
[0041]
[0042] 步骤d，计算待分析的正弦信号序列x(n)的频率、幅值和相位：
[0043] 所述方法，其特征在于步骤d中，采用加权后的频谱插值系数 来计算相应的待分析的正弦信号序列x(n)的频率 幅值A和相位θ，它们的表达式分别为：
[0044]
[0045]
[0046]
[0047] 传统的加窗插值FFT算法来计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用目前广泛应用的双谱线、三谱线插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。
[0048] 步骤e，将离散电压、离散电流序列U0(n)、I0(n)分别代入步骤b至步骤d，得到电压序列的幅值和相角分别记为Au和θu，电流序列的幅值和相角分别记为Ai和θi，动态功率为P＝AiAucos(θu－θi)。
[0049] 具体实施例：
[0050] 本实施例利用电压、电流互感器引入受到5次谐波干扰的普通居民生活用电，这里电压、电流互感器就是将居民所用的电压和电流按照一定比例，转换成对应的幅值模拟信号，以适合于A/D转换的电压值。本专利中所使用MAXIM公司推出的型号为MAX1166的16位逐次逼近型A/D转换器，具有体积小、功耗低、转换速度快、精度高等优点。串行接口的特点使其与CPU直接相连，大大简化了输人通道的设计。经过A/D转换成为数字信号后，送至计算机中进行分析。
[0051] 具体实过程步骤如下：
[0052] 第一步，以fs＝3200为采样频率对受到五次谐波干扰的实际居民生活用电的电压信号u(t)和电流信号i(t)进行离散采样，一个工频周期内采样点数N＝128，得到离散电压、电流序列U0(n)、I0(n)(n＝0,2,…,127)如图2、图3所示，电压U0(n)和电流I0(n)的值如表1所示；
[0053] 表1离散电压、电流序列的值
[0054]
[0055] 第二步，对待分析的正弦信号序列U0(n)、I0(n)做快速傅里叶变换，变换后的离散频谱序列U0(k)、I0(k)(k＝0,1,2,…,127)，得到的频谱如图4、图5所示。记U0(k)的峰值谱线的序号k为ku0，在该电压频谱序列中ku0＝3，
[0056] 则U0(ku0)＝a+jb＝U0(3)＝1239.2528-19769.0300i
[0057] U0(ku0+1)＝c+jd＝U0(4)＝-36.2848+471.7011i
[0058] U0(ku0－1)＝e+jf＝U0(2)＝15.9699-261.8091i
[0059] 将a＝1239.2528、b＝-19769.0300、c＝-36.2848、d＝471.7011、e＝15.9699、f＝-261.8091带入式(1)、(2)中计算电压序列的频谱插值系数δu+和δu－：
[0060]
[0061] 记I0(k)的峰值谱线的序号k为ki0，在该电流频谱序列中ki0＝3，
[0062] 则I0(ki0)＝a+jb＝I0(3)＝28.3751-449.8248i
[0063] I0(ki0+1)＝c+jd＝I0(4)＝-0.5687+9.8857i
[0064] I0(ki0－1)＝e+jf＝I0(2)＝0.5480-6.2068i
[0065] 将a＝28.3751、b＝-449.8248、c＝-0.5687、d＝9.8857、e＝0.5480、f＝-6.2068[0066] 带入式(3)、(4)中计算电流序列的频谱插值系数δi+和δi－：
[0067]
[0068] 第三步，根据峰值谱线邻近的两根谱线幅度，分别电压序列加权后的频谱插值系数 和电流序列加权后的频谱插值系数
[0069]
[0070] 第四步，计算电压序列的频率 幅值Au和相位θu，它们的表达式分别为：
[0071]
[0072] 第五步，计算电流序列的频率 幅值Ai和相位θi，它们的表达式分别为：
[0073]
[0074] 第六步，将求得的离散电压序列U0(n)的电压幅值和电压相角，离散电流序列I0(n)的电流幅值和电流相角带入公式(13)中，计算动态功率P：
[0075]
[0076] 该地区给居民供应的标准电压的基波幅值为 电流的基波幅值为受干扰时的实际基波的频率为50.5Hz。采用本专利所公开的方法计算得到电压信号的幅值为310.2965V，与所供应标准电压幅值 的误差仅为0.2669％，计算所得电压信号的频率为50.4948Hz，与实际频率50.5Hz的误差仅为0.0103％；采用本专利所公开的方法计算得到电流信号的幅值为7.0585A，与所供应标准电流幅值 的误差仅为
0.1770％；计算所得电流信号的频率为50.4652Hz，与实际频率50.5Hz的误差仅为
0.0690％。
0.1770％；计算所得电流信号的频率为50.4652Hz，与实际频率50.5Hz的误差仅为
0.0690％。
[0077] 本发明公开的改进频域插值计算动态功率的方法准确度较高，使用该方法来计算正弦信号的频率，误差仅有万分之几。使用该方法来计算正弦信号的幅值，误差仅有千分之几。
[0078] 传统的加窗插值FFT算法来计算频率时，存在因非同步采样或数据的非整数周期截断引起的误差，即使采用目前广泛应用的双谱线、三谱线插值算法，也存在忽略负频率部分影响，导致插值结果精度不高的问题。
[0079] 本专利所公开的方法在步骤b中，巧妙地运用三根相邻的幅度最大谱线X(k0－1)、X(k0)和X(k0+1)的实部和虚部，构建两个方向的插值多项式δ+和δ－，通过引入负频率部分的参数，使频域插值结果的准确度得到明显提升。
[0080] 本发明所提出的电功率计量方法中所采用的，改进FFT插值频率估计算法比传统插值算法数学模型简单，计算量小、抗干扰能力强，特别当峰值谱线邻近的两个峰值谱线收到严重干扰，产生错位时，由于本专利所公开的方法在步骤c中运用采用加权插值方法，根据峰值谱线邻近的两根谱线幅度来加权后的频谱插值系数 仍能够得到较为准确的频率、幅值和相角计算结果，从而确保动态功率计算的准确性。
[0081] 本发明所提出的电功率计量方法数学模型简单，用简单的MATLAB语句即可完成计算，便于在工程应用中实现。也可用C语言编程本专利所提出的电功率计量方法数据处理程序，将程序下载到移动设备中，便于生产生活中动态电功率的计量。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种内容可扩展的直线位置控制教学实验系统和方法 | 2020-09-09 | 1 |
| 等离子体处理装置及其基片直流偏置电压测量方法 | 2021-03-29 | 1 |
| 一种基于决策树模型的车辆强制换道决策方法 | 2022-06-08 | 1 |
| 频率估计的一种多段信号融合方法 | 2022-11-10 | 1 |
| 一种铣削加工刀具破损监测方法 | 2021-10-18 | 2 |
| 架空电缆电流数据分析方法、健康状态监测方法及系统 | 2021-01-01 | 0 |
| 一种基于FPGA和减法电路的正弦信号相位差测量方法 | 2020-06-18 | 2 |
| 电路、接收器以及用于处理无线信号的方法 | 2023-04-30 | 1 |
| 一种数字化变电站电能质量监测系统及方法 | 2020-12-02 | 2 |
| RADIO FREQUENCY ANTENNA DEVICE FOR GENERATING A DIGITAL MAGNETIC RESONANCE INFORMATION SIGNAL | 2021-10-12 | 2 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈