t=0:1:20 len=length(t) for k=1:len if (t(k)<=3) s(k)=0.5*0.0347*t(k).^2 else if ((3<t(k))&(t(k)<=17)) s(k)=0.104*t(k)-0.5*0.104*3 else if ((17<t(k))&(t(k)<=20)) s(k)=-0.5*0.0347*((t(k)-17).^2)+0.104*(t(k)-1.5) end end
t=0:1:20; len=length(t); s=zeros(len); for k=1:len if (t(k)<=3) s(k)=0.5*0.0347*t(k).^2; else if ((3<t(k))&&(t(k)<=17)) s(k)=0.104*t(k)-0.5*0.104*3; else %((17<t(k))&(t(k)<=20)) s(k)=-0.5*0.0347*((t(k)-17).^2)+0.104*(t(k)-1.5); end end end plot(t,s) 按上面语句试试