北理工_數(shù)據(jù)分析_實(shí)驗(yàn)5_數(shù)據(jù)擬合

1、.(封皮頁(yè))北京理工大學(xué)現(xiàn)代數(shù)據(jù)分析實(shí)驗(yàn)4實(shí)驗(yàn)報(bào)告主講：李明學(xué)生：李經(jīng)2012/11/7;2(目錄頁(yè))實(shí)驗(yàn)5數(shù)據(jù)擬合15.1實(shí)驗(yàn)?zāi)康?5.2實(shí)驗(yàn)內(nèi)容15.2.1數(shù)據(jù)分析15.2.2彈丸穿靶速度數(shù)據(jù)分析15.3實(shí)驗(yàn)代碼及結(jié)果15.3.1數(shù)據(jù)分析15.3.2彈丸穿靶速度數(shù)據(jù)分析2 實(shí)驗(yàn)5 數(shù)據(jù)擬合5.1 實(shí)驗(yàn)?zāi)康某醪绞煜atlab的數(shù)據(jù)處理基本功能。掌握數(shù)據(jù)擬合的原理與工程使用方法。5.2 實(shí)驗(yàn)內(nèi)容5.2.1 數(shù)據(jù)分析氣體在容器中被吸引的比率Y與氣體的溫度X1和吸收液體的蒸汽壓力X2有關(guān)，其數(shù)學(xué)模型為Y=AB1X1B2X2，測(cè)得試驗(yàn)數(shù)據(jù)為：表1 氣體被吸引比率數(shù)據(jù)求Y關(guān)于X1、X2的二元線性回歸

2、方程5.2.2 彈丸穿靶速度數(shù)據(jù)分析根據(jù)初步研究認(rèn)為穿透速度與彈丸直徑、靶板厚度之間有如下關(guān)系： ，其中a和n為待求擬合參數(shù)。表2 鎢球侵徹靶板彈道極限數(shù)據(jù)彈丸直徑d/mm6.108.126.18.126.18.127.567.56靶板厚度h/mm6.756.758.458.4510.410.413.38.45穿透速度v/ms-172352097363513077981302725利用建立的關(guān)系式對(duì)表中最后兩欄的彈靶情況對(duì)穿透速度進(jìn)行預(yù)測(cè)。5.3 實(shí)驗(yàn)代碼及結(jié)果5.3.1 數(shù)據(jù)分析代碼：x1=78 113.5 130 154 169 187 206 214'x2=1 3.2 4.8 8

3、.4 12 18.5 27.5 32'y=1.5 6 10 20 30 50 80 100'x=ones(8,1),x1,x2;b=regress(y,x)結(jié)果：b = 9.4397 -0.13843.6796所以：A=9.4397 B1=-0.1384 B2=3.6796 所求方程為：Y=9.4397-0.1384X1+3.6796X2；5.3.2 彈丸穿靶速度數(shù)據(jù)分析首先對(duì)方程取對(duì)數(shù)，得：ln(V)=ln(a)+n*ln(d/h)代碼：d=6.108.12 6.1 8.12 6.1 8.12;h=6.756.75 8.45 8.45 10.410.4;v=723 520 973 635 1307 798;x=d./h;y=v;x1=log(x);y1=log(y);p=polyfit(x1,y1,1)結(jié)果：p = -1.2891 6.4431所以：n=-1.2891 ln(a)= 6.4431所以：a=628.3517所求公式為：V=628.3517*(d/h) (-1.2891)預(yù)測(cè)值1=628.3517*(7.56/13.3) (-1.2891)= 1.30