function?H_total=ostu（threshold）
%?fxy=fstatus（xy）;
%?適應(yīng)值函數(shù)?xy為多維數(shù)組
%?對(duì)約束化問題，這里用懲罰函數(shù)的外點(diǎn)法，
%?懲罰函數(shù)法的基本思想是借助懲罰函數(shù)把約束化問題轉(zhuǎn)化為列約束化問題．
%?外點(diǎn)法　對(duì)違反約束的點(diǎn)在目標(biāo)函數(shù)中加入相應(yīng)的懲罰，而對(duì)可行點(diǎn)不懲罰．

%?求:max?f（xy）

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

w1=0;%定義參數(shù)
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;????%計(jì)算總的像素?cái)?shù)
%end
%for?a=1:254
?%???p（a）=f1（a）/s;
?%end???????????%計(jì)算各灰度值的概率
thr