y=-2:0.01:5;f=(0.7*(y+1).^2).*(1+(0.011+y*0.02)*8).*tanh(0.5./y*8).*(y>0)+0*(y<=0);df=diff(f);plot(y(2:end),df)按上面语句试试
