%SOM神經網絡，20051210
%實驗類內距離和類內離差陣的關系
%======每5類DOA相差9，檢查對DOA的適應情況=DOA變化從3倍sita變化=== 識別率為   100 ========================

[m0,m1,q]=som3(aa,bb)
p=q;
m=m0;
=m1;
grpcount=aa;
n=3;%向量的維數
N=aa*bb;
for i=1:grpcount
    m0{i}=mean(p{i});
end

Jw=0;
for i=1:grpcount
    cc=size(p{i})
    for j=1:cc(1)
        p0=p{i}
     Jw=Jw+(p0(j,:)-m0{i})*(p0(j,:)-m0{i})';   
    end
end
JB=0;
for i=1:grpcount
    JB=JB+(m0{i}-m)*(m0{i}-m)'
end
J1=JB/Jw
%%%%%%
Sw=zeros(n,n)；
for i=1:grpcount
    cc=size(p{i})
    for j=1:cc(1)
        p0=p{i}
     Sw=Sw+(p0(j,:)-m0{i})'*(p0(j,:)-m0{i});   
    end
end
Sw=1/N*Sw;
SB=zeros(n,n)；
for i=1:grpcount
     cc=size(p{i});
   SB=SB+cc(1)/N*(m0{i}-m)*(m0{i}-m)'
end
l=[1,0,0;0,1,0;0,0,1];
J2=sum(sum(l.*(SB*inv(Sw))))