一种两采样点无算法.docVIP

利用两采样点无功算法实现 张目 安徽一天电气技术有限公司 系统工程师 【摘要】 关键字 SVC、 MCR、 TSC、无功功率算法。 利用两采样点无功算法实现。使用两采样点即可计算出无功量的算法。此算法可用在SVG与MCR和TSC等无功补偿装置中应用。 一、几种无功算法公式的比较： 目前在电能质量的无功补偿算法中，采用补偿无功的算法有均方根公式算法、数字移相法、傅式算法（Fourior）等算法。 1．均方根公式算法：U=；I= ；P=； S=UI；Q=。很显然均方根算法采样点至少一个周波，算法复杂，占用资源多，效率低，实际中较少应用。 2．移相算法：Q=。即滞后90的电压与第K次电流的乘积，采样点至少一个周波，实际中很多使用这种算法。 3．傅式算法：U（K）=π；U（K）=π； I（K）=π；I（K）=π； Q=（ U（K）* I（K）- U（K）* I（K） ）； Q =。很繁琐，占用很多资源，有源滤波（APF）和无功发生器（SVG）有应用。 4．两采样点无功算法：Q=。公式简单明了，含义明显，算法简单。 通过以上列出的无功（Q）计算公式，可以看到两采样点无功算法计算简明扼要，且计算公式意思明了，电压（U）和电流（I）下坠K和K+1就是指电压和电流前后两次的采样点。 二、两采样点无功算法的推导： 假设电压和电流表达式为u=Usinωt; i=Isin(ωt-θ); 则第K次和第K+1次采样电压电流分别为u= Usinωt; i= Isin(ωt-θ); u= Usinωt= Usin(ωt+ωT) ; i= Isin(ωt-θ)= Isin(ωt+ωT-θ);然后将第K次采样电压与第K+1次采样电流相乘得： ui=[cos(θ-ωT)-cos(2ωt+ωT-θ)];将第K+1次采样电压与第K次采样电流相乘得：ui=[cos(θ+ωT)-cos(2ωt+ωT-θ)];将这两个相乘式 消去与t相关项，有： UIcosθ=（1式）; UIsinθ=（2式）; 通过以上变换，（1式）和（2式）是不是很眼熟，大家知道电网的单相有功功率是P= UIcosθ;无功功率是Q=UIsinθ.所以（1式）和（2式）可以写成： P=（3式）； Q=（4式）。 为了使计算更加简单，令ωT=π/2时，这样cosωT=0,sinωT=1，有功和无功计算公式即（3式）和（4式）就变成： P=； Q=（5式）。这样（5式）就是最终两采样点计算无功的算法公式。 三、两采样点无功算法的软件实现： 为了做到（5式）的计算方法，实际应用中只需将当前的采样点与1/4周期的采样点两点进行计算即可，这在实际中也容易实现。下面一种具体实现方法程序如下：定义一个存放采样数据的数组SampPoint[50]，这是24次取样，采样点是25个，存放2个周波的数据，采样点就是50个，所以数组长度是50。再定义两个指针一个是当前采样点指针CurrentPointPtr指向采样数组SampPiont[50]的中间段，一个是采样数组指针SampPointPtr指向采样数组。如图一所示 SampPointPtr CurrentPointPtr 1 。。。。。 。。 25 。。。。。。 。。。。。50 SampPoint[50]数组 图一 每次采样一个数据SampPointPtr和CurentPointPtr指针值增加一，采样24次以后指针值还原到原来的位置。当计算无功量（Q）时，利用当前的CurrentPointPtr指向的数组值，与CurrentPointPtr-6向前移动24/2=6个采用点指向的值进行计算。对于一个正弦周波用24 次采样，向前移动6个采样点正好是移相π/2（90）如图二所示： CurrentPointPtr-6 CurrentPointPtr 图二 下面就以一个三相电网中的A相单电源为例，采样24次，将采样的电压电流数据存入SampPoint数组中，假设先采样电压，再采样电流，以U表示电压，I表示电流，则数组存放24次采样数据如表所示： U0 I0 U1 I1 。。。。。。。 U24 I24 U0、I0表示第一次采样值，U1、I1表示第二次采样值 ….. U24、I24表示第24次采样值 按此数组存放数据格式实现的C代码如下： Unsigned int SampPoint[50],CurrentPointPtr,PointPtr;//定义采样数组长度指针 Unsign