Для дискретного сигнала и его спектра в Scilab как раз используется функция plot2d3 из вашего списка справки — она строит вертикальные отрезки (импульсы/отсчеты) от нуля до значения точки.
Обычный plot соединяет точки плавной линией, что больше подходит для аналоговых сигналов, а plot2d3 наглядно показывает именно дискретную природу отсчетов.
Код построения всех 6 графиков с использованием plot2d3:
// 1. Вектор частот для первой половины спектра (от 0 до fs/2)
f = (0 : N/2 - 1) * (fs / N);
// 2. Очистка и разбиение окна (3 строки, 2 столбца)
clf();
// График 1: Сигнал x1(t)
subplot(3, 2, 1);
plot2d3(t, x1);
xtitle("Дискретный сигнал x1(t)", "t, с", "x1");
xgrid();
// График 2: Спектр (АЧХ) S1(f)
subplot(3, 2, 2);
plot2d3(f, S1(1:N/2));
xtitle("Спектр x1", "f, Гц", "|X1(f)|");
xgrid();
// График 3: Сигнал x2(t)
subplot(3, 2, 3);
plot2d3(t, x2);
xtitle("Дискретный сигнал x2(t)", "t, с", "x2");
xgrid();
// График 4: Спектр (АЧХ) S2(f)
subplot(3, 2, 4);
plot2d3(f, S2(1:N/2));
xtitle("Спектр x2", "f, Гц", "|X2(f)|");
xgrid();
// График 5: Суммарный сигнал x(t)
subplot(3, 2, 5);
plot2d3(t, x);
xtitle("Суммарный сигнал x(t)", "t, с", "x");
xgrid();
// График 6: Спектр (АЧХ) суммарного сигнала
subplot(3, 2, 6);
plot2d3(f, S(1:N/2));
xtitle("Спектр суммарного сигнала", "f, Гц", "|X(f)|");
xgrid();
Полезный нюанс:
Если на концах палочек нужны точки-маркеры (классический вид дискретного сигнала «леденцы на палочке»), сразу после строки с plot2d3 можно добавить вызов точечного графика, например:
plot2d3(t, x1);
plot(t, x1, "r."); // красные точки на вершинах отсчетов