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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%調試用語句
% popsize=16;
% n=100000;
% pop1_init=randperm(n); %產生1到32的整數
% pop2_init=randperm(n);
% pop(1,:)=pop1_init(1:popsize);  %將前幾個整數作為初始群體
% pop(2,:)=pop2_init(1:popsize);;  %將前幾個整數作為初始群體
% 
% 
% 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(i)
            newpop_con(i,j)=up_range(i);
        elseif newpop_con(i,j)<down_range(i)
            newpop_con(i,j)=down_range(i);
        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;