数理统计课程设计

1、二氧化碳吸附量与活性炭孔隙结构的线性回归分析摘要：本文搜集了不同孔径下不同孔容的活性炭与 CO2吸附量的实验数据。分别以同一孔径下的不同孔容作为自变量，CO2吸附量作为因变量，作出散点图。选取分布大致呈直线的一组数据为拟合的样本数据。对样本数据利用最小二乘法进行回归分析，参数确定，并对分析结果进行显著性检验。同时利用matlab的regress函数进行直线拟合。结果表明：孔径在3. 03. 5 nm之间的孔容和CO2 吸附量之间存在较好的线性关系。关键字：活性炭 孔容CO2吸附量matlab一、问题分析1.1. 数据的收集和处理CO2吸附量本文主要研究同一孔径的孔容的活性炭和 co2吸附量之间

2、的线性关系，有关 实验数据是借鉴双全，罗雪岭等人的研究成果1。以太西无烟煤为原料、硝酸钾 为添加剂，将煤粉、添加剂和煤焦油经过充分混合后挤压成条状，在600C下炭化 15 min,然后用水蒸气分别在920E和860E下活化一定时间得到2组活性炭，测 定了 CO2吸附等温线，探讨了 2组不同工艺制备的活性炭的 CO2吸附量和孔容的 关系.数据如下表所示：1 1编号孔容 /( 10 L g )17.1816.224.475.270961156426.5914.418.453.75085.69155.134.541118.9716578.39153.745.1313.429.910.39076122

3、53.754.1610.518.983.87880.511361.764.9212.123.481.672569953.675.0812.623.893.58677.812265.585.291325.188.46966.410757.797.4716.926.946.47893.210758.2105.441321.444.19198.613776.6111.8164.618.353.111411014275121.2427.739.512611498.618398.71/(mL g )0.50.8 nm 0.8T.2 nm 1.21.8 nm 1.82.2 nm 2.22.2 nm 2.5

4、3.0 nm3.03.5 nm表1:孔分布与CO2吸附值编号112是在不同添加剂量，温度，活化时间处理下的对照组。因为处理方式不同得到不同结果是互不影响的，可以看出CO2的吸附量的值是互相独立的。我们将不同孔径下的孔容分为17组。作出不同孔径下与CO2吸附量的散点图如下:量附吸2OC100-I180 -A 0o60 - :°0 (XTU40':量附吸2OC1008060处,J4002468孑容10203040506070孑溶量附吸2OC100806040100806040152025303540孑容量附80'吸2 60OC40:050100150孑溶0F.LI)Q0

5、c c0A 0门 BnU UV-Pr量附吸2OC608010086孑容10040120量80 -附吸2 60OC40''丄5060708090100 110孑溶806040量附吸2OC10080604080100 120 140 160 180 200孑容图1:不同孔容与CO2吸附量的散点图图1中从左往右依次是第1到第7组孔容，从图中可以看出第五、六、七组 的点大致分散在一条直线附近，说明两个变量之间有一定的线性相关关系。 且自 变量的变化导致因变量CO2的浓度变化，因变量变化具有独立性。我们就选取第 七组的数据进行回归分析。、问题假设1. 假设误差分布服从正态分布。2. 为了

6、简化模型，便于回归分析，我们不考虑实验中各种因素对活性炭吸附 的影响，考虑孔容与co2吸附量的数据之间的线性关系。三、模型建立3.1. 回归参数的引进回归函数y f(x) E(Y|X x)是线性函数的回归分析称为线性回归,当可控制变量只有一个时，即回归函数为 y f(x) o !x，那么(1)Yoi x:N (0,2 )归直线，0、1称为回归系数，常数 0称为一元线性回归模型，上式称为丫对x的一元线性回归方程或者一元线性回 1、2均未知。3.2回归方程的构建1x中的参数 道，需要通过样本值对它们进行估计，得到估计值 ？， 方程丫 ？， ？x，此样本方程可用作总体回归方程 y由于总体回归方程y

7、f(x) 00、1在实际中并不知?，从而得到样本回归f (x) E(Y|X x)的估计。通常可用最小二乘法估计得到公式由于总体回归方程y f(x) 01X中的参数道，需要通过样本值对它们进行估计，得到估计值 ？， 方程丫 ？) ？x，此样本方程可用作总体回归方程 y0、1在实际中并不知?，从而得到样本回归f (x) E(Y|X x)的估计。通常可用最小二乘法估计得到公式n(Xi x)(yi y)? ±j1n(人 x)2i 1y ?x其x1nn1门x，y-y，记i 1ni 11212121 xyXi :y 12 x ylxx2 2xi12 xlyyyi2 12y2i 1=i 1i 1?

