% 2.6 變異
%遺傳算法子程序
%Name: mutation.m
%變異
function [newpop]=complete3_mutation(popsize,pop,pm,i,fun_num,var_num)
%函數(shù)說明
%入口參數(shù)：pop         當前處理的種群,
%          popsize     種群大小,
%          pm          變異概率
%          i           當前處理代數(shù)
%          fitvalue    當前代的適應度數(shù)組
%          fun_num     函數(shù)序列號,
%          up_range    自變量的相對上限,
%          down_range  自變量的相對下限
%出口參數(shù)：newpop      準子代
%功能說明：在相對定義域的個體上加地位的高斯擾動，實現(xiàn)小范圍內的變異。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%調試用語句
% popsize=16;
% n=100000;
% pop1_init=randperm(n); %產生1到32的整數(shù)
% pop2_init=randperm(n);
% pop(1,:)=pop1_init(1:popsize);  %將前幾個整數(shù)作為初始群體
% pop(2,:)=pop2_init(1:popsize);;  %將前幾個整數(shù)作為初始群體
% 
% 
% pop;
% i=5;

% popsize=10;
% pop=[1023,980,2133,726,1946,657,1702,1224,1522,807;
%      2025,2565,1773,2340,2354,1144,2061,240,1573,825]
% chromlength=12;
% pm=0.7;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% dai_t=i*10^i;
dai_t = i;
% up_range=5;
% down_range=-5;

% max_gauss=1/(sqrt(2*pi*1/(dai_t)));
pop_con = pop;
switch fun_num
    case 1
%          pop_con=(pop-10000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=100;
         down_range=-100;

    case 2
%          pop_con=(pop-1000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=10;
         down_range=-10;
    
    case 3
%          pop_con=(pop-10000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=100;
         down_range=-100;
         
    case 4
%          pop_con=(pop-10000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=100;
         down_range=-100;
         
    case 5
%          pop_con=(pop-30000)./1000;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=30;
         down_range=-30;
         
    case 6
%          pop_con=(pop-10000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=100;
         down_range=-100;
         
    case 7
%          pop_con=(pop-12800)./10000;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=1.28;
         down_range=-1.28;
    
    case 8
%          pop_con=(pop-50000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=500;
         down_range=-500;
         
    case 9
%          pop_con=(pop-51200)./10000;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=5.12;
         down_range=-5.12;
         
    case 10
%          pop_con=(pop-32000)./1000;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=32;
         down_range=-32;
         
    case 11
%          pop_con=(pop-60000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=600;
         down_range=-600;
         
    case 12
%          pop_con=(pop-10000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=100;
         down_range=-100;

    case 13
%          pop_con=(pop-10000)./100;  %將大定義域的值轉變到 實際的定義域中，便于進行高斯擾動變異
         up_range=100;
         down_range=-100;
    otherwise
%         break
          return;
end
    
% var_num = 10;
for i=1:var_num
    for j=1:popsize
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 對個體的后位進行高斯擾動變異
        if rand<pm
           newpop_con(i,j)=pop_con(i,j);
        else
%            newpop_con(i,j)=pop_con(i,j)+awgn(0,dai_t);
           newpop_con(i,j)=pop_con(i,j)+normrnd(0,1/dai_t);
        end
       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 判斷新生成的個體是否超界，如果超界則加減高斯擾動的最大值
        
        if newpop_con(i,j)>up_range
            newpop_con(i,j)=up_range;
        elseif newpop_con(i,j)<down_range
            newpop_con(i,j)=down_range;
        end
        newpop = newpop_con;
%         switch fun_num
%               case 1
%                    newpop(i,j)=newpop_con(i,j)*100+10000;
%               case 2
%                    newpop(i,j)=newpop_con(i,j)*100+1000;
%               case 3
%                    newpop(i,j)=newpop_con(i,j)*100+10000;
%               case 4
%                    newpop(i,j)=newpop_con(i,j)*100+10000;
%               case 5
%                    newpop(i,j)=newpop_con(i,j)*1000+30000;
%               case 6
%                    newpop(i,j)=newpop_con(i,j)*100+10000;
%               case 7
%                    newpop(i,j)=newpop_con(i,j)*10000+12800;
%               case 8
%                    newpop(i,j)=newpop_con(i,j)*100+50000;
%               case 9
%                    newpop(i,j)=newpop_con(i,j)*10000+51200;
%               case 10
%                    newpop(i,j)=newpop_con(i,j)*1000+32000;
%               case 11
%                    newpop(i,j)=newpop_con(i,j)*100+60000;
%               case 12
%                    newpop(i,j)=newpop_con(i,j)*100+10000;
%               case 13
%                    newpop(i,j)=newpop_con(i,j)*100+10000;
%               otherwise
%                    break
%          end
       
    
    end

end



% newpop;