clear,clc
%這是退火算法的主程序，它需要調用的函數是
%函數(1)nonLinearSumError1:計算非線性方程組總誤差的函數
%函數(2)newSolver1:在一組解的鄰域產生另一組解
%函數(3)isSolution:驗證方程是否得解
%設置初始值
i=0;T=10001;j=0;%i:同一溫度下狀態轉移次數;T:溫度；j:下降溫度
precision=0.1;
x1Group=1;%x1Group:可能解的組數
x1N=4;%非線性方程組的元數
x1=round((-0.5+rand(x1Group,x1N))*20);%隨機生成-10~10之間的初解
errorHold=Inf;
xHold=0;
%x1=[-7 5 1 -3];
i=0;
while i<200
i=i+1;
j=0;
T=T-50;%退火
while j<200
j=j+1;
functionError1=nonLinearSumError1(x1);%計算x1的誤差
x2=newSolver1(x1,functionError1,-10,1,10);%在x1的鄰域生成新一組解x2
functionError2=nonLinearSumError1(x2);%計算x2的誤差
%檢查方程是否得解
[solution1,minError1,isTrue1]=isSolution(x1,functionError1,precision);
[solution2,minError2,isTrue2]=isSolution(x2,functionError2,precision);
if isTrue1==1
%'方程得解'
functionError1
solutiourn
i,j
return
elseif isTrue2==1
%'方程得解'
solution2
functionError2
i,j
return
end
%x1
%x2
if functionError2-functionError1<0
x1=x2;%x2比x1好，用x2取代x1
elseif errorHold-functionError2<0
%x1=xHold;
else
p_x2x1=exp(-log(functionError2-functionError1)/T);
%狀態轉移概率，注意：誤差取對數，因為要解的非線性方程組比較復雜，
%可能解的一點偏差會引起方程很大的變化。所以通過取對數縮小差距。
if rand(1)<p_x2x1 %狀態轉移
xHold=x1;%hHold:把比較好的解保留下來
errorHold=functionError1;%比較好的解對應的誤差
x1=x2;
end
end

end
end
solution1
functionError1
solution2
functionError2