试验设计与数据处理第二全书

（优选）试验设计与数据处理第二版全书第一页，共三百二十二页。引言第二页，共三百二十二页。0.1试验设计与数据处理的发展概况20世纪20年代，英国生物统计学家及数学家费歇（R．A．Fisher）提出了方差分析

20世纪50年代，日本统计学家田口玄一将试验设计中应用最广的正交设计表格化数学家华罗庚教授也在国内积极倡导和普及的“优选法”我国数学家王元和方开泰于1978年首先提出了均匀设计

第三页，共三百二十二页。0.2试验设计与数据处理的意义0.2.1试验设计的目的:合理地安排试验,力求用较少的试验次数获得较好结果

例：某试验研究了3个影响因素：

A：A1，A2，A3B：B1，B2，B3C：C1，C2，C3

全面试验：27次正交试验：9次第四页，共三百二十二页。0.2.2数据处理的目的通过误差分析，评判试验数据的可靠性；确定影响试验结果的因素主次，抓住主要矛盾，提高试验效率；确定试验因素与试验结果之间存在的近似函数关系，并能对试验结果进行预测和优化；试验因素对试验结果的影响规律，为控制试验提供思路；确定最优试验方案或配方。第五页，共三百二十二页。第1章试验数据的误差分析第六页，共三百二十二页。误差分析（erroranalysis）：对原始数据的可靠性进行客观的评定误差（error）：试验中获得的试验值与它的客观真实值在数值上的不一致试验结果都具有误差，误差自始至终存在于一切科学实验过程中客观真实值——真值第七页，共三百二十二页。1.1真值与平均值1.1.1真值（truevalue）真值：在某一时刻和某一状态下，某量的客观值或实际值

真值一般是未知的相对的意义上来说，真值又是已知的平面三角形三内角之和恒为180°国家标准样品的标称值国际上公认的计量值

高精度仪器所测之值多次试验值的平均值第八页，共三百二十二页。1.1.2平均值（mean）（1）算术平均值（arithmeticmean）

等精度试验值适合：

试验值服从正态分布第九页，共三百二十二页。（2）加权平均值(weightedmean)适合不同试验值的精度或可靠性不一致时wi——权重加权和第十页，共三百二十二页。（3）对数平均值（logarithmicmean）说明：若数据的分布具有对数特性，则宜使用对数平均值对数平均值≤算术平均值如果1/2≤x1/x2≤2时，可用算术平均值代替设两个数：x1＞0，x2

＞0，则第十一页，共三百二十二页。（4）几何平均值（geometricmean）当一组试验值取对数后所得数据的分布曲线更加对称时，宜采用几何平均值。几何平均值≤算术平均值设有n个正试验值：x1，x2，…，xn，则第十二页，共三百二十二页。（5）调和平均值（harmonicmean）常用在涉及到与一些量的倒数有关的场合调和平均值≤几何平均值≤算术平均值设有n个正试验值：x1，x2，…，xn，则：第十三页，共三百二十二页。1.2误差的基本概念1.2.1绝对误差（absoluteerror）（1）定义

绝对误差＝试验值－真值或（2）说明真值未知，绝对误差也未知

可以估计出绝对误差的范围：绝对误差限或绝对误差上界或第十四页，共三百二十二页。绝对误差估算方法：最小刻度的一半为绝对误差；最小刻度为最大绝对误差；根据仪表精度等级计算：绝对误差=量程×精度等级%第十五页，共三百二十二页。1.2.2相对误差（relativeerror）（1）定义：或

或（2）说明：真值未知，常将Δx与试验值或平均值之比作为相对误差：或第十六页，共三百二十二页。可以估计出相对误差的大小范围：相对误差限或相对误差上界相对误差常常表示为百分数（%）或千分数（‰）∴第十七页，共三百二十二页。1.2.3算术平均误差(averagediscrepancy)定义式：可以反映一组试验数据的误差大小试验值与算术平均值之间的偏差——第十八页，共三百二十二页。1.2.4标准误差(standarderror)当试验次数n无穷大时，总体标准差：

试验次数为有限次时，样本标准差：表示试验值的精密度，标准差↓，试验数据精密度↑第十九页，共三百二十二页。

（1）定义：以不可预知的规律变化着的误差，绝对误差时正时负，时大时小（2）产生的原因：偶然因素（3）特点：具有统计规律小误差比大误差出现机会多正、负误差出现的次数近似相等当试验次数足够多时，误差的平均值趋向于零可以通过增加试验次数减小随机误差随机误差不可完全避免的

1.3.1随机误差（randomerror）1.3试验数据误差的来源及分类第二十页，共三百二十二页。1.3.2系统误差（systematicerror）

（1）定义：一定试验条件下，由某个或某些因素按照某一确定的规律起作用而形成的误差（2）产生的原因：多方面（3）特点：系统误差大小及其符号在同一试验中是恒定的它不能通过多次试验被发现，也不能通过取多次试验值的平均值而减小只要对系统误差产生的原因有了充分的认识，才能对它进行校正，或设法消除。

第二十一页，共三百二十二页。1.3.3过失误差（mistake）（1）定义：

一种显然与事实不符的误差（2）产生的原因：

实验人员粗心大意造成

（3）特点：可以完全避免没有一定的规律

第二十二页，共三百二十二页。1.4.1精密度（precision）（1）含义：反映了随机误差大小的程度在一定的试验条件下，多次试验值的彼此符合程度

例：甲：11.45，11.46，11.45，11.44

乙：11.39，11.45，11.48，11.50（2）说明：可以通过增加试验次数而达到提高数据精密度的目的试验数据的精密度是建立在数据用途基础之上的试验过程足够精密，则只需少量几次试验就能满足要求1.4试验数据的精准度

第二十三页，共三百二十二页。（3）精密度判断①极差（range）②标准差（standarderror）R↓，精密度↑标准差↓，精密度↑第二十四页，共三百二十二页。③方差（variance）

标准差的平方：样本方差（s2

）总体方差（σ2

）方差↓，精密度↑第二十五页，共三百二十二页。1.4.2正确度（correctness）

（1）含义：反映系统误差的大小（2）正确度与精密度的关系：

精密度不好，但当试验次数相当多时，有时也会得到好的正确度

精密度高并不意味着正确度也高

