1、误差理论与数据处理实验报告误差理论与数据处理实 验 报 告仪器与电子学院23杨松实验一 熟悉 MATLAB软件在误差处理中的应用(验证型)1、实验数据序i1234567891di2、代码di= m=mean(di)%m为所求的算术平均值v=di-m%v为所求的残差a=sum(v(:)%求残差的和 af=v.2b=sum(f(:) %残差的平方和 bc=sqrt(b/9)%单次测量的标准偏差d=c/sqrt(10)%算术平均值的标准偏差x=1:10plot(x,v, . );%残余误差的分布曲线3、结果算术平均值 d = 残余误差 vi di d =( 010vii110浮点数规则,实际为 0)
2、 vi2 =i110vi2 单次测量的标准偏差: i 1n1 标准偏差 d n极限误差 limd 3 d =圆柱直径的测量结果:ddlim d =4、利用 MATLAB画出残余误差 vi 分布曲线5、利用 MATLAB的标准差函数求出单次测量的标准偏差。s=std(di) %;用标准差函数 std 求单次测量的标准偏差s =实验二 利用 MATLAB对测试数据进行线性回归分析(设计型)1、求出某测试系统输出电压 (U) 与标准压力计读数 (P) 的回归方程;序号123456p/MPa020406080100U/V由 matlab 利用矩阵法可得 U= +2、对所求回归方程进行方差分析及显著性检
3、验;方差分析表源来和 方平度 由自差方F度 著显余残 归回)差 误算计(4140+计总5所得的回归方程式在 =水平上显著,可信赖程度为 99%以上,高 度显著。3、根据回归方程画出拟合曲线;MATLAB函4、利用 MATLAB的最小二乘法函数画出拟合曲线,体会用 数的优越性。5、代码U=u=mean(U) %各电压的平均值x=1 0;1 20;1 40;1 60;1 80;1 100 %第二列为标准压力计读数b=(x.*x)(-1)*(x.)*u %矩阵形式解回归系数%方差分析及显著性检验yba=mean(u) %观测值的算术平均值p=0 20 40 60 80 100y=+*p %回归点a=
4、(u-yba).2S=sum(a(:) %总的离差平方和 S c=(y-yba).2V=sum(c(:) %回归平方和 Ud=(u-y).2Q=sum(d(:) %残余平方和 QA=Q/4 %残余方差F=V/A %统计量 Ft=0:100y1=+*tfigure(1)plot(t,y1)title( 由方程画出的拟合曲线 ) xlabel( 标准压力计读数 p/Mpa ) ylabel( 输出电压 U/V ) %由最小二乘法函数画拟合曲线 f=polyfit(p,u,1) y2=polyval(f,t) figure(2) plot(p,u, *r ,t,y2, -b ) title( 由最小二乘法函数画出的拟合曲线 )xlabel( 标准压力计读数 p/Mpa )ylabel( 输出电压 U/V )