然後這個循環：
for
i=N+1:N+length(n),
y(i)
=
-a1*y(i-N:i-1)'
+
b1*x(i-N:i-N+M)';
end
其實是因為：
y[n]
+
a1*y[n-1]
+
a2*y[n-2]...
+
an*y[n-N]
=
b0*x[n]
+
b1*x[n-1]
+
...
+
bm*x[n-M]
所以：
y[n]
=
-（a1*y[n-1]
+
a2*y[n-2]...
+
an*y[n-N]
）+
b0*x[n]
+
b1*x[n-1]
+
...
+
bm*x[n-M]
具體來說，就是：
我們已知了y1、y2、y3。。。yN，然後通過循環依次求得yN+1、yN+2等等。。。