%SOM神經(jīng)網(wǎng)絡，20051210
%實驗類內距離和類內離差陣的關系
%======每5類DOA相差9，檢查對DOA的適應情況=DOA變化從3倍sita變化=== 識別率為   100 ========================
function yy=som33(bb)
aa=20;bb=bb;
[mm,m1,q]=som3(PP1,bb)
p=q;
m=mm;
grpcount=bb;
n=3;%向量的維數(shù)
N=aa*5;
% m0=[]
% for i=1:grpcount
%     mk=mean(p{i});
%     m0=[m0;mk]
% end

m0=m1
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
yy=J1
% %%%%%%
% 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))))
% J1
% 
%