z2=w*sin(pi/3)+w_n.*z1;                   %   <6>

subplot(1,3,1),plot(t,y,':r'),ylabel('y')

subplot(1,3,2),plot(t,z1,':r'),axis([0 10 -1 1])

subplot(1,3,3),plot(t,z2,'-b'),axis([0 10 -1 1])  

图2.13-2 

%处理方法二：

z=(y>=0).*y;                %   <11>

a=sin(pi/3);

z=(y>=a)*a+(y<a).*z;     %   <13>