function?H_total=ostu（threshold）
%?fxy=fstatus（xy）;
%?適應值函數?xy為多維數組
%?對約束化問題，這里用懲罰函數的外點法，
%?懲罰函數法的基本思想是借助懲罰函數把約束化問題轉化為列約束化問題．
%?外點法　對違反約束的點在目標函數中加入相應的懲罰，而對可行點不懲罰．

%?求:max?f（xy）

global?G_min?G_max
fvarsize=size（threshold）;
H_total=zeros（fvarsize（1）1）;
global?p

w1=0;%定義參數
w2=0;
u1=0;
u2=0;
u=0;
m=1;

%f1=csvread（‘d:\f1.txt‘）;
%p=f1;
%for?a=1:254
%s=f1（a）+s;????%計算總的像素數
%end
%for?a=1:254
?%???p（a）=f1（a）/s;
?%end???????????%計算各灰度值的概率
thr