电压波动与闪变的测量方法研究彭博，等 电压波动与闪变的测量方法研究 Ｒｅｓｅａｒｃｈ ｏｎ ｔｈｅ Ｍｅａｓｕｒｉｎｇ Ｍｅｔｈｏｄ ｆｏｒ Ｖｏｌｔａｇｅ Ｆｌｕｃｔｕａｔｉｏｎ ａｎｄ Ｆｌｉｃｋｅｒ 彭 傅 周 勇 史三省 （郑州大学电气工程学院 ，河南郑州４５０００１；河南电力试验研究院 ，河南郑州４５００５２） 摘要：为了简化电压波动与闪变的测量过程，避免国际电工委员会（ＩＥＣ）推荐的测量方法中的多次滤波运算，利用ＬａｂＷｉｎｄｏｗｓ／ ＣＶＩ虚拟仪器开发平台，实现了基于快速傅里叶变换的电压波动与闪变测量。给出了测量软件的流程图和用户操作界面，并对测量结 果进行了误差分析。通过与国际电工委员会推荐的测量方法对比，证实了此测量方法是一种准确、简便、实用的离散化计算方法。 关键词：快速傅里叶变换 电压波动 闪变虚拟仪器滤波器 中图分类号：ＴＨ７０１ 文献标志码：Ａ Ａｂｓｔｒａｃｔ：Ｔｏ ｓｉｍｐｌｉｆｙ ｔｈｅ ｍｅａｓｕｒｉｎｇ ｐｒｏｃｅｓｓ ｆｏｒ ｖｏｌｔａｇｅ ｆｌｕｃｔｕａｔｉｏｎ ａｎｄ ｆｌｉｃｋｅｒ，ａｎｄ ｔｏ ａｖｏｉｄ ｍｕｌｔｉｐｌｅ ｆｉｌｔｅｒｉｎｇ ｏｐｅｒａｔｉｏｎｓ ｉｎ ｔｈｅ ｍｅｔｈｏｄ ｒｅｃｏｍｍｅｎｄｅｄ ｂｙ ＩＥＣ，ｂｙ ａｄｏｐｔｉｎｇ ＬａｂＷｉｎｄｏｗｓ／ＣＶＩ ｖｉｒｔｕａｌ ｉｎｓｔｒｕｍｅｎｔ ｄｅｖｅｌｏｐｍｅｎｔ ｐｌａｔｆｏｒｍ，ｔｈｅ ｍｅａｓｕｒｅｍｅｎｔ ｂａｓｅｄ ｏｎ ｆａｓｔ Ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ （ＦＦｒ）ｆｏｒ ｖｏｌｔａｇｅ ｆｌｕｃｔｕａｔｉｏｎ ａｎｄ ｆｌｉｃｋｅｒ ｉｓ ｉｍｐｌｅｍｅｎｔｅｄ，ａｎｄ ｔｈｅ ｅｒｒｏｒｓ ｏｆ ｔｈｅ ｍｅａｓｕｒｅｍｅｎｔ ｒｅｓｕｌｔｓ ａｒｅ ａｎａｌｙｚｅｄ．Ｔｈｅ ｆｌｏｗｃｈａｒｔ ｏｆ ｔｈｅ ｍｅａｓｕｒｉｎｇ ｓｏｆｔｗａｒｅ ａｎｄ ｔｈｅ ｕｓｅｒ ｏｐｅｒａｔｉｏｎ ｉｎｔｅｒｆａｃｅ ａｒｅ ｇｉｖｅｎ．Ｔｈｒｏｕｇｈ ｉｎｔｅｒ—ｃｏｍｐａｒｉｓｏｎ ｗｉｔｈ ｔｈｅ ｍｅａｓｕｒｉｎｇ ｍｅｔｈｏｄ ｒｅｃｏｍｍｅｎｄｅｄ ｂｙ ＩＥＣ，ｉｔ ｉｓ ｖｅｒｉｉｆｅｄ ｔｈａｔ ｔｈｉｓ ｍｅｔｈｏｄ ｉｓ ａ ｐｒｅｃｉｓｅ，ｓｉｍｐｌｅ，ａｎｄ ｐｒａｃｔｉｃａｌ ｄｉｓｃｒｅｔｅ ｃａｌｃｕｌａｔｉｏｎ ｍｅｔｈｏｄ． Ｋｅｙｗｏｒｄｓ：Ｆａｓｔ Ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｉｌｌ（ＦＦＩ’） Ｖｏｌｔａｇｅ ｆｌｕｃｔｕａｔｉｏｎ Ｆｌｉｃｋｅｒ Ｖｉｒｔｕａｌ ｉｎｓｔｒｕｍｅｎｔ Ｆｉｌｔｅｒ ０ 引言 加权滤波器完成。 ＩＥＣ推荐的电压波动与闪变测量原理框图如图１ 为解决电网电压波动与闪变的问题，国际电工委 所示。 员会（ＩＥＣ）给出了电压波动与闪变的测量方法” ， 但这种方法需要设计多个滤波器并进行多次滤波运 框１ 增益调整与输入适配 算，实现过程较复杂 。而基于快速傅里叶变换（ｆａｓｔ ● 框２ 平方解调 Ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ，ＦＦＴ）的电压波动与闪变测量方法省 ● 略了复杂的滤波器设计，简化了计算过程 。故本文 框３ Ｏ．５－－３５ Ｈｚ带通滤波器 利用此方法实现了基于Ｌａｂｗｉｎｄｏｗｓ／ＣＶＩ虚拟仪器的 视感度加权滤波器 （ 电压波动与闪变的测量，并对测量结果进行了误差分 平方 析。通过与国际电工委员会推荐的测量方法对比，证 框４ 一阶低通滤波器 实了基于快速傅里叶变换的测量方法是一种准确、简 ０删 便、实用的离散化计算方法。 框５ 闪变的统计评定 １ ＩＥＣ推荐的测量方法简介 图１ ＩＥＣ推荐的电压波动与闪变测量原理框图 在ＩＥＣ推荐的测量方法中，首先要将输入的被测 Ｆｉｇ．１ Ｔｈｅ ｐｒｉｎｃｉｐｌｅ ｏｆ ｍｅａｓｕｒｉｎｇ ｖｏｌｔａｇｅ ｆｌｕｃｔｕａｔｉｏｎ ａｎｄ 电压信号ｕ（￡）调整为适合仪器内部参照水平的电压 ｌｆｉｃｋｅｒ ｒｅｃｏｍｍｅｎｄｅｄ ｂｙ ＩＥＣ 数值，再用平方检测法从被测电压信号中调解出反映 带通滤波器分别由截止频率为０．０５ Ｈｚ的高通滤 电压波动的调幅波；然后用带通滤波器滤除直流分量 波器和截止频率为３５ Ｈｚ的低通滤波器构成。其中， 和工频及以上的频率分量，要求对２倍工频的衰减须 ０．０５ Ｈｚ一阶高通滤波器的传递函数为： 在９０ ｄＢ以上。这个衰减作用由带通滤波器和视感度 ㈠＝ （１） 式中：ｓ为拉氏变换的算子；ＯＪ ＝２，ｒ Ｘ０．０５，为截止角 修改稿收到日期：２０１１—０４～１４。 第一作者彭博（１９８６一），男，现为郑州大学电力系统及其自动化专 频率。 业在读硕士研究生；主要从事电能质量分析与控制方面的研究。 ３５ Ｈｚ六阶巴特沃斯低通滤波器的传递函数为： 《自动化仪表》第３３卷第６期２０１２年６月 ６７ 电压波动与闪变的测量方法研究彭博。等 式中：Ｐ 为第ｉ个短时间闪变值； 为Ｐ。 测量时间段 内所包含Ｐ 的个数。 ＢＷ（ｓ）＝［１＋ｂ１（ ／∞ ）＋ｂ２（ｓ／ｗ ） ＋ ｂ３（ｓ／ｗ ）。＋ｂ４（ｓ／ｗ ） ＋ ｂ５（ｓ／∞ ） ＋６６（ｓ／（ｔ， ） ］ （２） 式中： 为截止角频率，∞ ＝２，ｒｒ Ｘ ３５；ｂ 为滤波系数 （ｉ＝１，２…．，６），ｂ１＝ｂ５＝３．８６４、ｂ２＝６４＝７．４６４、ｂ ＝ ９．１４１、ｂ６＝１．０。 ２基于ＦＦＴ的测量方法 ＩＥＣ推荐的闪变测量方法计算过程较复杂，在数 字式检测装置中较难实现。基于快速傅里叶变换 （ＦＦＴ）的闪变测量方法则无需设计滤波器，计算过程 视感度加权滤波器可模拟人眼视觉系统在白炽灯 受到正弦电压波动影响下的频率响应，其ＩＥＣ推荐的 传递函数为： ）＝ × （３） 较为简捷，其计算步骤如下 】。 ①对电压信号“（ｔ）进行每半个周期内ｍ点采样， 得到采样信号ｕ（ｎ），然后利用式（８）对ｕ（ｎ）每半个周 期求出一个电压均方根值，得到一段时间内的电压均方 式中： 为比例系数；Ａ为滤波系数；ｔＯ （ｉ＝１，２，…，４）为 各典型环节传递函数对应的截止角频率。这些系数是 以５个典型环节来逼近视感度频率特性曲线，并根据均 根值序列Ｕ（Ｎ），形成离散的电压均方根值曲线，即： 几— —■一Ｖ ｍ ｎ １ （Ⅳ）＝ ／ ∑Ｍ （ｎ） 电压波动即为： ｆ Ｔ —Ｔｉ （８） 方（偏）差值最小来确定的，其中，ｋ＝１．７４８ ０２，Ａ＝２竹× ４．０５９ ８１，∞１：２－ｒｒ×９．１５４ ９４、 ２＝２１Ｔ×２．２７９ ７９、∞３＝ ２盯×１．２２５ ３５、 ＝２ ×２１．９。 ｄ＝二 ＿二 ×１００％ Ⅳ （９） 视感度加权滤波器的输出为电压波动信号 （ｔ）。 为反映光照度暂态变化的能见度，还需加入平方和积 分两个滤波功能，以模拟人脑神经对视觉反映的非线 性响应和记忆效应。 式中： 、 分别为电压均方根值的两个相邻极值； （ 一 ）为对应的电压波动的幅值； 为额定电压。 ②对序列 （Ⅳ）进行快速傅里叶变换（Ｆ丌），求 出其离散频谱。将各频率对应的频谱幅值的２倍（该 频率下正弦电压均方根值曲线波动的峰峰值）除以分 解出的直流分量幅值，即可得到该频率下的电压波动 ｄ（ｉ），将ｄ（ｉ）代入式（１０），可得该频率下的瞬时闪变 值为： ｓ ：［ ｄ］ Ｉ 积分功能由一阶低通滤波器实现，其传递函数为： ＬＰ（ｓ） （４） 由波动信号 （ｔ），可利用式（５）求得瞬时闪变 值，即： ｓ（ｆ）＝Ｌ－１（ ３．３３ｅ ”‘ ） （ｆ）＝ （５） （１０） 式中：ｄ。为在该频率上产生１个单位瞬时闪变值所需 的电压波动。该值可从文献［５］中查得。 ③因为一段波形均方根值的平方等于组成该波 形的各频率分量均方根值的平方和，所以电压均方根 值序列 （Ⅳ）对应的瞬时闪变值Ｓ就等于其频谱上各 频率对应的瞬时闪变值ｓ 之和，即： （ ） 式中： 为卷积符号； （）为反拉普拉斯变换； （ｔ）为 波动信号。 对瞬时闪变值Ｓ（￡）进行等间隔采样并分级计时， 可计算出累积概率函数（ｃｕｍｕｌａｔｉｖｅ ｐｒｏｂａｂｉｌｉｔｙ ｆｕｎｃｔｉｏｎ， ＣＰＦ），最后，可得到短时间闪变值Ｐ 和长时间闪变 值ＰｌＩ。 Ｐ 的计算公式为 ： Ｐ ＝（０．０３１ ４Ｐｏ．１＋０．０５２ ５Ｐｌ＋ Ｓ：＝∑ ＝∑［ Ｓ ：至［ 掣］］ ｚ ｌ ｌ Ｕ （（）１１） 以上３个步骤构成了瞬时闪变值的离散化计算过 程。利用式（６）和式（７）可计算短时间闪变值Ｐ 和长时 间闪变值尸ｌ 。需要说明的是，通过计算ＣＰＦ来确定 式（６）用到的５个规定值很繁琐，不利于软件编程。可采 用排序法求出这５个规定值，它们分别为该段时间内ｓ序 列中的９９．９％ ９９％ ９ｒ７％ ９０％、５ｏ％概率大值。 一。 ｓｔ一 ０．０６５ ７Ｐ３＋０．２８Ｐ１ｏ＋０．０８Ｐ５０）。（６） 式中：Ｐ０ 、Ｐ。、Ｐ，、Ｐ 。、Ｐ 。分别为１０ ｒａｉｎ内ＣＰＦ曲线上 等于０．１％、１％、３％、１０％和５０％时间的ｓ（ ）值。 闪 。、 量时间段内所包含的短 ３ 基于ＦＦＴ的测量算法的实现 ……一………… 于ＦＦＴ的电压波动与闪变测量方法， 时间闪变值Ｐ 计算获得，即： Ｐ ，一１ ￣ ｐ３ｔ ．（７） ３．１软件设计６８ ＰＲＯＣＥＳＳ ＡＵＴＯＭＡＴＩＯＮ ＩＮＳＴＲＵＭＥＮＴＡＴＩＯＮ ＶｏＬ ３３ Ｎｏ．６ Ｊｕｎｅ ２０１２ 采用美国国家仪器（ＮＩ）公司的ＬａｂＷｉｎｄｏｗｓ／ｃＶＩ虚拟 仪器开发平台 ，设计了电压波动与闪变测量软件，其 流程图如图２所示。 ｌ 初始化 Ｉ ＇’ Ｉ 数据采集 Ｉ ● Ｉ对这１２０个数从大到小排序Ｉ● 每半个周期计算 ● 一得到 个平方根值，序列，Ｎ＝５１２ ｌＩ 并计算短时闪变值 ｔ求出５个规定值 ｌｌ ● 对 呐系列进行ＦＦＴ ． Ｉ ＋＋ Ｉ ● 求一个瞬时闪变值 并存储 ＜ ｔ Ｓ—＿＋ Ｉ 计算长时闪变值Ｐｌｔ Ｉ ＜ 盏≥ ｌ否 图２测量软件流程图 Ｆｉｇ．２ Ｆｌｏｗｃｈａｒｔ ｏｆ ｔｈｅ ｍｅａｓｕｒｉｎｇ ｓｏｆｔｗａｒｅ 电压波动与闪变测量软件可以实现参数设置、信 号采集、电压均方根值计算及其频谱分析、电压波动和 闪变值测量、结果显示以及数据存盘等功能。 在对电网电压信号采样时，需要注意以下事项。 ①每半个周期内的采样点数取为ｍ＝６４。 ②Ｆ丌要求Ⅳ值必须为２的整数次幂。 ③根据加权滤波频谱特性，ＦＦＴ频谱分辨率应介 于０．１～１ Ｈｚ之间，即计算１个瞬时闪变值５对应的 时间应在１～１０ Ｓ之间。 如果每半个周期（０．０１ Ｓ）求得一个电压均方根 值，Ⅳ取５１２，则每５．１２ ｓ（０．０１ Ｓ×５１２＝５．１２ Ｓ）可计 算得到一个瞬时闪变值Ｓ，１０ ｍｉｎ内大约可计算得到 １２０个瞬时闪变值Ｓ（５．１２ Ｓ×１２０ １０ ｍｉｎ）。 （Ⅳ）的频谱分辨率为１／５．１２＝０．１９５ ３（Ｈｚ）， （Ⅳ）的频谱分布在０～５０ Ｈｚ（５１２×０．１９５ ３／２ Ｈｚ＝ ５０ Ｈｚ）范围内，满足闪变的最大觉察频率范围， ｌ￣ｌｆｏ．０５～３５ Ｈｚ。 ３．２算法验证及误差分析 为简化运算，假设电网电压信号为： Ｍ（ｔ）＝５［１＋ （ｔ）］ｃｏｓ（２７ｒ Ｘ５０×ｔ） （１２） 其中波动信号为： Ｖ（ｔ）＝０．１７６ｃｏｓ（２叮Ｔ×８．８×ｔ） （１３） 由测量软件，根据采样值绘制的电压信号波形和 电压均方根值序列的频谱分别如图３和图４所示。 由图４可知，电压频谱的直流分量为３．５３７，该值对 应于电压信号中工频分量的有效值；对应于４５点的频率 为８．８ Ｈｚ（４５ ｘ０．１９５ ３ Ｈｚ），该点的幅值约为０．６２，该值的 《自动化仪表》第３３卷第６期２０１２年６月 电压波动与闪变的测量方法研究彭博，等 ２倍对应于电压波动的幅值（２×５ ｘ０．１７６／１．４１４）；其余频 率对应的幅值均几乎为０。最后计算的电压波动为 ３４．４７５％，短时间闪变值为９７．９６９。 图３电压信号的波形 Ｆｉｇ．３ Ｗａｖｅｆｏｒｍ ｏｆｔｈｅ ｖｏｌｔａｇｅ ｓｉｇｎａｌ 喜 图４电压均方根值序列的频谱 Ｆｉｇ．４ Ｓｐｅｃｔｒｕｍ ｏｆ ｖｏｌｔａｇｅ ＲＭＳ ｓｅｑｕｅｎｃｅ 查询ＩＥＣ给出的单位瞬时闪变值Ｓ＝１时不同频率 的正弦电压波动表可知，当Ｓ＝１时，８．８ Ｈｚ正弦电压波动 为０．２５％。 当电压波动波形、频率固定时，式（６）中的５个规定值 相等，且设为Ｓ（ｔ），代人式（６）则有： Ｐ ＝／０．５０９ ６Ｓ（ ）＝０．７１４４￣－（ｔ） （１４） 将Ｓ＝１代入式（１４），可得Ｐ ＝０．７１４。由于 ３４．４７５％／９７．９６９ ０．２５％／０．７１４，符合电压波动波 形、频率固定时电压波动ｄ与短时间闪变值Ｐ 成正比 的规律，从而验证了该软件的计算结果基本上是准确 的。计算结果之所以会有一些误差，原因分析如下。 ①采集误差 对于时域离散系统，有限的采样点数会导致还原 出来的波形失真、频谱特性与原始的时域连续信号不 同。因此，采样本身存在采集误差。 ②频谱混叠 当采样频率 低于原始信号最高频率 的２倍 时，原始信号中高于ｆＪ２的频谱分量会在低于 ／２的 频谱中再现，即出现频谱的混叠。这会使频谱分析出 现误差。 ③频谱泄漏 将无限长非周期离散的均方根值序列进行截断， 限定为有限长的信号进行分析，会引起频谱的泄漏，从 而导致频谱分析出现误差。 （下转第７２页） ６９ 嵌入式系统雷达电子干扰模拟器的设计２．１ 系统软件总体结构 邹向阳，等 系统软件采用多任务方式实现控制，其总体结构 如图４所示。 样式任务ｌ干扰ｌ Ｉ干扰ｌＩ干扰ｌｌ干扰ｌ Ｉ Ｉ干扰 频率任务ｌＩ强度任务ｌｌ带宽任务ｆＩ速度任务 ３结束语 雷达电子干扰模拟器采用ＡＲＭ７系列高性能微处 理器ＬＰＣ２１３８和嵌入式操作系统ｂ￣Ｃ／ＯＳ．ＩＩ，具有实时 性好、处理多任务能力强、人机对话更人性化等特点， 适合在复杂恶劣的现场环境中使用。该模拟器自推广 到部队后，始终担负着模拟逼真电子干扰环境的任务， 满足了防空兵部队电子对抗训练的实际需要，解决了 茬妻Ｈ 。ｓ－Ⅱ嵌入式操作系统 ＡＲＭ７位ＰＣ２１３８）硬件系统 显示 任务 某型雷达缺少干扰训练器材的难题。 图４ 系统软件总体结构图 Ｆｉｇ．４ Ｔｈｅ ｏｖｅｒａｌｌ ｓｏｆｔｗａｒｅ ｓｔｒｕｃｔｕｒｅ 参考文献 ［１］邹向阳．雷达干扰模拟及训练仿真［Ｍ］．桂林：桂林空军学院出 版社，２００６． ２．２软件功能设计 ［２］杨震伦，熊茂华．嵌入式操作系统及编程［Ｍ］．北京：清华大学 出版社，２００９． 系统软件采用多进程的方式实现多任务功能，其 功能结构如图５所示。其中，键盘输入管理进程完成 ［３］周立功．ＡＲＭ嵌入式系统基础教程［Ｍ］．北京：北京航空航天 大学出版社，２００６． 对键盘按键的监视和按键信息的接收；系统状态管理 进程根据按键信息完成不同系统状态间的转换和控 制；显示输出管理进程将显示堆栈中的显示信息输出 到ＬＣＤ液晶显示屏；干扰输出进程负责将干扰堆栈中 ［４］周立功．深入浅出ＡＲＭ７－－ＬＰＣ２１３ｘ／２１４ｘ［Ｍ］．北京：北京航空 航天大学出版社，２００５． ［５］ＭＡＸＩＭ Ｃｏ．Ｌｔｄ．Ｈｉｇｈ—ｆｒｅｑｕｅｎｃｙ ｗａｖｅｆｏｒｍ ｇｅｎｅｒａｔｏｒ［ＥＢ／ＯＬ］． ［２０１１一Ｏ４—２０］．ｈｔｔｐ：／／ｗｗｗ．ｍａｘｉｍ—ｉｃ．ｃｏｎ． ［６］韩叶祥，朱兆优，张斌．基于ＭＡＸ０３８的程控函数信号发生器的 设计［Ｊ］．电子元器件应用，２００８，１０（５）：２２—２５． ［７］胡卫华．串行Ｄ／Ａ转换器ＭＡＸ５３８接口技术及应用［Ｊ］．工业 仪表与自动化装置，２００２（３）：４６—４８． 的干扰控制命令输出到干扰设备。 键盘输入管理进程Ｉ Ｉ干扰输出管理进程 按键堆栈接口 ｌ ｌ 干扰堆栈接口 Ｉ ｌ显示堆栈 按键堆栈Ｉ Ｊ Ｉ干扰堆栈ｌ ｌ Ｉ显示堆栈接口 ［８］温鑫．雷达噪声干扰模拟器的设计［Ｄ］．哈尔滨：哈尔滨工程大 学，２００７． 系统状态管理进程 ｌｌ显示输出管理进程 ［９］任哲．嵌入式实时操作系统ＩｘＣ／ＯＳ－ＩＩ原理及应用［Ｍ］．北京： 北京航空航天大学出版社，２００５：１—１５． ［１０］王忠凯，赵磊．ＩｘＣ／ＯＳ—ＩＩ任务调度研究［Ｊ］．山东理工大学学 报：自然科学版，２００９，２３（２）：３０—３５． 图５ 系统软件功能结构图 Ｆｉｇ．５ Ｔｈｅ ｆｕｎｃｔｉｏｎａｌ ｓｔｒｕｃｔｕｒｅ ｏｆ ｓｙｓｔｅｍ ｓｏｆｔｗａｒｅ （上接第６９页） 证了其测量结果的准确性，从而为进一步研制全数字式 闪变仪奠定了基础。 参考文献 ［１］肖湘宁．电能质量分析与控制［Ｍ］．北京：中国电力出版社， ２００４：８５—１２１． ④计算误差 方面，在求取式（６）中用到的５个规定值（Ｐ 、 一Ｐ。、Ｐ，、Ｐ Ｐ 。）时要用到插值拟合的方法 ，这必然 引入一些误差。 另一方面，在闪变计算过程中，电压波动与瞬时闪 变值的关系均视为理想的线性关系，然而实际情况并 非如此，这样也会产生误差。 ［２］高师湃，李群湛，贺建闽．基于ＩＥＣ标准的闪变测试系统研究［Ｊ］．电 力自动化设备，２００３，２３（７）：２７—３０． ［３］袁世英，曹晖．基于ＤＳＰ的电压闪变仪的开发［Ｊ］．继电器， ２００５，３３（４）：５７—７Ｏ． ４结束语 电压波动与闪变测量方法的实质在于提取出电压 ［４］赵刚，施围，林海雪．闪变值计算方法的研究［Ｊ］．电网技术， ２００１，２５（１１）：１５—１８． ［５］ＩＥＣ．ＩＥＣ ６０ ８６８－０－１９９１ Ｆｌｉｃｋｅｒｍｅｔｅｒ ｐａｒｔ ０：Ｅｖａｌｕａｔｉｏｎ ｏｆ ｆｌｉｃｋｅｒ ｓｅｖｅｒｉｔｙ［Ｓ］．１９９１． 波动分量，然后对其进行视感度加权处理，得到瞬时闪 变值Ｓ（ｔ）曲线，再对该曲线进行统计，计算出短时间 闪变值Ｐ 和长时间闪变值Ｐ 本文采用基于ＦＦＴｒ的闪变测量方法，在ＬａｂＷｉｎｄｏｗｓ／ ＣＶＩ虚拟仪器环境下实现了电压波动与闪变的测量，并验 ７２ ［６］孙晓云，郭立炜，孙会琴．基于ＬａｂＷｉｎｄ０ｗ ｃＶ１的虚拟仪器设 计与应用［Ｍ］．北京：电子工业出版社，２００５． ［７］赵海翔，陈默子，戴慧珠．闪变值离散化计算的误差分析及其应 用［Ｊ］．电网技术，２００４，２８（１３）：８４—８７． ＰＲＯＣＥＳＳ ＡＵＴＯＭＡＴＩＯＮ ＩＮＳＴＲＵＭＥＮＴＡＴＩＯＮ Ｖｏ１．３３ Ｎｏ．６ Ｊｕｎｅ ２０１２