系统函数是1+5*z^(-1)+5*z^(-2)+z^(-3) ，给出Matlab中求该系统频率响应的代码

假设系统函数为H(z)，则其频率响应为H(w)，可以通过以下代码求解：

syms z w
H = 1 + 5*z^(-1) + 5*z^(-2) + z^(-3);  % 定义系统函数
Hw = subs(H, z, exp(1i*w));  % 将z用e^(jw)代替
Hw = simplify(Hw);  % 化简
absHw = abs(Hw);  % 求幅度响应
angleHw = angle(Hw);  % 求相位响应

其中，`simplify`函数用于化简表达式，`abs`函数用于求绝对值，`angle`函数用于求相位。