function?[XminGmin]=DIEDAI（x0max_step）
%該函數(shù)用于進行單純性替換，需要輸入初始點坐標x0和終止限zz
n=length（x0）;%問題的維度
[VGF]=CSDCX（x02）;%得到單純形矩陣與單純形頂點函數(shù)值矩陣棱長l為2
[GF?index]=sort（GF）;%對單純形頂點函數(shù)值排序??goalfunction
V=V（:index）;%按照函數(shù)值排序的順序對單純形矩陣排序
%下面進行迭代準備，求出初始單純形的最好，最差，次差及中心點
for?k=1:max_step
????Vbest=V（:1）;%最好點
????????GFbest=GF（1）;
????Vsoso=V（:n）;%次差點
????????GFsoso=GF（n）;
????Vworst=V（:n+1）;%最差點
????????GFworst=GF（n+1）;
????center?=?sum（V（:1:n）2）./?n;??%單純形的中心點
????%下面進行反射，通過中心點center反射最差點Vworst得到反射點Vr（refelct）
????Vr=2*center-Vworst;%得到反射點向量
????????GFr=feval（‘ZHIBIAO‘Vr）;%反射點函數(shù)值
????????if?GFr????????????Ve=center+2*（Vr-center）;%得到延伸點Ve（extend）
????????????????GFe=feval（‘ZHIBIAO‘Ve）;%得到延伸點函數(shù)值
????????????if?GFe<=GFbest%如果延伸點比最好電還要好，則以延伸點Ve代替原最好點Vbest；否則，以反射點Vr代替原最好點Vbest
????????????????V（:n+1）?=Ve;
????????????????GF（n+1）?=GFe;
????????????????%替換后需要重新形成單純形
????????????????????[GF?index]=sort（GF）;%排序
????????????????????V=V（:index）;%排序
????????????????%單純形更新完成
????????????else
????????????????V（:n+1）?=Vr;
????????????????GF（n+1）?=GFr;
?????????????????%替換后需要重新形成單純形
????????????????????[GF?index]=sort（GF）;%排序
????????????????????V=V（:index）;%排序
?????????????????%單純形更新完成
????????????end%延伸點判斷
????????else
????????????if?GFr????????????????V（:n+1）?=Vr;%用Vr代替Vworst
????????????????GF（n+1）?=GFr;
?????????????????%替換后需要重新形成單純形
????????????????????[GF?index]=sort（GF）;%排序
????????????????????V=V（:index）;%排序
?????????????????%單純形更新完成
????????????else%如果反射點不比最差點好
????????????????if?GFr>GFworst%如果反射點不比最差點好，則收縮，收縮因子為0.5，得到收縮點Vs（shink）
????????????????????Vs=center+0.5*（Vworst