誤差理論與數(shù)據(jù)處理實(shí)驗(yàn)報(bào)告

1、誤差理論與數(shù)據(jù)處理實(shí) 驗(yàn) 報(bào) 告儀器與電子學(xué)院23楊松實(shí)驗(yàn)一 熟悉 MATLAB軟件在誤差處理中的應(yīng)用（驗(yàn)證型）1、實(shí)驗(yàn)數(shù)據(jù)序i1234567891di2、代碼di= m=mean(di)%m為所求的算術(shù)平均值v=di-m%v為所求的殘差a=sum(v(:)%求殘差的和 af=v.2b=sum(f(:) %殘差的平方和 bc=sqrt(b/9)%單次測量的標(biāo)準(zhǔn)偏差d=c/sqrt(10)%算術(shù)平均值的標(biāo)準(zhǔn)偏差x=1:10plot(x,v, '.' );%殘余誤差的分布曲線3、結(jié)果算術(shù)平均值 d = 殘余誤差 vi di d =( 010vii110浮點(diǎn)數(shù)規(guī)則，實(shí)際為 0) vi

2、2 =i110vi2 單次測量的標(biāo)準(zhǔn)偏差： i 1n1 標(biāo)準(zhǔn)偏差 d n極限誤差 limd ±3 d =±圓柱直徑的測量結(jié)果：dd±lim d =±4、利用 MATLAB畫出殘余誤差 vi 分布曲線5、利用 MATLAB的標(biāo)準(zhǔn)差函數(shù)求出單次測量的標(biāo)準(zhǔn)偏差。s=std(di)%；用標(biāo)準(zhǔn)差函數(shù) std 求單次測量的標(biāo)準(zhǔn)偏差s =實(shí)驗(yàn)二 利用 MATLAB對測試數(shù)據(jù)進(jìn)行線性回歸分析(設(shè)計(jì)型)1、求出某測試系統(tǒng)輸出電壓 (U) 與標(biāo)準(zhǔn)壓力計(jì)讀數(shù) (P) 的回歸方程；序號123456p/MPa020406080100U/V由 matlab 利用矩陣法可得 U= +

3、2、對所求回歸方程進(jìn)行方差分析及顯著性檢驗(yàn)；方差分析表源來和 方平度 由自差方F度 著顯余殘 歸回）差 誤算計(jì)（4140+計(jì)總5所得的回歸方程式在 =水平上顯著，可信賴程度為 99%以上，高 度顯著。3、根據(jù)回歸方程畫出擬合曲線；MATLAB函4、利用 MATLAB的最小二乘法函數(shù)畫出擬合曲線，體會用 數(shù)的優(yōu)越性。5、代碼U=u=mean(U')' %各電壓的平均值x=1 0;1 20;1 40;1 60;1 80;1 100%第二列為標(biāo)準(zhǔn)壓力計(jì)讀數(shù)b=(x.'*x)(-1)*(x.')*u%矩陣形式解回歸系數(shù)%方差分析及顯著性檢驗(yàn)yba=mean(u)%觀測值

4、的算術(shù)平均值p=0 20 40 60 80 100y=+*p%回歸點(diǎn)a=(u-yba).2S=sum(a(:)%總的離差平方和 S c=(y-yba).2V=sum(c(:)%回歸平方和 Ud=(u'-y).2Q=sum(d(:)%殘余平方和QA=Q/4%殘余方差F=V/A%統(tǒng)計(jì)量 Ft=0:100y1=+*tfigure(1)plot(t,y1)title( ' 由方程畫出的擬合曲線 ' ) xlabel( ' 標(biāo)準(zhǔn)壓力計(jì)讀數(shù) p/Mpa' ) ylabel( ' 輸出電壓 U/V' ) %由最小二乘法函數(shù)畫擬合曲線 f=polyfit(p,u',1) y2=polyval(f,t) figure(2) plot(p,u, '*r' ,t,y2, '-b' ) title（ ' 由最小二乘法函數(shù)畫出的擬合曲線 &