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

R1：系统函数：
H(z)=Y(z)/X(z) =(∑_(m=0)^M▒〖b_m z^(-m) 〗)/(1+∑_(k=1)^N▒〖a_k z^(-k) 〗)
R2：频率响应：
H(e^jω )=Y(e^jω )/X(e^jω ) =(∑_(m=0)^M▒〖b_m e^(-jωm) 〗)/(1+∑_(k=1)^N▒〖a_k z^(-jωk) 〗)
R3：在Matlab中的表示：
在有关Matlab的系统函数中，分子和分母的系数被定义为向量，即
b=[b_0,b_1,⋯,b_M ]
a=[a_0,a_1,⋯,a_M ]
且a_0=1。如果a_0≠1，则程序将自动地将其归一化为1。

求LSI系统h(n)={▁1,4,3,6,9,7}的频率响应（例如，在实验报告中，此处是需要提供理论说明的），画出它的幅频特性及相频特性。
matlab代码如下：

clc
clear
close all
b = [1,1,1