投石器實(shí)驗(yàn)報(bào)告數(shù)據(jù)處理部分

1、1.1運(yùn)用doe部分析因設(shè)計(jì)生成實(shí)驗(yàn)表格進(jìn)行實(shí)驗(yàn)并記錄數(shù)據(jù)。數(shù)據(jù)如下：運(yùn)行序點(diǎn)類型區(qū)組ABC標(biāo)準(zhǔn)序1距離（cm）101000121221111124903-111.681793003*401000423050100052456-11001.6817936*7010007234811-11-1811291111-1929610111-1-1101801101000112521211-111128913-110-1.68179013*1411-1-1-1144115-11-1.681790015*1611-1-111612817010001722918-1101.681793018*19111-1

2、11927620-1100-1.6817920*101000121721111124703-111.681793003*4010004*5010005*6-11001.6817936*7010007*811-11-18*91111-1919610111-1-1101771101000112371211-1111210113-110-1.68179013*1411-1-1-114*15-11-1.681790015*1611-1-1116*17010001722718-1101.681793018*19111-1119*20-1100-1.6817920*1.2我們運(yùn)用響應(yīng)曲面回歸分析:距離（c

3、m）與A,B,C。其中，響應(yīng)曲面回歸:距離（cm）與A,B,C.其中，不能估計(jì)B*B,C*C并且已經(jīng)刪除。分析是使用已編碼單位進(jìn)行的。我們得到如下回歸方程：距離（cm）的估計(jì)回歸系數(shù)系數(shù)標(biāo)項(xiàng)系數(shù)準(zhǔn)誤TP常量233.6721.0311.1120.000A109.0018.215.9850.000B45.2518.212.4850.035C44.2518.212.4300.038A*A-32.1727.82-1.1560.277S=51.5109PRESS=93039.8R-Sq（調(diào)整）=77.68%R-Sq=84.55%R-Sq（預(yù)測）=39.79%距離（cm）的方差分析來源自由度SeqSSAd

4、jSSAdjMSFP回歸413064113064132660.112.310.001線性312709312709342364.315.970.001A1950489504895048.035.820.000B1163801638016380.56.170.035C1156651566515664.55.900.038平方1354835483547.51.340.277A*A1354835483547.51.340.277殘差誤差923880238802653.4失擬422911229115727.729.540.001純誤差5969969193.9合計(jì)13154521距離（cm）的異常觀測值擬

5、合值標(biāo)準(zhǔn)化觀測值標(biāo)準(zhǔn)序距離（cm）擬合值標(biāo)準(zhǔn)誤殘差殘差2289.000182.00036.424-93.000-2.55R1414490.000400.00036.42490.0002.47RR表示此觀測值含有大的標(biāo)準(zhǔn)化殘差距離（cm）的估計(jì)回歸系數(shù)，使用未編碼單位的數(shù)據(jù)項(xiàng)系數(shù)常量233.667A109.000B45.2500C44.2500A*A-32.1667其中殘差圖如下：結(jié)論：響應(yīng)曲面分析得到的結(jié)果非常不理想!首先擬合方程的系數(shù)只有84.55%，說明這個(gè)擬合方程本身就非常不準(zhǔn)確；其次，B*B、C*C的交互影響不知道什么原因估計(jì)不了。1.2我們再次使用廣義回歸來分析數(shù)據(jù)。得出結(jié)果如下：廣

6、義回歸分析：距離（cm）與A,B,C*注*無法估計(jì)B*B,它已被刪除。*注*無法估計(jì)C*C,它已被刪除?；貧w方程距離（cm）=231.444+100.563A+38.6875B+50.1875C-36.8819A*A+29.1875A*B+32.6875A*C+4.0625B*C+30.0625A*B*C使用了21個(gè)案例，19個(gè)案例包含缺失值系數(shù)項(xiàng)系數(shù)系數(shù)標(biāo)準(zhǔn)誤TP常量231.4447.778629.75400.000A100.5637.145114.07430.000B38.6887.14515.41450.000C50.1877.14517.02400.000A*A-36.88210.56

