誤差理論與數(shù)據(jù)處理 實驗報告.doc

誤差理論與數(shù)據(jù)處理實 驗 指 導(dǎo) 書 姓名 學(xué)號 機(jī)械工程學(xué)院2016年05月實驗一 誤差的基本性質(zhì)與處理一、實驗內(nèi)容1對某一軸徑等精度測量8次，得到下表數(shù)據(jù)，求測量結(jié)果。序號（10-4）1234567824.67424.67524.67324.67624.67124.67824.67224.674-0.0001 0.0009 -0.0011 0.0019 -0.0031 0.0039 -0.0021 -0.00010.0002 0.0077 0.0127 0.0352 0.0977 0.1502 0.0452 0.0002Matlab程序：l=24.674,24.675,24.673,24.676,24.671,24.678,24.672,24.674;%已知測量值x1=mean(l);%用mean函數(shù)求算數(shù)平均值disp(1.算術(shù)平均值為： ,num2str(x1);v=l-x1;%求解殘余誤差disp(2.殘余誤差為： ,num2str(v);a=sum(v);%求殘差和ah=abs(a);%用abs函數(shù)求解殘差和絕對值 bh=ah-(8/2)*0.001;%校核算術(shù)平均值及其殘余誤差,殘差和絕對值小于n/2*A,bh0，故以上計算正確if bh0 disp(3.經(jīng)校核算術(shù)平均值及計算正確);else disp(算術(shù)平均值及誤差計算有誤);endxt=sum(v(1:4)-sum(v(5:8);%判斷系統(tǒng)誤差（算得差值較小，故不存在系統(tǒng)誤差）if xt0.1 disp(4.用殘余誤差法校核，差值為：,num2str(x1),較小，故不存在系統(tǒng)誤差);else disp(存在系統(tǒng)誤差);endbz=sqrt(sum(v.2)/7);%單次測量的標(biāo)準(zhǔn)差disp(5.單次測量的標(biāo)準(zhǔn)差,num2str(bz);p=sort(l);%用格羅布斯準(zhǔn)則判斷粗大誤差，先將測量值按大小順序重新排列g(shù)0=2.03;%查表g(8,0.05)的值g1=(x1-p(1)/bz;g8=(p(8)-x1)/bz;%將g1與g8與g0值比較，g1和g8都小于g0，故判斷暫不存在粗大誤差if g1g0&g8=max(Falpha) table2,7=高度顯著;elseif (Fb=min(Falpha) table2,7=顯著;else table2,7=不顯著;endtable3、程序x=12 13 14 15 16 18 20 22 24 26;y=52.0 55.0 58.0 61.0 65.0 70.0 75.0 80.0 85.0 91.0;plot(x,y,*k)title(散點圖);X=ones(size(x), x;b= regress(y,X