基于LabVIEWDSP模块的手持式电力谐波分析仪研究

1基于LabVIEWDSP模块的手持式电力谐波分析仪研究谢勃1，卫志农1，陈斌2(1.河海大学电气工程学院，江苏省南京市210098；2.江苏省电力科学研究院，江苏省南京市210036)摘要：本文提出一种基于全功能可视化DSP设计工具LabVIEWDSP模块的手持式电力谐波分析仪开发方案。方案以32位浮点DSP芯片为核心对准同步采样的数据进行数字滤波、加窗和FFT等算法运算。通过充分利用LabVIEWDSP模块的交互性、易用性和集成资源以及DSP芯片强大的数据处理能力，研究开发具有数据采集、实时分析、结果显示、数据保存及与其它设备通信等功能的谐波分析仪，达到缩短开发周期、降低开发难度、减小仪器体积、提高仪器性能、方便功能扩展等目的。同时，本文还指出日益发展的虚拟仪器技术，在电力系统中的将有更广阔的应用前景。关键词：谐波分析；LabVIEWDSP模块；DSP；准同步采样；FFT0引言随着电力电子技术的广泛应用，非线性负荷的容量日趋增大、数量日趋增多和工作时间不确定形成了多谐波源的运行特点[1]，使谐波问题日显突出，电能质量下降，给电网中设备的安全可靠运行带来了极大的危害，同时也给传统测量方法及常用工频仪表的使用带来了一系列的影响。因此有必要对电网谐波进行治理。治理谐波的先决条件是从电网中方便、实时和准确地对电网谐波含量进行测量与分析。过去对于不同目的的谐波测量，往往需要采用不同的测量分析仪器，有时甚至在同一项目中测量和分析需要分别采用不同的仪器。这样不但测试设备繁多、利用率低、测试成本高，而且也给测量和分析工作带来诸多不便[2]。为此，本文引入虚拟仪器技术的思想和设计方法，对具有强大数据处理能力的32位浮点DSP芯片直接进行图形化编程，建立一个主要由软件完成各项测试功能的多功能手持式电力谐波分析仪，以满足多方位、多需求的电力谐波测量和分析需要。1LabVIEW和LabVIEWDSP模块近年来，计算机技术的发展为虚拟仪器的发展和应用奠定了基础，使传统的以硬件为主体的检测仪器迅速向虚拟仪器方向发展。目前公认的代表虚拟仪器技术最高水平的是美国NationalInstrumentCorporation生产的系列产品，主要有LabVIEW、LabWindows/CVI及MeasurementStudio开发平台。LabVIEW是创建虚拟仪器系统的理想工具，它为用户提供了一个简单易用的程序开发环境，尤其适合于从事科研开发的科学家、工程技术人员，被誉为工程师和科学家的语言[3]。新版LabVIEW(如7.1和8.0等版本)及其附加模块，不仅将广泛的数据采集、分析、显示、仿真和网络通讯等功能集中在了同一个环境中，让用户可以在自己的平台上无缝地集成一套完整的应用方案，还为用户提供了包括台式操作系统、实时操作系统、便携式PDAs(个人电子助理)、FPGA(现场可编程门阵列)和DSP等广泛的应用平台[4]。LabVIEWDSP模块是一个全功能的DSP算法和系统的设计、实现与分析的工具，基于LabVIEW可视化程序环境，生成专为数字信号处理器优化的应用程序。与LabVIEW集成，功能相当强大，是DSP开发的新思路。特点包括[4]：(1)将LabVIEW的核心性能扩展至可直接利用信号处理函数(包括数百个内置分析和开发函数)和编程组件(如频谱分析和滤波)对TI的TMS320C6711和C6713DSK以及NISPEEDY-33进行编程，无需单独的DSP编译器。且不必重写代码即可在支持的DSP目标之间切换，具有多DSP目标的代码移植性。(2)全面支持ExpressVI功能，具有图形化和系统级的选择来对DSP设备进行编程，具有面向未来设计的代码模块性和独立应用开发能力等。(3)可更快地减少反复的设计周期和总体的开发时间，并且在实时地对目前的DSP硬件作出改变的同时马上显示结果。(4)具有易于访问DSP目标的模拟和数字I/O线以及满足实时交互需要的交互式GUI，可以建立实时的应用程序。2仪器硬件设计方案仪器结构如图1所示，主要由电流和电压传2感器、采样保持和模数转换单元、数据处理单元(DSP)、存储单元、显示单元及扩展接口等几个部分组成。2.1电流和电压传感器传统的电压、电流传感器常采用电磁型互感器，但其动态响应慢，传递频带窄[5]，会产生相移(0.1°～0.3°)、磁滞和饱和等现象；霍尔传感器虽然具有快速的动态响应和传递频带宽等特性，但输出的温度漂移很大、精度不高、价格较高，所以都不宜采用。经过多方面比较，本文最终选择包含信号调理电路的电流传感器SCT354和电压传感器SPT304。两者都是精密传感器，分别将电流信号和电压信号隔离变换成适当的电压信号，后面均可不处理直接接A/D采样电路，响应带宽均为20～4000Hz，且都具有精度高、相移小(不大于5’)、体积小等特点[6]。2.2采样保持及模数转换单元将模拟信号转化为数字化的系列电平信号将会引入A/D转换的量化误差e：REFb1Ve21(1)其中，VREF是输入电压满量程值，b是A/D转换器输出二进制数的位数。为减小量化误差，应尽量选用位数多的A/D转换器。因需测量较高次数的谐波，为避免出现混叠现象，保证数据采样精度，A/D转换器还要具有较高的采样频率。此外，还要考虑通道数、转换时间等指标。本文选用了TI公司提供的适合电力系统电气量测量的六路、高速、低功耗、十六位模数转换器ADS8364。其六路输入通道分成A、B、C三组，每组两个通道，分别由HOLDA、HOLDB、和HOLDC来启动转换，三组正好分别对应输入的三相电压和电流模拟量。ADS8364内部集成了六个ADC并且每个通道均具有采样保持器，能以250kHz的采样频率同时对六通道信号采样并进行转换，且六个通道的数据输出接口电压介于2.7V到5.5V之间，便于与DSP直接接口，省去了中间的电平转换[7]。2.3数据处理单元(DSP)影响仪器实时性的主要因素是FFT的数据分析计算时间，因此选择具有高速的实时信号处理能力的DSP芯片是提高仪器实时性的关键之一。本文采用了TI公司的32位浮点型DSP芯片TMS320C6711作为核心处理器件。该芯片是TMS320C6000系列中性价比较好的一款，200MHz时钟时完成1024点的FFT算法的时间只要66µs，比传统的DSP芯片要快一个数量级；该芯片采用TI公司专利的超长指令字结构，片内具有两级Cache；具有8/16/32bit的高性能外部存储器接口(EMIF)，提供了与SDRAM、SBSRAM和SRAM等同步/异步存储器的直接接口；片内集成多种外设，如多通道EDMA控制器、多通道缓冲串口(McBSP)、32bit通用计数器和可访问DSP的整个存储空间的主机口(HPI)[7]。本仪器采用TMS320C6711芯片作为主处理器，正是充分利用其强大的信号分析处理能力，以满足实时性的设计要求。3仪器软件开发软件采用模块化结构设计，不仅为开发带来方便，也为日后仪器的调试和升级带来方便。而且，利用LabVIEWDSP模块能将图形化开发环境下开发的程序直接下载至DSP硬件设备，无需单独的DSP编译器，操作非常方便。在软件的设计及实现过程中，VI的概念得以体现和应用，灵活、高效的开发特点也得以显现。3.1数据采集模块数据采集模块使用一组开发软件内置的高级模拟输入函数，采用循环缓冲区的方式，实现对信号的连续采集，实时分析。为防止缓冲区中未被读取的数据被提前覆盖掉，需协调好采样频率、缓存区大小和每次读取数据量三者间的关系。而为防止在运算中频域出现混叠现象，对连续信号的采样频率fs(每一通道独立采样频率)需满足奈奎斯特采样理论，即：s1maxf2fn(2)而工程应用中，为保证数据采样精度，一般留有一定的裕度[8]：s1maxf(4~6)fn(3)其中，f1为基波频率(50Hz)，nmax为最高谐波次数。A相电压A相电流B相电压B相电流C相电压C相电流存储器扩展接口采样保持和模数转换单元DSP电压传感器电压传感器电压传感器电流传感器电流传感器电流传感器显示单元图1仪器硬件结构框图Fig.1BlockDiagramofHardwareCircuit33.2准同步采样模块目前，谐波分析中通常采用FFT算法，得到准确结果的前提是同步采样。由于实际电网中基波频率总是波动的，而且用于测量的仪器不可能不存在误差，就目前技术即使采用复杂的硬件同步电路也很难实现严格的同步采样，且采用硬件实现开销较高，还不利于减小仪器的体积和重量。因而，本文选择了经典的准同步采样算法来实现减少误差的目标。准同步采样测量的最大特点是去掉了同步环节，利用增加每周期的采样点数和增加采样周期，较好地解决了同步误差对测量准确度的影响问题。文献[9]较详细介绍了其基本原理，这里只给出该算法的递推公式：x2kk1x1F(x)F(x)dx2(4)式中，Fk(x)是周期信号f(x)一周期内平均值的第k次递推估计结果，△是同步偏差；积分上下标中x，与Fk(x)中的x意义相同，而与Fk-1(x)及dx中x意义不同。为了便于用计算机处理，用求和代替积分，使得算法离散化。最终递推估计结果是一系列连续等间隔采样值的加权平均：wwk1iiiiwki0i0ii011FffN(5)其中，k是采样周期数，N是采样点数，ηi为权系数，w由积分方法定，采用复化梯形时w=kN，采用复化矩形时w=k(N－1)。准同步采样采样周期不要求与信号周期严格同步，即准同步采样对第一次采样点无任何要求；此方法允许信号周期有小范围的波动，即△不大，N个采样点落在区间2π+△上，但只要满足条件：N(2π+△)n/2π，其中n为高次谐波次数，该方法仍有效，且收敛速度非常快，在电网谐波测量中，正好符合网频在小范围内波动的实际情况。3.3数字滤波模块现场谐波测量过程中存在不同的干扰信号，所以除了采用良好的屏蔽措施和硬件滤波外，还必须使用数字滤波技术对采集到的电压、电流信号进行处理。数字滤波能够在程序中控制，具有高效、稳定和可预知等特性，无温度和湿度漂移，无需精度补偿，有很高的性价比，适合于实时处理[10]。LabVIEWDSP模块集成了最新的NILabVIEW数字滤波器设计工具包，使滤波器的设计更加轻松。本文采用FIR滤波器对高频带干扰信号进行滤波。设计FIR滤波器最简单的方法是窗函数法，但设计出的滤波器效率不高。另一种设计FIR滤波器的方法是使用Parks-McClellan算法[11]，Parks-McClellan算法将加权后的纹波均匀分配到通带和阻带中去，并且频率响应拥有陡峭的过渡带，此种方法设计出的滤波器较优。虚拟仪器开发工具提供了Parks-McClellan滤波器函数，在软件中可直接调用，设置好相关参数即可。但需注意，滤波器的设计部分和执行部分是分开的，因为滤波器的设计很费时间，而滤波过程则很快。一般在含有循环过程的程序中，将滤波器设计放在循环外，将设计好的滤波系数传递到循环中，在循环内只进行滤波处理，以提高程序的运行效率。Parks-McClellan.vi只计算数字FIR滤波器的系数，并不执行滤波功能，需与Convolution.vi联用滤波，即将输出系数与信号进行卷积。3.4信号加窗处理模块由于在谐波分析中使用快速傅立叶变换FFT存在频率泄漏问题，影响了谐波测量的精度。为此，在输入信号经数字滤波后通过选择适当的窗函数抑制长范围泄漏。也就是将采样信号与窗函数相乘，窗函数的幅度变化很平滑，并且在边缘处逐渐趋近于零。虚拟仪器开发工具提供了海明窗、汉宁窗、三角窗、布拉克曼窗等十多种窗函数。本文选择了汉宁窗，其减少泄露的效果比较好，计算量小，应用比较广泛。汉宁窗的定义为[11]：iiy0.5x(1cos),i0,1,,N1(6)其中，ω＝2πi/N，N为信号样本量。3.5谐波测量与分析模块谐波测量与分析是软件部分的核心。求模拟信号连续频谱的一般方法是对它作傅立叶变换。而任何一个满足狄里赫利条件的非正弦周期函数f(t)均可分解为傅立叶级数[12