累计流量

累计流量累计流量累计流量累计流量(西门子程序西门子程序西门子程序西门子程序) FUNCTION "FLUXbase" : VOID V AR_IN_PUT PV : REAL ; //瞬时流量CYC : REAL ; //采样时间(MS) PLL : REAL ; //最小瞬时流量END_V AR V AR_IN_OUT OV_ER1 : REAL ; //累计偏差OV_ER2 : REAL ; //累计偏差PV_ALT : REAL ; //上周期瞬时流量OV : REAL ; //流量累积输出END_V AR V AR_TEMP tmpReal1 : REAL ; tmpReal2 :REAL ; OValt : REAL ; END_V AR BEGIN OV alt := OV ; IF PV >= PLL THEN tmpReal1 := PV * (CYC / 3.6E+6 ) + OV_ER1 ; tmpReal2 := tmpReal1 + OV_ER2 ; OV := tmpReal2 + OValt ; OV_ER1 := OV_ER2 - tmpReal2 + tmpReal1 ; OV_ER2 := OValt - OV + tmpReal2 ; ELSE OV := OValt ; END_IF ; IF OV >= 1.000000e+009 THEN OV := 0.000000e+000; END_IF ; END_FUNCTION 这是我的累积计算，程序简单，精度也比较高OValt = 1.0e+006 tmpReal2 = 0.1 OV = OValt + tmpReal2 = 1.0e+006 + 0.1 = 1.0e+006 OV_ER2 := OValt - OV + tmp Real2 = 1.0e+006 - 1.0e+006 + 0.1 = 0.1 qlhcco 2008-5-13 18:36:36 我把我的家底都抖给你了，这是浮点数在执行加减法时产生的误差，我把每次运算的误差再加的本次的积分里，不管再小的数，累加定量以后，他就可以加到OV里了，所以我的这个算法不会因为OV值大了而不能累积。

你可以模拟一下，有偏差变量和没有偏差变量V AR_INPUTFUNCTION "FLUXbase" : VOIDVAR_INPUTPV : REAL ; //瞬时流量CYC : REAL ; //采样时间(MS)PLL : REAL ; //最小瞬时流量END_VARVAR_IN_OUTOV_ER1 : REAL ; //累计偏差OV_ER2 : REAL ; //累计偏差PV_ALT : REAL ; //上周期瞬时流量OV : REAL ; //流量累积输出END_VARVAR_TEMPtmpReal1 : REAL ;tmpReal2 :REAL ;OValt : REAL ;END_VARBEGINOValt := OV ;IF PV >= PLL THENtmpReal1 := PV * (CYC / 3.6E+6 ) + OV_ER1 ;tmpReal2 := tmpReal1 + OV_ER2 ;OV := tmpReal2 + OValt ;OV_ER1 := OV_ER2 - tmpReal2 + tmpReal1 ;OV_ER2 := OValt - OV + tmpReal2 ;ELSEOV := OValt ;END_IF ;IF OV >= 1.000000e+009 THENOV := 0.000000e+000;END_IF ;END_FUNCTION这是我的累积计算，程序简单，精度也比较高OValt = 1.0e+006tmpReal2 = 0.1OV = OValt + tmpReal2 = 1.0e+006 + 0.1 = 1.0e+006OV_ER2 := OValt - OV + tmp Real2 = 1.0e+006 - 1.0e+006 + 0.1 = 0.1qlhcco2008-5-13 18:36:36我把我的家底都抖给你了，这是浮点数在执行加减法时产生的误差，我把每次运算的误差再加的本次的积分里，不管再小的数，累加定量以后，他就可以加到OV里了，所以我的这个算法不会因为OV值大了而不能累积。

你可以模拟一下，有偏差变量和没有偏差变量的区别。

你会发现其中的奥妙了。

yzdatou2008-5-14 7:11:04模拟了一下确实是对的。

谢谢ls指教。

你的思路很独特，受教了。

qlhcco2008-5-14 9:56:14我这个思路也不是完全是我自己的，其实是我在研究SETP7 FB41的PID功能是发现的，最后我就应用在累积计算的算法里了yzdatou2008-5-14 19:20:34很佩服你的钻研精神，自愧弗如。

另外，我思考再三，仍然要给你来个鸡蛋里面挑骨头，呵呵。

假定累积量越来越大，最终达到8位数，而实际有效数字仅7位，那么这个累计量的个位数将不再变化了。

