基于虚拟仪器用于闪变计算的补偿迭代算法.docVIP

　　基于虚拟仪器用于闪变计算的补偿迭代算法|第1 内容加载中... 摘　要：电压波动和闪变是评定电能质量的重要指标，文章通过对IEC设计规范的闪变测试系统的分析，提出一种可准确计算瞬时闪变视感度的补偿迭代算法，从而计算出短时间闪变值和长时间闪变值，并基于虚拟仪器开发平台Labvieouseg(this) 电压波动和闪变已成为衡量电能质量的一个重要指标。IEC制定的有关标准对电压波动和闪变已作出定义并规定了限值，其中电压波动的计算较为简单，而闪变的概念较为复杂，特别是对国际上通用的衡量闪变强弱的量值Pst的计算方法未作明确说明，仅给出了电压信号的处理框图，但运用该框图无法进行任何实用的Pst数值计算。3 采用补偿迭代算法测量闪变 3.1 调幅波的检波方法 　　检测闪变需先检测出调幅波。常用的调幅波检波方法[12]有：平方检波法、整流检波法、有效值检波法。这些方法对于数字测量来说都太复杂，且精度不高。 　　本文采用补偿迭代的方法检测调幅波信号。先采用数字滤波器滤除采样数据中的高次谐波和直流分量，再检测出信号中幅值最大的信号（即基波）的频率、幅值、相角，通过信号发生器产生一反向补偿信号，抵消掉基波分量，剩余信号就是调幅波信号。 3.2 闪变值的获得 　　闪变是电网电压的波动所引起的灯光闪烁对人眼视觉产生的刺激，是根据人眼视感度曲线对电压波动信号进行相应处理后获得的，目前已有多种不同的计算方法： 　　（1）模拟加权滤波器法 　　采用一权值符合视感度曲线的加权模拟滤波器，将各频率的电压波动分量折算为8.8Hz标称频率上的波动量，从而可采用标称频率波动的大小来反映闪变的大小。该过程就是通过非线性响应来模拟灯–人眼–大脑对闪变的感觉。最后求出折算后的1min内电压波动的均方根值来表示该段时间内的闪变值。 　　（2）数字加权滤波器法 　　采用数字加权滤波器将各频率的电压波动分量按人眼视感度系数进行等效折算，其作用与模拟加权滤波器完全相同，不同之处是通过数字量的运算和处理来实现这一功能，在信号进入数字加权滤波器之前采用模/数（A/D）转换器将待处理的波动信号数字化。 　　（3）频谱分析法 　　采用A/D转换器对待测信号进行数字化采样，然后按大约为0.1Hz的频率间隔进行FFT变换，可获得0.1~30Hz之间各频率谱线的幅值，再按视感度曲线上对应各个频率分量上的视感度系数进行等效折算，从而得到所需的闪变值。 　　由于频谱分析法变换后的频谱产生了泄漏，数据不易处理，且误差较大。本文对频谱分析法进行了改进，通过补偿迭代的方法测量闪变。 3.3 补偿迭代算法测量闪变的具体实现 　　（1）采集数据的前期处理 　　为精确得到调幅波信号，必须消除直流分量、基波分量和谐波分量对采样数据的影响。先用截止频率为0.05Hz的1阶巴特沃斯高通数字滤波器滤除直流分量，再用截止频率为35Hz的6阶巴特沃斯低通数字滤波器滤除基波和高次谐波分量，而滤波过程有一从初始状态过渡到稳态的过程，舍去前面一小部分数据，只取滤波后的稳态数据进行后续的分析。由于基波频率与低通滤波器的截止频率接近，并不能完全滤除基波分量，需增加一能检测出幅值最大信号（即剩余的基波分量）的频率、幅值和相角的补偿环节，采用信号发生器产生一反向补偿信号，抵消掉基波分量，剩余信号就是调幅波信号，再采用补偿迭代算法即可得到S(t)。 　　（2）补偿迭代算法的数学描述 　　任意时刻的闪变波形在短时间内可近似视为周期性非正弦函数，采用傅立叶级数的三角级数形式展开为 　　500)this.style.ouseg(this) 式中 V(t)为闪变信号；C0为直流分量；Ci、ouseg(this) 4 算例分析 　　为验证补偿迭代算法的正确性，仿真中输入电压信号 　　500)this.style.ouseg(this) 式中 m为调幅波峰值与额定电压峰值之比；ouseg(this) 为调幅波角频率；r(t)为谐波分量和直流分量。由于电压波动d等于调幅波的峰谷差值除以载波电压额定有效值U，即 　　500)this.style.ouseg(this) 当调幅波为两个频率的叠加时，如频率8Hz、电压波动d = 0.256%的信号和频率12Hz、电压波动d = 0.312%的信号叠加时，测得的S(t)=1.00080＋1.00000=2.00080，与理论值S(t)=2相差很小，这表明本文方法对于存在多个调幅波信号的测量也是可行的。算例的波形如图3所示。仿真计算得到不同的调幅波频率f和电压波动d对应的S(t)值如表1所示。 500)this.style.ouseg(this) 500)this.style.ouseg(this)5 闪变的统计评定 　　将测量得到的S(t)从最小值到最大值进行分级