使用impz函数编制程序求解下列系统的单位冲激响应并绘出其图形。 y(n)-0.1y(n-1)+0.2y(n-2)-0.3y(n-3) =2x(n)+0.3x(n-1)+2x(n-2) +3x(n-3)

你可以使用MATLAB的`impz`函数结合`lsim`函数来求解给定线性差分方程（LDE）的单位冲激响应，并通过`plot`函数绘制其图形。这里给出一个简单的步骤说明：

首先，将你的线性差分方程转换成传递函数的形式。对于给定的方程：
\[ y(n) - 0.1y(n-1) + 0.2y(n-2) - 0.3y(n-3) = 2x(n) + 0.3x(n-1) + 2x(n-2) + 3x(n-3) \]

它的Z变换对应于传递函数为：
\[ G(z) = \frac{Y(z)}{X(z)} = \frac{2 + 0.3z^{-1} + 2z^{-2} + 3z^{-3}}{1 - 0.1z^{-1} + 0.2z^{-2} - 0.3z^{-3}} \]

接下来，使用`tf`函数创建传递函数模型：

num = [2 0.3 2 3];   % 分子多项式系数
den = [1 -0.1 0.2 -0.3];   % 分母多项式系数
sys = tf(num, den);

然后，计算单位冲激响应（即 impulse response）：

ir = lsim(sys, ones(1, 100));  % 选择足够长的时间序列以观察完整响应

最后，画出响应图形：

plot(ir)
xlabel('Time Index')
ylabel('Response')
title('Unit Impulse Response of the System')

这将显示系统的单位冲激响应随时间变化的情况。