matlab数据分析方法

第四讲 数据分析方法 第一节、数据拟合 问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 ：给定一批数据点（输入变量与输出变量的数据），需确定满足特定要求的曲线或 曲面。如果输入变量和输出变量都只有一个，则属于一元函数的拟合和插值；而若输入变量 有多个，则为多元函数的拟合和插值（有点回归分析的意思） 解决 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 ： （1） 若要求所求曲线（面）通过所给所有数据点，就是插值问题； （2） 若不要求曲线（面）通过所有数据点，而是要求它反映对象整体的变化趋势，这就 是数据拟合，又称曲线拟合或曲面拟合。 注意：插值和拟合都是要根据一组数据构造一个函数作为近似，由于近似的要求不同， 二者的数学方法上是完全不同的。而面对一个实际问题，究竟应该用插值还是拟合，有时 容易确定，有时则并不明显。 例 1：下面数据是某次实验所得，希望得到 X 和 f 之间的关系？ x 1 2 4 7 9 12 13 15 17 f 1.5 3.9 6.6 11.7 15.6 18.8 19.6 20.6 21.1 曲线拟合问题最常用的解法——最小二乘法的基本思路 第一步：确定拟合的函数类型 ),,,;( 21 maaaxfy "= ，其中 为待定系数。 （函数类型的确定可以根据内在的规律确定，如果无现成的规则，则可以通过散点图，联系 曲线的形状进行分析） maaa ,,, 21 " 第二步：确定 的最小二乘准则：要求 个已知点 与曲线 的距离 的平方和 最小 。 maaa ,,, 21 " n ),( ii yx )(xfy = id ∑ = − n i ii xfy 1 2))(( 用 MATLAB 作拟合 1．多项式拟合。作多项式 拟合，可利用 mmm axaxay +++= − "110 a=polyfit(x,y,m)—其中 x,y 为给出的数据，m 为多项式的次数。 多项式在 x 处的值 y 可用以下命令计算： y=polyval（a,x） 2．用 MATLAB 作非线性最小二乘拟合 Matlab 的提供了两个求非线性最小二乘拟合的函数：lsqcurvefit 和 lsqnonlin。两个命令 都要先建立 M-文件 fun.m，在其中定义函数 f(x)。 （1）x = lsqcurvefit (‘fun’,x0,xdata,ydata); （2）x =lsqcurvefit (‘fun’,x0,xdata,ydata,options); （3）x = lsqcurvefit (‘fun’,x0,xdata,ydata,options,’grad’); （4）[x, options] = lsqcurvefit (‘fun’,x0,xdata,ydata,…); （5）[x, options,funval] = lsqcurvefit (‘fun’,x0,xdata,ydata,…); 1 （6）[x, options,funval, Jacob] = lsqcurvefit (‘fun’,x0,xdata,ydata,…); （7）x=lsqnonlin（‘fun’，x0）； （8）x= lsqnonlin （‘fun’，x0，options）； （9）x= lsqnonlin （‘fun’，x0，options，‘grad’）； （10）[x，options]= lsqnonlin （‘fun’，x0，…）； （11）[x，options，funval]= lsqnonlin （‘fun’，x0，…）； 在使用 lsqcurvefit 与 lsqnonlin 命令时，共同的问题是要先知道函数的类型，而拟合其实是决 定函数中的待定系数。 第二节 插值 插值的基本问题：给出 个数对，(n , ( )), 1, 2, ,i iP f P i n= " ，求点 处对应的函数值P ( )f P 。 一、一维插值 已知 个节点 ，求任意点 处的函数值 。常用的插值方法 有拉格朗日多项式插值、牛顿插值、分段线性插值、Hermite 插值和三次样条插值。 1+n niyx ii ,,1,0),,( "= *x *y 9 分段线性插值：将各数据点用折线连接起来 9 多项式插值：求一个多项式通过所有数据点，可以假设出多项式的系数，最后通过 求解方程得到每个系数（拉格朗日插值，用 次多项式描述n 1+n 个点） 9 样条插值：分段多项式的光滑连接（三次样条插值） 9 牛顿插值：利用节点之间的各阶差商和差分构造多项式 9 Hermite 插值：对插值函数，不仅要求它在节点处与函数同值，而且要求它与函数有 相同的一阶、二阶甚至更高阶的导数值 （1）MATLAB 命令：y=interp1(x0,y0,x,'method') method 指定插值的方法，默认为线性插值。其值可为： 'nearest' 最近项插值 'linear' 线性插值 'spline' 立方样条插值 'cubic' 立方插值。 所有的插值方法要求 x0 是单调的。 当 x0 为等距时可以用快速插值法，使用快速插值法的格式为'*nearest'、'*linear'、'*spline'、 '*cubic'。 （2）三次样条插值在 Matlab 中的实现 在 Matlab 中数据点称之为断点。如果三次样条插值没有边界条件，最常用的方法，就 是采用非扭结（not-a-knot）条件。这个条件强迫第 1 个和第 2 个三次多项式的三阶导数相 等。对最后一个和倒数第 2 个三次多项式也做同样地处理。 Matlab 中三次样条插值也有现成的函数： y=interp1(x0,y0,x,'spline')； y=spline(x0,y0,x)； pp=csape(x0,y0,conds), pp=csape(x0,y0,conds,valconds)，y=ppval(pp,x)。 其中 x0,y0 是已知数据点，x 是插值点，y 是插值点的函数值。 对于三次样条插值，我们提倡使用函数 csape，csape 的返回值是 pp 形式，要求插值点 2 的函数值，必须调用函数 ppval。 pp=csape(x0,y0)：使用默认的边界条件，即 Lagrange 边界条件。 pp=csape(x0,y0,conds,valconds)中的 conds 指定插值的边界条件，其值可为： 'complete' 边界为一阶导数，一阶导数的值在 valconds 参数中给出，若忽略 valconds 参数，则按缺省情况处理。 'not-a-knot' 非扭结条件 'periodic' 周期条件 'second' 边界为二阶导数，二阶导数的值在 valconds 参数中给出，若忽略 valconds 参数，二阶导数的缺省值为[0, 0]。 'variational' 设置边界的二阶导数值为[0,0]。 对于一些特殊的边界条件，可以通过 conds 的一个 21× 矩阵来 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示，conds 元素的取值 为 0，1，2。 conds(i)=j 的含义是给定端点 的i j 阶导数，即 conds 的第一个元素表示左边界的条件， 第二个元素表示右边界的条件，conds=[2,1]表示左边界是二阶导数，右边界是一阶导数，对 应的值由 valconds 给出。 例 2 机床加工 待加工零件的外形根据工艺要求由一组数据 给出（在平面情况下），用程控铣床 加工时每一刀只能沿 ),( yx x方向和 方向走非常小的一步，这就需要从已知数据得到加工所要求 的步长很小的 坐标。 y ),( yx 表中给出的 yx, 数据位于机翼断面的下轮廓线上，假设需要得到 x坐标每改变 0.1 时的 坐标。试完成加工所需数据，画出曲线，并求出y 0=x 处的曲线斜率和 范围 内 的最小值。 1513 ≤≤ x y x 0 3 5 7 9 11 12 13 14 15 y 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 要求用分段线性和三次样条计算。 解 编写以下程序： x0=[0 3 5 7 9 11 12 13 14 15]; y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]; x=0:0.1:15; y2=interp1(x0,y0,x); y3=interp1(x0,y0,x,'spline'); pp1=csape(x0,y0); y4=ppval(pp1,x); pp2=csape(x0,y0,'second'); y5=ppval(pp2,x); [x',y1',y2',y3',y4',y5'] subplot(2,2,1) plot(x0,y0,'+',x,y2) title('Piecewise linear') subplot(2,2,2) plot(x0,y0,'+',x,y3) 3 title('Spline1') subplot(2,2,3) plot(x0,y0,'+',x,y4) title('Spline2') dx=diff(x); dy=diff(y3); dy_dx=dy./dx; dy_dx0=dy_dx(1) ytemp=y3(131:151); ymin=min(ytemp); index=find(y3==ymin); xmin=x(index); [xmin,ymin] 二、二维插值： 前面讲述的都是一维插值，即节点为一维变量，插值函数是一元函数（曲线）。若节点 是二维的，插值函数就是二元函数，即曲面。如在某区域测量了若干点（节点）的高程（节 点值），为了画出较精确的等高线图，就要先插入更多的点（插值点），计算这些点的高程（插 值）。 （1） 网格节点：已知 个节点 nm× njmizyx ijji ,,2,1;,,2,1),,,( "" == ，构造 一个二元函数 通过全部已知节点，即),( yxfz = ),( jiij yxfz = ，再利用 插值。 ),( yxfz = （2） 散乱节点：已知 n 个节点 nizyx iii ,,2,1),,,( "= ，构造一个二元函数 通过全部已知节点，即),( yxfz = ),( iii yxfz = ，再利用 插值。 ),( yxfz = 二维插值方法： （1）最邻近插值：二维或高维情形的最邻近插值，与被插值点最邻近的节点的函数值 即为所求。 （2）分片线性插值：将四个插值点（矩形的四个顶点）处的函数值依次简记为： f (xi, yj)=f1，f (xi+1, yj)=f2，f (xi+1, yj+1)=f3，f (xi, yj+1)=f4 分两片的函数表达式如下： 第一片(下三角形区域)： ))(())((),( 23121 ji yyffxxfffyxf −−+−−+= 第二片(上三角形区域)： ))(())((),( 43141 ij xxffyyfffyxf −−+−−+= （3）双线性插值：双线性插值是一片一片的空间二次曲面构成。双线性插值函数的形 式如下： 。其中有四个待定系数，利用该函数在矩形的四个顶点 （插值节点）的函数值，得到四个代数方程，正好确定四个系数。 ))((),( dcybaxyxf ++= 用 MATLAB 作网格节点数据的插值： 4 z=interp2(x0,y0,z0,x,y,’method’) —method 可以选择 ‘nearest’：最邻近插值；‘linear’： 双线性插值；‘cubic’：双三次插值；缺省时, 双线性插值。 如果是三次样条插值，可以使用命令 pp=csape({x0,y0},z0,conds,valconds)，z=fnval(pp,{x,y}) 其中 x0,y0 分别为 维和n维向量，z0 为m nm× 维矩阵， 为矩阵，它的行数为z x的维数， 列数为 的维数，表示得到的插值，具体使用方法同一维插值。 y 例3 在一丘陵地带测量高程， x和 方向每隔100米测一个点，得高程如下表，试插值 一曲面，确定合适的模型，并由此找出最高点和该点的高程。 y x y 100 200 300 400 500 100 636 697 624 478 450 200 698 712 630 478 420 300 680 674 598 412 400 400 662 626 552 334 310 解 编写程序如下： x=100:100:500; y=100:100:400; z=[636 697 624 478 450; 698 712 630 478 420; 680 674 598 412 400; 662 626 552 334 310]; pp=csape({x,y},z') xi=100:10:500;yi=100:10:400 cz1=fnval(pp,{xi,yi}) cz2=interp2(x,y,z,xi,yi','spline') [i,j]=find(cz1==max(max(cz1))) x=xi(i),y=yi(j),zmax=cz1(i,j) 用 MATLAB 作散点数据的插值计算 cz =griddata（x,y,z,cx,cy,‘method’） —method 可以选择‘nearest’：最邻近插值；‘linear’ ： 双线性插值；‘cubic’：双三次插值；'v4'：Matlab 提供的插值方法；缺省时, 双线性插值 例 4 在某海域测得一些点（x,y）处的水深 z 由下表给出，在矩形区域（75,200）×(-50,150) 内画出海底曲面的图形。 x 129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5 y 7.5 141.5 23 147 22.5 137.5 85.5 –6.5 -81 3 56.5 –66.5 84 -33.5 z 4 8 6 8 6 8 8 9 9 8 8 9 4 9 解 编写程序如下： x=[129,140,103.5,88,185.5,195,105,157.5,107.5,77,81,162,162,117.5]; y=[7.5,141.5,23,147,22.5,137.5,85.5,-6.5,-81,3,56.5,-66.5,84,-33.5]; z=[4,8,6,8,6,8,8,9,9,8,8,9,4,9]; xi=75:1:200; 5 yi=-50:1:150; zi=griddata(x,y,z,xi,yi','cubic') subplot(1,2,1) plot(x,y,'*') subplot(1,2,2) mesh(xi,yi,zi) 第三节、回归分析方法 “回归”问题最早来源于生物界，英国生物学家兼统计学家高尔顿（Galton,1822-1911） 发现同一种族中儿子的平均高度介于其父亲的高度与种族平均高度之间。儿子的身高有返归 于种族平均身高的趋势，即回归于种族的平均身高。 回归分析是指对具有相关关系的现象，根据其关系形态，选择一个合适的数学模型，用 来近似地表示变量间的平均变化关系的一种统计方法。 回归分析的分类：按照回归模型中变量个数分（一元回归，多元回归）；按照回归曲线 的形态分（线性回归，非线性回归）；按照是否要求总体分布类型已知分（参数回归，非参 数回归） 一元线性回归特点 ¾ 两个变量中，一个是自变量，一个是因变量 ¾ 回归方程不是抽象的数学模型，而随机方程，可以进行实证 ¾ 因果关系不明显时，应同时作两个回归方程 ¾ 回归系数具有较强的经济含义 ¾ 作为回归模型的因变量是随机变量，而自变量是确定性变量，即可控变量 1．一元线性回归模型 一般地，称由 εββ ++= xy 10 确定的模型为一元线性回归模型，记为 ⎩⎨ ⎧ == ++= 2 10 ,0 σεε εββ DE xy 其中固定的未知参数 0β ， 1β 称为回归系数，自变量 x 也称为回归变量， xY 10 ββ += 称 为 y 对 x的回归直线方程。 一元线性回归分析的主要任务 1．用试验值（样本值）对 0β ， 1β 和σ 作点估计； 2．对回归系数 0β ， 1β 作假设检验； 3．在 处对 作预测，对 作区间估计。 0xx = y y 回归系数的最小二乘估计 有 组独立观测值，n niyx ii ,,2,1),,( "= 。 设 ⎩⎨ ⎧ == =++= 相互独立且 ， nii ii DE nixy εεεσεε εββ ..., ,0 ,...,2,1, 21 2 10 6 记 。最小二乘法就是选择(∑∑ == −−=== n i ii n i i xyQQ 1 2 10 1 2 10 ),( ββεββ ) 0β 和 1β 的估 计 ， 使得 。计算得到 0βˆ 1βˆ ),(min)ˆ,ˆ( 10,10 10 ββββ ββ QQ = ⎪⎩ ⎪⎨ ⎧ − −= −= 221 10 ˆ ˆˆ xx yxxy xy β ββ 其中 ∑∑ == == n i i n i i yn yx n x 11 1,1 ， ∑∑ == == n i ii n i i yxn xyx n x 11 22 1,1 。（经验）回归方程为： )(ˆˆˆˆ 110 xxyxy −+=+= βββ 。 记 ，称 为残差平方和或剩 余平方和。 的无偏估计为 ( )∑ ∑ = = −=−−== n i n i iiiie yyxyQQ 1 1 22 1010 )ˆ(ˆˆ)ˆ,ˆ( ββββ eQ 2σ )2(ˆ 2 −= nQeeσ 。 回归方程的显著性检验 对回归方程 xY 10 ββ += 的显著性检验，归结为对假设 0:;0: 1110 ≠= ββ HH 进行 检验。 假设 0: 10 =βH 被拒绝，则回归显著，认为 与y x存在线性关系，所求的线性回归方 程有意义；否则回归不显著， 与y x的关系不能用一元线性回归模型来描述，所得的回归方 程也无意义。 F 检验法：当 成立时， 0H )2,1(~)2/( −−= nFnQ UF e ，其中 ( )∑ = −= n i i yyU 1 2ˆ （回 归平方和）。若 ，拒绝 ，否则就接受 。 )2,1(1 −> − nFF α 0H 0H 回归系数的置信区间 0β 和 1β 置信水平为 α−1 的置信区间分别为 ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ +−++−− −− xx e xx e L x n nt L x n nt 2 2 10 2 2 10 1ˆ)2(ˆ,1ˆ)2(ˆ σβσβ αα 和 ⎥⎦ ⎤⎢⎣ ⎡ −+−− −− xxexxe LntLnt /ˆ)2(ˆ,/ˆ)2(ˆ 2 11 2 11 σβσβ αα ； 7 2σ 的置信水平为 α−1 的置信区间为 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ −− − )2( , )2( 2 2 2 2 1 n Q n Q ee αα χχ 。 预测 用 的回归值 作为 的预测值。 的置信水平为0y 0100 ˆˆˆ xy ββ += 0y 0y α−1 的预测区间为 [ ])(ˆ),(ˆ 0000 xyxy δδ +− ，其中 ( ) xx e L xx n ntx 2 0 2 10 11)2(ˆ)( −++−= −α σδ 。特别，当 很 大 且 在 n 0x x 附 近 取 值 时 , 的 置 信 水 平 为y α−1 的 预 测 区 间 近 似 为 ⎥⎦ ⎤⎢⎣ ⎡ +− −− 2 1 2 1 ˆˆ,ˆˆ αα σσ uyuy ee 。 2．可线性化的一元非线性回归（需要配曲线） 先对两个变量 x和 作n次试验观察得y niyx ii ,...,2,1),,( = 画出散点图，根据散点图确 定须配曲线的类型.然后由n对试验数据确定每一类曲线的未知参数 和 。采用的方法是通 过变量代换把非线性回归化成线性回归，即采用非线性回归线性化的方法。 a b 通常选择的六类曲线如下 （1）双曲线 x ba y +=1 （2）幂函数曲线 ，其中 baxy = 0,0 >> ax （3）指数曲线 其中参数 bxaey = 0>a （4）倒指数曲线 其中 xbaey /= 0>a （5）对数曲线 0,ln >+= xxbay （6）S 型曲线 xbeay −+= 1 3．多元线性回归 一般称 为高斯—马尔柯夫线性模型(k 元线性回归模 型)，并简记为 。 ⎩⎨ ⎧ == += nICOVE XY 2),(,0)( σεεε εβ ),,( 2 nIXY σβ 8 其 中 ， ， ， ， ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = ny y Y ... ... 1 ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = nknn k k xxx xxx xxx X ...1 ............... ...1 ...1 21 22221 11211 ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = kβ β β β ... 1 0 ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = nε ε ε ε ... 2 1 kk xxy βββ +++= ...110 称为回归平面方程。 线性模型 考虑的主要问题是： ),,( 2 nIXY σβ （1）用试验值（样本值）对未知参数 β 和 作点估计和假设检验，从而建立 与 之间的数量关系； 2σ y kxxx ,...,, 21 （2）在 处对 的值作预测与控制，即对 作区间估计。 kk xxxxxx 0022011 ,...,, === y y 对 iβ 和 作估计，用最小二乘法求2σ kββ ,...,0 的估计量：作离差平方和 ( )∑ = −−−−= n i ikkii xxyQ 1 2 110 ... βββ 选择 kββ ,...,0 使 达到最小。解得估计值Q ( ) ( )YXXX TT 1ˆ −=β ，得到的 代入回归平面 方程得 ，称为经验回归平面方程. 称为经验回归系数 iβˆ kk xxy βββ ˆ...ˆˆ 110 +++= iβˆ 4．多元非线性回归 多项式回归 设变量 的回归模型为 Yx, εββββ +++++= pp xxxY ...2210 其中 是已知的，p ),,2,1( pii "=β 是未知参数， ε 服从正态分布 。 称为回归多项式，上面的回归模型称为多项式回归。 ),0( 2σN k k xxxY ββββ ++++= ...2210 令 ，i=1，2，…，k 多项式回归模型变为多元线性回归模型. ii xx = 多元线性回归中的检验与预测 假设 0...: 100 ==== kH βββ 。 当 成立时，0H )1,(~)1/( / −−−−= knkFknQ kU F e 。如果 ， 则拒绝 ，认为 )1,(1 −−> − knkFF α 0H y 与 之间显著地有线性关系；否则就接受 ，认为kxx ,,1 " 0H y 与 之间线性关系不显著。 kxx ,,1 " 9 求 出 回 归 方 程 ， 对 于 给 定 自 变 量 的 值 ， 用 来预测 .称 为 的点预测。 kk xxy βββ ˆ...ˆˆˆ 110 +++= ∗kxx ,...,*1 ** 110 * ˆ...ˆˆˆ kk xxy βββ +++= εβββ ++++=∗ **110 ... kk xxy *yˆ *y 5．逐步回归分析 “最优”的回归方程就是包含所有对Y 有影响的变量, 而不包含对Y 影响不显著的变量 回归方程。选择“最优”的回归方程有以下几种方法： （1）从所有可能的因子（变量）组合的回归方程中选择最优者； （2）从包含全部变量的回归方程中逐次剔除不显著因子； （3）从一个变量开始，把变量逐个引入方程； （4）“有进有出”的逐步回归分析 逐步回归分析法的思想： （1）从一个自变量开始，视自变量Y 作用的显著程度，从大到小地依次逐个引入回归方程。 当引入的自变量由于后面变量的引入而变得不显著时，要将其剔除掉。 （2）引入一个自变量或从回归方程中剔除一个自变量，为逐步回归的一步。 （3）对于每一步都要进行Y 值检验，以确保每次引入新的显著性变量前回归方程中只包含 对Y 作用显著的变量。 （4）这个过程反复进行，直至既无不显著的变量从回归方程中剔除，又无显著变量可引入 回归方程时为止。 6．回归 MATLAB 命令： （1）多元线性回归 确定回归系数的点估计值：b=regress( Y, X ) ⎥⎥ ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢ ⎣ ⎡ = p b β β β ˆ ... ˆ ˆ 1 0 ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = nY Y Y Y ... 2 1 ⎥⎥ ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢ ⎣ ⎡ = npnn p p xxx xxx xxx X ...1 ............... ...1 ...1 21 22221 11211 求回归系数的点估计和区间估计、并检验回归模型： [b, bint,r,rint,stats]=regress(Y,X,alpha) Stats给出用于检验回归模型的统计量，有三个数值：相关系数r2、F值、与F对应的概率p。 判别规则：相关系数r2越接近 1，说明回归方程越显著； )1,(1 −−> − knkFF α 时拒绝H0， F越大，说明回归方程越显著；与F对应的概率p α< 时拒绝H0，回归模型成立。 （2）多元二项式回归 rstool（x，y，’model’, alpha） ‘model’：由下列 4 个模型中选择 1 个（用字符串输入，缺省时为线性模型） linear（线性）： mm xxy βββ +++= "110 ； 10 purequadratic（纯二次）： ∑ = ++++= n j jjjmm xxxy 1 2 110 ββββ " interaction（交叉）： ∑ ≤≠≤ ++++= mkj kjjkmm xxxxy 1 110 ββββ " quadratic（完全二次）： ∑ ≤≤ ++++= mkj kjjkmm xxxxy ,1 110 ββββ " （3）非线性回归 确定回归系数的命令： [beta，r，J]=nlinfit（x，y，’model’, beta0）：r-返回残差。 nlintool（x，y，’model’, beta0，alpha） [Y，DELTA]=nlpredci（’model’, x，beta，r，J） 求 nlinfit 或 nlintool 所得的回归函数在 x 处的预测值 Y 及预测值的显著性为 1-alpha 的置信 区间。 11