用matlab求离散时间系统 y[n] +4y[n -1] + 2y[n - 2] + y[n - 3] = x[n] 的单位脉冲响应 h[k]，给出程序代码。

可以使用 MATLAB 的 `filter` 函数求解该离散时间系统的单位脉冲响应 h[k]。具体的程序代码如下：

% 系统的差分方程系数
a = [1 4 2 1];
b = 1;

% 求解单位脉冲响应
h = impz(b, a);

% 绘制单位脉冲响应的图像
stem(h);
xlabel('k');
ylabel('h[k]');
title('The Impulse Response of the System');

运行上述程序，即可得到该离散时间系统的单位脉冲响应 h[k]，并绘制出其图像。