（a）（b）（c）第二十六页，共三百二十二页。1.4.3准确度（accuracy）（1）含义：反映了系统误差和随机误差的综合表示了试验结果与真值的一致程度（2）三者关系无系统误差的试验精密度：A＞B＞C正确度：A＝B＝C准确度：A＞B＞C第二十七页，共三百二十二页。有系统误差的试验精密度：A＇＞B＇＞C＇准确度：A＇＞B＇＞C＇，A＇＞B，C第二十八页，共三百二十二页。1.5.1随机误差的检验

1.5试验数据误差的统计假设检验

1.5.1.1检验（

-test）

（1）目的：对试验数据的随机误差或精密度进行检验。在试验数据的总体方差已知的情况下，（2）检验步骤：若试验数据服从正态分布，则①计算统计量第二十九页，共三百二十二页。②查临界值

服从自由度为的分布显著性水平——一般取0.01或0.05，表示有显著差异的概率双侧（尾）检验(two-sided/tailedtest)：③检验若则判断两方差无显著差异，否则有显著差异第三十页，共三百二十二页。单侧（尾）检验(one-sided/tailedtest)：左侧（尾）检验：则判断该方差与原总体方差无显著减小，否则有显著减小右侧（尾）检验则判断该方差与原总体方差无显著增大，否则有显著增大若若（3）Excel在检验中的应用

第三十一页，共三百二十二页。1.5.1.2F检验(F-test)

（1）目的：

对两组具有正态分布的试验数据之间的精密度进行比较

（2）检验步骤①计算统计量设有两组试验数据：都服从正态分布，样本方差分别为和和，则第一自由度为第二自由度为服从F分布，第三十二页，共三百二十二页。②查临界值给定的显著水平α查F分布表临界值双侧（尾）检验(two-sided/tailedtest)：③检验若则判断两方差无显著差异，否则有显著差异第三十三页，共三百二十二页。单侧（尾）检验(one-sided/tailedtest)：左侧（尾）检验：则判断该判断方差1比方差2无显著减小，否则有显著减小

右侧（尾）检验则判断该方差1比方差2无显著增大，否则有显著增大

若若（3）Excel在F检验中的应用

第三十四页，共三百二十二页。1.5.2系统误差的检验1.5.2.1t检验法（1）平均值与给定值比较①目的：检验服从正态分布数据的算术平均值是否与给定值有显著差异②检验步骤：计算统计量：服从自由度的t分布(t-distribution)——给定值（可以是真值、期望值或标准值）第三十五页，共三百二十二页。双侧检验：若则可判断该平均值与给定值无显著差异，否则就有显著差异单侧检验左侧检验若且则判断该平均值与给定值无显著减小，否则有显著减小右侧检验若且则判断该平均值与给定值无显著增大，否则有显著增大第三十六页，共三百二十二页。（2）两个平均值的比较目的：判断两组服从正态分布数据的算术平均值有无显著差异①计算统计量：两组数据的方差无显著差异时服从自由度的t分布s——合并标准差：第三十七页，共三百二十二页。两组数据的精密度或方差有显著差异时服从t分布，其自由度为：②t检验第三十八页，共三百二十二页。双侧检验：若则可判断两平均值无显著差异，否则就有显著差异单侧检验左侧检验若且则判断该平均值1较平均值2无显著减小，否则有显著减小右侧检验若且则判断该平均值1较平均值2无显著增大，否则有显著增大第三十九页，共三百二十二页。（3）成对数据的比较目的：试验数据是成对出现，判断两种方法、两种仪器或两分析人员的测定结果之间是否存在系统误差①计算统计量：

——成对测定值之差的算术平均值：——零或其他指定值——n对试验值之差值的样本标准差：服从自由度为的t分布第四十页，共三百二十二页。②t检验若否则两组数据之间存在显著的系统误差，则成对数据之间不存在显著的系统误差，（4）Excel在t检验中的应用

第四十一页，共三百二十二页。1.5.2.2秩和检验法（ranksumtest）（1）目的：两组数据或两种试验方法之间是否存在系统误差、两种方法是否等效等，不要求数据具有正态分布（2）内容：设有两组试验数据，相互独立，n1，n2分别是两组数据的个数，总假定n1≤n2；将这个试验数据混在一起，按从小到大的次序排列每个试验值在序列中的次序叫作该值的秩（rank）将属于第1组数据的秩相加，其和记为R1

R1——第1组数据的秩和（ranksum）如果两组数据之间无显著差异，则R1就不应该太大或太小第四十二页，共三百二十二页。查秩和临界值表：根据显著性水平和n1，n2，可查得R1的上下限T2和T1

检验：如果R1＞T2

或R1

＜T1，则认为两组数据有显著差异，另一组数据有系统误差如果T1＜R1＜T2，则两组数据无显著差异，另一组数据也无系统误差

第四十三页，共三百二十二页。（3）例：

设甲、乙两组测定值为：

甲：8.6，10.0，9.9，8.8，9.1，9.1

乙：8.7，8.4，9.2，8.9，7.4，8.0，7.3，8.1，6.8

已知甲组数据无系统误差，试用秩和检验法检验乙组测定值是否有系统误差。（＝0.05）解:（1）排序：秩1234567891011.511.5131415甲8.68.89.19.19.910.0乙6.87.37.48.08.18.48.78.99.2第四十四页，共三百二十二页。（2）求秩和R1

R1=7＋9＋11.5＋11.5＋14＋15＝68（3）查秩和临界值表对于＝0.05，n1=6，n2=9得T1=33，T2＝63，∴R1＞T2

故：两组数据有显著差异，乙组测定值有系统误差

第四十五页，共三百二十二页。1.5.3异常值的检验

可疑数据、离群值、异常值

一般处理原则为：在试验过程中，若发现异常数据，应停止试验，分析原因，及时纠正错误试验结束后，在分析试验结果时，如发现异常数据，则应先找出产生差异的原因，再对其进行取舍在分析试验结果时，如不清楚产生异常值的确切原因，则应对数据进行统计处理；若数据较少，则可重做一组数据对于舍去的数据，在试验报告中应注明舍去的原因或所选用的统计方法第四十六页，共三百二十二页。1.5.3.1拉依达（）检验法①内容：可疑数据xp

，若则应将该试验值剔除。②说明：计算平均值及标准偏差s时，应包括可疑值在内3s相当于显著水平＝0.01，2s相当于显著水平＝0.05第四十七页，共三百二十二页。可疑数据应逐一检验，不能同时检验多个数据

首先检验偏差最大的数

