基于改进s变换的电能扰动分析仪的制作方法
【专利摘要】本实用新型涉及电能扰动信号处理技术领域，是一种基于改进S变换的电能扰动分析仪，包括电源模块、模数转换器和数字信号处理器，电源模块的输出端分别与模数转换器的输入端和数字信号处理器的输入端电连接，模数转换器的输入端电连接有信号调理电路，模数转换器的输出端与数字信号处理器的输入端电连接。本实用新型的电能扰动分析仪结构简单，布局合理；本实用新型通过提取关键特征频率点，只对特征频率点进行加DPSSW窗运算，计算量少，效率高且实用性强。改进S变换所得时频矩阵结果能够直观、有效的表示各种扰动的时频特征，分类速度快，不存在误分、拒分区域，适合嵌入式实现。
【专利说明】
基于改进S变换的电能扰动分析仪
技术领域
[0001] 本实用新型设及电力系统电能扰动信号处理技术领域，是一种基于改进S变换的 电能扰动分析仪。
【背景技术】
[0002] 电能是幅值和频率均稳定的标准工频正弦电压，但在实际的输送过程中，电能受 到各种因素的影响，到达用户的电能会偏离正弦波形而发生崎变，出现电能扰动问题。近年 来，电能扰动问题对电网和配电系统造成的直接危害W及给人类生产生活造成的损失越来 越大。为保证电力系统的稳定运行，必须对其进行控制和处理。目前，国内外一些常用的电 能扰动信号检测方法主要包括：
[0003] (I)FFT及其改进算法:FFT在电能检测领域有着广泛的应用。但是在非同步采样 时，FFT算法会产生频谱泄漏和栅栏效应，此外，FFT在信号分析中还面临着时域和频域的局 部化矛盾，使其只适合分析平稳信号，无法满足对具有暂态、突变等特性的电能扰动非平稳 信号进行分析的要求。
[0004] (2)短时傅立叶变换：为解决FFT时域和频域的局部化矛盾，Dennis Gabor于1946 年引入了短时傅立叶变换(STFT)，将非平稳信号看作是一系列短时平稳信号的叠加，而短 时性则通过时间窗来获得。STFT虽然在一定程度上克服了 FFT及其改进算法不具有局部分 析能力的缺陷，但它也存在自身不可克服的缺陷：1)窗函数难W选择。对于同一信号，选择 的窗函数不同其STFT分析结果相差很远。2)窗函数确定后，无法改变信号分析的时域、频域 分辨率，若需调节时间分辨率与频率分辨率，则必须重新选择窗函数。
[0005] (3)小波变换:小波变换是傅里叶变换研究思想的发展和延拓。小波变换(WT)是一 种具有多分辨率特性的时频局部化分析方法，它克服了FFT及其改进算法的一些缺点，特别 适合于突变的非平稳信号的分析，但是它的局限性妨碍了它在电能扰动分析领域的应用： 1 )WT中的频域分辨率很粗糖，各频带间可能存在严重的频率混叠现象，远不能达至IjFFT及其 改进算法的程度。2)不同尺度的小波函数在频域上相互干扰，且易受噪声影响，不能很好地 分离频率较近的谐波和间谐波。3)利用WT在突变点的特性，可对电能扰动信号进行定性分 析，而由WT直接检测信号幅值或谐波分量存在一定困难。4)WT不同程度地存在计算量大、不 利于实时计算、难W用嵌入式系统实现等不足。
[0006] (4化i化ert-化ang变换化HT):皿T是近年来应用于非平稳信号分析的一种新方 法，它依赖于经验模态分解(EMD)，将扰动信号分解为不同频带上的固有模态分量(IMF),但 其的不足之处在于：1)在电压骤降问题中，HHT方法不能精确地在线分析，抗噪性能差。2)用 高次样条插值的EMD算法能显著提高EMD算法的精度，但所耗费的计算时间增加。3化MD并不 彻底，容易导致频带混叠，且IMF的物理意义并不明确，故实际应用仍缺乏可靠。
[0007] (5)S变换：S变换是由Stockwell等人于1996年提出，是由连续WT和STFT结合的一 种新型时频分析方法，继承和发展了STFT和WT的局部化思想，并可由运两种变换导出，克服 了STFT窗口高度和宽度固定的缺陷。因此，通过S变换不仅可W获得某一时刻的频率信息， 还可获得在某一频率上的信号幅值信息。与连续WT、STFT等时间-频率域分析方法相比，S变 换有其独特优点：1)信号S变换的分辨率与频率（即尺度)有关;2)信号的S变换结果与其傅 里叶变换保持直接的联系；3)基本小波不必满足容许性条件。因此，S变换非常适合进行电 能扰动信号分析。但是，S变换及其改进算法的时频矩阵信息量非常大，计算繁琐，难W嵌入 式实现。
[000引国内外一些常用的电能扰动信号分类方法主要包括：
[0009] (1)人工神经网络:神经网络具有简单的结构和很强的问题求解能力且可较好地 处理噪声数据是分类识别的重要方法，但是它自身有几个比较大的缺陷：算法存在局部最 优问题、收敛性较差、训练时间较长、可靠性有限。
[0010] (2)模糊技术:模糊技术的优点是通过简单明了的IF-T肥N形式的知识规则形成判 断，识别效率较高，但是同时运限制了模糊技术运用，因为许多电能扰动，例如电压波动W 及谐波，很难建立IF-T肥N运样的明显知识规则。
[OOW (3)支持向量机理论（SVMs): SVMs有效地解决了小样本高维数、非线性等问题，具 有很好的泛化能力，并克服了人工神经网络存在局部最优和训练时间较长等缺点，已广泛 应用于故障分类、时间序列预测和非线性系统建模和辨识等领域并取得了良好的效果，但 其在分类时存在不可分的现象。
【发明内容】
[0012] 本实用新型提供了基于改进S变换的电能扰动分析仪，克服了上述现有技术之不 足，其能有效解决现有技术中数据量大，计算繁琐，难W嵌入式实现W及电能扰动信号分类 算法存在的训练时间长，可靠性有限，收敛性差的问题。
[0013] 本实用新型的技术方案之一是通过W下措施来实现的：一种基于改进S变换的电 能扰动分析仪，包括电源模块、模数转换器和数字信号处理器，所述的电源模块的输出端分 别与模数转换器的输入端和数字信号处理器的输入端电连接，模数转换器的输入端电连接 有信号调理电路，模数转换器的输出端与数字信号处理器的输入端电连接，数字信号处理 器的输出端口分别电连接有复位模块、仿真调试接口、同步动态随机存储器和化ASH存储 器。
[0014] 下面是对上述发明技术方案的进一步优化或/和改进：
[0015] 上述信号调理电路输入端连接有被测电路，信号调理电路用于将被测电路信号进 行低通滤波处理，根据被测信号的时频分辨率要求，确定模拟低通滤波器的电阻电容的取 值及模数转换器的采样率和位数。
[0016] 本实用新型的电能扰动分析仪结构简单，布局合理;本实用新型通过提取关键特 征频率点，只对特征频率点进行加DPSSW窗运算，计算量少，效率高且实用性强。改进S变换 所得时频矩阵结果能够直观、有效的表示各种扰动的时频特征;其中，时域特性曲线由改进 S变换模矩阵的行向量包络线得到，它反映了被测扰动信号的幅值信息;频域特性曲线由改 进S变换的列向量包络线得到，它反映了被测扰动信号的频谱信息；改进S变换具有从复杂 崎变信号中准确提取电能扰动信号的能力;而采用DAG-SVMs分类方法，简单易行，只需要使 用n-1个决策函数即可得到分类结果，分类速度快，不存在误分、拒分区域，适合嵌入式实 现。
【附图说明】
[0017] 附图1为本发明实施例1的原理图。
[0018] 附图2为本发明实施例2的DPSSW窗函数自适应调整示意图；
[0019] 附图3为本发明实施例2方法流程图；
[0020] 附图4为本发明实施例2中改进S变换的运算示意图；
[0021 ]附图5为本发明实施例2中改进S变换流程示意图；
[0022] 附图6为本发明实例2中四个分类问题的DAG-SVMs结构示意图；
[0023] 附图7为本发明实施例2中正常电压、电压骤升、电压骤降信号时频分析图：
[0024] (a)(b)正常电压信号时域波形和其经过改进S变换后的时频=维分布；
[0025] (c)(d)电压骤升信号时域波形和其经过改进S变换后的时频=维分布；
[0026] (e)(f)电压骤降信号时域波形和其经过改进S变换后的时频=维分布；
[0027] 附图8为本发明实施例2中电压中断、谐波、电压骤降+谐波信号时频分析图：
[0028] (a)(b)电压中断信号时域波形和其经过改进S变换后的时频=维分布；
[0029] (c)(d)谐波信号时域波形和其经过改进S变换后的时频=维分布；
[0030] (e)(f)电压骤降+谐波信号时域波形和其经过改进S变换后的时频=维分布；
[0031] 附图9为本发明实施例2中电压骤降+谐波、瞬时震荡、电压闪变信号时频分析图：
[0032] (a) (b)电压骤降+谐波信号时域波形和其经过改进S变换后的时频=维分布；
[0033] (c)(d)瞬时震荡信号时域波形和其经过改进S变换后的时频=维分布；
[0034] (e)(f)电压闪变信号时域波形和其经过改进S变换后的时频=维分布。
【具体实施方式】
[0035] 本实用新型不受下述实施例的限制，可根据本实用新型的技术方案与实际情况来 确定具体的实施方式。
[0036] 在本实用新型中，为了便于描述，各部件的相对位置关系的描述均是根据说明书 附图1的布图方式来进行描述的，如:前、后、上、下、左、右等的位置关系是依据说明书附图 的布图方向来确定的。
[0037] 下面结合实施例及附图对本实用新型作进一步描述：
[0038] 实施例1:如附图1所示，一种基于改进S变换的电能扰动分析仪，其特征在于包括 电源模块、模数转换器和数字信号处理器，所述的电源模块的输出端分别与模数转换器的 输入端和数字信号处理器的输入端电连接，模数转换器的输入端电连接有信号调理电路， 模数转换器的输出端与数字信号处理器的输入端电连接，数字信号处理器的输出端口分别 电连接有复位模块、仿真调试接口、同步动态随机存储器和FLASH存储器。在实际工作过程 中，仿真调试可使用JTAG接口与数字信号处理器连接，进行仿真调试;模数转换器可W采用 TI公司生产的ADS8556进行模拟信号与数字信号的转换;数字信号处理器可W采用TI公司 生产的TMS320VC6745，其主要的参数选择为:采样频率化SPS，FFT运算数据长度N= 1024，进 行采样运算。
[0039] W上给出的模数转换器和数字信号处理器，及其参数的选取是给出了一个典型实 例，参数只要满足采样定理和实时处理的要求，就可W按W上所述方法进行电能扰动信号 的识别。
[0040] 可根据实际需要，对上述基于改进S变换的电能扰动分析作进一步优化或/和改 进：
[0041] 如附图1所示，信号调理电路输入端连接被测电路，信号调理电路用于将被测电路 信号进行低通滤波处理，根据被测信号的时频分辨率要求，确定模拟低通滤波器的电阻电 容的取值及模数转换器的采样率和位数。
[0042] 实施例2:如图2所示，基于S变换的电能扰动分析仪的使用方法，包括W下步骤：
[0043] 第一步，低通滤波:将被测电路信号输入信号调理电路进行低通滤波处理，根据被 测信号时频分辨率要求，确定模拟低通滤波器的电阻、电容取值W及模数转换器的采样率 和位数；
[0044] 第二步，模数转换:将经过滤波后的电压模拟信号经模数转换器转换为数字信号；
[0045] 第=步，改进S变换:将经过转换后的数字信号送入数字信号处理器内，在数字信 号处理器中完成构建离散长球序列窗并进行改进S变换，对改进S变换运算结果所得时频矩 阵进行分析处理，得到电能扰动信号的特征量；
[0046] 第四步，基于DAG-SVMs的电能扰动分类:将得到的电能扰动信号的特征向量输入 到训练好的DAG-SVMs进行电能扰动信号的分类识别。
[0047] 如附图2、3、6、7、8、9所示，第S步中改进S变换包括W下步骤：
[0048] 步骤1，构建离散长球序列窗函数化PSSW):被测信号x(t)经模数转换器(ADC)转换 后的离散数字序列为x(n)，n = 0，l，. . .，N-1，其采样时间间隔为Ts,对应的离散时间傅立叶 变换(DTFT)为
[0049]
[0050] 假设采样间隔为1，则其对应的频率间隔为[-1/2,1/2]，在频率a为(〇《a《l/2)的 频率范围内，其所占整个信号能量的比例可表示为
[0化1 ]
[0052] 其中，A表示信号在频域的能量聚集度。
[0053] 离散长球序列(DPSS)能够保证其构造的离散序列能够获得最大的能量聚集度，离 散长球序列窗(DPSSW)是由时间受限的离散长球序列的零阶特征向量构成，是如下方程的 解：
[0化4]
[0055] 其中n，i = 0，l，2，...，N-l，N是窗函数的长度，(6 = k/N代表在频域[0，l/2]之间所 期望最大主瓣的截止频率，k代表DPSSW的DFT变换的主瓣内最大的整数值，受香农采样定理 的约束，k必须满足k《N/2条件，紀/?表示离散长球序列，而AW则表示对应离散长球 序列的能量聚集度，其也可看作上述方程解得的每个特征向量所对应的特征值，此外，该方 程的特征向量必须满足如下条件
[0化6]
[0057] 解得方程最大特征值aW所对应的的零阶特征向量wp!紙巧，巧特征向量又可改 写为w(n，k)，其即是我们需要的DPSSW窗函数；
[0058] 步骤2,特征频率点的选取:结合电能扰动信号的特点，采用只计算关键频率点信 息、剔除无关信息的方法选取特征频率点。在运里，关键频率点的选取根据实际计算的实时 性来选取，一般选取基波到40谐波共40个频率点及运些点附近的2到4个点。
[0059] 步骤3,窗函数自适应调整：步骤1中w(n，k)为离散长球序列窗函数(DPSSW)，其宽 度可随频率自适应调整；
[0060] 步骤4，改进S变换：由步骤1、步骤2、步骤3构建的离散改进S变换可表示为
[0061]
[0062]
[0063]
[0064] 其中m = 0，l，. . .，N-l，m和n代表时间分量，k和1代表频率分量，N代表采样点数，X a+k)和W(l，k)分别代表妍马拜和w(n，k)的快速傅立叶(FFT)变换； / .、
[0065] 步骤5，FFT运算:将步骤1采样所得序列x(n)进行FFT变换得到X[l];计算步骤1中 的DPSSW窗函数w(m，k)的FFT变换W(l，k)，如果k>N/2,则《(m，k)的FFT变换为W(l，N/2);将 X[l]按照步骤2进行频率点提取，将满足条件的特定频率点X[l]平移到X[l+k];计算X[l+k] 乘 WW[l，k]得到 Y[l，k];
[0066] 步骤6，IFFT运算:将步骤5相乘得到的频谱Y[l，k]进行IFFT运算，即可得到改进S 变换的时频矩阵;改进S变换所得时频矩阵为复数矩阵，可表示为极坐标形式
[0067] SdpsswU，m)= I SdpsswU，m) I eWn,m)
[006引其中，I SDPssw(n，m) I为改进S变换模值矩阵，珍如縱）-治?織絲從知说)为改进S变 换相位矩阵；
[0069] 步骤7,时频矩阵特征向量的提取:将步骤6运算的复数时频矩阵进行分析处理得 到电能扰动自动识别所需的九个特征向量，包括时域幅值超过标准值10%的持续时间tup、 时域幅值低于标准值90%的持续时间tdDwn、时域幅值低于标准值10%的持续时间tinter、频 域特征量为最大3个波峰值的极大值Pl，P2，P3及其对应的频率f 1，f 2，f 3 ;
[0070] 如附图6至8中所示，时域图是为了说明信号时域的特点，而频域图是为了说明专 利方法提取出来电能扰动信号特征，从而更好的说明专利所述方法的有效性，也能更容易 得出电能扰动的原因，方便工作人员进一步进行处理。
[0071] 如附图5所示，第四步中所述的基于DAG-SVMs的电能扰动分类方法包括W下步骤：
[0072] 步骤1 ,DAG-SVMs参数的选取:在训练阶段，采用"一对一"的方法，即在n类训练样 本中，构造所有可能的二分分类器，每类仅在n类中的两类训练样本上训练，因此，共需构造 n(n-l)/2个SVM子分类器。在构造类i和类j的SVM子分类器时，在样本数据集选取属于类i和 类j的样本数据作为训练样本数据，并将属于类i的数据标记为正，将属于类j的数据标记为 负，即得到构建DAG-SVMs模型所需的相关参数；
[0073] 步骤2，DAG-SVMs模型的构建:在分类阶段，将所有分类器构成一个两向有向无环 图，包括n(n-l)/2个节点和n个叶，其中每个节点为一个分类器，并与下一层的两个节点相 连。
[0074] 如附图5所示，第四步的步骤2中，构造DAG-SVMs时，将分类准确率最高的分类器作 为根节点，次高的作为第二层节点，其他分类器依照W上分类方式依次进行排布。附图5给 出了一个四分类的典型实例，n分类的模型和运个类似，具体哪个出于根节点和下层节点是 根据在训练阶段，其分类准确率的大小确定。
[0075] 对于一个N点的S变换来说，其运算复杂度为0(妒），运算量极大。若充分利用基于 基2的FFT算法，则总的需要(N+l)Nlog2N次复数加法和(N+l)N/21og2N+N2次复数乘法，则N点 的S变换计算量将为O(N2IogsN)D W德州仪器(TI)型号为6745的DSP(主频375MHz)为例，执行 一次长度为1024个点的S变换分别需60ms，计算时间已经超过工频电压周期（20ms)，不利于 信号的实时快速分析。而计算包含h次谐波的改进S变换，其只需要(3h+l)Nlog2N次复数加 法和（（3hパ）N/21og2N+3hN，整个运算复杂度降为0(hNlog2N)，同样Wl024个点和 115320￥〔6745为例，其计算时间为3.51113，其远小于工频电网基波额定周期(2〇1113)便于嵌入 式系统实现和在线实时检测应用。
[0076] W上技术特征构成了本实用新型的实施例，其具有较强的适应性和实施效果，可 根据实际需要增减非必要的技术特征，来满足不同情况的需求。
【主权项】
1. 一种基于改进S变换的电能扰动分析仪，其特征在于包括电源模块、模数转换器和数 字信号处理器，所述的电源模块的输出端分别与模数转换器的输入端和数字信号处理器的 输入端电连接，模数转换器的输入端电连接有信号调理电路，模数转换器的输出端与数字 信号处理器的输入端电连接，数字信号处理器的输出端口分别电连接有复位模块、仿真调 试接口、同步动态随机存储器和FLASH存储器。2. 根据权利要求1所述的基于改进S变换的电能扰动分析仪，其特征在于所述的信号调 理电路输入端连接有被测电路，信号调理电路用于将被测电路信号进行低通滤波处理，根 据被测电路信号的时频分辨率要求，确定模拟低通滤波器的电阻和电容的取值及模数转换 器的采样率和位数。
【文档编号】G06Q50/06GK205721844SQ201620458799
【公开日】2016年11月23日
【申请日】2016年5月18日
【发明人】孙谊媊, 李宁, 张科, 王璐, 杨冰, 卫程, 王刚, 罗庆
【申请人】国网新疆电力公司电力科学研究院, 国家电网公司