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;