8、l '/ l1xyxx?0 y ?xSeSTSl?21Rxxxx2n 2可得?l /lj'xyxx?-0y ?x2.3求一定孔容下的C02勺吸附量的回归直线方程 利用matlab对数据进行计算，结果如下表所示头验编号孔容XiCO吸附量yi2x2Yi1111564132254096736029155.182813036.015014.139153.782812883.694886.7412253.7148842883.696551.4511361.7127693806.896972.169953.698012872.965306.4712265.5148844290.257991

9、810757.7114493329.296173.9910758.2114493387.246227.41013776.6187695867.5610494.21114275201645625106501218398.7334899741.6918062.11429773.517744551820.2795689.3表2：孔容与C02吸附度的回归计算讲结果代入上上述公式可得下列计算表:x =1429.00X =119.08X2 =177445.00nx2=2129340.00i =7274.92lXXSe=201.66n=12X yi =95689.30n X y=1148271.60lxy

10、=3578.34y =773.50y =64.46y2 =51820.27 ny2 =621843.24lyy =1961.752=63.77l /1 =0.491 xy xx0 y x ? =5.88(4)表3:回归参数的计算表 由此可得线性回归方程为：y 0.49 x 5.88四、回归方程的显著性检验对回归方程是否有意义做判断就是对如下的检验问题做出判断:Ho：10vsH1 :10(5)拒绝域H。表示回归方程是显著的。利用F检验对参数进行检验。经计算有St lyy 63.77fT 11sr Jxx 48.42fR 1Se ST SR 15.35fe 10(6)(7)(8)4.1F值检验取显

11、著水平a =0.05，其拒绝域为:5 0 5 09 9 8 85 O7 7查表可得拒绝域的值为：f 4.96计算得FSr87.28，远远大于F的临界值，说明拒绝原假设，原Se/ (n 2)假设不成立，自变量和因变量有着显著的线性关系。4.2. p值检验将(6)(8)中的各平方和和自由度移入方差分析表，继续进行计算可得：来源平方和自由度均方F比P值回归1760.11760.87.2820.000b残差201.6561020.166总计1961.74911这里p值很小，因此，在显著性水平 0.01下回归方程是显著的五、计算方法的涉及和计算机的实现4.1用matlab拟合直线：先将数据以txt格式保

12、存，再用dlmread读取ASCII码文件。调用matlab 中的regress多元线性回归函数(代码见附录)，对 12个样本数据进行拟合， 作出散点图和直线拟合图在一图上如下：图2 :孔容和CO2吸附量的直线拟合O5 0 5 096 6 5 540容1孔100从图中可以看出样本点大致分布在直线附近，拟合效果比较好4.2直线参数的估计值的置信区间以及三种检验利用regess函数求出参数的估计值和置信区间以及参数的检验统计量（设 置a =0.05 ）如下：匚omnnand Windovj蚩数估计值置信区间P- G5. S843 -6. 3805 20. 1495 10.49L9 0. 3746

13、. 6092R*2=0.3972F=E7. 2SL3 p<29541-06s"2=201656A»|图3:用matlab计算的参数值和检验值。其中，RA2=0.8972指因变量（CO2吸附度）有89.7%可由模型确定，F的值 远远超过F的临界值。P远小于a,因而模型从整体上看是可用的。六、主要的结论孔容和CO2吸附量之间存在线性关系，经过显著性检验，线性方程回归效果较好,即线性方程能基本描述孔径围3. 03. 5 nm的活性炭孔容和CO2吸附量七、参考文献1双全，罗雪岭，郭哲，董明建，岳晓明.CO2吸附量与活性炭孔隙结构线性关系的研究J.中国矿业大学学报.2008（0

14、4）附录Matlab制作散点图：M=dlmread( 'co2.txt' ); %读取 ASCII 码文件 for i=1:1:7subplot(4,2,i)x1=M(:,i); y=M(:,8);plot(x1,y, 'bo');xlabel('孔容'),ylabel( 'CO2 吸附量'); endMatlab直线拟合：clc; format short g;M=dlmread( 'co2.txt' ); %读取 ASCII 码文件x1=M(:,7); y=M(:,8);plot(x1,y, 'bo&#

15、39;);b=regress(y,ones(size(x1),x1);% b= p 0 p 1'，列向量x1=sort(x1);%按升序排序，用于画图y=ones(size(x1),x1*b;%使用矩阵乘法hold on;plot(x1,y, '-r' );title('图2:孔容和C02吸附量的直线拟合') xlabel( ' 孔容 ' );ylabel( 'CO2 吸附量 ' );hold off ;Matlab 参数估计: clc; format compact ; format short g;M=dlmread( 'co2.txt' ); %读取 ASCII 码文件 x1=M(:,7); y=M(:,8);b,bint,r,rint,stats=regress(y,ones(size(x1),x1,0.05);fprintf( '%2s%5s%11sn' , ' 参数