% 定义omega的值
omega = 2 * pi; % 假设omega是2*pi，表示周期为1秒

% 幅度图数据点
freqs_amp = [-omega, omega]; % 频率
amps = [0.5, 0.5]; % 幅度

% 相位图数据点
freqs_phase = [-omega, omega]; % 频率
phases = [0, 0]; % 相位

% 绘制幅度图
subplot(2, 1, 1); % 在上方的子图位置
stem(freqs_amp, amps, 'b', 'filled'); % 'filled'使得标记实心
title('幅度图');
xlabel('频率 (rad/s)');
ylabel('幅度');
xlim([-3*omega, 3*omega]);
ylim([0, max(amps)+0.1]);
grid on;

% 绘制相位图
subplot(2, 1, 2); % 在下方的子图位置
stem(freqs_phase, phases, 'b', 'filled');
title('相位图');
xlabel('频率 (rad/s)');
ylabel('相位 (弧度)');
xlim([-3*omega, 3*omega]);
ylim([-pi/2, pi]);
grid on;

% 定义omega的值
omega = 2 * pi; % 假设omega是2*pi，表示周期为1秒

% 幅度图数据点
freqs_amp = [-omega, omega]; % 频率
amps = [0.5, 0.5]; % 幅度

% 相位图数据点
freqs_phase = [-omega, omega]; % 频率
phases = [pi/2, -pi/2]; % 相位

% 绘制幅度图
subplot(2, 1, 1); % 在上方的子图位置
stem(freqs_amp, amps, 'b', 'filled'); % 'filled'使得标记实心
title('幅度图');
xlabel('频率 (rad/s)');
ylabel('幅度');
xlim([-3*omega, 3*omega]);
ylim([0, max(amps)+0.1]);
grid on;

% 绘制相位图
subplot(2, 1, 2); % 在下方的子图位置
stem(freqs_phase, phases, 'b', 'filled');
title('相位图');
xlabel('频率 (rad/s)');
ylabel('相位 (弧度)');
xlim([-3*omega, 3*omega]);
ylim([-3*pi/2, 3*pi/2]); % 或者可以用 ylim([-pi, pi]); 但这将使得π/2在上限边缘
grid on;

% 设置图形大小
set(gcf, 'Position', [100, 100, 700, 500]);