（1）用MATLAB计算差分方程 y(n)+0.7y(n-1)-0.45y(n-2)-0.6y(n-3) =0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3) 当输入序列为x(n)=δ(n)时的输出结果y(n)，0≤ n≤40

在MATLAB中，你可以使用`ode45`函数来求解差分方程。该函数用于数值积分，特别适合处理微分和常微分方程。对于给定的线性非齐次差分方程，首先需要编写一个自定义函数描述方程，然后传递给`ode45`。

这里是一个简单的步骤：

1. 定义差分方程函数 `dydt.m` (假设输入为`t`, `y`, 和 `u`):

function dy = dydt(t,y,u)
    dy = [y(2); -0.7*y(1) + 0.45*y(2) + 0.6*y(3) - 0.8*u(1) + 0.44*u(2) - 0.36*u(3) + 0.02*u(4)];
end

2. 初始化时间向量 `tspan` 和状态向量 `y0`:

tspan = 0:0.1:40;
y0 = zeros(1, 4); % 因为有四个状态变量 y(n), y(n-1), y(n-2), y(n-3)

3. 输入序列 `x(n)` 为单位脉冲序列，可以表示为 `u` 向量：

u = dirac(tspan, 1); % 或者直接设置为 u = ones(size(tspan));

4. 使用 `ode45` 求解方程并得到输出序列 `y(n)`:

[t, y] = ode45(@dydt, tspan, y0, u);

5. 现在，`y(:,1)` 将包含从 `n=0` 到 `n=40` 的输出序列 `y(n)`。

注意：由于实际运行这段代码需要MATLAB环境，如果你是在文本环境中，我会展示大致思路，但无法在这里直接生成结果。如果你想要查看具体的结果，建议你在MATLAB环境下运行上述代码。