剔除一个数后，如果还要检验下一个数，应重新计算平均值及标准偏差方法简单，无须查表该检验法适用于试验次数较多或要求不高时3s为界时，要求n＞102s为界时，要求n＞5第四十八页，共三百二十二页。

有一组分析测试数据：0.128，0.129，0.131，0.133，0.135，0.138，0.141，0.142，0.145，0.148，0.167，问其中偏差较大的0.167这一数据是否应被舍去?（＝0.01）解：（1）计算③例：（2）计算偏差（3）比较3s＝3×0.01116＝0.0335＞0.027故按拉依达准则，当＝0.01时，0.167这一可疑值不应舍去第四十九页，共三百二十二页。（2）格拉布斯（Grubbs）检验法

①内容：可疑数据xp

，若

则应将该值剔除。——Grubbs检验临界值第五十页，共三百二十二页。格拉布斯（Grubbs）检验临界值G(,n)表第五十一页，共三百二十二页。②说明：计算平均值及标准偏差s时，应包括可疑值在内可疑数据应逐一检验，不能同时检验多个数据首先检验偏差最大的数

剔除一个数后，如果还要检验下一个数，应重新计算平均值及标准偏差能适用于试验数据较少时格拉布斯准则也可以用于检验两个数据偏小，或两个数据偏大的情况③例：例1-13第五十二页，共三百二十二页。（3）狄克逊（Dixon）检验法

①单侧情形将n个试验数据按从小到大的顺序排列：

x1≤x2≤…≤xn-1≤xn

如果有异常值存在，必然出现在两端，即x1

或xn计算出统计量D或D′查单侧临界值检验xn时，当

时，可剔除xn检验检验x1时，当

时，可剔除x1第五十三页，共三百二十二页。②双侧情形计算D和D′查双侧临界值检验当，判断为异常值当，判断为异常值第五十四页，共三百二十二页。③说明适用于试验数据较少时的检验，计算量较小单侧检验时，可疑数据应逐一检验，不能同时检验多个数据剔除一个数后，如果还要检验下一个数，应重新排序④例：例1-14

第五十五页，共三百二十二页。1.6.1有效数字（significancefigure）

能够代表一定物理量的数字有效数字的位数可反映试验或试验仪表的精度数据中小数点的位置不影响有效数字的位数例如：50㎜，0.050m，5.0×104μm第一个非0数前的数字都不是有效数字，而第一个非0数后的数字都是有效数字例如：29㎜和29.00㎜第一位数字等于或大于8，则可以多计一位例如：9.99

1.6有效数字和试验结果的表示第五十六页，共三百二十二页。1.6.2有效数字的运算（1）加、减运算：与其中小数点后位数最少的相同（2）乘、除运算以各乘、除数中有效数字位数最少的为准（3）乘方、开方运算：与其底数的相同：例如：2.42=5.8（4）对数运算：与其真数的相同

例如ln6.84＝1.92；lg0.00004＝－4第五十七页，共三百二十二页。（5）在4个以上数的平均值计算中，平均值的有效数字可增加一位（6）所有取自手册上的数据，其有效数字位数按实际需要取，但原始数据如有限制，则应服从原始数据。（7）一些常数的有效数字的位数可以认为是无限制的

例如，圆周率π、重力加速度g、、1/3等（8）一般在工程计算中，取2～3位有效数字第五十八页，共三百二十二页。1.6.3有效数字的修约规则≤4：舍去≥5，且其后跟有非零数字

，进1位例如：3.14159→3.142＝5，其右无数字或皆为0时，“尾留双”：若所保留的末位数字为奇数则进1若所保留的末位数字为偶数则舍弃例如：3.1415→3.1421.3665→1.366第五十九页，共三百二十二页。1.7误差的传递误差的传递：根据直接测量值的误差来计算间接测量值的误差1.7.1误差传递基本公式间接测量值y与直接测量值xi之间函数关系：全微分第六十页，共三百二十二页。函数或间接测量值的绝对误差为：相对误差为：——误差传递系数——直接测量值的绝对误差；——间接测量值的绝对误差或称函数的绝对误差。第六十一页，共三百二十二页。函数标准误差传递公式：1.7.2常用函数的误差传递公式表1-4第六十二页，共三百二十二页。1.7.3误差传递公式的应用（1）根据各分误差的大小，来判断间接测量或函数误差的主要来源：

例1-16（2）选择合适的测量仪器或方法：

例1-17第六十三页，共三百二十二页。秩和临界值表第六十四页，共三百二十二页。n检验高端异常值检验低端异常值3～78～1011～1314～30统计量D计算公式第六十五页，共三百二十二页。第2章试验数据的表图表示法第六十六页，共三百二十二页。2.1列表法将试验数据列成表格，将各变量的数值依照一定的形式和顺序一一对应起来（1）试验数据表①记录表试验记录和试验数据初步整理的表格表中数据可分为三类：原始数据中间数据最终计算结果数据第六十七页，共三百二十二页。②结果表示表表达试验结论应简明扼要第六十八页，共三百二十二页。（2）说明：三部分：表名、表头、数据资料

必要时，在表格的下方加上表外附加

表名应放在表的上方，主要用于说明表的主要内容，为了引用的方便，还应包含表号

表头常放在第一行或第一列，也称为行标题或列标题，它主要是表示所研究问题的类别名称和指标名称数据资料：表格的主要部分，应根据表头按一定的规律排列表外附加通常放在表格的下方，主要是一些不便列在表内的内容，如指标注释、资料来源、不变的试验数据等第六十九页，共三百二十二页。（3）注意：表格设计应简明合理、层次清晰，以便阅读和使用；数据表的表头要列出变量的名称、符号和单位；要注意有效数字位数；试验数据较大或较小时，要用科学记数法来表示，并记入表头，注意表头中的与表中的数据应服从下式：数据的实际值×10±n＝表中数据；数据表格记录要正规，原始数据要书写得清楚整齐，要记录各种试验条件，并妥为保管。第七十页，共三百二十二页。2.2.1常用数据图（1）线图（linegraph/chart）表示因变量随自变量的变化情况

