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

1、.封皮頁)現(xiàn)代數(shù)據(jù)分析實(shí)驗(yàn)4 實(shí)驗(yàn)報(bào)告主講：李明學(xué)生：李經(jīng)2012/11/7.(目錄頁)實(shí)驗(yàn) 5數(shù)據(jù)擬合.15.1 實(shí)驗(yàn)?zāi)康?15.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)康某醪绞煜?Matlab 的數(shù)據(jù)處理基本功能。掌握數(shù)據(jù)擬合的原理與工程使用方法。5.2 實(shí)驗(yàn)內(nèi)容5.2.1 數(shù)據(jù)分析氣體在容器中被吸引的比率與氣體的溫度X1和吸收液體的蒸汽壓力X2有關(guān)，其數(shù)學(xué)模型為Y=A，測(cè)得試驗(yàn)數(shù)據(jù)為：表1 氣體被吸引比率數(shù)據(jù)求 Y關(guān)于 X1、X2的二元

2、線性回歸方程5.2.2 彈丸穿靶速度數(shù)據(jù)分析根據(jù)初步研究認(rèn)為穿透速度與彈丸直徑、靶板厚度之間有如下關(guān)系： ndV a* ，其中 a和 n 為待求擬合參數(shù)。h表2 彈丸直徑 d/mm 6.10 8.12 6.1靶板厚度 h/mm 6.75 6.75 8.45 8.45 10.4穿透速度 v/ms-1 723 520 973 6357.568.451307 7981302 725利用建立的關(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.4 12 1

3、8.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.10 8.12 6.1 8.12 6.1 8.12;h=6.75 6.75 8.45 8.45 10.4 10.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.28916.4431所以：n=-1.2891所以：a=628.3517ln(a)= 6.4431所求公式為：V=628.3517*(d/h) (-1.2891)預(yù)測(cè)值 1=628.3517*(7.56/13.3) (-1.2891)= 1.3015