exm03132_2.m

matlab程序

%exm03132_2.m
t=linspace(0,3*pi,500);y=sin(t); %产生正弦波 
%处理方法一，从自变量入手进行逐段处理
z1=((t<pi)|(t>2*pi)).*y;		 %获得整流半波
w=(t>pi/3&t<2*pi/3)+(t>7*pi/3&t<8*pi/3);%关系逻辑运算和数值运算
w_n=~w;															
z2=w*sin(pi/3)+w_n.*z1;					%获得削顶整流半波
shg,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])  
disp('按任意键，继续！'),pause
%处理方法二，从函数值入手进行逐段处理
z=(y>=0).*y;							%正弦整流半波
a=sin(pi/3);
z=(y>=a)*a+(y<a).*z;					%削顶的正弦整流半波

clf
plot(t,y,':r');hold on;plot(t,z,'-b')
xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数')
legend('y=sin(t)','z=f(t)',4),hold off 
disp('结束')