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

1、精選優(yōu)質文檔-傾情為你奉上誤差理論與數(shù)據(jù)處理 實驗報告 儀器與電子學院楊松 專心-專注-專業(yè)實驗一 熟悉MATLAB軟件在誤差處理中的應用（驗證型）1、 實驗數(shù)據(jù)序號12345678910di24.23424.23824.23124.23024.23224.23724.23324.23524.23424.2362、 代碼di=24.234 24.238 24.231 24.230 24.232 24.237 24.233 24.235 24.234 24.236m=mean(di) %m為所求的算術平均值v=di-m %v為所求的殘差a=sum(v(:) %求殘差的和af=v.2b=sum(f

2、(:) %殘差的平方和bc=sqrt(b/9) %單次測量的標準偏差d=c/sqrt(10) %算術平均值的標準偏差x=1:10plot(x,v,'.'); %殘余誤差的分布曲線s=std(di) %；用標準差函數(shù)std求單次測量的標準偏差3、 結果算術平均值= 24.2340 殘余誤差=（0 0.0040 -0.0030 -0.0040 -0.0020 0.0030 -0.0010 0.0010 0 0.0020）-1.4211e-14（浮點數(shù)規(guī)則，實際為0）= 6.0000e-05 單次測量的標準偏差：= 0.0026 標準偏差= 8.1650e-04 極限誤差±

3、3=±0.0024圓柱直徑的測量結果：±=(24.2340±0.0024)4、 利用MATLAB畫出殘余誤差vi分布曲線。5、 利用MATLAB的標準差函數(shù)求出單次測量的標準偏差。s=std(di) %；用標準差函數(shù)std求單次測量的標準偏差s = 0.0026實驗二 利用MATLAB對測試數(shù)據(jù)進行線性回歸分析（設計型）1、 求出某測試系統(tǒng)輸出電壓(U)與標準壓力計讀數(shù)(P)的回歸方程；序號123456p/MPa020406080100U/V-0.044-0.014-0.0140.0190.0190.0360.0360.0153.2603.3553.3103.38

4、03.3373.3983.3423.3826.6756.7506.6876.7906.7226.8076.6926.77810.08510.18910.12210.21710.15110.21810.20110.23113.57513.64413.62013.58113.64013.68513.58913.61317.12017.12017.16717.16717.16417.16417.20117.201由matlab利用矩陣法可得U= -0.0663+ 0.1715p2、 對所求回歸方程進行方差分析及顯著性檢驗；方差分析表來源平方和自由度方差F顯著度回歸殘余205.8858（計算誤差）0.

5、018514 0.00464.4433e+040.01總計205.86875所得的回歸方程式在=0.01水平上顯著，可信賴程度為99%以上，高度顯著。3、 根據(jù)回歸方程畫出擬合曲線；4、 利用MATLAB的最小二乘法函數(shù)畫出擬合曲線，體會用MATLAB函數(shù)的優(yōu)越性。5、代碼U=-0.044 -0.014 -0.014 0.019 0.019 0.036 0.036 0.015; 3.260 3.355 3.310 3.380 3.337 3.398 3.342 3.382; 6.675 6.750 6.687 6.790 6.722 6.807 6.692 6.778; 10.085 10.1

6、89 10.122 10.217 10.151 10.218 10.201 10.231; 13.575 13.644 13.620 13.581 13.640 13.685 13.589 13.613; 17.120 17.120 17.167 17.167 17.164 17.164 17.201 17.201u=mean(U')' %各電壓的平均值x=1 0; 1 20; 1 40; 1 60; 1 80; 1 100 %第二列為標準壓力計讀數(shù)b=(x.'*x)(-1)*(x.')*u %矩陣形式解回歸系數(shù)%方差分析及顯著性檢驗yba=mean(u) %觀

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