基于克隆选择算法与改进快速s变换获取谐波参数的方法
【专利摘要】本发明公开了一种基于克隆选择算法与改进快速S变换获取谐波参数的方法，包括：对信号进行改进快速S变换；依据改进快速S变换的结果矩阵得到被测信号中谐波成分个数及其频率、幅值、相位参数；对改进快速S变换的结果矩阵进行线性分解；用逆快速S变换重构各谐波成分的时域波形；对检测到的各谐波成分的频率值进行修正；依据前述检测信息构建被测信号中各谐波成分的时域模型；利用克隆选择算法对检测到的谐波成分的幅值、相位信息进行校正。本发明采用克隆选择算法与改进快速S变换算法，具有检测精度高、收敛速度快以及搜索能力强等优点，可以在复杂畸变信号中提取谐波成分的各类参数。
【专利说明】基于克隆选择算法与改进快速S变换获取谐波参数的方法
【技术领域】
[0001] 本发明属于电力系统领域，涉及一种获取谐波参数的方法，具体涉及一种基于克 隆选择算法与改进快速S变换获取谐波参数的方法。
【背景技术】
[0002] 谐波是电能质量问题的主要根源，对其参数的准确估计是研究电网谐波问题的重 要内容，是实现谐波治理、潮流计算、设备状态维护的重要基础，对保证电网运行、经济安全 具有重要意义。
[0003] 目前谐波分析方法有：傅里叶变换、短时傅里叶变换、小波变换、Hilbert-Huang 变换、S变换。
[0004] 1)傅里叶变换
[0005] FFT对稳态信号能够实现有效地分析。但是FFT无法实现谐波时域、频域参数的全 面分析，无法满足对非平稳信号的检测要求。
[0006] 2)短时傅里叶变换
[0007] STFT是Gabor于1946年提出，其基本思想是在对信号加一个滑动窗的基础上进行 傅里叶变换。由于函数窗的宽度足够小的，故可以看成是将非平稳信号分解成可在时间轴 上拼接起来的一系列短时平稳信号。因此STFT得到的是可以反映信号局部频谱特征的时 频矩阵。STFT在继承了傅里叶变换优点的同时通过滑动窗实现了对时间信息的保留。
[0008] 但是STFT的滑动窗一旦确定下来就无法改变，因此它只有单一的时频分析能力。 依据海森堡测不准原理，STFT无法兼顾低频与高频信号的分析要求。因此，STFT只适合于 分析准平稳信号。
[0009] 3)小波变换
[0010] WT是Morlet于1980年在进行地震数据分析时创造的，其基本思想是用一簇函数 来表示或逼近一个信号。这一簇函数被称为小波函数系，是由母小波函数通过不同尺度的 平移与伸缩变换得到的函数组成。WT的函数窗是一个面积固定但形状可随频率变化的窗 口，因此它具有多分辨率特性，在非平稳信号检测上具有FFT与STFT所没有的优点。但是， WT结果的好坏与母小波函数的选择有关，且它存在易受噪声信号影响、变换结果的相位难 以理解以及变换尺度与频率的关系难以明确等缺点。
[0011] 4) Hilbert-Huang 变换
[0012] HHT变换是由Norden E. Huang于1998年提出，主要包括经验模态分解（Empirical Mode Decomposition，EMD)与Hilbert分析两个部分。HHT是一种时频分析方法，具有局部 适应力与分辨能力强以及易于精确分析的优点，因此其适合于平稳信号、非平稳信号的分 析。目前，HHT在谐波检测、电压暂升与暂降的时间点定位方面都有重要应用，但是由于其存 在模态混叠、端点效应、实时性不好、受噪声影响等缺点，使其目前的研究仍处于初级阶段。
[0013] 5) S 变换方法（Stockwell Transform，ST)
[0014] ST是Stockwell于1996年为了分析地震波，在STFT与WT的基础上提出来的一种 时频可逆分析方法。ST采用宽度可随着频率变化的高斯窗作为窗函数，从而保证了其具有 与WT相同的多分辨率分析的优点。但是S变换存在结果矩阵信息量大，运算时间长、耗费 的存储空间大等缺点，限制了应用场合。
【发明内容】
[0015] 本发明的目的在于克服上述现有技术的缺点，提供了一种基于克隆选择算法与改 进快速S变换获取谐波参数的方法，该方法可以准确的获取各类稳态信号及非稳态信号情 况下的谐波参数，并且能够在消耗较小空间与时间成本的情况下，快速准确的得到待测信 号的谐波参数。
[0016] 为达到上述目的，本发明所述的基于克隆选择算法与改进快速S变换获取谐波参 数的方法包括以下步骤：
[0017] 1)对待测信号进行改进快速S变换，得结果矩阵；
[0018] 2)依据步骤⑴得到的结果矩阵绘制幅值频率曲线与相位频率曲线，并根据结果 矩阵得待测信号中谐波分量的数量以及各谐波分量的频率、幅值及相位；
[0019] 3)利用改进快速S变换的线性特性采用线性分解的方法分解结果矩阵，得到各频 率分量的改进快速S变换结果矩阵；
[0020] 4)对步骤3)得到各频率分量的改进快速S变换结果矩阵根据逆变无损性分别进 行逆变改进快速S变换，重构各频率分量的时域信号；
[0021] 5)所述待测信号的FFT频率分辨率为A f = fs/N，其中，fs及N分别为时域信号 对应的采样频率及采样点数，当所述谐波分量的频率为待测信号的FFT频率分辨率△ f的 整数倍时，则步骤2)测得的待测信号中谐波分量的数量及各谐波分量的频率、幅值及相位 为最终值；
[0022] 当谐波分量的频率不为待测信号的FFT频率分辨率A f的整数倍时，则采用具有 偏移校正技术的Jacobsen频率修正方法对步骤2)得到的各谐波分量的频率进行修正，并 根据修正后的结果构建各谐波分量的时域模型，然后再采用克隆选择算法对各谐波分量的 时域模型进行修正，得待测信号中谐波分量的数量以及各谐波分量频率、幅值及相位的最 终值。
[0023] 所述改进快速S变换为
【权利要求】
1. 一种基于克隆选择算法与改进快速S变换获取谐波参数的方法，其特征在于，包括 以下步骤： 1) 对待测信号进行改进快速S变换，得结果矩阵； 2) 依据步骤（1)得到的结果矩阵绘制幅值频率曲线与相位频率曲线，并根据结果矩阵 得待测信号中谐波分量的数量以及各谐波分量的频率、幅值及相位； 3) 利用改进快速S变换的线性特性采用线性分解的方法分解结果矩阵，得到各频率分 量的改进快速S变换结果矩阵； 4) 对步骤3)得到各频率分量的改进快速S变换结果矩阵根据逆变无损性分别进行逆 变改进快速S变换，重构各频率分量的时域信号； 5) 所述待测信号的FFT频率分辨率为Λf=fs/N，其中，fs及N分别为时域信号对应 的采样频率及采样点数，当所述谐波分量的频率为待测信号的FFT频率分辨率△f的整数 倍时，则步骤2)测得的待测信号中谐波分量的数量及各谐波分量的频率、幅值及相位为最 终值； 当谐波分量的频率不为待测信号的FFT频率分辨率△f的整数倍时，则采用具有偏移 校正技术的Jacobsen频率修正方法对步骤2)得到的各谐波分量的频率进行修正，并根据 修正后的结果构建各谐波分量的时域模型，然后再采用克隆选择算法对各谐波分量的时域 模型进行修正，得待测信号中谐波分量的数量以及各谐波分量频率、幅值及相位的最终值。
2. 根据权利要求1所述的基于克隆选择算法与改进快速S变换获取谐波参数的方法， 其特征在于，所述改进快速S变换为：-S'(/./)=fm办，其中，S(t，f)为 J-OO 变换后的结果，11(1：)为待分析信号,>^-?~，/)=1'^|_￡242为广义高斯窗函数，t为时 移变量，f为频率，τ是时间，k为用户预设的参量。
3. 根据权利要求1所述的基于克隆选择算法与改进快速S变换获取谐波参数的方法， 其特征在于，步骤2)中通过幅值频率曲线中波峰的个数确定谐波分量的数量，并根据所述 波峰的数值确定各谐波分量的幅值，通过相位频率曲线上各频率点的数值确定各谐波分量 的相位。
【文档编号】G06F19/00GK104239702SQ201410446002
【公开日】2014年12月24日 申请日期:2014年9月3日 优先权日:2014年9月3日
【发明者】丁忠安, 雷万钧, 袁露, 刘博
申请人:西安交通大学