7、22-3.49190.004A*B29.1877.14514.08500.002A*C32.6887.14514.57480.001B*C4.0627.14510.56860.580A*B*C30.0627.14514.20740.001R-Sq（調(diào)整）=95.35%模型匯總S=23.3358R-Sq=97.21%PRESS=*R-Sq（預(yù)測）=*%方差分析來源自由度SeqSSAdjSSAdjMSFP回歸82277132277132846452.2700.000000A1116349107870107870198.0870.000000B123711159651596529.3170.0001

8、56C138708268672686749.3370.000014A*A195586640664012.1930.004448A*B1130249087908716.6870.001513A*C115820113971139720.9290.000638B*C19031761760.3230.580132A*B*C196409640964017.7020.001216誤差1265356535545合計(jì)20234248異常觀測值的擬合和診斷擬合值觀測值距離（cm）擬合值標(biāo)準(zhǔn)誤殘差標(biāo)準(zhǔn)化殘差811211223.3358-0*X929624616.5009503.03014R14414123.335

9、8-0*X1612812823.33580*X1927627623.33580*X2919624616.5009-50-3.03014RR表示此觀測值含有大的標(biāo)準(zhǔn)化殘差X表示受X值影響很大的觀測值。距離（cm）殘差圖率頻4t4-1亠+*V*1+擬合值與順序結(jié)論：廣義回歸的相關(guān)系數(shù)已經(jīng)達(dá)到了97%多，擬合度較好，殘差圖較之之前也理想多了。然后卻有較多異常值，我們分析異常出現(xiàn)的原因是測量時(shí)候的操作和測量失誤。操作員在放手的時(shí)候，可能對投石器臂施加了一點(diǎn)壓力，而這點(diǎn)在試驗(yàn)中很難避免；其次測量員在測量球砸中的距離時(shí)可能產(chǎn)生較大偏差，因?yàn)榍虻乃俣葮O快，落腳點(diǎn)又較小，很容易產(chǎn)生偏差最后，也是我們實(shí)驗(yàn)中的重

10、大失誤，沒有很好固定住投石器，導(dǎo)致每次試驗(yàn)中投石器會出現(xiàn)一定的晃動(dòng)，這對球投出去的距離有非常大的影響。1.4剔除異常值，進(jìn)行廣義回歸廣義回歸分析：距離（cm）與A,B,C回歸方程cm）=232.207-16.2436A-18.8547B+8.84825C+0.25827A*A+0.501459A*B-2.60973A*C+0.592714B*B-4.36696B*C-31.1512C*C+0.508654A*B*C系數(shù)項(xiàng)系數(shù)系數(shù)標(biāo)準(zhǔn)誤TP常量232.20724.90509.323710.000A-16.2447.0580-2.301440.026B-18.8557.1736-2.628330.

11、012C8.84825.31640.349510.728A*A0.2580.39660.651180.518B*B0.5930.31061.908200.063C*C-31.15132.1245-0.969700.338A*B0.5010.43781.145500.258A*C-2.6104.2169-0.618880.539B*C-4.3676.1007-0.715810.478A*B*C0.5090.64750.785550.437模型匯總S=65.0323R-Sq=72.03%R-Sq（調(diào)整）=65.37%PRESS=291448R-Sq（預(yù)測）=54.10%方差分析來源自由度SeqSS

12、AdjSSAdjMSFP回歸1045734445734445734.410.8140.000000A12142882240022400.45.2970.026397B11242112921629215.76.9080.011933C118185517516.60.1220.728456A*A12631817931793.30.4240.518476B*B1416781540015399.53.6410.063216C*C1897739773976.80.9400.337746A*B12055255495549.41.3120.258486A*C152516201619.80.3830.539340B*C1121672167.00.5120.478069A*B*C1261026102609.80.6170.436542誤差421776261776264229.2失擬1717665717665710391.6268.0080.000000純誤差2596996938.8合計(jì)52634970異常觀測值的擬合和診斷擬合值觀測值距離