shiwangzhibiao.m

这是数模中电梯概率模型源码

function y=shiwangzhibiao(u,t)
[m n]=size(u);

for i=1:n
    x=u(i);
    if x<=t/2
        y(i)=4*x^2/t-4*x^3/t+2*x^2-2*x+1-x/t;
    elseif x<=t
        y(i)=t*x+2*x^2/t-2*x^2+1-t/2-x/t;
    elseif x<=0.5
        y(i)=(-8*x^3+4*t*x^2-6*x+t^2+12*x^2-2*x*t^2-4*t*x+2)/(2-2*t);
    elseif x<=(0.5+t/2)
        y(i)=(-2*t^2*x+4*t*x^2+t^2-4*t*x+1)/(2-2*t);
    else 
        y(i)=(8*x^3-4*t*x^2+6*x+4*t*x-12*x^2-2*t)/(2-2*t);
    end
end