kkxchangenew.m

来自「利用MATLAB编写」· M 代码 · 共 19 行

function newbujian=kkxchangenew(bujian,jiage,kekaoxing,r0)
%bujian  1*7
r1=0.001;
while (r1<r0)
    for j=1:5,
            if rand>0.5
                        bujian(j)=bujian(j)+1;
            else
                        bujian(j)=bujian(j)-1;
            end
            if bujian(j)<0  bujian(j)=0;  end  
    end
    r1=1;
    for j=1:5,
        r1=r1*(1-(1-kekaoxing(j))^bujian(j));
    end
end
bujian(6)=r1;bujian(7)=bujian(1:5)*jiage;
newbujian=bujian;