test2_2.m

来自「两点间的头镇试验源码」· M 代码 · 共 17 行

clear
t=linspace(0,3*pi,100);
y=sin(t);
y1=y;
y1(find(y1<-1/2))=-1/2;
y1(find(y1>sqrt(2)/2))=sqrt(2)/2;
plot(t,y,':k')
hold on
plot(t,y1,'-b')
legend('y=sin(t)','z=f(t)',0)
axis([0 10 -1 1])
set(gca,'xtick',[0:1:10],'ytick',[-1:.25:1])
xlabel('t'),ylabel('z=f(t)'),
title('限幅函数')
set(gcf,'color','w')
hold off