给出第2题对应的完整代码

对于题目中的第二题，涉及到的是一个连续系统的分析，这里提供一个MATLAB代码框架供参考：

% 定义符号变量
syms s t;

% 给定系统的传递函数 H(s)
num = [2 1]; % 分子多项式系数
den = [1 5 6]; % 分母多项式系数
H = tf(num, den);

% 初始条件
y_0 = 2;
yd_0 = 1;

% 输入信号 x(t)
x_t = exp(- 3);

% 频率范围设定
w = logspace(-2, 1, 500); 

% 绘制频率响应
bode(H, w);
title('系统的频响');
grid on;

% 使用 ilaplace 函数求解微分方程初值问题，找出状态变量初始条件
Y = ilaplace(laplace(y_0 + yd_0*t + int(x_t, t)), t, 0);

% 零输入响应 (仅考虑初始条件的影响)
zi_response = initial(H, y_0, yd_0, t);

% 零状态响应 (假设初始条件下为零，只考虑输入对输出的影响)
zs_response = lsim(H, x_t, t);

figure; % 创建一个新的图形窗口
subplot(2, 1, 1);
plot(t, zi_response);
title('零输入响应');
xlabel('时间(s)');
ylabel('幅值');

subplot(2, 1, 2);
plot(t, zs_response);
title('零状态响应');
xlabel('时间(s)');
ylabel('幅值');

% 计算全响应并比较两种方法的结果
full_response_method1 = zi_response + zs_response;
full_response_method2 = lsim(H, x_t, t, y_0, yd_0);

figure;
hold on;
plot(t, full_response_method1, 'r', 'DisplayName', 'Method 1');
plot(t, full_response_method2, 'g', 'DisplayName', 'Method 2');
legend show;
title('全响应');
xlabel('时间(s)');
ylabel('幅值');
grid on;

注意:
- `tf` 是用来定义传递函数的；
- `bode` 命令用于绘制系统的Bode图（频率响应）；
- `initial` 和 `lsim` 分别用于模拟零输入响应和零状态响应；
- `heaviside` 表示阶跃函数；
- 这段代码是示意性的，某些命令或步骤可能需要依据具体情况进行调整。

请确保在实际使用前对上述代码进行适当的修改和完善，使其符合您的特定需求及环境配置。