bound.m

来自「various matlab programs to slove various」· M 代码 · 共 26 行

function [alfa,beta]=bound(dd,bb)
n=length(dd); 
alfa=dd(1)-abs(bb(1)); 
temp=dd(n)-abs(bb(n-1));
if (temp < alfa)
   alfa=temp; 
end
for i=2:(n-1)
    temp=dd(i)-abs(bb(i-1))-abs(bb(i));
    if (temp < alfa)
      alfa=temp; 
    end
end
beta=dd(1)+abs(bb(1)); 
temp=dd(n)+abs(bb(n-1));
if (temp > beta)
   beta=temp; 
end
for i=2:(n-1)
    temp=dd(i)+abs(bb(i-1))+abs(bb(i));
    if (temp > beta)
      beta=temp; 
    end
end
return