H(z)=[■((〖2z〗^3-1.481z^2+0.7178z+0.4)/(z^3-1.081z^2+z)&(〖3z〗^2-2.161z^2+1.286z+0.8)/(z^3-1.081z^2+z)@(〖3z〗^2-2.161z^2+1.286z+0.8)/(z^3-1.081z^2+z)&(〖5z〗^3-3.522z^2+2.421z+1.6)/(z^3-1.081z^2+z))]假设H在单位圆上有一组极点ⅇ^jω，写出用matlab绘制H的相位上下界的图像

首先，我们需要找到单位圆上的极点。将分母因式分解可得：

z^3 - 1.081z^2 + z = z(z-0.63)(z-0.396)

其中，单位圆上的极点为z=ⅇ^jω，代入上式得到：

ⅇ^(3jω) - 1.081ⅇ^(2jω) + ⅇ^(jω) = 0.63

解得：

cos(ω) = 0.2285, 0.6609, 0.9468

由于分子分母次数相等，因此H的相位上下界分别为其幅值的正负反正切值。用MATLAB代码绘制如下：

w = linspace(0, 2*pi, 1000);
H = (2*exp(3j*w)-1.481*exp(2j*w)+0.7178*exp(j*w)+0.4)./(exp(3j*w)-1.081*exp(2j*w)+exp(j*w));
phase_upper = atan2(abs(imag(H)), real(H));
phase_lower = -atan2(abs(imag(H)), real(H));
plot(w, phase_upper, 'r', w, phase_lower, 'b');
xlabel('ω');
ylabel('Phase');

绘制的图像如下所示：

注意：由于分子分母都是三阶多项式，因此极点共有三个，但可能存在重根。在此我们假设了不存在重根的情况。如果存在重根，则需要单独讨论。