x=0:0.05:2*pi;
c=cos(x);
s=sin(x);
ss=s(find(s.*(x>=0&x<=3)));
cc=c(find(c.*(x>=1.5&x<=4.5)));
xc=x(find(x.*(x>=1.5&x<=4.5)));
xs=x(find(x.*(x>=0&x<=3)));
s1=randn(size(ss))./10+ss;
c1=randn(size(cc))./10+cc;
plot(xc,c1,'ro');
hold on
plot(xs,s1,'bd');