线图分类：单式线图：表示某一种事物或现象的动态复式线图：在同一图中表示两种或两种以上事物或现象的动态，可用于不同事物或现象的比较2.2图示法第七十一页，共三百二十二页。图1高吸水性树脂保水率与时间和温度的关系第七十二页，共三百二十二页。图2某离心泵特性曲线第七十三页，共三百二十二页。（2）XY散点图（scatterdiagram）表示两个变量间的相互关系散点图可以看出变量关系的统计规律图3散点图第七十四页，共三百二十二页。（3）条形图和柱形图用等宽长条的长短或高低来表示数据的大小，以反映各数据点的差异两个坐标轴的性质不同数值轴：表示数量性因素或变量分类轴：表示的是属性因素或非数量性变量

图4不同提取方法提取率比较第七十五页，共三百二十二页。分类：单式：只涉及一个事物或现象复式：涉及到两个或两个以上的事物或现象

图5不同提取方法对两种原料有效成分提取率效果比较第七十六页，共三百二十二页。（4）圆形图和环形图①圆形图（circlechart）也称为饼图（piegraph）表示总体中各组成部分所占的比例只适合于包含一个数据系列的情况饼图的总面积看成100%，每3.6°圆心角所对应的面积为1%，以扇形面积的大小来分别表示各项的比例图6全球天然维生素E消费比例第七十七页，共三百二十二页。②环形图（circulardiagram）每一部分的比例用环中的一段表示

可显示多个总体各部分所占的相应比例，有利于比较图7全球合成、天然维生素E消费比例比较第七十八页，共三百二十二页。（5）三角形图（ternary）常用于表示三元混合物各组分含量或浓度之间的关系

三角形：等腰Rt△、等边△、不等腰Rt△等顶点：纯物质边：二元混合物三角形内：三元混合物MABS●xAxSxB＝1－xA－xS●图8等腰直角三角形坐标图第七十九页，共三百二十二页。ABCxCxBxA●xAxAxCxCxBxBMEF图9等边三角形坐标图第八十页，共三百二十二页。（6）三维表面图（3Dsurfacegraph）

三元函数Z=f(X,Y)对应的曲面图，根据曲面图可以看出因变量Z值随自变量X和Y值的变化情况

图10三维表面图

第八十一页，共三百二十二页。（7）三维等高线图（contourplot）

三维表面图上Z值相等的点连成的曲线在水平面上的投影

图11三维等高线图

第八十二页，共三百二十二页。绘制图形时应注意：（1）在绘制线图时，要求曲线光滑,并使曲线尽可能通过较多的实验点，或者使曲线以外的点尽可能位于曲线附近，并使曲线两侧的点数大致相等；（2）定量的坐标轴，其分度不一定自零起；（3）定量绘制的坐标图，其坐标轴上必须标明该坐标所代表的变量名称、符号及所用的单位，一般用纵轴代表因变量；（4）坐标轴的分度应与试验数据的有效数字位数相匹配；（5）图必须有图号和图题（图名），以便于引用，必要时还应有图注。第八十三页，共三百二十二页。2.2.2坐标系的选择坐标系（coordinatesystem）笛卡尔坐标系（又称普通直角坐标系）、半对数坐标系、对数坐标系、极坐标系、概率坐标系、三角形坐标系…...对数坐标系（semi-logarithmiccoordinatesystem）半对数坐标系双对数坐标系

第八十四页，共三百二十二页。（1）选用坐标系的基本原则：①根据数据间的函数关系线性函数：普通直角坐标系幂函数：双对数坐标系指数函数：半对数坐标②根据数据的变化情况两个变量的变化幅度都不大，选用普通直角坐标系；有一个变量的最小值与最大值之间数量级相差太大时，可以选用半对数坐标；两个变量在数值上均变化了几个数量级，可选用双对数坐标；在自变量由零开始逐渐增大的初始阶段，当自变量的少许变化引起因变量极大变化时，此时采用半对数坐标系或双对数坐标系，可使图形轮廓清楚第八十五页，共三百二十二页。例：x10204060801001000200030004000y24146080100177181188200图12普通直角坐标系第八十六页，共三百二十二页。图13对数坐标系第八十七页，共三百二十二页。（2）坐标比例尺的确定

①在变量x和y的误差Δx，Δy已知时，比例尺的取法应使试验“点”的边长为2Δx，2Δy，而且使2Δx＝2Δy＝1～2㎜，若2Δy＝2㎜，则y轴的比例尺My应为：②推荐坐标轴的比例常数M＝（1、2、5）×10±n

（n为正整数），而3、6、7、8等的比例常数绝不可用；③纵横坐标之间的比例不一定取得一致，应根据具体情况选择，使曲线的坡度介于30°～60°之间第八十八页，共三百二十二页。例2：研究pH值对某溶液吸光度A的影响，已知pH值的测量误差ΔpH＝0.1，吸光度A的测量误差ΔA＝0.01。在一定波长下，测得pH值与吸光度A的关系数据如表所示。试在普通直角坐标系中画出两者间的关系曲线。pH8.09.010.011.0吸光度A1.341.361.451.36设2ΔpH＝2ΔA＝2mm解：∵ΔpH＝0.1，ΔA＝0.01∴横轴的比例尺为纵轴的比例尺为第八十九页，共三百二十二页。图14坐标比例尺对图形形状的影响第九十页，共三百二十二页。2.3.1Excel在图表绘制中的应用（1）利用Excel生成图表的基本方法（2）对数坐标的绘制（3）双Y轴（X轴）复式线图的绘制（4）图表的编辑和修改2.3.2Origin在图形绘制中的应用

（1）简单二维图绘制的基本方法（2）三角形坐标图的绘制（3）三维图的绘制2.3计算机绘图软件在图表绘制中应用第九十一页，共三百二十二页。表2-1离心泵特性曲线测定实验的数据记录表序号流量计读数/(L/h)真空表读数/MPa压力表读数/MPa功率表读数/W12附：泵入口管径：

__________mm；泵出口管径：_______mm；真空表与压力表垂直距离：______mm；水温：_____________℃；电动机转速

r/min。第九十二页，共三百二十二页。第九十三页，共三百二十二页。第3章试验的方差分析

第九十四页，共三百二十二页。方差分析（analysisofvariance，简称ANOVA）检验试验中有关因素对试验结果影响的显著性试验指标（experimentalindex）衡量或考核试验效果的参数

因素（experimentalfactor）影响试验指标的条件

可控因素（controllablefactor）水平（leveloffactor）因素的不同状态或内容

第九十五页，共三百二十二页。3.1单因素试验的方差分析

（one-way

analysis

of

