orthexpanalysis3.m

Matlab在化学工程中的应用

function OrthExpAnalysis3
% 正交试验的方差分析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: 2002/08/18 $
%
%   [Ref] 胡上序和陈德钊编著，实验数据的分析和处理，杭州：浙江大学出版社，
%   1996，p.273
%
%   DATA: 正交试验方案及结果---最后一列为试验指标（收率）,其余是以水平表示的正交试验方案
%   NF: 因素数（列数）
%   NL: 水平数
%   r:  每个水平有r次试验

clear all
clc

DATA = ...              % 最后一列为试验指标（收率）,其余是以水平表示的正交试验方案
   [1  1  1  1  31;
    1  2  2  2  54;
    1  3  3  3  38;
    2  1  2  3  53;
    2  2  3  1  49;
    2  3  1  2  42;
    3  1  3  2  57;
    3  2  1  3  62;
    3  3  2  1  64];
[m,n] = size(DATA);
NF = n-1;              % 因素数（列数）
NL = max(DATA(:,1))    % 水平数
r = 0;
for i=1:m;
    if DATA(i,1)==1
        r = r+1;        % 每个水平有r次试验
    end
end

A = DATA(:,1:NF);
Y = DATA(:,NF+1);
% 计算每个因素（列）相同水平的指标总和K
for i=1:NL              % i:水平
    for j=1:NF          % j:因素（列号）
        K(i,j) = sum(Y(find(A(:,j)==i)))
    end
end

xa = sum(Y)/length(Y)
xT = sum(Y)
P = 1/(NL*r)*xT^2
j=1:NF;
Q(j) = sum(K(:,j).^2)/r
S = Q - P               % 各因素的平方和
ST = sum(S)             % 总平方和
Sav = S/(NL-1)          % 各因素的均方
% 显著性检验
F = Sav(1:NF-1)/Sav(end)