function OrthExpAnalysis1
% 正交試驗(yàn)的極差分析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] 胡上序和陳德釗編著，實(shí)驗(yàn)數(shù)據(jù)的分析和處理，杭州：浙江大學(xué)出版社，
%   1996，p.273

NL = 3;             % 水平數(shù)
NF = 3;             % 因素?cái)?shù)
r = 3;              % 每個(gè)水平有r次試驗(yàn)
DATA = ...          % DATA矩陣的前三列是正交試驗(yàn)方案，第四列為試驗(yàn)指標(biāo)--收率
    [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:因素（列號(hào)）
        K(i,j) = sum(Y(find(A(:,j)==i)))
    end
end
k = K/r
% 極差
R = max(k(:,:)) - min(k(:,:))