variance）3.1.1单因素试验方差分析基本问题（1）目的：检验一个因素对试验结果的影响是否显著性（2）基本命题：设某单因素A有r种水平：A1，A2，…，Ar，在每种水平下的试验结果服从正态分布在各水平下分别做了ni（i＝1，2，…，r）次试验判断因素A对试验结果是否有显著影响

第九十六页，共三百二十二页。（3）单因素试验数据表试验次数A1A2…Ai…Ar1x11x21…xi1…xr12x12x22…xi2…xr2…………………jx1jx2j…xij…xrj…………………nix1n1x2n2…xini…xrnr第九十七页，共三百二十二页。3.1.2单因素试验方差分析基本步骤（1）计算平均值组内平均值：总平均：第九十八页，共三百二十二页。（2）计算离差平方和①总离差平方和SST（sumofsquaresfortotal）表示了各试验值与总平均值的偏差的平方和反映了试验结果之间存在的总差异②组间离差平方和SSA（sumofsquareforfactorA）反映了各组内平均值之间的差异程度由于因素A不同水平的不同作用造成的第九十九页，共三百二十二页。③组内离差平方和SSe（sumofsquareforerror）反映了在各水平内，各试验值之间的差异程度由于随机误差的作用产生三种离差平方和之间关系：第一百页，共三百二十二页。（3）计算自由度（degreeoffreedom）总自由度：dfT＝n－1组间自由度：dfA

＝r－1组内自由度：dfe

＝n－r

三者关系：dfT＝dfA＋dfe（4）计算平均平方均方＝离差平方和除以对应的自由度MSA——组间均方MSe——组内均方/误差的均方第一百零一页，共三百二十二页。（5）F检验服从自由度为（dfA,dfe）的F分布（Fdistribution）对于给定的显著性水平，从F分布表查得临界值F(dfA，dfe)

如果FA

＞F(dfA，dfe)

，则认为因素A对试验结果有显著影响否则认为因素A对试验结果没有显著影响第一百零二页，共三百二十二页。（6）方差分析表若FA

＞F0.01(dfA，dfe)

，称因素A对试验结果有非常显著的影响，用“**”号表示；若F0.05(dfA，dfe)

＜

FA

＜F0.01(dfA，dfe)

，则因素A对试验结果有显著的影响，用“*”号表示；若FA

＜F0.05(dfA，dfe)

，则因素A对试验结果的影响不显著单因素试验的方差分析表差异源SSdfMSF显著性组间（因素A）SSAr－1MSA＝SSA／（r－1）MSA／MSe组内（误差）SSen－rMSe＝SSe／（n－r）总和SSTn－1第一百零三页，共三百二十二页。3.1.3Excel在单因素试验方差分析中的应用利用Excel“分析工具库”中的“单因素方差分析”工具第一百零四页，共三百二十二页。3.2双因素试验的方差分析讨论两个因素对试验结果影响的显著性，又称“二元方差分析”3.2.1双因素无重复试验的方差分析（1）双因素无重复试验B1B2…BsA1x11x12…x1sA2x21x22…x2s……………Arxr1xr2…xrs第一百零五页，共三百二十二页。（2）双因素无重复试验方差分析的基本步骤①计算平均值总平均：

Ai水平时：

Bj水平时：第一百零六页，共三百二十二页。②计算离差平方和总离差平方和：因素A引起离差的平方和：因素B引起离差的平方和：误差平方和：第一百零七页，共三百二十二页。③计算自由度SSA的自由度：dfA

＝r－1SSB的自由度：dfB＝s－1SSe的自由度：dfe＝（r－1）（s－1）SST的自由度：dfT＝n－1＝rs－1dfT＝dfA＋dfB＋

dfe④计算均方

第一百零八页，共三百二十二页。⑤F检验FA服从自由度为（dfA,dfe)的F分布；FB服从自由度为（dfB,dfe)的F分布；对于给定的显著性水平，查F分布表：

F（dfA,dfe)，F（dfB,dfe)若FA＞F

（dfA,dfe)，则因素A对试验结果有显著影响，否则无显著影响；若FB＞F

（dfB,dfe)，则因素B对试验结果有显著影响，否则无显著影响；第一百零九页，共三百二十二页。差异源SSdfMSF显著性因素ASSAr－1因素BSSBs－1误差SSe总和SSTrs－1⑥无重复试验双因素方差分析表无重复试验双因素方差分析表第一百一十页，共三百二十二页。因素B1B2…BsA1…A2………………Ar…3.2.2双因素重复试验的方差分析（1）双因素重复试验方差分析试验表双因素重复试验方差分析试验表第一百一十一页，共三百二十二页。（2）双因素重复试验方差分析的基本步骤①计算平均值总平均：任一组合水平（Ai，Bj）上：Ai水平时：Bj水平时：第一百一十二页，共三百二十二页。②计算离差平方和总离差平方和：因素A引起离差的平方和：因素B引起离差的平方和：交互作用A×B引起离差的平方和：误差平方和：第一百一十三页，共三百二十二页。③计算自由度SSA的自由度：dfA

＝r－1SSB的自由度：dfB＝s－1SSA×B的自由度：dfA×B

＝(r－1)(s－1)SSe的自由度：dfe＝rs(c－1)SST的自由度：dfT＝n－1＝rsc－1dfT＝dfA＋dfB＋

dfA×B＋

dfe第一百一十四页，共三百二十二页。④计算均方第一百一十五页，共三百二十二页。⑤F检验若FA＞F

（dfA,dfe)，则认为因素A对试验结果有显著影响，否则无显著影响；若FB＞F

（dfB,dfe)，则认为因素B对试验结果有显著影响，否则无显著影响；若FA×B＞F

（dfA×B,dfe)，则认为交互作用A×B对试验结果有显著影响，否则无显著影响。第一百一十六页，共三百二十二页。⑥重复试验双因素方差分析表第一百一十七页，共三百二十二页。3.2.3Excel在双因素方差分析中的应用（1）双因素无重复试验方差分析利用“分析工具库”中的“无重复双因素方差分析”工具（2）双因素重复试验方差分析利用“分析工具库”中的“重复双因素方差分析”工具第一百一十八页，共三百二十二页。第4章试验数据的回归分析第一百一十九页，共三百二十二页。4.1基本概念（1）相互关系①确定性关系：变量之间存在着严格的函数关系②相关关系：变量之间近似存在某种函数关系（2）回归分析（regressionanalysis）

