计算物理学pptchp2讲述.ppt

计算物理 (Computational Physics) 2.2 平均值 方差 标准偏差 目录 计算机数和误差 实验数据的统计处理 正态分布密度函数 ?为数学期望, ?为标准误差 如?=0, ?=1 x=-5:0.01:5; y=exp(-0.5*x.^2)/sqrt(2*pi); plot(x,y) 计算物理 (Computational Physics) 2.2 平均值 方差 标准偏差 目录 计算机数和误差 实验数据的统计处理 对于物理量x的n次测量值，通常用其算术平均值作为观测量的真值的估计值。 算术平均值： 测量值对其平均值的均方偏差： 其数学期望： 故令测量值的方差为： 其数学期望为： 测量值的标准偏差： 计算物理 (Computational Physics) 2.2 平均值 方差 标准偏差 目录 计算机数和误差 实验数据的统计处理 相关MATLAB命令： 标准差:std(), 方差:var(), 最大值:max()，最小值:min(), 平均值:mean(), 求和:sum(),求积: prod(); 例： x=rand(1,100); mean(x) ans = 0.5286 std(x) ans = 0.2818 var(x) ans = 0.0794 max(x) ans = 0.9883 x=1:4 x = 1 2 3 4 sum(x) ans = 10 prod(x) ans = 24 A=[1 2 3;4 5 6] A = 1 2 3 4 5 6 sum(A) ans = 5 7 9 prod(A) ans = 4 10 18 计算物理 (Computational Physics) 2.3 错误值的剔除 目录 计算机数和误差 实验数据的统计处理 错误值：测量数据在采集、传输与记录的过程中由于很多原因可能造成数据丢失或个别数据产生不切合实际的偏差，这种数据称为错误值（或奇异项，坏值）。 拉依达方法: 当某测量值与平均值之差大于标准偏差的3倍，则剔除。 则剔除 xi 肖维勒方法: 当 其中 则剔除 xi 注：肖维勒方法和拉依达方法均属事后处理方法，不适用于实时处理过程。但他们共同的优点是简单易行。 计算物理 (Computational Physics) 2.3 错误值的剔除 目录 计算机数和误差 实验数据的统计处理 一阶差分法: 预估值： 比较判别：如果 则保留 xn否则,用 代替 xn , W 为误差窗 function y = lyd(x) % 拉依达法 sx = std(x); mx = mean(x); y = [];%空向量 for i=1:length(x) d = abs(x(i)-mx); if d = 3*sx y = [y,x(i)];%连接 end end x=randn(1,10000); hist(x,20) y=lyd(x); hist(y,20) 计算物理 (Computational Physics) 2.3 错误值的剔除 目录 计算机数和误差 实验数据的统计处理 function y = cf(x,w) % 差分法剔除错误值 % w为误差窗 n = length(x); y(1) = x(1); y(2) = x(2); for i=1:n-2 xx = 2*x(i+1)-x(i); d = abs(x(i+2)-xx); if d = w y(i+2) = xx; else y(i+2) = x(i+2); end end 例： x=rand(1,1000); hist(x,20) y=cf(x,0.1); hist(y,20) 计算物理 (Computational Physics) 思考题 目录 计算机数和误差 实验数据的统计处理 1.已知一组某物理量的测量数据， 求 x 的平均值, 方差, 标准偏差, 并用拉依达方法剔除所测数据的错误值。 2.用肖维勒方法剔除x=randn(1,10000)中的错误值，并用直方图表示。 3.用差分法剔除x=rand(1,10000)中的错误值，并用直方图表示。 肖维勒方法 function y = xwl(x,n) % 肖维勒法 sx = std