数字信号处理-LSI系统的Z域频率响应

称为系统函数：H(z)=\frac{Y(z)}{X(z)}=\frac{\sum_{m=0}^{M} b_m z^{-m}}{1+\sum_{k=1}^{N} a_k z^{-k}}
频率响应函数：H(e^{j\omega})=\frac{Y(e^{j\omega})}{X(e^{j\omega})}=\frac{\sum_{m=0}^{M} b_m e^{-j\omega m}}{1+\sum_{k=1}^{N} a_k e^{-j\omega k}}
在Matlab环境中，在涉及系统函数的文献中，默认情况下分子与分母的各项系数被编码为向量序列：
即b = [b_0, b_1, \dots, b_M]，
a = [a_0, a_1, \dots, a_N]，
并且a_0 = 1。
如果a_0 \neq 1，
程序会自动将其归一化为a_0 = 1。

请对LSI系统h(n)={1,4,3,6,9,7}进行频率响应计算（例如，在实验报告中此处需要提供理论说明），并绘制其幅频特性和相频特性曲线

    clc
    clear
    close all
    b = [1,1,1];
    %b= [1,4,3,6,9,7];
    a = [1,0.5,-0.25];