处理变量之间相关关系的统计方法确定回归方程：变量之间近似的函数关系式检验回归方程的显著性

试验结果预测第一百二十页，共三百二十二页。4.2一元线性回归分析4.2.1一元线性回归方程的建立

（1）最小二乘原理设有一组试验数据（如表），若x，y符合线性关系xx1x2……xnyy1y2……yn第一百二十一页，共三百二十二页。计算值与试验值yi不一定相等

与yi之间的偏差称为残差：a，b——回归系数（regression

coefficient）——回归值/拟合值，由xi代入回归方程计算出的y值。一元线性回归方程：第一百二十二页，共三百二十二页。残差平方和：残差平方和最小时，回归方程与试验值的拟合程度最好求残差平方和极小值：第一百二十三页，共三百二十二页。正规方程组（normal

equation）：解正规方程组：第一百二十四页，共三百二十二页。简算法：第一百二十五页，共三百二十二页。4.2.2一元线性回归效果的检验（1）相关系数检验法①相关系数（correlation

coefficient）：描述变量x与y的线性相关程度定义式：第一百二十六页，共三百二十二页。②相关系数特点：－1≤r≤1r＝±1：x与y有精确的线性关系第一百二十七页，共三百二十二页。r＜0：x与y负线性相关（negativelinearcorrelation）r＞0：x与y正线性相关（positive

linearcorrelation）第一百二十八页，共三百二十二页。r≈0时，x与y没有线性关系，但可能存在其它类型关系相关系数r越接近1，x与y的线性相关程度越高试验次数越少，r越接近1第一百二十九页，共三百二十二页。当，说明x与y之间存在显著的线性关系对于给定的显著性水平α，查相关系数临界值rmin③相关系数检验第一百三十页，共三百二十二页。（2）F检验①离差平方和总离差平方和：回归平方和（regression

sum

of

square）：残差平方和：三者关系：第一百三十一页，共三百二十二页。②自由度SST的自由度：dfT＝n－1SSR的自由度：dfR＝1SSe的自由度：dfe＝n－2三者关系：dfT＝dfR＋dfe③均方

第一百三十二页，共三百二十二页。④F检验F服从自由度为（1，n－2）的F分布给定的显著性水平α下，查得临界值：Fα（1，n－2)若F＞Fα（1，n－2)，则认为x与y有明显的线性关系，所建立的线形回归方程有意义第一百三十三页，共三百二十二页。⑤方差分析表第一百三十四页，共三百二十二页。4.3多元线性回归分析（1）多元线性回归形式试验指标（因变量）y与m个试验因素（自变量）xj（j=1,2,…,m）多元线性回归方程：4.3.1多元线性回归方程的建立偏回归系数：第一百三十五页，共三百二十二页。（2）回归系数的确定根据最小二乘法原理：求偏差平方和最小时的回归系数偏差平方和：根据：

得到正规方程组，正规方程组的解即为回归系数。第一百三十六页，共三百二十二页。4.3.2多元线性回归方程显著性检验（1）F检验法总平方和：回归平方和：残差平方和：第一百三十七页，共三百二十二页。F服从自由度为（m，n－m－1）的分布给定的显著性水平α下，若F＞Fα(m，n－m－1)，则y与x1，x2，…，xm间有显著的线性关系方差分析表：

第一百三十八页，共三百二十二页。（2）相关系数检验法复相关系数（multiplecorrelationcoefficient）R：

反映了一个变量y与多个变量（x1，x2，…，xm

）之间线性相关程度

计算式：R＝1时，y与变量x1，x2，…，xm之间存在严格的线性关系R≈0时，y与变量x1，x2，…，xm之间不存在线性相关关系当0＜R＜1时，变量之间存在一定程度的线性相关关系R＞Rmin时，y与x1，x2，…，xm之间存在密切的线性关系

R一般取正值，0≤R≤1

第一百三十九页，共三百二十二页。4.3.3因素主次的判断（1）偏回归系数的标准化设偏回归系数bj的标准化回归系数为Pj：

Pj越大，则对应的因素（xj）越重要第一百四十页，共三百二十二页。（2）偏回归系数的显著性检验计算每个偏回归系数的偏回归平方和SSj

：

SSj＝bjLjy

SSj的大小表示了因素xj对试验指标y影响程度，对应的自由度dfj＝1

服从自由度为（1，n－m－1）的F分布

如果若F＜Fα(1，n－m－1)，，则说明xj对y的影响是不显著的，这时可将它从回归方程中去掉，变成（m－1）元线性方程第一百四十一页，共三百二十二页。（3）偏回归系数的t检验计算偏回归系数的标准差：

t值的计算：单侧t分布表检验：→如果说明xj对y的影响显著，否则影响不显著，第一百四十二页，共三百二十二页。4.4.1一元非线性回归分析通过线性变换，将其转化为一元线性回归问题：直角坐标中画出散点图；推测y与x之间的函数关系；线性变换；用线性回归方法求出线性回归方程；返回到原来的函数关系，得到要求的回归方程

4.4非线性回归分析第一百四十三页，共三百二十二页。4.4.2一元多项式回归任何复杂的一元连续函数都可用高阶多项式近似表达：可以转化为多元线性方程：4.4.3多元非线性回归如果试验指标y与多个试验因素xj之间存在非线性关系，如二次回归模型：第一百四十四页，共三百二十二页。4.5Excel在回归分析中的应用4.5.1“规划求解”在回归分析中应用解方程组最优化

4.5.2Excel内置函数在回归分析中应用4.5.3Excel图表功能在回归分析中的应用4.5.4分析工具库在回归分析中应用第一百四十五页，共三百二十二页。第5章优选法第一百四十六页，共三百二十二页。优选法：根据生产和科研中的不同问题，利用数学原理，合理地安排试验点，减少试验次数，以求迅速地找到最佳点的一类科学方法。适用于：试验指标与因素间不能用数学形式表达表达式很复杂第一百四十七页，共三百二十二页。x1x2bx35.1单因素优选法基本命题试验指标f(x)是定义区间(a，b)的单峰函数用尽量少的试验次数，来确定f(x)的最大值的近似位置

5.1.1来回调试方法

