Using stem plot for two plots

I have very basic question regarding drawing stem plot for this example, it is supposed to have 21 data values.

T=4e-3; n=21;

x=linspace(0,T/2,2001);

x1=linspace(T/2,T,2001);

y=sin(pi*x/T);

y1=(1+cos(pi*x1/T));

plot(x,y,'b',x1,y1);

grid on; zoom on; hold on

Answer 1

[Voss]

T=4e-3; n=21;
% x=linspace(0,T/2,2001);
% x1=linspace(T/2,T,2001);
x=linspace(0,T/2,(n+1)/2);
x1=linspace(T/2,T,(n+1)/2);
y=sin(pi*x/T);
y1=(1+cos(pi*x1/T));
% plot(x,y,'b',x1,y1);
stem(x,y,'b');
hold on
stem(x1,y1);
grid on; zoom on;

[Voss]

Or

T=4e-3; n=21;
% x=linspace(0,T/2,2001);
% x1=linspace(T/2,T,2001);
x=linspace(0,T/2,(n+1)/2);
x1=linspace(T/2,T,(n+1)/2);
y=sin(pi*x/T);
y1=(1+cos(pi*x1/T));
% plot(x,y,'b',x1,y1);
stem(x1,y1,'Color',[0.85 0.325 0.098]);
hold on
stem(x,y,'b');
grid on; zoom on;

[Arkadius882]

This one works which I'm very thankful of, but the plot should also be visible. I've been trying to connect these two but the plot must be drawn few thousand times (academic task) and because of that vectors are not the same length, at least it sounds like it to me.

T=4e-3; n=21;

x=linspace(0,T/2,2001);

xd=linspace(0,T/2,(n+1)/2);

x1=linspace(T/2,T,2001);

xd1=linspace(T/2,T,(n+1)/2);

y=sin(pi*x/T);

y1=(1+cos(pi*x1/T));

plot(x,y,'b',x1,y1);

grid on; zoom on; hold on;

stem(xd,y,'r*');

stem(xd1,y1,'r*');

[Voss]

T=4e-3; n=21;
x=linspace(0,T/2,2001);
xd=linspace(0,T/2,(n+1)/2);
x1=linspace(T/2,T,2001);
xd1=linspace(T/2,T,(n+1)/2);
y=sin(pi*x/T);
yd=sin(pi*xd/T);
y1=(1+cos(pi*x1/T));
yd1=(1+cos(pi*xd1/T));
plot(x,y,'b',x1,y1);
grid on; zoom on; hold on;
stem(xd,yd,'r*');
stem(xd1,yd1,'r*');

[Arkadius882]

Thank you!

[Voss]

You're welcome!

Answer 2

[Steven Lord]

This looks like a reasonable plot to me.

T=4e-3; n=21;
x=linspace(0,T/2,2001);
x1=linspace(T/2,T,2001);
y=sin(pi*x/T);
y1=(1+cos(pi*x1/T));
plot(x,y,'b',x1,y1);
grid on; zoom on; hold on

Let's plot it as a stem plot instead.

figure
stem([x,x1], [y,y1]);

There are too many stems to show each one as an individual line. Let's plot every 100 stems and see how that looks.

figure
stem([x(1:100:end), x1(1:100:end)], [y(1:100:end), y1(1:100:end)])

That looks nice to me.