利用matab编写程序，能够接受从键盘输入的系统函数的分子分母多项式系数向量，并绘制出系统的零极点图、系统的单位冲激响应、系统的幅度频率响应和相位频率相应的图形。

以下是一个利用matlab实现的程序，可以接受从键盘输入的系统函数的分子分母多项式系数向量，并绘制出系统的零极点图、系统的单位冲激响应、系统的幅度频率响应和相位频率响应的图形。

% 输入系统函数的分子分母多项式系数向量
num = input('Enter the numerator coefficients: ');
den = input('Enter the denominator coefficients: ');

% 绘制系统的零极点图
figure;
zplane(num, den);
title('Pole-Zero Plot');

% 绘制系统的单位冲激响应
figure;
impz(num, den);
title('Unit Impulse Response');

% 绘制系统的幅度频率响应
figure;
freqz(num, den, 'whole');
title('Magnitude Frequency Response');

% 绘制系统的相位频率响应
figure;
[h,w] = freqz(num, den, 'whole');
plot(w, unwrap(angle(h)));
title('Phase Frequency Response');

注意：运行程序时，需要先输入系统函数的分子分母多项式系数向量。例如，如果系统函数为 s^2 + 2s + 1 / s^2 + s，则输入的分子分母多项式系数向量为 [1 2 1] 和 [1 1]。