误差分布函数matlab,基于MATLAB的函数误差分布模拟与计算

基于 MATLAB的函数误差分布模拟与计算 陈 　光 　任志良 　孙海柱 (海军工程大学兵器工程系 ,武汉 430033) 摘 　要 　函数误差是间接测量误差的一种 ,本文介绍了一种用于求解函数误差分布的计算机随机模拟法。结合具体算例编制了 MATLAB 仿真程序 ,仿真结果表明该方法准确、可行 ,对复杂测量模型的误差分析尤为有效。 关键词 　函数误差 ;间接测量 ;标准差 ;随机模拟法 ;MATLAB 0 　引言 在测量活动中 ,有时因被测量对象的特点或为了追求更好的测量效果 ,而采用间接测量。间接测量误差通常由直接测量与被测量之间的函数关系来计算 ,因此又称为函数误差。研究函数误差的分布问题 ,传统的做法为 :根据概率论的知识 ,设间接测 量 x1 , x2 , ⋯, x n 的分布密度函数为 pi ( x) , i = 1 ,2 , ⋯, n。则函数 y = f ( x1 , x2 , ⋯, x n) 的概率分布函 数为 : F( y) = P( f ( x1 , x2 , ⋯, x n) ≤y) =∫f ( x 1 , x 2 , ⋯, x n ) ≤y p ( x1 , x2 , ⋯, x n ) d x1d x2 ⋯d x n (1) p ( x1 , x2 , ⋯, x n) 为 x1 , x2 , ⋯, x n 的联合概率分布密度函数。通过解析法由式(1) 即可求得函数 y 的概率分布密度函数 p ( y) = F′( y) 。这种方法在间接测量很多的情况下 ,求解是十分困难的 ,且精度也比较差。 为解决上述问题 ,本文介绍一种计算机随机模拟法 ,这种方法使用简单、计算精度高 ,不仅可以模拟出函数误差的分布 , 还可计算标准差等数字特征。 1 　计算机随机模拟法 计算机模拟是指用计算机及其软件模拟某个实际系统的部分特征或状态。对测量误差问题 ,模拟的对象为测量系统 ,所关心的问题是测量诸因素的变化对测量值的影响。计算机模拟测量系统的原理框图如图 1 所示。用随机模拟法求解函数随机误差 图 1 　计算机模拟测量系统原理图 分布的一般步骤如下 : ①输入各间接测量 x1 , x2 , ⋯, x n 及其算术平 均值x1 ,x2 , ⋯,gx n 和标准偏差σ 1 ,σ 2 , ⋯,σ n 。 ②产生所需误差分布如正态分布或均匀分布等的大样本伪随机数 ,并绘制描述各输入直接量误差分布的统计直方图。 ③按函数测量模型计算间接量 y ,并绘制该函 数误差分布的统计直方图。 ④统计并输出该间接量的最佳估计值gy ,标准差 s 与gsy 及误差分布区间半宽度Δ。 2 　算例 用两块标称长度为 50mm 的标准块规校准某 块规 ,输出两者的长度差如下 : y = f ( x1 , x2 , x3 , x4 , x5 , x6) = x1 + x2 - x1 ( x4 x5 + x3 x6) (2) 式中 , x1 为被校验的块规在 20 ℃时的校准长度 50. 000623mm ,服从正态分布 ,σ x1 = 25nm ; x2 为两标准块规在 20 ℃时的长度差 215nm , 服从正态分布 ,σ x2 = 9. 7nm ; x3 为标准块规的热膨胀系数 11. 5 ×10- 6 ℃- 1 , 　误差与数据处理 　 ·54 　 · 计量技术 20061No 5 服从均匀分布 ,σ x3 = 1. 2 ×10 - 6 ℃- 1 ; x4 为试验温度偏离标准温度 - 011 ℃,服从反正弦分布 , σ x4 = 0. 41 ℃; x5 为两块规的热