orthexpanalysis1.m

Matlab在化学工程中的应用

function OrthExpAnalysis1
% 正交试验的极差分析Variance Analysis of Orthogonal experiment
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/07/12 $
%
%   [Ref] 胡上序和陈德钊编著，实验数据的分析和处理，杭州：浙江大学出版社，
%   1996，p.273

NL = 3;             % 水平数
NF = 3;             % 因素数
r = 3;              % 每个水平有r次试验
DATA = ...          % DATA矩阵的前三列是正交试验方案，第四列为试验指标--收率
    [1  1  1  31;
    1  2  2  54;
    1  3  3  38;
    2  1  2  53;
    2  2  3  49;
    2  3  1  42;
    3  1  3  57;
    3  2  1  62;
    3  3  2  64];
A = DATA(:,1:NF);
Y = DATA(:,NF+1);

for i=1:NL          % i:水平
    for j=1:NF      % j:因素（列号）
        K(i,j) = sum(Y(find(A(:,j)==i)))
    end
end
k = K/r
% 极差
R = max(k(:,:)) - min(k(:,:))