matlab求差分方程通解,【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法...

例题：

一个3阶低通滤波器由下面差分方程描述：

y(n) = 0.0181 x(n) + 0.0543 x(n-1) + 0.0543 x(n-2) + 0.0181 x(n-3) + 1.76 y(n-1) - 1.1829 y(n-2) + 0.2781 y(n-3)

画出这个滤波器的幅度和相位响应，并验证它是一个低通滤波器。

第一种方法是博文里给出的：

第二个案例中，类比，如果我知道一个LTI系统的脉冲响应h(n)，那么我也能求出它的频率响应：

k = [0:M];

n = [n1:n2];

X = x * (exp(-j * pi/M)).^(n'*k);

博文里面由具体的推荐，上述程序中X就是频率响应，也就是x的DTFT，如果x换成h，则X可以换成H。

第二种方法是，通过差分方程直接求出系统的频率响应，求解的方法是通过向量化的方法。

设某一LTI系统的差分方程表示为：

可以用一种简单的矩阵向量乘法来完成。如果在