mpbfx2.m

MATLAB科学计算与工程分析源代码源程序3

function [r1,r2,alpha,r]=mpbfx2(X1,X2,X)
X1=X1';
X2=X2';
miu1=mean(X1,2);
miu2=mean(X2,2);
[m,n1]=size(X1);
[m,n2]=size(X2);
[m,n]=size(X);
 for i=1:m
     ss11(i,:)=X1(i,:)-miu1(i); 
     ss12(i,:)=X1(i,:)-miu2(i); 
     ss22(i,:)=X2(i,:)-miu2(i);  
     ss21(i,:)=X2(i,:)-miu1(i); 
     ss2(i,:)=X(i,:)-miu2(i);  
     ss1(i,:)=X(i,:)-miu1(i); 
 end 
s1=ss11*ss11';
s2=ss22*ss22';
V1=(s1)/(n1-1);
V2=(s2)/(n2-1);
for j=1:n1
     r1(j)=ss12(:,j)'*inv(V2)*ss12(:,j)-ss11(:,j)'*inv(V1)*ss11(:,j);
end
for k=1:n2
     r2(k)=ss22(:,k)'*inv(V2)*ss22(:,k)-ss21(:,k)'*inv(V1)*ss21(:,k);
end
r1(r1>=0)=1;
r1(r1<0)=2;
r2(r2>=0)=1;
r2(r2<0)=2;
num1=n1-length(find(r1==1));
num2=n2-length(find(r2==2));
alpha=(num1+num2)/(n1+n2);
for l=1:n
    r(l)=ss2(:,k)'*inv(V2)*ss2(:,k)-ss1(:,k)'*inv(V1)*ss1(:,k);
end
r(r>0)=1;
r(r<0)=2;