模拟退火法【matlab】

function JG=MF(G,M,b,Mx)
% G    系数矩阵n*m
% M    反演矩阵m*1
% Mx   搜索范围【最小 最大】
% b    G*M=b   b为n*1
%  结果中最后一行为最终的平方和平均误差
% error1为直接计算结果
%example:
%  G=[1 -1;-1 1;2 2;-3 1];
%  b= [1; 2; 3; 4];
%  Mx=[-2 4];
%  M=[1;1];
%  MF(G,M,b,Mx)
%  error1=sum((G*(G\b)-b).^2)/length(b)
format long
L=length(M);
T0=1000;
Tend=0.10;
i=1;ZZ=[];%初始化最值记录矩阵

    while T0>Tend
        ML=[];%初始化记录
        for j=1:100
            M1=min(Mx)+(max(Mx)-min(Mx))*rand(L,1);
            f2=F(G,M1,b);f1=F(G,M,b);
            dc=(f2-f1);%求最小值
            if dc<0
                M=M1;%f2小，f2对应的（x1，y1）换到最小值对应的（x,y）
            elseif (exp(-1000*dc/T0))>=rand
                M=M1;
            end
            ML(:,j)=M;% 记录M链
            ML1(j)=F(G,M,b);% 记录M链的值
        end
        [a,b1]=min(ML1);
        M=ML(:,b1);
        error(i)=a;
        ZZ=[ZZ M];%记录最值及坐标
        T0=0.99*T0;
        i=i+1;
    end
    JG=[ZZ(:,end);error(end)];
plot(error,'g^')

 

function f=F(G,M,b)

f=sum((G*M-b).^2)/(length(b));

% 

% f1=4*x.^3+4*x.*y+2*y.^2-42*x-14;
% f2=4*y.^3+4*x.*y+2*x.^2-26*x-22;
% f=f1.^2+f2.^2;