905正交实验设计与结果分析 - read.doc

905正交实验设计与结果分析 ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 一． 功能，利用正交表安排实验称为正交实验设计，该程序是对正交实验的结果，运用方差分析表和F检验个因素对实验结果是否有显著影响。程序是针对不考虑个因素的交互作用的正交设计编写的。 二. 算法简介[16] 正交表一般用Ln(S)表示，其中L表示正交表，n表示安排试验的次数，S代表每一个因子的水平数；r代表试验中安排因子的最大的个数。正交表较多，如L(2) ,L(2),L(3) 选择正交表的原则：首先要求正交表水平数S与试验中每一个水平数一致；其次，要求正交表中因子数r大于或等于试验的实际因子数r1，第三，在满足前两条的前提下，适当选择试验次数n较小的正交表。 例如，考察影响化工产品转化率的因素选取了三个有关因素；反应因素（A），反应时间（B），用碱量（c），而每个因素取三种水平，如表9-3所示。 水平 因子 1 2 3 温度（A） 80℃(A1) 85℃(A2) 90℃(A3) 时间 (B) 90分(B1) 120分（B2） 150分（B3） 用碱量（C） 5%(C1) 6%(C2) 7%(C3) 应当选择正交表L9-3安排实验，该表中最多可安排4个因子。现为三个，三个因子可安排在正交表L9-3的四个列中的任意三个列，如前三列。得到表9-4。 按着表9-4可作方差分析。、 首先将总离差（平方和）分解 Q= 表9-4 1 2 3 (A) (B) (C) 组合水平 试验值 1 1 1 1 A1 B1 C1 Y1 2 1 2 2 A1 B2 C2 Y2 3 1 3 3 A1 B3 C3 Y3 4 2 1 2 A2 B1 C2 Y4 5 2 2 3 A2 B2 C3 Y5 6 2 3 1 A2 B3 C1 Y6 7 3 1 3 A3 B1 C3 Y7 8 3 2 1 A3 B2 C1 Y8 9 3 3 2 A3 B3 C2 Y9 其中 Y=1/n 可用下面公式计算Q ,Q,Q及Q.令 K= P=k W= (2) 则有 U= U= U= (3) 其中K =Y+Y+Y,即表9-4中因子A取水平时三个组合水平的试验值之和，K=Y+Y+Y,即表9-4中因子C取3水平时三个组合水平的试验值互和，余类推。利用（2）式之和（3）式，有 Q Ln(S’)安排试验，且实际因子数为r1（r1≤r）,各因子引起的离差（平方和）除以6都服从自由度为S-1的X分布，Qe/6n-1-r,(S-1)的X分布Q1/6n-1的X分布。 S=Qa(s-1), S=Qb/(s-1),S=Qc/(s-1) (5) S=Qc/[n-1-r1(s-1)] (5) 计算F值 Fa= S/ S,Fb= S/S,Fc=S/ S (6) Fa,Fb,Fc,都服从自由度为（S-1,n-1-r1(s-1)）的F分布。 =0.05）,作F检验，表查得Fα(s-1,n-1-r1(s-1))的值，将Fa、Fb、Fc，分别与Fα的比较，若 Fa≥Fα(s-1,n-1-r1(s-1)) Hb1,即认为因子A对试验结果有显著影响；否则，接受假设Hb1, 即认为因子A对试验结果没有显著影响。余类推。 F检验的结果，可安排第二批试验，以寻求更好的方案。如确认因素A对试验结果有显著影响，因素B与C无显著影响，则可根据A的原来的水平与试验结果的关系，另选A的几个水平，而保持B、C的某个水平继续作几次试