流量稳定性计算

按照JJG164-2000中的流量稳定性检定介绍章节中的说明，装置的流量稳定性计算公式主要有两种计算公式，分别如下：
1、累计时间内流量稳定性检定
2、各累积时间之间流量稳定性检定。
此文只是针对第一种检定方法（累计时间内流量稳定性检定）

计算公式

由于JJG 164-2000中关于此点的计算公式错误，应该是如下公式：

Python2 程序

需要python2，以及指定目录下的“Data.xlsx“ 流量数据文件。
Data.xlsx 文件有3个sheet，第二个sheet和第三个sheet分别记录装置小流量数据和大流量数据。记录的数据均来源于仪表读数，数据超过60个。此程序只取自记录的前70个数据。

from openpyxl import load_workbook
def stable_flow(Data):sum=0for S_Data in Data:sum=sum+S_DataE_average=sum*1.0/len(Data)  # calc average of flowrateE_diff=[]n=len(Data) # calc num of flowratefor i in range(n):E_diff.append((Data[i]-E_average)*100.0/E_average) #E_diff is percent for each flowrate from averageprint "E_diff is %s"%E_diff,len(E_diff)R_diff=[]for j in range(n):s_sum=0for i in range(n-j):s_sum=s_sum+E_diff[i]*E_diff[i+j]R_diff.append(s_sum*1.0/(n-j)) # R_diff is the relation parameterprint "R_diff is %s"%R_diff,len(R_diff)for j_min in range(n):Rate=R_diff[j_min]*1.0/R_diff[0]   #to find the min Jprint Rateif abs(Rate)<=0.1:breakprint j_minR_sum=0for j in range(j_min+1):R_sum=R_sum+abs(R_diff[j])k=2Eq=((2.0*R_sum/n)**0.5)*kprint "Eq is %.2f %%"%Eqreturn Eqworkbook=load_workbook('''D:\\MASS\\Stability Flow\\Data.xlsx''') #use the address
sheets=workbook.get_sheet_names()
print sheets
booksheet_1=workbook.get_sheet_by_name(sheets[1]) # 第二个sheet “Rig136 small flow rate”
booksheet_2=workbook.get_sheet_by_name(sheets[2]) # 第三个sheet “Rig136 big flow rate”
booksheet=[booksheet_1,booksheet_2]
flowData=[]
smallflowData=[]
bigflowData=[]
Num=72 # need 70 pcs data
for i in range(Num-2):
#Data in column9 beginning from row 3smallflowData.append(booksheet_1.cell(row=i+3,column=9).value)
for j in range(Num-2):bigflowData.append(booksheet_2.cell(row=j+3,column=9).value)print len(smallflowData),len(bigflowData)
Eq_small=stable_flow(smallflowData)
Eq_big=stable_flow(bigflowData)
if Eq_small> Eq_big:flow_Stable=Eq_small
else:flow_Stable=Eq_big
print "Flow stable is %.2f"%flow_Stable

“Data.xlsx”文件格式如下

使用Python计算水流量装置的流量稳定性相关推荐

python ks值计算_利用Python计算KS的实例详解
在金融领域中,我们的y值和预测得到的违约概率刚好是两个分布未知的两个分布.好的信用风控模型一般从准确性.稳定性和可解释性来评估模型.sOf免费资源网一般来说.好人样本的分布同坏人样本的分布应该是有很 ...
使用OpenCV和Python计算图像的“彩色度”
使用OpenCV和Python计算图像"彩色度" 1. 效果图 2. 炫彩度量方法是什么? 3. 源代码参考你是否尝试过计算每个图像的炫彩值,并根据炫彩值对自己的图像数据集进行 ...
Python计算训练数据集（测试集）中某个分类变量阴性(阳性)标签样本的不同水平（level）或者分类值的统计个数以及比例
Python计算训练数据集(测试集)中某个分类变量阴性(阳性)标签样本的不同水平(level)或者分类值的统计个数以及比例目录
Python计算两个numpy数组的交集（Intersection）实战：两个输入数组的交集并排序、获取交集元素及其索引、如果输入数组不是一维的，它们将被展平(flatten)，然后计算交集
Python计算两个numpy数组的交集(Intersection)实战:两个输入数组的交集并排序.获取交集元素及其索引.如果输入数组不是一维的,它们将被展平(flatten),然后计算交集目录
Python使用datetime中的timedelta模块实现时间增减：python计算100天后是哪年那月那日？
Python使用datetime中的timedelta模块实现时间增减:python计算100天后是哪年那月那日? 目录
python计算特征的统计值并文本输出
python计算特征的统计值并文本输出 # 输出统计分位数 df.describe(percentiles=[0.05, 0.25, 0.5, 0.75, 0.95]) featname:A => ...
python计算特征与目标的相关性并可视化
python计算特征与目标的相关性并可视化 pandas计算相关性 # 相关性 tips.corr() Out[2]: total_bill tip size total_bill 1.000000 ...
python计算时间差
python计算时间差 # 使用datetime 和 timedelta import datetime old_time = datetime.datetime.now() print(old_ti ...
python计算两个字典的相同点(从key、value、item角度)
python计算两个字典的相同点(从key.value.item角度) 字典是由一个一个的 key-value 构成的,key 是找到数据的关键,Python 对字典的操作都是通过 key 来完成的. ...

使用Python计算水流量装置的流量稳定性