som33.m

本程序研究在单个脉冲发现概率一定的条件下

%SOM神经网络，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;%向量的维数
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
% 
%