部分代码： function Population1=GA_tubian(Population,pe_tubian) %遗传算法突变算子 %pe为突变概率 Population1=Population; n=length(Population(:,1)); m=length(Population(1,:)); for i=1:n for j=1:m test=rand; if test<pe_tubian Population1(i,j)=1-Popula