orthexpanalysis2.m

matlab在化学工程中的应用实例程序,对大家研究学习有指导作用。

function OrthExpAnalysis2
% 正交试验的极差分析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] 何为和陈际达编，优化试验设计法及其在化学中的应用，成都：
%   电子科技大学出版社，1994(p.20)

NL = 2;            % 水平数
NF = 6;            % 因素数
r = 4;             % 每个水平有r次试验
DATA = ...         % DATA矩阵的前三列是正交试验方案，第四列为试验指标--收率
    [1  1  1  1  1  1  65;
    1  1  1  2  2  2  74;
    1  2  2  1  1  2  71;
    1  2  2  2  2  1  73;
    2  1  2  1  2  2  70;
    2  1  2  2  1  1  73;
    2  2  1  1  2  1  62;
    2  2  1  2  1  2  67];
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(:,:))