给出第一题对应的matlab代码

对于第一题中的 MATLAB 编程任务，可以按照以下步骤实现：

1. **定义状态空间模型** - 利用 `ss` 函数；
2. **绘制单位脉冲响应** - 使用 `impulse` 函数；
3. **使用 Routh-Hurwitz 准则评估稳定性** - 可通过 `routh` 或者手动构建劳斯表并检查其来完成；
4. **绘制系统的零极点图** - 利用 `pzmap` 函数。

以下是可能的一个解决方案的代码片段：

% 定义A, B, C矩阵
A = [0 1; 0 1];
B = [0; 1];
C = [1 0];
D = 0;

% 创建状态空间对象
sys = ss(A,B,C,D);

% 绘制单位脉冲响应曲线
figure;
impulse(sys);
title('Unit Impulse Response');
grid on;

% 检查稳定性（此处简化展示，未直接调用Routh相关工具）
p = eig(A); % 获取特征值
is_stable = all(real(p)<0); % 如果所有特征值的实部都小于0，则稳定
disp(['System is stable: ', num2str(is_stable)]);

% 绘制零极点分布图
figure;
pzmap(sys);
title('Pole-Zero Map');
grid on;

注意上述代码假设了给定的状态空间模型参数，实际应用中需要替换为具体的问题数据。对于稳定性检验部分，这里仅给出了通过特征值得到稳定性的简单判断逻辑。如果要严格依照Routh-Hurwitz准则操作，还需要进一步编写或查找相关的MATLAB函数。