julei.m

数模中的钻井问题的matlab源码

function y=julei(x,chidu)
[m n]=size(x);
x0=ones(1,n);
for m=1:n
    chushijulizhen=chushijulizhen(x);
    s=min(min(chushijulizhen));
    if s<chidu
        [a b]=find((chushijulizhen-s)==0);
        a=a(1);
        b=b(1);
        x(a)=(x(a)*x0(a)+x(b)*x0(b))/(x0(a)+x0(b));
        x(b)=[];
        x0(a)=x0(a)+x0(b);
        x0(b)=[];
    end
end
y=[x;x0];