pu.m

来自「有兴趣学系matlab的朋友」· M 代码 · 共 16 行

function y=pu(u,t)
[m n]=size(u);
for i=1:n
    x=u(i);
    if x<=t/2
        y(i)=x/t-2*x^2/t+x;
    elseif x<=t
        y(i)=x/t-x+t/2;
    elseif x<=0.5
        y(i)=(4*x-2*t-4*x^2-2*t^2+2*t*x+t^2)/(2-2*t);
    elseif x<=(0.5+t/2)
        y(i)=(-2*t-2*t^2+2*t*x+1+t^2)/(2-2*t);
    else 
        y(i)=(4*x^2-t^2+2-4*x+t^2-2*t*x)/(2-2*t);
    end
end