x1x2ab若f(x1)<f(x2)若f(x2)<f(x3)x3x1x2x4……第一百四十八页，共三百二十二页。x35.1.2黄金分割法（0.618法）黄金分割：优选步骤：x20.6180.382x1ab0.6180.382x2x1b……第一百四十九页，共三百二十二页。5.1.3分数法菲波那契数列：F0＝1，F1＝1，Fn＝Fn-1＋Fn-2(n≥2)1，1，2，3，5，8，13，21，34，55，89，144，…分数：第一百五十页，共三百二十二页。x42/5x3分数法优选方法：适用于：试验值只能取整数的情况试验次数有限时x1x25/83/8x1x23/5x1x32/31/3第一百五十一页，共三百二十二页。分数法试验次数：第一百五十二页，共三百二十二页。B(无电)甲（有电）乙（无电）A(有电)5.1.4对分法特点：每次只做1次试验每次试验区间可以缩小一半适用条件：要有一个标准（或具体指标）要预知该因素对指标的影响规律优选方法：第一百五十三页，共三百二十二页。5.1.5抛物线法在三个试验点x1，x2，x3，且x1＜x2＜x3，分别得试验值y1，y2，y3，根据Lagrange插值法可以得到一个二次函数：设二次函数在x4取得最大值：第一百五十四页，共三百二十二页。在x

＝x4处做试验，得试验结果y4假定y1，y2，y3，y4中的最大值是由xi’给出除xi’之外，在x1，x2，x3和x4中取较靠近xi’的左右两点，将这三点记为x1’，x2’，x3’此处x1’＜x2’＜x3,

，若在处的函数值分别为y1’，y2’，y3’,……第一百五十五页，共三百二十二页。5.1.6分批试验法（1）均分法每批做2n个试验

先把试验范围等分为（2n+1）段，在2n个分点上作第一批试验比较结果，留下较好的点，及其左右一段然后把这两段都等分为（n+1）段分点处做第二批试验**第一百五十六页，共三百二十二页。（2）比例分割法每一批做2n＋1个试验把试验范围划分为2n+2段，相邻两段长度为a和b(a＞b)在（2n＋1）个分点上做第一批试验，比较结果，在好试验点左右留下一长一短把a分成2n＋2段，相邻两段为a1，b1(a1＞b1)，且a1＝b第一百五十七页，共三百二十二页。长短段的比例：当n=0时，λ=0.618第一百五十八页，共三百二十二页。5.1.7逐步提高法（爬山法）方法：找一个起点寻找方向

注意：起点步距：“两头小，中间大”AB＜AC＞AD＞CE＜DF第一百五十九页，共三百二十二页。5.1.8多峰情况（1）不论“单峰”还是“多峰”，按前述方法优选（2）先做一批分布得比较均匀、疏松的试验，看是否有“多峰”现象，分别找出这些“峰”第一百六十页，共三百二十二页。5.2双因素优选法命题迅速地找到二元函数z＝f(x，y)的最大值，及其对应的（x，y）点的问题假定是单峰问题双因素优选法的几何意义第一百六十一页，共三百二十二页。Q5.2.1对开法优选范围：a＜x＜b，c＜y＜d优选方法：abdcPbQR第一百六十二页，共三百二十二页。P2P15.2.2旋升法（从好点出发法）优选范围：

a＜x＜b，c＜y＜d优选方法：abdcbP2P3第一百六十三页，共三百二十二页。RPQ5.2.3平行线法两个因素：一个易调整，另一个不易调整时优选范围：

a＜x＜b，c＜y＜d优选方法：（设：x易调整，y不易调整）abdc0.3820.618第一百六十四页，共三百二十二页。5.2.4按格上升法将试验区域画上格子将分数法与上述方法结合起来第一百六十五页，共三百二十二页。5.2.5翻筋斗法ACBDEFGF′G′第一百六十六页，共三百二十二页。优选法在因素主次判断中的应用：在因素的试验范围内做两个试验（可选0.618和0.382两点）如果这两点的效果差别显著，则为主要因素如果这两点效果差别不大在（0.382～0.618）、（0～0.382）和（0.618～1）三段的中点分别再做一次试验如果仍然差别不大，则此因素为非主要因素可将该因素固定在0.382～0.618间的任一点当对某因素做了五点以上试验后，如果各点效果差别不明显，则该因素为次要因素第一百六十七页，共三百二十二页。第6章正交试验设计第一百六十八页，共三百二十二页。6.1概述适合多因素试验全面试验：每个因素的每个水平都相互搭配进行试验例：3因素4水平的全面试验次数≥43=64次正交试验设计(orthogonaldesign)：利用正交表科学地安排与分析多因素试验的方法例：3因素4水平的正交试验次数：16第一百六十九页，共三百二十二页。6.1.1正交表(orthogonaltable)（1）等水平正交表：各因素水平数相等的正交表①记号：Ln(rm)

L——正交表代号n——正交表横行数（试验次数）r——因素水平数m——正交表纵列数(最多能安排的因数个数)第一百七十页，共三百二十二页。第一百七十一页，共三百二十二页。第一百七十二页，共三百二十二页。②等水平正交表特点表中任一列，不同的数字出现的次数相同表中任意两列，各种同行数字对（或称水平搭配）出现的次数相同两性质合称为“正交性”：使试验点在试验范围内排列整齐、规律，也使试验点在试验范围内散布均匀第一百七十三页，共三百二十二页。（2）混合水平正交表

各因素的水平数不完全相同的正交表第一百七十四页，共三百二十二页。混合水平正交表性质：（1）表中任一列，不同数字出现次数相同（2）每两列，同行两个数字组成的各种不同的水平搭配出现的次数是相同的，但不同的两列间所组成的水平搭配种类及出现次数是不完全相同第一百七十五页，共三百二十二页。6.1.2正交试验设计的优点能均匀地挑选出代表性强的少数试验方案由少数试验结果，可以推出较优的方案可以得到试验结果之外的更多信息第一百七十六页，共三百二十二页。6.2.1单指标正交试验设计及其结果的直观分析例：单指标：乳化能力因素水平：3因素3水平（假定因素间无交互作用）6.2正交试验设计结果的直观分析法第一百七十七页，共三百二十二页。（1）选正交表要求：因素数≤正交表列数因素水平数与正交表对应的水平数一致选较小的表选L9(34)第一百七十八页，共三百二十二页。（2）表头设计将试验因素安排到所选正交表相应的列中因不考虑因素间的交互作用，一个因素占有一列（可以随机排列）空白列（空列）：最好留有至少一个空白列第一百七十九页，共三百二十二页。（3）明确试验方案第一百八十页，共三百二十二页。（4）按规定的方案做试验，得出试验结果注意：按照规定的方案完成每一号试验试验次序可随机决定试验条件要严格控制第一百八十一页，共三百二十二页。（5）计算极差，确定因素的主次顺序三个符号：Ki：表示任一列上水平号为i时，所对应的试验结果之和。ki

