本发明专利技术公开了一种基于滑窗频谱分离的基波参数测量方法，首先建立电力信号模型，然后通过短时傅里叶变换推导出频谱分离方法，针对谱峰分离出的多个频率成分，设定真实频率成分判定条件，最后计算并修正基波参数。本发明专利技术方法能够大大降低负频率成分和间谐波成分对基波的干扰，准确地提取基波成分，获得良好的参数测量精度；对含噪声信号的分析也表明了算法具有良好的抗噪性，适用于电网工频动态测量。
【技术实现步骤摘要】
本专利技术属于电力系统领域，具体涉及一种基于滑窗频谱分离的基波参数测量方法。
技术介绍
随着大量同步相量测量单元和频率扰动记录仪等测量设备在电网中的投入使用，电力系统的动态监控成为了电网可靠运行和控制的重要保证。其中，电力系统工频参数是电力系统稳定、高效、安全运行的重要指标，尤其是电力系统频率，对评估电网和分布电机同步运行至关重要。因此，快速准确地估计工频参数具有重要的工程实用价值。现有的基波参数估计方法，大多只单纯考虑了单独基波信号的情况，常见的有短时傅里叶变换、卡尔曼滤波、过零点检测等。但由于变频器等大量电力电子设备的使用会给电力系统信号带来间谐波干扰，给工频参数的动态监测带来了新的困难。而上述方法在短时窗下，信号存在间谐波且接近基波时，估计精度较低，其实际应用受到了限制。因此，在电力系统信号含有间谐波干扰的条件下，提供准确的基波参数估计结果，是目前基波动态监测领域亟待解决的一个重要课题，对电力系统运行控制和电能质量分析具有重要意义。
技术实现思路
本专利技术所要解决的技术问题是提供一种基于滑窗频谱分离的基波参数测量方法，利用电力信号模型和频谱分离算法，更准确地估计基波参数。为解决上述技术问题，本专利技术采用的技术方案是：一种基于滑窗频谱分离的基波参数测量方法，包括以下步骤：步骤1：采集输配电网的电压信号V或电流信号I；步骤2：电力信号建模，包括：步骤2.1：电力系统电压或电流信号的单频复指数模型为：其中，Ap，fp和分别为幅值、频率和初相角；步骤2.2：设定采样频率为fs，则采样间隔Δt＝1/fs，连续信号xp(t)离散后表示为：其中，n＝0,1,2,…,N-1，N为采样点数；步骤2.3：离散信号xp[n]的短时傅里叶变换表示为：其中，l为移动的采样点数，h[n]为矩形窗，且M为矩形窗的窗长；步骤2.4：由ω＝2πk/M，得xp[n]的离散短时傅里叶变换为：其中，Δf＝1/(NΔt)为频率分辨率，k为频域谱线位置；步骤2.5：采样频率fs远高于分析的频率成分fp，则Xp(l,kΔf)对应的求和运算近似为积分计算：其中，为滑窗移动l个采样点后的初相角，记作θp；步骤2.6：由Xp(l,k)得滑窗移动l个采样点时的频谱，则在频域k位置对应的离散傅里叶变换的值Xp(l,k)表示为：其中，βp＝fp/Δf；步骤3：频谱分离算法，包括：步骤3.1：考虑某次滑窗移动l个采样点后的频谱X(l,k)，将X(l,k)简写为X(k)，若信号的频谱由m个频率成分叠加而成，则信号在频域位置k时离散傅里叶变换的值表示为：步骤3.2：通过频谱搜索到频谱峰值，得到相应的最高和次高谱线及其位置，其中靠近零点的谱线位置由kr表示，根据X(k)表达式，结合谱峰附近的2m条谱线得方程组为：其中，D＝diag(X(k1),X(k2),…,X(k2m))，Δ1＝[δ1δ2…δm]T，Δ2＝[δm+1δm+2…δ2m]T，ki＝(kr-m+i)，δi由αj和βj构成，且i＝1,2,…,2m；j＝1,2,…,m；步骤3.3：向量Δ1和Δ2由参数α和β表示为：步骤3.4：向量Δ1和Δ2又由δ1,δ2,…,δ2m表示，根据计算参数δ1,δ2,…,δ2m的值为：步骤3.5：取信号离散傅里叶变换后的负频率部分，即取双边谱并修正X(k)的值为：步骤3.6：求得δ1,δ2,…,δ2m后，根据Δ2列出一元m次线性方程为：βm-δ2mβm-1+…+(-1)m-1δm+2β+(-1)mδm+1＝0；步骤3.7：根据一元m次线性方程构造矩阵为Ψ：得到矩阵Ψ后，通过特征值分解计算该矩阵的特征值，即一元m次方程的根β1,β2,…,βm；步骤3.8：观察向量Δ1的结构，分解为矩阵和A相乘：其中，A＝[α1α2…αm]T，步骤3.9：计算出对应β1,β2,…,βm的α1,α2,…,αm的值为：其中，是矩阵的逆矩阵；步骤4：真实频率成分判别及基波参数计算，包括：步骤4.1：得到全部α和β参数值后，计算频域谱峰处的m个频率成分自k1至k2m的离散傅里叶变换的值为：X＝[X1X2…Xj…Xm]T其中，Xj＝[Xj(k1),Xj(k2),…,Xj(k2m)]；步骤4.2：X表示某个谱峰处的频谱由m个频率成分构成，设定两个判断条件，用于提取出真实频率成分，条件一是：识别出的频率成分Xj中的谱线的幅值随着k不是单调变化的，此性质由矩形窗函数的幅频响应特性决定，条件二是Xj满足不等式：max(abs(Xj(ki)))＞μ·AF其中，AF为步骤4.3中基波幅值的计算值，参数μ由用户或噪声等级决定；步骤4.3：指定l下信号的真实频率成分的参数表示为与现有技术相比，本专利技术的有益效果是：1)建立电力信号模型，视每个频率分量为独立成分，通过频谱分离算法提取主要频率成分，避免了负频率成分、间谐波等对基波成分的影响，提高了算法的抗干扰性。2)通过识别真实频率成分，以及根据正弦或余弦信号修正计算结果，提高了在电网含有间谐波等干扰成分时的基波参数估计精度。3)通过频谱分离技术不仅减小负频率的影响，还大大降低间谐波或谐波对基波参数估计偏差，同时还具有计算效率高，分析时窗短的特点，适用于电网动态变化工况下的基波参数估计。附图说明图1为仿真的电网信号波形图。图2为本专利技术的实施例计算频率结果对比图。图3为本专利技术的实施例计算幅值结果对比图。具体实施方式下面结合附图和具体实施方式对本专利技术做进一步的详细说明。图1示出为仿真的电网电压信号，含有变化的基波成分和间谐波成分，采用的步骤为：1、信号采集采集输配电网的电压信号V或电流信号I，供离线或在线分析之用，即采用本专利技术方法既可以对实时采集的信号进行分析，也可以对某一时间段的历史信号进行动态分析。2、电力信号建模1)、电力系统电压或电流信号的单频复指数模型为：其中Ap，fp和分别为幅值、频率和初相角；2)、设定采样频率为fs，则采样间隔Δt＝1/fs，那么连续信号xp(t)离散后表示为：其中n＝0,1,2,…,N-1，N为采样点数；3)、离散信号xp[n]的短时傅里叶变换表示为：l为移动的采样点数，可以看作时间，h[n]为矩形窗，且其中M为矩形窗的窗长；4)、由ω＝2πk/M，得xp[n]的离散短时傅里叶变换为：其中，Δf＝1/(NΔt)为频率分辨率，k为频域谱线位置；5)、通常采样频率fs要远高于分析的频率成分fp，因此Xp(l,kΔf)对应的求和运算近似为积分计算：其中，为滑窗移动l个采样点后的初相角，记作θp；6)、由Xp(l,k)可得滑窗移动l个采样点时的频谱，则在频域k位置对应的离散傅里叶变换的值Xp(l,k)又表示为：其中，3、频谱分离算法1)、如果只考虑某次滑窗移动l个采样点后的频谱X(l,k)，将X(l,k)简写为X(k)，由于傅里叶变换的线性可叠加性，如果信号的频谱主要由m个频率成分叠加而成，因此信号在频域位置k时离散傅里叶变换的值表示为：2)、通过频谱可搜索到频谱峰值，也可得到相应的最高和次高谱线及其位置，其中靠近零点的谱线位置由kr表示，根据X(k)表达式，那么结合谱峰附近的2m条谱线得方程组为：其中，Η＝D·M，D＝diag(X(k1),X(k2),…,X(k2m))，Δ1＝[δ1δ2…δm]T，Δ2＝[δm+1δm+2…δ2m]本文档来自技高网...
【技术保护点】
一种基于滑窗频谱分离的基波参数测量方法，其特征在于，包括以下步骤：步骤1：采集输配电网的电压信号V或电流信号I；步骤2：电力信号建模，包括：步骤2.1：电力系统电压或电流信号的单频复指数模型为：其中，Ap，fp和分别为幅值、频率和初相角；步骤2.2：设定采样频率为fs，则采样间隔Δt＝1/fs，连续信号xp(t)离散后表示为：其中，n＝0,1,2,…,N‑1，N为采样点数；步骤2.3：离散信号xp[n]的短时傅里叶变换表示为：X(l,ω)=Σn=-∞∞xp[n+l]h[n]e-jωn]]>其中，l为移动的采样点数，h[n]为矩形窗，且M为矩形窗的窗长；步骤2.4：由ω＝2πk/M，得xp[n]的离散短时傅里叶变换为：其中，Δf＝1/(NΔt)为频率分辨率，k为频域谱线位置；步骤2.5：采样频率fs远高于分析的频率成分fp，则Xp(l,kΔf)对应的求和运算近似为积分计算：Xp(l,k)=Apejθp(ej2π(fp/Δf)-1)2πj(fp/Δf-k)]]>其中，为滑窗移动l个采样点后的初相角，记作θp；步骤2.6：由Xp(l,k)得滑窗移动l个采样点时的频谱，则在频域k位置对应的离散傅里叶变换的值Xp(l,k)表示为：Xp(l,k)=αpβp-k]]>其中，βp＝fp/Δf；步骤3：频谱分离算法，包括：步骤3.1：考虑某次滑窗移动l个采样点后的频谱X(l,k)，将X(l,k)简写为X(k)，若信号的频谱由m个频率成分叠加而成，则信号在频域位置k时离散傅里叶变换的值表示为：X(k)=X1(k)+X2(k)+...+Xm(k)=α1β1-k+α2β2-k+...+αmβm-k;]]>步骤3.2：通过频谱搜索到频谱峰值，得到相应的最高和次高谱线及其位置，其中靠近零点的谱线位置由kr表示，根据X(k)表达式，结合谱峰附近的2m条谱线得方程组为：[K-H]·Δ=X^]]>其中，Η＝D·M，D＝diag(X(k1),X(k2),…,X(k2m))，Δ1＝[δ1 δ2 … δm]T，Δ2＝[δm+1 δm+2 … δ2m]T，ki＝(kr‑m+i)，δi由αj和βj构成，且i＝1,2,…,2m；j＝1,2,…,m；步骤3.3：向量Δ1和Δ2由参数α和β表示为：Δ1=Σi=1m(αiΠj=1,j≠imβj)Σi=1m(αiΣj=1,j≠im(Πl=1,l≠i,l≠jmβl))...Σi=1mαiT,]]>Δ2=Πi=1mβiΣi=1m(Πj=1,j≠imβj)...Σi=1mβiT;]]>步骤3.4：向量Δ1和Δ2又由δ1,δ2,…,δ2m表示，根据计算参数δ1,δ2,…,δ2m的值为：Δ=K-H-1·X^;]]>步骤3.5：取信号离散傅里叶变换后的负频率部分，即取双边谱并修正X(k)的值为：X(k)=X(k)k≥0X(M+k)k<0;]]>步骤3.6：求得δ1,δ2,…,δ2m后，根据Δ2列出一元m次线性方程为：βm‑δ2mβm‑1+…+(‑1)m‑1δm+2β+(‑1)mδm+1＝0；步骤3.7：根据一元m次线性方程构造矩阵为Ψ：Ψ=δ2m-δ2m-1...-(-1)m-1δm+2-(-1)mδm+110...0001...00...............00...10]]>得到矩阵Ψ后，通过特征值分解计算该矩阵的特征值，即一元m次方程的根β1,β2,…,βm；步骤3.8：观察向量Δ1的结构，分解为矩阵和A相乘：Δ1=B^·A]]>其中，A＝[α1 α2 … αm]T，B^=Πj=1,j≠1mβjΠj=1,j≠2mβj...Πj=1,j≠mmβjΣj=1,j≠1m(Πl=1,l≠1,l≠jmβl)Σj=1,j&NotEq...
【技术特征摘要】
1.一种基于滑窗频谱分离的基波参数测量方法，其特征在于，包括以下步骤：步骤1：采集输配电网的电压信号V或电流信号I；步骤2：电力信号建模，包括：步骤2.1：电力系统电压或电流信号的单频复指数模型为：其中，Ap，fp和分别为幅值、频率和初相角；步骤2.2：设定采样频率为fs，则采样间隔Δt＝1/fs，连续信号xp(t)离散后表示为：其中，n＝0,1,2,…,N-1，N为采样点数；步骤2.3：离散信号xp[n]的短时傅里叶变换表示为：X(l,ω)=Σn=-∞∞xp[n+l]h[n]e-jωn]]>其中，l为移动的采样点数，h[n]为矩形窗，且M为矩形窗的窗长；步骤2.4：由ω＝2πk/M，得xp[n]的离散短时傅里叶变换为：其中，Δf＝1/(NΔt)为频率分辨率，k为频域谱线位置；步骤2.5：采样频率fs远高于分析的频率成分fp，则Xp(l,kΔf)对应的求和运算近似为积分计算：Xp(l,k)=Apejθp(ej2π(fp/Δf)-1)2πj(fp/Δf-k)]]>其中，为滑窗移动l个采样点后的初相角，记作θp；步骤2.6：由Xp(l,k)得滑窗移动l个采样点时的频谱，则在频域k位置对应的离散傅里叶变换的值Xp(l,k)表示为：Xp(l,k)=αpβp-k]]>其中，βp＝fp/Δf；步骤3：频谱分离算法，包括：步骤3.1：考虑某次滑窗移动l个采样点后的频谱X(l,k)，将X(l,k)简写为X(k)，若信号的频谱由m个频率成分叠加而成，则信号在频域位置k时离散傅里叶变换的值表示为：X(k)=X1(k)+X2(k)+...+Xm(k)=α1β1-k+α2β2-k+...+αmβm-k;]]>步骤3.2：通过频谱搜索到频谱峰值，得到相应的最高和次高谱线及其位置，其中靠近零点的谱线位置由kr表示，根据X(k)表达式，结合谱峰附近的2m条谱线得方程组为：[K-H]·Δ=X^]]>其中，Η＝D·M，D＝diag(X(k1),X(k2),…,X(k2m))，Δ1＝[δ1δ2…δm]T，Δ2＝[δm+1δm+2…δ2m]T，ki＝(kr-m+i)，δi由αj和βj构成，且i＝1,2,…,2m；j＝1,2,…,m；步骤3.3：向量Δ1和Δ2由参数α和β表示为：Δ1=Σi=1m(αiΠj=1,j≠imβj)Σi=1m(αiΣj=1,j≠im(Πl=1,l≠i,l≠jmβl))...Σi=1mαiT,]]>Δ2=Πi=1mβiΣi=1m(Πj=1,j&NotEqual...
【专利技术属性】
技术研发人员：臧天磊，孙仲民，何正友，钱清泉，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。