MATLAB编写程序，能够接受从键盘输入的系统函数H（s）=s/((s+1)(s+2))的分子分母多项式系数向量，并绘制出系统的零极点图、系统的单位冲激响应、系统的幅度频率响应和相位频率相应的图形。分析从图形中可以看出该系统的零点和极点位置是多少？从时域和零极点分布特征两个方面说明该系统是否是稳定的系统？从频率响应特性上看，该系统具有何种滤波特性？

以下是MATLAB代码实现：

% 从键盘输入系统函数H(s)分子分母多项式系数向量
num = input('请输入H(s)的分子多项式系数向量：');
den = input('请输入H(s)的分母多项式系数向量：');

% 绘制系统的零极点图
figure(1)
zplane(num, den)
title('系统的零极点图')

% 计算系统的单位冲激响应
impulseResponse = impz(num, den);
n = length(impulseResponse);

% 绘制系统的单位冲激响应图形
figure(2)
stem(0:n-1, impulseResponse)
title('系统的单位冲激响应')

% 计算系统的幅度频率响应和相位频率响应
[w, mag, phase] = bode(tf(num, den));

% 将幅度频率响应从对数坐标转换为线性坐标
mag = 20*log10(mag);

% 绘制系统的幅度频率响应图形
figure(3)
semilogx(w, mag)
title('系统的幅度频率响应')
xlabel('频率(rad/s)')
ylabel('幅度(dB)')

% 绘制系统的相位频率响应图形
figure(4)
semilogx(w, phase)
title('系统的相位频率响应')
xlabel('频率(rad/s)')
ylabel('相位(deg)')

% 分析从图形中可以看出该系统的零点和极点位置是多少？
% 该系统的零点为0，极点为-1和-2

% 从时域和零极点分布特征两个方面说明该系统是否是稳定的系统？
% 该系统是稳定的系统，因为其所有极点都位于左半平面

% 从频率响应特性上看，该系统具有何种滤波特性？
% 该系统是一个带通滤波器，其通带范围在-1rad/s和-2rad/s之间，阻带范围在0rad/s和-3rad/s之外