：ki=Ki/s，其中s为任一列上各水平出现的次数R（极差）：在任一列上

R=max｛K1,K2,K3｝－min｛K1,K2,K3｝，或R=max｛k1,k2,k3｝－min｛k1,k2,k3｝第一百八十二页，共三百二十二页。第一百八十三页，共三百二十二页。R越大，因素越重要若空列R较大，可能原因：漏掉某重要因素因素之间可能存在不可忽略的交互作用

第一百八十四页，共三百二十二页。（6）优方案的确定优方案：在所做的试验范围内，各因素较优的水平组合若指标越大越好，应选取使指标大的水平若指标越小越好，应选取使指标小的水平还应考虑：降低消耗、提高效率等第一百八十五页，共三百二十二页。（7）进行验证试验，作进一步的分析优方案往往不包含在正交实验方案中，应验证优方案是在给定的因素和水平的条件下得到的，若不限定给定的水平，有可能得到更好的试验方案对所选的因素和水平进行适当的调整，以找到新的更优方案趋势图第一百八十六页，共三百二十二页。正交试验设计的基本步骤：(1)明确试验目的，确定评价指标(2)挑选因素(包括交互作用)，确定水平(3)选正交表，进行表头设计(4)明确试验方案，进行试验，得到结果(5)对试验结果进行统计分析(6)进行验证试验，作进一步分析第一百八十七页，共三百二十二页。6.2.2多指标正交试验设计及其结果的直观分析两种分析方法：综合平衡法综合评分法第一百八十八页，共三百二十二页。（1）综合平衡法先对每个指标分别进行单指标的直观分析对各指标的分析结果进行综合比较和分析，得出较优方案第一百八十九页，共三百二十二页。②例三个指标：提取物得率总黄酮含量葛根素含量

三个指标都是越大越好第一百九十页，共三百二十二页。

对三个指标分别进行直观分析：提取物得率：因素主次：CAB优方案：C3A2B2

或C3A2B3

总黄酮含量：因素主次：ACB优方案：A3C3B3

葛根素含量：因素主次：CAB优方案：C3A3B2

综合平衡：A3B2C3

第一百九十一页，共三百二十二页。③综合平衡原则：次服从主（首先满足主要指标或因素）少数服从多数降低消耗、提高效率

④综合平衡特点：计算量大信息量大有时综合平衡难第一百九十二页，共三百二十二页。（2）综合评分法①综合评分法：根据各个指标的重要程度，对得出的试验结果进行分析，给每一个试验评出一个分数，作为这个试验的总指标进行单指标试验结果的直观分析法第一百九十三页，共三百二十二页。②评分方法：直接给出每一号试验结果的综合分数对每号试验的每个指标分别评分，再求综合分若各指标重要性相同：各指标的分数总和若各指标重要性不相同：各指标的分数加权和第一百九十四页，共三百二十二页。③如何对每个指标评出分数非数量性指标：依靠经验和专业知识给出分数有时指标值本身就可以作为分数，如回收率、纯度等用“隶属度”来表示分数：第一百九十五页，共三百二十二页。④例两个指标：取代度、酯化率两个指标重要程度不同综合分数＝取代度隶属度×0.4＋酯化率隶属度×0.6第一百九十六页，共三百二十二页。⑤综合评分法特点将多指标的问题，转换成了单指标的问题，计算量小准确评分难第一百九十七页，共三百二十二页。6.2.3有交互作用的正交试验设计（1）交互作用的判断设有两个因素A和B，各取两水平在每个组合水平上做试验，根据试验结果判断第一百九十八页，共三百二十二页。A1A2B12535B23040A1A2B12535B23015第一百九十九页，共三百二十二页。（2）有交互作用的正交试验设计及其结果的直观分析例：3因素2水平交互作用：A×B、A×C指标：吸光度，越大越好第二百页，共三百二十二页。①选表应将交互作用看成因素按5因素2水平选表：L8（27）②表头设计

交互作用应该占有相应的列——交互作用列交互作用列是不能随意安排表头设计两种方法：查交互作用表查表头设计表

第二百零一页，共三百二十二页。③明确试验方案、进行试验、得到试验结果第二百零二页，共三百二十二页。④计算极差、确定因素主次注意：排因素主次顺序时，应该包括交互作用⑤优方案的确定

如果不考虑因素间的交互作用，优方案：A2B2C1

交互作用A×C比因素C对试验指标的影响更大因素A，C水平搭配表第二百零三页，共三百二十二页。因素A，C水平搭配表A1A2C1(y1+y3)/2=(0.484+0.532)/2=0.508(y5+y7)/2=(0.472+0.554)/2=0.513C2(y2+y4)/2=(0.448+0.516)/2=0.482(y6+y8)/2=(0.480+0.552)/2=0.516第二百零四页，共三百二十二页。说明：表头设计中的“混杂”现象（一列安排多个因素或交互作用）高级交互作用，如A×B×C，一般不考虑r水平两因素间的交互作用要占r－1列，当r＞2时，不宜用直观分析法即使不考虑交互作用，最好仍与有交互作用时一样，按规定进行表头设计第二百零五页，共三百二十二页。6.2.4混合水平的正交试验设计两种方法：直接利用混合水平的正交表拟水平法：将混合水平的问题转化成等水平问题来处理第二百零六页，共三百二十二页。6.2.5Excel在直观分析中应用函数SUMIF绘制趋势图第二百零七页，共三百二十二页。（1）直接利用混合水平的正交表例注意：不同列Ki与ki的计算计算极差时，按ki计算混合水平正交表也可以安排交互作用

第二百零八页，共三百二十二页。（2）拟水平法例拟水平：将现有较好的水平重复一次注意：有拟水平的列，Ki，ki计算计算极差时，按ki计算有拟水平的因素确定优水平时，应按ki确定可以对多个因素虚拟水平

第二百零九页，共三百二十二页。6.3正交试验设计结果的方差分析法

能估计误差的大小能精确地估计各因素的试验结果影响的重要程度第二百一十页，共三百二十二页。6.3.1方差分析的基本步骤与格式设：用正交表Ln(rm)来安排试验试验结果为yi（i=1,