exm0212_4.m

《matlab教程》配套光盘

%exm0212_4.m
t=linspace(0,3*pi,500);y=sin(t);%产生正弦波  
%从自变量着手进行逐段处理。
z1=((t<pi)|(t>2*pi)).*y;					%获得整流半波				<3>
w=(t>pi/3&t<2*pi/3)+(t>7*pi/3&t<8*pi/3);%关系逻辑运算和数值运算	<4>
w_n=~w;																%	<5>
z2=w*sin(pi/3)+w_n.*z1;						%获得削顶整流半波			<6>
clf
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])  
shg