qlhcco2008-5-15 0:07:57确实是这样的，这是由于浮点数的问题，这个问题不是不可以解决，dint型数据可以显示2147483647；可以将OV结果转换为DINT类型数据，OV_ER2也转换成DINT类型数据，OV（dint）+ OV_ER2( dint )= OUT(显示用)该输出不用于计算只是用于显示。

你认为如何？yzdatou2008-5-15 16:26:57你说的没错。

我就是用dint类型变量来存储累积量，再用一个real变量存储小于1的累积量到1就进位。

FUNCTION_BLOCK FB618//TITLE = 'LEIJI'//know_how_protectVERSION: '6.0'AUTHOR: nyjNAME: nyjFAMILY: china{S7_m_c:='true'}// Block ParametersVAR_INPUT// Input ParametersTongDaoDiZhi{S7_m_c:='true'}:WORD:=0;// AI adressXunHuanShiJian{S7_m_c:='true'}:REAL:=100;// unit:msLiangcheng{S7_m_c:='true'}:REAL:=0.0;// m3/h(max) , t/h(max)QingLing{S7_m_c:='true';S7_string_0:='false';S7_string_1:='true'}:BOOL:=0;// back to zeroEND_VARVAR_OUTPUT// Output ParametersShunShiZhi{S7_m_c:='true'}:REAL;//m3/h , t/hLeiJiZhi{S7_m_c:='true'}:DINT;//m3 , tEND_VARVAR_IN_OUTLeiJiHH:DINT:=0;//(1e--??ee)LeiJiH:DINT:=0;//(1--1e)LeiJiL:REAL:=0.0;//(0.0--1.0)END_VARBEGINIF WORD_TO_INT(TongDaoDiZhi)< 0 OR WORD_TO_INT(TongDaoDiZhi)= 32767 OR WORD_TO_INT(TongDaoDiZhi)= -32767 THENTongDaoDiZhi:=0;END_IF;IF WORD_TO_INT(TongDaoDiZhi)> 27648 AND WORD_TO_INT(TongDaoDiZhi)<> 32767 AND WORD_TO_INT(TongDaoDiZhi)<> -32767 THENTongDaoDiZhi:=27648;END_IF;ShunShiZhi:=WORD_TO_INT(TongDaoDiZhi)*Liangcheng/27648.0;LeiJiL:=LeiJiL+ShunShiZhi*XunHuanShiJian/3600000.0;IF LeiJiL >= 1.0 AND LeiJiL < 100000000.0 THENIF REAL_TO_DINT(LeiJiL)<=LeiJiL THENLeiJiH:= LeiJiH + REAL_TO_DINT(LeiJiL / 1.0);LeiJiL:= LeiJiL- REAL_TO_DINT(LeiJiL / 1.0);ELSELeiJiH:= LeiJiH + REAL_TO_DINT(LeiJiL / 1.0)-1;LeiJiL:= LeiJiL - REAL_TO_DINT(LeiJiL / 1.0)+1;END_IF;END_IF;IF LeiJiL >= 100000000.0 THENIF REAL_TO_DINT(LeiJiL/100000000.0)<=(LeiJiL/100000000.0) THENLeiJiHH:= LeiJiHH + REAL_TO_DINT(LeiJiL / 100000000.0);LeiJiL:= LeiJiL- REAL_TO_DINT(LeiJiL / 100000000.0)*100000000.0;ELSELeiJiHH:= LeiJiHH + REAL_TO_DINT(LeiJiL / 100000000.0)-1;LeiJiL:= LeiJiL - (REAL_TO_DINT(LeiJiL / 100000000.0)-1.0)*100000000.0; END_IF;END_IF;IF LeiJiH >= 100000000.0 THENLeiJiHH:= LeiJiHH+(LeiJiH DIV 100000000);LeiJiH:= LeiJiH-(LeiJiH DIV 100000000)*100000000;END_IF;IF LeiJiL < 0.0 THENLeiJiL:=0.0;END_IF;IF LeiJiH < 0 THENLeiJiH:=0;END_IF;IF LeiJiHH < 0 THENLeiJiHH:=0;END_IF;LeiJiZhi:=REAL_TO_DINT(LeiJiHH*100000000.0+LeiJiH+LeiJiL);IF QingLing = true THENQingLing := false ;LeiJiHH:= 0 ;LeiJiH:= 0 ;LeiJiL:= 0.0 ;END_IF;END_FUNCTION_BLOCK―――――――――――――――――――――――――――――――――粗粗看了下，程序先进行量程转换，得出瞬时流量，然后进行累加。