制药工程正交分析设计.doc

专业班级：姓 名：学 号：一 、变量说明：i,j，k 是随机变量result9 输入的数组Qa , Qb ,Qc ,Qe 是A、B、C、e的差方和p是因素水平数r是每个水平数在列中重复次数f自由度Fx3 F值F3 临界值Sa , Sb ,Sc,Se是方差二、 程序清单：#include float maxFind(float a ,float b ,float c)float max; if(ab) max = a; else max = b;if(max c) max = c ;return max;float minFind(float a ,float b,float c) float min ;if(ab) min = a;else min = b;if(cmin) min = c;return min;void main()float result9;float K01,K02,K03,K11,K12,K13,K21,K22,K23,K31,K32,K33;float x01,x02,x03,x11,x12,x13,x21,x22,x23,x31,x32,x33;float R1,R2,R3,R4;int r,p;cout 请从表中查出水平重复次数 r 和水平数 p 进行输入 r;cin p;cout请输入实验中记录的实验结果endl;for(int i = 0 ; i resulti;K01 = result0 + result1 + result2;K02 = result3 + result4 + result5;K03 = result6 + result7 + result8;K11 = result0 + result3 + result6;K12 = result1 + result4 + result7;K13 = result2 + result5 + result8;K21 = result0 + result5 + result7;K22 = result1 + result3 + result8;K23 = result2 + result4 + result6; K31 = result0 + result4 + result8;K32 = result1 + result5 + result6;K33 = result2 + result3 + result7;x01 = K01/3; x02 = K02/3; x03 = K03/3;x11 = K11/3; x12 = K12/3; x13 = K13/3;x21 = K21/3; x22 = K22/3; x23 = K23/3;x31 = K31/3; x32 = K32/3; x33 = K33/3;R1 = maxFind(x01,x02,x03)-minFind(x01,x02,x03);R2 = maxFind(x11,x12,x13)-minFind(x11,x12,x13);R3 = maxFind(x21,x22,x23)-minFind(x21,x22,x23);R4 = maxFind(x31,x32,x33)-minFind(x31,x32,x33); cout 各个影响因素的极差依次为(极差大的表明本实验受该因素影响大) n; cout A : R1 n; cout B : R2 n;cout C : R3 n;cout e: R4 n;float Ka ,Kb,Kc,Ke;float Qa , Qb ,Qc ,Qe;float Fa , Fb , Fc ;int f = p-1;float Sa , Sb ,Sc,Se;float F3;Ka = K01 + K02 + K03; Kb = K11 + K12 + K13;Kc = K21 + K22 + K23;Ke = K31 + K32 + K33;Qa = (float)1/r*(K01*K01+K02*K02+K03*K03) - (float)1/(p*r)*Ka*Ka;Qb = (float)1/r*(K11*K11+K12*K12+K13*K13) - (float)1/(p*r)*Kb*Kb;Qc = (float)1/r*(K21*K21+K22*K22+K23*K23) - (float)1/(p*r)*Kc*Kc;Qe = (float)1/r*(K31*K31+K32*K32+K33*K33) - (float)1/(p*r)*Ke*Ke; Sa = Qa/f; Sb = Qb/f ; Sc = Qc/f ; Se = Qe/f;if(Sa Se) Sa = (Sa+Se)/2;if(Sb Se) Sb = (Sb+Se)/2;if(Sc Se) Sc = (Sc+Se)/2;Fa = Sa/Se ; Fb = Sb/Se; Fc = Sc/Se; float Fx3 = Fa,Fb,Fc ;cout 请从表中查出F统计量的三个临界值依次输入：(从大到小的输入)endl;for(int j = 0;j Fj;for(int k = 0 ; k F0) cout 第 k+1 个因素影响特别显著，标注为 * n; else if(Fxk F1) cou