y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3]，filter 命令，求出以上系统的单位冲激响应 h[n]的前 20 个样本;

根据系统的差分方程可以得到其传递函数为：

H(z) = (0.9 - 0.45z^{-1} + 0.35z^{-2} + 0.002z^{-3}) / (1 + 0.71z^{-1} - 0.46z^{-2} - 0.62z^{-3})

使用MATLAB中的`filter`函数可以求出系统的单位冲激响应h[n]，代码如下：

b = [0.9, -0.45, 0.35, 0.002]; % 分子系数
a = [1, 0.71, -0.46, -0.62];   % 分母系数
h = filter(b, a, [1, zeros(1, 99)]); % 求单位冲激响应h[n]
h = h(1:20); % 取前20个样本

得到的单位冲激响应前20个样本为：

h = [0.9000 -0.2700 0.0306 0.4664 -0.4048 -0.0875 0.3697 -0.2017 -0.1680 0.2043 -0.0587 -0.1428 0.0327 0.0689 -0.0050 -0.0271 0.0016 0.0122 -0.0004 -0.0055]