%?clear??
%?clc??
%?
%?%生成初始解??
%?sol_new2=1;%（1）解空間（初始解）??
%?sol_new1=2-sol_new2^2;??
%?sol_current1?=?sol_new1;???
%?sol_best1?=?sol_new1;??
%?sol_current2?=?sol_new2;???
%?sol_best2?=?sol_new2;??
%?E_current?=?inf;??
%?E_best?=?inf;??
%?
%?rand（‘state‘sum（clock））;?%初始化隨機數發生器??
%?t=90;?%初始溫度??
%?tf=89.9;?%結束溫度??
%?a?=?0.99;?%溫度下降比例??
%?
%?while?t>=tf%（7）結束條件??
%?????for?r=1:100?%退火次數??
%?
%?????????%產生隨機擾動（3）新解的產生??
%?????????sol_new2=sol_new2+rand*0.2;??
%?????????sol_new1=2-sol_new2^2;??
%?
%?????????%檢查是否滿足約束??
%?????????if?sol_new1^2-sol_new2>=0?&&?-sol_new1-sol_new2^2+2==0?&&?sol_new1>=0?&&sol_new2>=0??
%?????????else??
%?????????????sol_new2=rand*2;??
%?????????????sol_new1=2-sol_new2^2;??
%?????????????continue;??
%?????????end??
%?
%?????????%退火過程??
%?????????E_new=sol_new1^2+sol_new2^2+8;%（2）目標函數??
%?????????if?E_new%?????????????????E_current=E_new;??
%?????????????????sol_current1=sol_new1;??
%?????????????????sol_current2=sol_new2;??
%?????????????????if?E_new%?????????????????????%把冷卻過程中最好的解保存下來??
%?????????????????????E_best=E_new;??
%?????????????????????sol_best1=sol_new1;??
%?????????????????????sol_best2=sol_new2;??
%?????????????????end??
%?????????else??
%?????????????????if?rand%?????????????????????E_current=E_new;??
%?????????????????????sol_current1=sol_new1;??
%?????????????????????sol_current2=sol_new2;??
%?????????????????else??
%?????????????????????sol_new1=sol_current1;??
%?????????????????????sol_new2=sol_current2;??
%?????????????????end??
%?????????end??
%?????????plot（rE_b