主题：【原创】［胖记原创］不确定度的计算程序。

注：在FOXBASE 或FOXPRO里都能运行，是按5 个标准平行和5 个样品平行加一个B类标准不确定度合并算的，用时直接修改里面的数据再 RUN 即可。。



编给自己用的，所以也就没搞什么界面，你们凑合使使吧。。。







** 修改以下5 个标准峰面积**

a01=5307

a02=5315

a03=5312

a04=5311

a05=5315



***以下不修改***

a0p=(a01+a02+a03+a04+a05)/5

** AOP:标准峰面积平均值**

? '标准平均值A0'+str(a0p)

UA01=(A01-A0P)*(A01-A0P)+(A02-A0P)*(A02-A0P)+(A03-A0P)*(A03-A0P)+(A04-A0P)*(A04-A0P)+(A05-A0P)*(A05-A0P)

UA0=SQRT(UA01/20)

? '标准不确度度U(A0)'+str(UA0,20,10)

? 'U(A)/A0)='+STR(UA0/A0P,20,10)

?

?

** 修改以下5 样品峰面积**

aX1=6128

aX2=6124

aX3=6119

aX4=6085

aX5=6075

***以下不修改***

aXp=(aX1+aX2+aX3+aX4+aX5)/5

** AXP:样品峰面积平均值**

? '样品平均值AX='+str(aXp)

UAX1=(AX1-AXP)*(AX1-AXP)+(AX2-AXP)*(AX2-AXP)+(AX3-AXP)*(AX3-AXP)+(AX4-AXP)*(AX4-AXP)+(AX5-AXP)*(AX5-AXP)

UAX=SQRT(UAX1/20)

? '样品不确度度U(AX)='+str(UAX,20,10)

? 'U(AX/AX)='+STR(UAX/AXP,20,10)

?

** 修改以下标准溶液的相对不确度**



UC0C0=0.0035

***以下不修改***



zh1=UC0C0*UC0C0+(UA0/A0P)*(UA0/A0P)+(UAX/AXP)*(UAX/AXP)

zh=sqrt(zh1)

? 'U='+str(zh,20,10)

? 'u*3='+str(zh*3,20,10)