作者简介:于华楠(1981-),女,副教授,博士.研究方向:压缩传感理论,无线通信信道估计.E-mail:yhn810117@163.com
首先,选取傅立叶基作为压缩感知中的过完备字典、高斯随机矩阵作为观测矩阵、正交匹配追踪(OMP)算法作为重构算法,对电能质量扰动信号进行压缩采样。然后,采用灰度共生矩阵纹理特征中的能量特征值、灰度值出现概率两种方法对压缩感知重构信号进行分类检测。试验结果表明:本文方法可以同时实现对三相电能质量信号的压缩重构,也可以实现对信号的准确分类,并减少了信号压缩采样过程中的数据量。
First, in order to reconstruct three-phase power quality disturbance signals, the Fourier basis is chosen as the over-complete dictionary; Gaussian random matrices as the observation matrix; and orthogonal pursuit algorithm as the reconstruction algorithm. Then, the energy eigenvalue of gray level cooccurrence matrix texture features and the probability of the occurrence of grey value are used to classify the disturbance signals. The experiment results show that the proposed method can not only realize simultaneous compression and reconstruction of three-phase electric signals and accurate classification of such signals, but also reconstruct the original signal samples, reducing the data volume of signal compression and sampling.
近年来, 电能质量已成为影响用电设备安全可靠运行的重要因素之一。研究电能质量信号的压缩和分类具有重要意义。对信号压缩、检测和分类的方法主要有小波变换理论、S变换等, 它们都是以Nyquist采样定理为基础的传统处理方法, 采样数据量较大。Donoho等提出了一种新的理论— — 压缩感知(Compressed sensing, CS)理论[1, 2, 3], 它同时实现了信号采样和压缩, 突破了Nyquist采样定理的限制, 极大地减少了采样数据量, 具有广泛的应用前景[4, 5, 6, 7]。文献[8]首次将压缩感知应用于短时电能质量扰动信号压缩采样中, 研究了CS短时电能质量信号压缩采样的实现方法。文献[9]研究了应用于电能质量信号的CS稀疏基及其性能。文献[10]将CS应用于暂态短时电能质量扰动信号处理中, 并将一维时间信号转换为二维矩阵进行处理, 进一步减少了采样数据量。目前对于电能质量扰动信号的分类方法研究已取得一定进展。文献[11]提出了一种基于dq0变换和专家系统的扰动信号区分方法, 但是对信号的区分过程比较复杂。文献[12]利用小波神经网络分类器对电能质量扰动进行识别, 该方法需要进行大量的训练, 计算量非常大。文献[13]提出了电力信号短时扰动检测及自动分类方法, 分类时利用二进制序列, 不直观、不易于观测。文献[14]提出了一种基于广义S变换的暂态电能质量扰动定位及识别。
目前, 电能质量信号压缩检测方法几乎都是针对单相信号进行压缩检测。然而, 实际的电能质量监测系统往往是同时监测三相数据信息, 而且三相电信号之间存在很大的相关性。因此, 本文提出了一种同时处理三相电能质量扰动信号的方法。通过变换使三相电信号同时通过本文设计的压缩感知算法进行重构, 然后利用灰度共生矩阵纹理特征中的能量特征值、检测灰度值出现概率两种方法对三相电信号进行分类。从而同时实现对三相电能质量扰动信号的快速、准确压缩重构及分类。
本文所提出的三相电能质量短时扰动信号压缩及分类整体方案如图1所示。
首先将三相电能质量短时扰动信号(ua、ub、uc)转换为连续的一维时间信号u0, 便于在下一步同时对三相电信号进行压缩感知; 然后选择稀疏表示中的稀疏基、随机投影中的观测矩阵和重构算法, 对信号进行压缩感知; 在对信号进行压缩感知时, 为了减少信号之间的冗余, 采用了将一维时间信号转换为二维灰度图的方法, 因此, 重构的信号是以二维灰度图的形式表示的。然后, 通过以下两种方法对灰度图进行特征提取:①通过灰度共生矩阵对重构信号进行纹理特征提取:提取出各种扰动信号的纹理特征, 经过比较分析, 找到能够准确区分各种信号的特征量, 然后通过观测找到这种特征量的某种变换的值, 作为以后试验时对信号进行区分的标准; ②信号灰度图的各灰度值出现概率统计:绘制出各种扰动信号对应的灰度值出现概率图, 经过分析, 找到能够区分各种信号的规律。
当对某未知的每相至多含有一种短时扰动的三相电信号, 采用本文方法进行CS重构出二维灰度图后, 采用以上两种方法对灰度图进行特征提取, 提取时涉及到对灰度图的分块, 按照第①种分类方法将得到3个纹理特征量(ua1、ub1、uc1), 这3个量代表了三相电信号扰动时刻的特征量, 通过将这3个量与已经得到的区分各类扰动信号规律的量进行比较, 可以对信号进行快速分类。采用第②种方法可以得到三幅灰度值概率图, 根据这3幅图和得到的信号分类规律来对信号进行分类。综合比较两种方法对信号分类的准确性、复杂性等特征, 选择其中一种较好的方法。
压缩感知(CS)的核心思想是:若长度为N的信号x本身是稀疏的, 或者在某个正交基Ψ 上具有稀疏性, 那么可将原始信号x通过一个与正交基Ψ 不相关的测量矩阵
CS理论直接获取数据压缩后的表达, 省略了获取
式中:Θ =Ф Ψ 为M× N的线性测量矩阵; S中只有K个非零元素, 即S为N维K-稀疏向量; 观测量y的维数满足M=CKlogN≪N, C为与重建精度有关的过采样系数, C≥ 1。
CS的重构过程就是从测量向量y=Θ S中恢复N维信号x的过程。首先通过对式(1)的逆问题(见式(2))进行求解, 获取稀疏系数S然后将S代入
稀疏基在CS理论中的角色与传统采样理论中的Nyquist采样频率类似, 是实现压缩重构的必要条件之一。能否高精度地重构原始信号主要依赖于该信号的稀疏表示。对短时电能质量扰动信号进行压缩采样时, 需要选择合适的稀疏基对短时电能质量扰动信号进行稀疏表示。由于电能质量信号的基波是频率为50 Hz的正弦信号, 而傅立叶基也是正弦信号, 且为正交基。因此, 本文选择傅立叶基作为稀疏表示中的稀疏基。本文选择高斯随机矩阵作为测量矩阵, 为了提高测量矩阵的随机性, 将对测量矩阵进行近似QR分解及单位化处理。本文选择了运行效率和重构精度适中的正交匹配追踪(Orthogonal matching pursuit, OMP)算法作为重构算法, 再利用测量值
本文对信号进行CS处理过程中采用了文献[10]中的将一维时间信号转换为二维灰度图的方法。一维时间信号的相邻数据为2个, 只能表示信号在水平方向上的变化, 而二维灰度图像中的相邻数据为8个, 能表示信号在水平、垂直和对角线三个方向上的变化。所以, 对信号进行分析时将一维时间信号转换为二维灰度图进行分析具有分析角度灵活、直观等特点。由于二维图像比一维信号具有更大的冗余度, 因此可以通过去除图像各像素之间的相关性, 进一步提高压缩性能和重构质量。因此, 本文将时间信号转换为灰度图进行处理。
3.2.1算法实现步骤
灰度共生矩阵是通过计算图像中特定方向和特定距离的两像素之间灰度差值出现的次数, 得到图像在方向、间隔、幅度变化大小和速度的综合信息, 从而将图像的灰度值转化为纹理信息[15], 进而对图像进行分类等研究。
灰度共生矩阵是指灰度级为i的点到灰度级为j的点的概率, 其中灰度级为i的点是按照某个固定的位置关系d=(Dx, Dy)离开i点的。用Pd (i, j)(i, j=0, 1, 2, L-1)来表示灰度共生矩阵, 其中i, j为像素的灰度, L为图像共有多少个灰度级, d为两个像素的空间位置间的关系, d的不同代表了两个像素之间的距离和方向不同。灰度共生矩阵的生成方向用
式中:Pd 的每个元素代表一种灰度组合下出现的次数, 例如:Pd (0, 1)表示在图像上位置关系为d时, 两个像素的灰度为0和1的情况出现的次数, 本文中d取1, 则Pd 可记为P。
灰度共生矩阵的纹理特征有14种, 一般采用下面4个最常用的不相关的纹理特征:
式中:
本文将对三相电能质量短时扰动信号进行压缩重构和分类研究。短时电压波动的具体分类情况如表1所示。
由表1可以看出, IEEE制定的电压波动特征参数中, 0.5周波(10 ms)、30周波(0.6 s)、3 s、1 min是各类别的持续时间的分界点; 0.1、0.9、1.1、1.2、1.4、1.8 p.u.是各类别的幅值的分界点。由上文知当三相电为表1中的短时扰动时, 通过CS方法重构的信号以灰度图的形式存在, 本文利用灰度共生矩阵的4种纹理特征对电能质量扰动信号进行分类检测, 实现步骤为:
(1)通过CS方法对表1所示的电压中断、电压暂升、电压暂降及标准的电能质量信号进行重构, 得到重构的灰度图。
(2)对步骤(1)所得的4个灰度图进行灰度共生矩阵分析, 计算本文所提到的4个纹理特征。通过4个纹理特征, 对其进行分析, 找到能够区分信号的特征量。
(3)利用步骤(2)得到的特征量做进一步变换, 使之在对灰度图进行区分时直观明了。
3.2.2 纹理特征量的选取
利用CS方法对实现步骤(1)进行Matlab仿真试验。仿真的三相电信号基波频率
3.3.1 算法实现步骤
利用各种扰动的灰度图计算灰度图中每种灰度值出现的概率。利用灰度值出现概率检验扰动信号类型的步骤为:
(1)计算正弦信号、电压中断、电压暂升及电压暂降信号的灰度值的出现概率。此处, 扰动存在于整个扰动时间内。
(2)对步骤(1)所得的4个灰度值出现概率图进行分析, 找到区分4种信号的方法。
3.3.2 信号分类标准的选择
灰度图中黑色的灰度值为0, 白色的灰度值为255, 其他介于黑色和白色之间的灰度值为0~255。由图2中各图像可以通过Matlab得到对应的灰度值出现概率图, 如图4所示。
由图4可见:中断信号只在灰度值为0处出现概率不为0, 概率为1; 电压暂升信号在灰度值为0及255两处出现概率较大, 灰度值为0处出现概率与正弦信号值一样, 灰度值为255处出现概率大于正弦信号此处的出现概率; 电压暂降信号的灰度值的取值范围小于255, 在灰度值为0处出现概率与正弦信号值一样; 而正弦信号的灰度值为0~255, 且在灰度值为0处出现概率最大且只有一个最大值点。因此, 通过以上分析可知:当在灰度值为0处的出现概率大于0.5时, 存在中断扰动; 当灰度值为255处的出现概率大于0.025时, 存在电压暂升扰动; 当灰度值的取值范围小于255且在非零处有值时, 存在电压暂降扰动。所以, 通过灰度值的出现概率图可以对以上几种扰动信号进行分类。
对本文方法进行仿真试验。三相电信号的表达式为:
式中:基波频率
仿真的三相电信号为正弦信号, 奈奎斯特采样频率为6.4 kHz, 数据长度N为16384; 观测数M为12 160; x=a, b, c; 0.01 s≤ t2-t1≤ 3 s, 此处的t1=0.24 s, t2=0.64 s。
(1)当
(2)
(3)
三相电信号的Matlab仿真结果如图5所示。
构造出三相电信号后, 采用本文方法先将三相电信号转换为二维灰度图, 然后对其进行压缩重构, 原灰度图及压缩感知重构的灰度图如图6所示。
本文采用峰值信噪比PSNR来评价CS算法对信号重构的效果。图7为128× 128图像原灰度图和重构灰度图的PSNR对比。试验结果为重复20次后取平均值得到的。从图7中可知, 随着采样值的增加即压缩比的增大, PSNR呈上升趋势, 重构准确度提高。图6中原灰度图数据为16 384, 重构灰度图数据为12 160。此外, 采用本文方法当使用的采样数据比Nyquist采样数据少23%时, PSNR为35 dB, 可以很好地重构出原灰度图。
由图6可见, 其中有三块位置变化明显, 此三块灰度图为扰动发生时刻的灰度图。为了提取出扰动时刻的灰度图, 对重构的灰度图进行分块。因为已知重构的灰度图由三相电信号顺次组成, 所以, 首先将灰度图分为三块(见图8), 理论上每一块灰度图中都应该只包含一种扰动; 然后, 再根据本文的两种方法对灰度图进行扰动信号检测。当采用灰度共生矩阵检测方法时, 将图8得到的每一块灰度图均分为三小块灰度图, 并将这三小块灰度图作为一组, 计算每小块灰度图的4个角度的能量平均值, 如表3所示。
在表3中, 对每一组中的三个能量值进行分析, 找出与表2中接近的数值, 进而区分出扰动的类型及扰动存在的灰度图小块。当采用灰度值出现概率的方法对扰动信号进行分类时, 对图8得到的3幅灰度图进行Matlab仿真, 得到如图9所示的灰度值出现概率图, 然后对其进行分类。最后, 对两种方法的分类效果进行比较。
由表3可知, 图8中的(a)(b)(c)图各均分为三块后每一块在4个角度的能量的平均值。由于图8(a)的三个能量值中的第二个值是1, 且每个图只有一种扰动, 所以可以快速地判定出图8(a)为中断信号的重构结果, 即三相电信号中ua 为中断信号。依次进行判定, 可知ub 为电压暂升信号, uc为电压暂降信号。
图9中的(a)(b)(c)图对应图8中的(a)(b)(c)图, 图(d)为正弦信号的灰度值概率图。
由图9(a)可见, 灰度值为0处出现概率约为0.65, 远大于0.5, 所以图9(a)中存在电压中断信号; 图9(b)中的灰度值为255处的出现概率接近于0.06, 远大于0.025, 所以图9(b)中存在电压暂升扰动; 图9(c)中灰度值为255处的出现概率小于正弦信号的值, 且其他值处出现概率密集, 因此, 图9(c)为暂降信号。
综上可知, 两种方法都可以对扰动信号进行检测。灰度共生矩阵方法是对几个数据进行比较就可以得到结果, 比较简单方便; 灰度值出现概率的方法是对直方图之间进行观测, 这种方法的优点是当扰动时间较长时可以快速地分辨出各种扰动, 但是当扰动时间较短时各个灰度值出现概率的直方图就比较接近于正弦信号的直方图, 此时, 很难用肉眼观测出信号有没有扰动及为何种扰动, 而灰度共生矩阵的方法却可以对这种情况进行快速检测。所以, 由以上分析可知, 采用灰度共生矩阵的方法更加合理。
采用灰度共生矩阵的方法的Matlab仿真试验时间为2.1910 s; 采用对三相电信号的每一相分别进行仿真的时间为4.6890 s。由此可知:采用本文方法将三相电信号进行转换为一维时间信号再对信号进行处理的运行时间比对每一相信号分别处理的时间要短, 而且可以同时对三相电信号进行处理, 试验结果一次性得出, 还避免了重复进行试验的麻烦。所以, 本文方法是简便易行的。
采用本文灰度共生矩阵的方法对表1中的各类信号进行扰动分类试验100次, 结果如表4所示。
由表4可知, 采用灰度共生矩阵方法区分瞬时扰动信号的正确率较低, 因为当扰动时间较短时, 灰度图的能量值与正弦信号的能量值较接近, 所以容易造成区分错误; 区分暂时扰动信号时, 准确率有所提升, 因为扰动时间增加, 扰动部分在灰度图中占的面积变大, 能量值越来越接近表2中的值, 容易区分; 区分短时扰动信号时, 由于扰动时间很长, 所以扰动部分的灰度图更明显, 能量值与表2的值更为接近, 准确率更高。通过以上试验可知, 在信号扰动时长、幅度适当时, 本文方法可以快速、准确地区分出信号的扰动类型。
采用本文方法可以同时对三相电信号进行压缩, 且采样值比以Nyquist采样定理为基础的传统处理方法少20%时也能够很好地重构原灰度图。本文方法减少了信号压缩采样过程中的数据量, 突破了奈奎斯特采样限制。在对信号进行分类时, 利用灰度共生矩阵纹理特征量中的能量平均值作为区分扰动信号的标准, 可以简单、快速地对扰动信号进行准确区分。通过试验验证可知, 本文方法简单、容易实现, 是一种采样数据少、可以同时区分三相电扰动信号的有效方法。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|