function fit=Fitness_Function(Pos,x,y0)
    yp=fun1(Pos,x);%计算参数拟合下的值Pos为两个参数k1,k2的值,x为x1,x2,x3的值?
    %yp为在该粒子下计算得y值，记为yp与y同维度，因此计算残差平方和视为适应度，使其最小
    fit=sum((yp-y0).^2);
  

end