%  本程序是基于人工免疫原理的天氣預報系統的MATLAB實現程序
%  實現了對6月份的氣象預報
%  具體算法請參考該文獻 
%  易世華2006年3月15日
clc
clear
load data2 
load data3
n=size(agf,1);
result=zeros(n,6);
result(:,1:3)=agf(:,13:15);
agnor=1/sqrt(12)*agf(:,1:12);
for i=1:n                            %第i個抗原  
    for c=1:8
        for j=1:nm(c)
            M(j,16,c)=1-sqrt((M(j,1:12,c)-agnor(i,:))*(M(j,1:12,c)-agnor(i,:))');
        end        
    end 
    [Mmax1 num1]=max(M(:,16,:));    %找出刺激值最大的記憶細胞
    [Mmax2 num2]=max(Mmax1(1,1,:));
    M(num1(1,1,num2),16,num2)=0;
     
    [Mmax3 num3]=max(M(:,16,:));
    [Mmax4 num4]=max(Mmax3(1,1,:));
    M(num3(1,1,num4),16,num4)=0;
   
    [Mmax5 num5]=max(M(:,16,:));
    [Mmax6 num6]=max(Mmax5(1,1,:));
    M(num5(1,1,num6),16,num6)=0;
    
    fore1=zeros(1,8);
    fore1(num2)=0.45;
    
    fore2=zeros(1,8);    
    fore2(num4)=0.35;
    
    fore3=zeros(1,8);
    fore3(num6)=0.2;
    
    fore=fore1+fore2+fore3;         %3個刺激值最大的記憶細胞投票決定天氣情況的預報值             
    [maxfore c]=max(fore);    
    switch c
         case 1
              result(i,4:6)=[0 0 0];                 
         case 2
              result(i,4:6)=[0 0 1];  
         case 3
              result(i,4:6)=[0 1 0];
         case 4
              result(i,4:6)=[0 1 1];    
         case 5
              result(i,4:6)=[1 0 0];    
         case 6
              result(i,4:6)=[1 0 1];    
         case 7
              result(i,4:6)=[1 1 0];    
         case 8
              result(i,4:6)=[1 1 1];        
    end 
end
today=(n-sum(bitxor(result(:,1),result(:,4))))/n
tomorrow=(n-sum(bitxor(result(:,2),result(:,5))))/n
tomorrow2=(n-sum(bitxor(result(:,3),result(:,6))))/n