binsearch.m

可以方便的调用这个函数

%用折半搜索算法在单调递增值中寻找一个元素下标
function ia=binSearch(x,xhat)
n=length(x);
if xhat<x(1)|xhat>x(n)
    error(sprintf('Test value of %g is not in range of x' ,xhat'));
end
ia=1;ib=n;
while ib-ia>1
    im=fix((ia+ib)/2);
    if x(im)<xhat
        ia=im;
    else
        ib=im;
    end
end